まず,矩形でないデータはdata.frameにはできません。なのでNAを詰めて使うならば,
- list型として利用,あるいはそれぞれをベクトルで切り離して利用
- long data(縦型データ)へ変換して利用
のいずれかとなります。以降は後者を念頭においています:
library(tidyr)
# サンプル用のデータ生成
df <- data.frame(
A1 = sample(c(4:6, NA), 10, prob = c(4, 2, 2, 2), replace = TRUE),
A2 = rep(NA, 10),
A3 = sample(c(4:6, NA), 10, prob = c(2, 4, 2, 2), replace = TRUE),
A4 = sample(c(4:6, NA), 10, prob = c(2, 2, 4, 2), replace = TRUE)
)
# サンプル用データの構造を確認
str(df)
#> 'data.frame': 10 obs. of 4 variables:
#> $ A1: int 6 5 NA 4 6 5 5 4 5 NA
#> $ A2: logi NA NA NA NA NA NA ...
#> $ A3: int 5 6 5 5 6 5 NA NA 5 4
#> $ A4: int 6 6 NA 4 6 5 6 6 4 6
# long data(縦型データ)へ変換
# tidyr::gather関数を使ってます
df_long <- gather(df, key = fac, value = value, factor_key = TRUE)
str(df_long)
#> 'data.frame': 40 obs. of 2 variables:
#> $ fac : Factor w/ 4 levels "A1","A2","A3",..: 1 1 1 1 1 1 1 1 1 1 ...
#> $ value: int 6 5 NA 4 6 5 5 4 5 NA ...
# valueの列でNAがあるものを除去
# tidyr::drop_na関数を使ってます
df_long_rmNA <- drop_na(df_long, value)
str(df_long_rmNA)
#> 'data.frame': 25 obs. of 2 variables:
#> $ fac : Factor w/ 4 levels "A1","A2","A3",..: 1 1 1 1 1 1 1 1 3 3 ...
#> $ value: int 6 5 4 6 5 5 4 5 5 6 ...
# NAを別の値へ置換
# tidyr::replace_na関数を使ってます
df_long_repNA <- replace_na(df_long, list(value = 0))
str(df_long_repNA)
#> 'data.frame': 40 obs. of 2 variables:
#> $ fac : Factor w/ 4 levels "A1","A2","A3",..: 1 1 1 1 1 1 1 1 1 1 ...
#> $ value: num 6 5 0 4 6 5 5 4 5 0 ...
おそらくはこれで希望されたデータの加工ができるかと思います。なお,この縦型データのままでwilcox.exact
も実行可能です:
library(exactRankTests)
#> Package 'exactRankTests' is no longer under development.
#> Please consider using package 'coin' instead.
# wilcox.exact関数を実行(fomula形式で指定)
# 引数subsetで,用いるレコード(行)を絞る
# subsetの使い方についてはsubset関数を参照
wilcox.exact(value ~ fac, df_long,
subset = df_long$fac %in% c("A1", "A3"))
#>
#> Exact Wilcoxon rank sum test
#>
#> data: value by fac
#> W = 29, p-value = 0.9776
#> alternative hypothesis: true mu is not equal to 0
wilcox.exact(value ~ fac, df_long_rmNA,
subset = df_long_rmNA$fac %in% c("A1", "A3"))
#>
#> Exact Wilcoxon rank sum test
#>
#> data: value by fac
#> W = 29, p-value = 0.9776
#> alternative hypothesis: true mu is not equal to 0
wilcox.exact(value ~ fac, df_long_repNA,
subset = df_long_repNA$fac %in% c("A1", "A3"))
#>
#> Exact Wilcoxon rank sum test
#>
#> data: value by fac
#> W = 47, p-value = 0.9131
#> alternative hypothesis: true mu is not equal to 0
ただ,exactRankTestsパッケージはまもなく開発完了するのでcoinパッケージを使うように,とメッセージがでました。一応coinパッケージのwilcox_test関数を利用したバージョンも以下に示します:
library(coin)
#> Loading required package: survival
#>
#> Attaching package: 'coin'
#> The following objects are masked from 'package:exactRankTests':
#>
#> dperm, pperm, qperm, rperm
wilcox_test(value ~ fac, df_long,
subset = df_long$fac %in% c("A1", "A3"),
distribution = "exact")
#>
#> Exact Wilcoxon-Mann-Whitney Test
#>
#> data: value by fac (A1, A3)
#> Z = -0.35161, p-value = 0.9776
#> alternative hypothesis: true mu is not equal to 0
出力もつけたほうがわかりやすいかと思ったのですが,長文となり申し訳ありません。