假设您有一个data.frame
,其中包含多个不同级别的因素:
V1<-factor(sample(c(1:5,9),100,TRUE))
V2<-factor(sample(c(1:5,9),100,TRUE))
V3<-factor(sample(c(1:5),100,TRUE))
V4<-factor(sample(c(1:5),100,TRUE))
dat<-data.frame(V1,V2,V3,V4)
目标是估计两个因素的水平频率差异。但是,由于级别数不同,基于V1 / V2和V3 / V4的两个表中的数组不一致,例如:
table(dat$V1)-table(dat$V3)
Error in table(dat$V1) - table(dat$V3) : non-conformable arrays
目标是使V3和V4一致,以便操作有效。一种选择是:
dat$V3<-factor(dat$V3,levels=c('1','2','3','4','5','9')
然而,它需要为每个变量设置因子水平,这对许多变量V5,...,Vn来说是不切实际的。我想
dat[,3:4]<-apply(dat[,3:4],2,factor,levels=c('1','2','3','4','5','9'))
可能更通用,但is.factor(dat$V3)
则为FALSE。
编辑:此功能可能会完成SimonO101的答案:
correct_factors<-function(df_object,range){
if(is.data.frame(df_object)==FALSE){stop('Requires data.frame object')}
levs <- unique( unlist( lapply( df_object[,range[1]:range[2]] , levels ) ) )
df_object[,range[1]:range[2]] <-
data.frame( lapply( df_object[,range[1]:range[2]] , factor , levels = levs ) )
return(df_object)
}
答案 0 :(得分:4)
尝试这样来协调各个级别......
# Get vector of all levels that appear in the data.frame
levs <- unique( unlist( lapply( dat , levels ) ) )
# Set these as the levels for each column
dat2 <- data.frame( lapply( dat , factor , levels = levs ) )
table(dat2$V1)-table(dat2$V3)
# 1 2 3 4 5 9
#-15 -5 4 7 -5 14