1

我有一个包含一些缺失数据（随机缺失）的数据集（data1.csv），我正在从该数据集（d1）创建一个子集，以便在 d2 中仅保留完整的观察结果。我正在使用ftable， as.data.frame函数并创建一个表示每个组合百分比的列p 。

 d1 = read.csv("C:/Users/....../Data1.csv",header=T)

   y x1 x2 x3
1  1  2  1 NA
2  0  2 NA  0
3  1 NA  1  1
4  0  3  1  0
5  0  2  0 NA
6  0  1  0  1
7  1 NA NA  1
8  1  3  0  1
9  1  2  0  0
10 0  3  0 NA
11 1 NA  0  1
12 1  3 NA  1
13 0  3  1 NA
14 0 NA  1  0
15 1  1  0 NA
16 0  1 NA  1
17 1 NA  0 NA
18 0  3  1  1
19 0  1 NA  0
20 0 NA  0  1

 d2= d1[complete.cases(d1),]
 d2

   y x1 x2 x3
4  0  3  1  0
6  0  1  0  1
8  1  3  0  1
9  1  2  0  0
18 0  3  1  1


 d3= ftable(d2)
 d4=as.data.frame(d3)
 d4


 y x1 x2 x3 Freq   p
 1  0  1  0  0    0 0.0
 2  1  1  0  0    0 0.0
 3  0  2  0  0    0 0.0
 4  1  2  0  0    1 0.2
 5  0  3  0  0    0 0.0
 6  1  3  0  0    0 0.0
 7  0  1  1  0    0 0.0
 8  1  1  1  0    0 0.0
 9  0  2  1  0    0 0.0
 10 1  2  1  0    0 0.0
 11 0  3  1  0    1 0.2
 12 1  3  1  0    0 0.0
 13 0  1  0  1    1 0.2
 14 1  1  0  1    0 0.0
 15 0  2  0  1    0 0.0
 16 1  2  0  1    0 0.0
 17 0  3  0  1    0 0.0
 18 1  3  0  1    1 0.2
 19 0  1  1  1    0 0.0
 20 1  1  1  1    0 0.0
 21 0  2  1  1    0 0.0
 22 1  2  1  1    0 0.0
 23 0  3  1  1    1 0.2
 24 1  3  1  1    0 0.0

函数ftable， as.data.frame工作正常，但问题是我仍然在结果（d4）中看到所有这些观察结果，这很有意义。我以为我在完成时摆脱了这个。cases(d1)

所以需要帮助摆脱缺失值并只为这些完整的观察做频率表

   y x1 x2 x3
4  0  3  1  0
6  0  1  0  1
8  1  3  0  1
9  1  2  0  0
18 0  3  1  1
4

2 回答 2

2

查找 ?complete.cases()。它会做你想要的；您可以指定缺失值以及要为哪些列执行此操作。

于 2014-08-24T23:17:46.253 回答
0

这是您对输入数据进行采样

d1<-structure(list(y = c(1L, 0L, 1L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 
1L, 1L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L), x1 = c(2L, 2L, NA, 3L, 
2L, 1L, NA, 3L, 2L, 3L, NA, 3L, 3L, NA, 1L, 1L, NA, 3L, 1L, NA
), x2 = c(1L, NA, 1L, 1L, 0L, 0L, NA, 0L, 0L, 0L, 0L, NA, 1L, 
1L, 0L, NA, 0L, 1L, NA, 0L), x3 = c(NA, 0L, 1L, 0L, NA, 1L, 1L, 
1L, 0L, NA, 1L, 1L, NA, 0L, NA, 1L, NA, 1L, 0L, 1L)), .Names = c("y", 
"x1", "x2", "x3"), class = "data.frame", row.names = c("1", "2", 
"3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", 
"15", "16", "17", "18", "19", "20"))

这些命令

d2 <- d1[complete.cases(d1),]
d3 <- ftable(d2)
d4 <- as.data.frame(d3)
d4

返回下表。

   y x1 x2 x3 Freq
1  0  1  0  0    0
2  1  1  0  0    0
3  0  2  0  0    0
4  1  2  0  0    1
5  0  3  0  0    0
6  1  3  0  0    0
7  0  1  1  0    0
8  1  1  1  0    0
9  0  2  1  0    0
10 1  2  1  0    0
11 0  3  1  0    1
12 1  3  1  0    0
13 0  1  0  1    1
14 1  1  0  1    0
15 0  2  0  1    0
16 1  2  0  1    0
17 0  3  0  1    0
18 1  3  0  1    1
19 0  1  1  1    0
20 1  1  1  1    0
21 0  2  1  1    0
22 1  2  1  1    0
23 0  3  1  1    1
24 1  3  1  1    0

不完全是原始行。在您的完整案例中，这些是 、 、 和 的所有y可能x1组合x2。x3如果你只想要那些真正有观察的人，你可以这样做

 subset(d4, Freq>0)

但这基本上会让你回到d2

4  1  2  0  0    1
11 0  3  1  0    1
13 0  1  0  1    1
18 1  3  0  1    1
23 0  3  1  1    1

似乎不清楚您期望ftable做什么，或者您期望结果如何不同于您的专栏d2或您的p专栏来自哪里。代码完全按照我的预期工作，所以我不确定你认为结果会是什么。

于 2014-08-25T03:36:34.717 回答