我有与此主题相关的新问题 。 在新的情况下,变量x和x1具有不同的长度
x <- c(-10, 1:6, 50)
x1<- c(-20, 1:5, 60)
z<- c(1,2,3,4,5,6,7,8)
bx <- boxplot(x)
bx$out
bx1 <- boxplot(x1)
bx1$out
x<- x[!(x %in% bx$out)]
x1 <- x1[!(x1 %in% bx1$out)]
x_to_remove<-which(x %in% bx$out)
x <- x[!(x %in% bx$out)]
x1_to_remove<-which(x1 %in% bx1$out)
x1 <- x1[!(x1 %in% bx1$out)]
z<-z[-unique(c(x_to_remove,x1_to_remove))]
z
data.frame(cbind(x,x1,z))
然后我收到警告
Warning message:
In cbind(x, x1, z) :
number of rows of result is not a multiple of vector length (arg 2)
所以在新的数据框架中。 Z的不对应于x和x1。 我该如何判断这个问题? 这种解决对我没有帮助 deleting outlier in r with account of nominal var 或者我只是做错了什么。
x_to_remove<-which(x %in% bx$out)
x <- x[!(x %in% bx$out)]
x1_to_remove<-which(x1 %in% bx1$out)
x1 <- x1[!(x1 %in% bx1$out)]
z<-z[-unique(c(x_to_remove,x1_to_remove))]
z
d=data.frame(cbind(x,x1,z))
d
这是错的 警告信息:
In cbind(x, x1, z) :
number of rows of result is not a multiple of vector length (arg 2)
d
x x1 z
1 1 1 2
2 2 2 3
3 3 3 4
4 4 4 5
5 5 5 6
6 6 1 2
这3个columg如何获得此输出
Na Na Na
1 1 2
2 2 3
3 3 4
4 4 5
5 5 6
Na Na Na
Na Na Na
六行(d)是多余的
答案 0 :(得分:1)
原始x,x1和z列表中的差异长度是第一个问题,你怎么能说哪个z值与每个x和x1值有关?
x <- c(-10, 1:6, 50)
x1<- c(-20, 1:5, 60)
z<- c(1,2,3,4,5,6,7,8)
length(x)
[1] 8
length(x1)
[1] 7
length(z)
[1] 8
另一个问题是:
x<- x[!(x %in% bx$out)] #remove this
x1 <- x1[!(x1 %in% bx1$out)] #remove this
x_to_remove<-which(x %in% bx$out)
x <- x[!(x %in% bx$out)]
x1_to_remove<-which(x1 %in% bx1$out)
x1 <- x1[!(x1 %in% bx1$out)]
您在计算x
和x1
之前清除x_to_remove
和x1_to_remove
编辑: 要获得所需的输出,请尝试使用此代码(/ ode行添加在注释中签名):
x <- c(-10, 1:6, 50)
x1<- c(-20, 1:5, 60)
z<- c(1,2,3,4,5,6,7,8)
length_max<-min(length(x),length(x1),length(z)) #Added: identify max length before outlier detection
bx <- boxplot(x)
bx1 <- boxplot(x1)
x_to_remove<-which(x %in% bx$out)
x <- x[!(x %in% bx$out)]
x1_to_remove<-which(x1 %in% bx1$out)
x1 <- x1[!(x1 %in% bx1$out)]
z<-z[-unique(c(x_to_remove,x1_to_remove))]
length_min<-min(length(x),length(x1),length(z)) #Minimum length after outlier remove
d=data.frame(cbind(x[1:length_min],x1[1:length_min],z[1:length_min])) #Bind columns
colnames(d)<-c("x","x1","z")
d_NA<-as.data.frame(matrix(rep(NA,(length_max-length_min)*3),nrow=(length_max-length_min))) #Create NA rows
colnames(d_NA)<-c("x","x1","z")
d<-rbind(d,d_NA) #Your desired output
d
x x1 z
1 1 1 2
2 2 2 3
3 3 3 4
4 4 4 5
5 5 5 6
6 NA NA NA
7 NA NA NA