Question

我的数据如下：

ID  gamesAlone  gamesWithOthers  gamesRemotely  tvAlone  tvWithOthers  tvRemotely
1   1                                                    1
2                                1                       1
3                                1              1
4                                1              1
5                                1                       1
6                                1              1
7                                1              1
8               1                                        1
9   1                                                                   1

我希望代码可以执行以下两项操作：

首先，将其转换为一个整齐的列联表，如下所示：

        Alone   WithOthers   Remotely
games   2       1            6
tv      4       4            1

其次，使用卡方来观察这些活动（游戏与电视）的社交背景是否有所不同。

这是生成数据框的代码：

data<-data.frame(ID=c(1,2,3,4,5,6,7,8,9),
             gamesAlone=c(1,NA,NA,NA,NA,NA,NA,NA,1),
             gamesWithOthers=c(NA,NA,NA,NA,NA,NA,NA,1,NA),
             gamesRemotely=c(NA,1,1,1,1,1,1,NA,NA),
             tvAlone=c(NA,NA,1,1,NA,1,1,NA,NA),
             tvWithOthers=c(1,1,NA,NA,1,NA,NA,1,NA),
             tvRemotely=c(NA,NA,NA,NA,NA,NA,NA,NA,1))

Answer 1

忽略第一列 id （[-1]），然后取出每列的总和（colSums），同时删除 NA 值（{ {1}}），并将长度为6的结果向量放入一个有2行的矩阵中。如果需要，还可以相应地标记矩阵尺寸（na.rm=TRUE参数）：

dimnames

Answer 2

这将使您以您提供的形式进入列联表。建议：拨打您的数据框data1而不是data，以避免混淆。

library(dplyr)
library(tidyr)
data1_table <- data1 %>% 
  gather(key, value, -ID) %>% 
  mutate(activity = ifelse(grepl("^tv", key), substring(key, 1, 2), substring(key, 1, 5)), 
         context = ifelse(grepl("^tv", key), substring(key, 3), substring(key, 6))) %>% 
  group_by(activity, context) %>% 
  summarise(n = sum(value, na.rm = TRUE)) %>% 
  ungroup() %>% 
  spread(context, n)

# A tibble: 2 x 4
  activity Alone Remotely WithOthers
*    <chr> <dbl>    <dbl>      <dbl>
1    games     2        6          1
2       tv     4        1          4

对于卡方：它取决于您想要比较的内容，我会假设您的真实数据具有更高的计数。你可以把这一切都安排到这样的chisq.test，但我不认为这些内容非常丰富：

data1_table %>%
  select(2:4) %>%
  chisq.test()

对具有宽数据的数据帧进行卡方检验

2 个答案: