컷() 오류 - '나누기'

Question

내가 dataframe 다음 한 고유하지 않은 :이 dataframe에 다음 코드를 적용 할 때

a   
    ID a.1 b.1  a.2 b.2 
1 1 40.00 100.00 NA 88.89 
2 2 100.00 100.00 100 100.00 
3 3 50.00 100.00 75 100.00 
4 4 66.67 59.38 NA 59.38 
5 5 37.50 100.00 NA 100.00 
6 6 100.00 100.00 100 100.00 

가 :

이

temp <- do.call(rbind,strsplit(names(df)[-1],".",fixed=TRUE)) 
dup.temp <- temp[duplicated(temp[,1]),] 

res <- lapply(dup.temp[,1],function(i) { 
breaks <- c(-Inf,quantile(a[,paste(i,1,sep=".")], na.rm=T),Inf) 
cut(a[,paste(i,2,sep=".")],breaks) 
}) 

컷() 함수가 있습니다 오류 :

Error in cut.default(a[, paste(i, 2, sep = ".")], breaks) : 
'breaks' are not unique 

그러나 동일한 코드가 비슷한 dataframe에 완벽하게 잘 작동합니다

varnames<-c("ID", "a.1", "b.1", "c.1", "a.2", "b.2", "c.2") 

a <-matrix (c(1,2,3,4, 5, 6, 7), 2,7) 

colnames (a)<-varnames 

df<-as.data.frame (a) 


    ID a.1 b.1 c.1 a.2 b.2 c.2 
    1 1 3 5 7 2 4 6 
    2 2 4 6 1 3 5 7 

res <- lapply(dup.temp[,1],function(i) { 
breaks <- c(-Inf,quantile(a[,paste(i,1,sep=".")], na.rm=T),Inf) 
cut(a[,paste(i,2,sep=".")],breaks) 
}) 

res 
[[1]] 
[1] (-Inf,3] (-Inf,3] 
Levels: (-Inf,3] (3,3.25] (3.25,3.5] (3.5,3.75] (3.75,4] (4, Inf] 

[[2]] 
[1] (-Inf,5] (-Inf,5] 
Levels: (-Inf,5] (5,5.25] (5.25,5.5] (5.5,5.75] (5.75,6] (6, Inf] 

[[3]] 
[1] (5.5,7] (5.5,7] 
Levels: (-Inf,1] (1,2.5] (2.5,4] (4,5.5] (5.5,7] (7, Inf] 

왜이 오류가 발생 했습니까? 어떻게 고칠 수 있습니까? 감사합니다. 열 b.1, a.2 및 b.2에 대한 데이터 분위수 값이 어떤 수준에 대해 동일하기 때문에

Answer 1

는이 오류를 얻을, 그래서 그들은 직접 기능 cut()에서 휴식 값으로 사용할 수 없습니다. 이 문제를 해결하기

apply(a,2,quantile,na.rm=T) 
     ID  a.1 b.1 a.2  b.2 
0% 1.00 37.5000 59.38 75.0 59.3800 
25% 2.25 42.5000 100.00 87.5 91.6675 
50% 3.50 58.3350 100.00 100.0 100.0000 
75% 4.75 91.6675 100.00 100.0 100.0000 
100% 6.00 100.0000 100.00 100.0 100.0000 

한 가지 방법은 unique() 함수 내에서 quantile()를 넣어하는 것입니다 - 당신이 고유하지 않은 모든 분위수 값을 제거 할 수 있도록. quantiles가 고유하지 않은 경우 물론 이것은 물론 단점을 덜 만듭니다.

res <- lapply(dup.temp[,1],function(i) { 
    breaks <- c(-Inf,unique(quantile(a[,paste(i,1,sep=".")], na.rm=T)),Inf) 
    cut(a[,paste(i,2,sep=".")],breaks) 
}) 

[[1]] 
[1] <NA>  (91.7,100] (58.3,91.7] <NA>  <NA>  (91.7,100] 
Levels: (-Inf,37.5] (37.5,42.5] (42.5,58.3] (58.3,91.7] (91.7,100] (100, Inf] 

[[2]] 
[1] (59.4,100] (59.4,100] (59.4,100] (-Inf,59.4] (59.4,100] (59.4,100] 
Levels: (-Inf,59.4] (59.4,100] (100, Inf] 

Answer 2

오히려 quantile의 수를 유지하려면 다른 지터를 추가하는 것이 좋습니다.

breaks = c(-Inf,quantile(a[,paste(i,1,sep=".")], na.rm=T),Inf) 
breaks = breaks + seq_along(breaks) * .Machine$double.eps 

Answer 3

컷 대신, 고유하지 않은 끊김 벡터를 허용하는 .bincode를 사용할 수 있습니다. 당신이 분위, 분위수 등 당신이 실제로 인구의 10 % 또는 25 % 부분과 여기 분위/분위 버킷 아닌 실제 숫자 값을 의미하는 경우 말할 때

Answer 4

당신이 사용할 수있는 것입니다 :

a <- c(1,1,1,2,3,4,5,6,7,7,7,7,99,0.5,100,54,3,100,100,100,11,11,12,11,0) 
ar<-rank(a,ties.method = "first") 
decile <- cut(ar, quantile(ar, probs=0:10/10), include.lowest=TRUE, labels=FALSE)