【发布时间】:2014-12-22 03:48:07
【问题描述】:
我正在开展一个项目来模拟缺失数据并对采样数据进行回归。 这是我到目前为止所拥有的。 图书馆(MASS)
#specifying the covariance matrix
sigma <- matrix(c(1,.7,.49,.343,.2401,.7,1,.7,.49,.343,.49,.7,1,.7,.49,.343,.49,
.7,1,.7,.2401,.343,.49,.7,1),5,5,byrow=TRUE)
#generating the data
data <- mvrnorm(n=1000, c(5,5.25,5.5,5.75,6), sigma)
split(data,)
#specifying the missing data mechanism for MCAR
LogoddsratioMCAR <- -.5
OddsRatioMCAR <-exp(LogoddsratioMCAR)
OddsMCAR <- OddsRatioMCAR/(1+OddsRatioMCAR)
Probability2 <- 1-OddsMCAR
Probability3 <- Probability2 - OddsMCAR*(Probability2)
Probability4 <- Probability3 - OddsMCAR*(Probability3)
Probability5 <- Probability4 - OddsMCAR*(Probability4)
#sampling from each column
dataframe <- as.data.frame(data)
dataMCAR1 <- dataframe$V1
dataMCAR2 <- dataframe$V2[sample(1:nrow(data),Probability2*nrow(data))]
dataMCAR3 <- dataframe$V3[sample(1:nrow(data),Probability3*nrow(data))]
dataMCAR4 <- dataframe$V4[sample(1:nrow(data),Probability4*nrow(data))]
dataMCAR5 <- dataframe$V5[sample(1:nrow(data),Probability5*nrow(data))]
现在我需要将 NA 添加到 dataMCAR2-dataMCAR5 以使列表具有相同的长度。我想将它们组合成一个数据框并对它们进行回归。
如何将这些 NA 附加到列表中?
【问题讨论】:
-
this 对你有用吗?