htest クラスの結果の日本語化
をテンプレートにして作成
[
トップ
] [
新規
|
一覧
|
検索
|
最終更新
|
ヘルプ
]
開始行:
[[RjpWiki]]
ほとんど冗談でですが。。。
というか，どんな些細なことでも，日本語で書かれていれば。...
ということで，
以下のような関数を定義しておけば，base の統計的検定関数（...
すべての検定関数，全てのオプション（要するに出力される全...
大本の検定関数を日本語化するという方向性もあるのだけど。...
* 関数
print.htest <- function (x, digits = 4, quote = TRUE, pr...
{
conv1 <- function(x) ##### 検定の名称
{
# print(x)
if(any(grep("[0-9]+-sample test for equality of proport...
x <- "割合の一様性の検定（連続性の補正なし）"
}
else if(any(grep("[0-9]+-sample proportions test wi...
x <- "割合の検定（連続性の補正なし）"
}
else if (x == "Box-Pierce test") { # Box.test
x <- "Box-Pierce検定"
}
else if (x == "Box-Ljung test") {
x <- "Box-Ljung検定"
}
else if (x == "Phillips-Perron Unit Root Test") { #...
x <- "Phillips-Perronの単位根検定"
}
else if (x == "Ansari-Bradley test") { # ansari.test
x <- "アンサリ・ブラドレイ検定"
}
else if (x == "Bartlett test of homogeneity of var...
x <- "分散の一様性の検定（バートレット検定）"
}
else if (x == "Exact binomial test") { # binom.test
x <- "二項検定"
}
else if (x == "Pearson's Chi-squared test with Yate...
x <- "ピアソンのカイ二乗検定（イエーツの連続性補正...
}
else if (x == "Pearson's Chi-squared test") { # ch...
x <- "ピアソンのカイ二乗検定（連続性補正なし）"
}
else if (x == "Chi-squared test for given probabilit...
x <- "理論比が与えられたときのカイ二乗検定（適合度...
}
else if (x == "Pearson's product-moment correlation...
x <- "ピアソンの積率相関係数"
}
else if (x == "Spearman's rank correlation rho") { ...
x <- "スピアマンの順位相関係数"
}
else if (x == "Kendall's rank correlation tau") { #...
x <- "ケンドールの順位相関係数"
}
else if (x == "Fisher's Exact Test for Count Data")...
x <- "計数データにおけるフィッシャーの正確確率検定"
}
else if (x == "Fligner-Killeen test of homogeneity ...
x <- "分散の一様性の検定（Fligner-Killeen検定）"
}
else if (x == "Friedman rank sum test") { # fredman...
x <- "フリードマン検定"
}
else if (x == "Kruskal-Wallis rank sum test") { # k...
x <- "クラスカル・ウォリス検定"
}
else if (x == "Two-sample Kolmogorov-Smirnov test")...
x <- "二標本コルモゴロフ・スミルノフ検定"
}
else if (x == "One-sample Kolmogorov-Smirnov test")...
x <- "一標本コルモゴロフ・スミルノフ検定"
}
else if (x == "Mantel-Haenszel chi-squared test wit...
x <- "マンテル・ヘンツェルのカイ二乗検定（連続性の...
}
else if (x == "Exact conditional test of independen...
x <- "2 x 2 x k 分割表における条件付き独立性の正確...
}
else if (x == "Cochran-Mantel-Haenszel test") { # m...
x <- "コクラン・マンテル・ヘンツェル検定"
}
else if (x == "McNemar's Chi-squared test with cont...
x <- "マクネマー検定（連続性の補正）"
}
else if (x == "Mood two-sample test of scale") { # ...
x <- "尺度についての二標本ムード検定"
}
else if (x == "One-way analysis of means") { # one...
x <- "一元配置分散分析"
}
else if (x == "One-way analysis of means (not assu...
x <- "一元配置分散分析（等分散を仮定しない場合）"
}
else if (x == "Chi-squared Test for Trend in Propo...
x <- "割合の傾向についてのカイ二乗検定（傾向検定）"
}
else if (x == "Quade test") { # quade.test
x <- "Quade検定"
}
else if (x == "Shapiro-Wilk normality test") { # s...
x <- "シャピロ・ウィルクの正規性検定"
}
else if (x == "Welch Two Sample t-test") { # t.test
x <- "二標本t検定（Welchの方法）"
}
else if (x ==" Two Sample t-test") { # t.test
x <- "二標本t検定（分散が等しいと仮定できるとき）"
}
else if (x == "One Sample t-test") { # t.test
x <- "一標本t検定（母平均の検定）"
}
else if (x == "Paired t-test") { # t.test
x <- "対応のある場合のt検定"
}
else if (x == "F test to compare two variances") { ...
x <- "二群の等分散性の検定"
}
else if (x == "Wilcoxon signed rank test") { # wilc...
x <- "ウィルコクソンの符号付順位和検定"
}
else if (x == "Wilcoxon rank sum test") { # wilcox....
x <- "ウィルコクソンの順位和検定（マン・ホイットニ...
}
else if (x == "Wilcoxon rank sum test with continui...
x <- "ウィルコクソンの順位和検定（連続性の補正）"
}
else if (x == "Wilcoxon signed rank test with conti...
x <- "ウィルコクソンの符号付順位和検定（連続性の補...
}
return(x)
}
conv2 <- function(x)
{
if (length(x) == 2) {
if (x[1] == "num df") {
x[1] <- "第1自由度"
}
if (x[2] == "denom df") {
x[2] <- "第2自由度"
}
}
else if (x == "df") {
x <- "自由度"
}
else if (x == "Truncation lag parameter") {
x <- "切り捨てラグ・パラメータ"
}
else if (x == "number of trials") {
x <- "試行数"
}
return(x)
}
conv3 <- function(x) ##### 検定対象の名前
{
if (x =="difference in means") {
x <- "母平均の差"
}
else if (x == "mean") {
x <- "母平均"
}
else if (x == "mu") {
x <- "母平均"
}
else if (x == "correlation") {
x <- "母相関"
}
else if (x == "rho") {
x <- "母相関（ρ）"
}
else if (x == "tau") {
x <- "母相関（τ）"
}
else if (x == "probability of success") {
x <- "成功確率（母比率）"
}
else if (x == "ratio of scales") {
x <- "尺度の比"
}
else if (x == "common odds ratio") {
x <- "共通オッズ比"
}
else if (x == "odds ratio") {
x <- "オッズ比"
}
else if (x == "ratio of variances") {
x <- "分散比"
}
return(x)
}
conv4 <- function(x) ##### 推定値の名前
{
names(x) <- sub("mean of the differences", "差の平...
names(x) <- sub("mean of ", "平均値", names(x))
names(x) <- sub("prop ", "割合", names(x))
if (any(grep("mean in group [0-9]+", names(x)[1]))...
names(x) <- paste("グループ", 1:length(x), "の平均...
}
names(x) <- sub("cor", "相関係数", names(x))
names(x) <- sub("rho", "ρ", names(x))
names(x) <- sub("tau", "τ", names(x))
names(x) <- sub("probability of success", "成功確率...
names(x) <- sub("ratio of scales", "尺度の比", name...
names(x) <- sub("common odds ratio", "共通オッズ比"...
names(x) <- sub("odds ratio", "オッズ比", names(x))
names(x) <- sub("ratio of variances", "分散比", nam...
names(x) <- sub("difference in location", "位置母数...
}
conv5 <- function(x)
{
if (any(grep(" and ", x))) {
return(gsub(" and ", " と ", x))
}
else if (any(grep(" by ", x))) {
y <- unlist(strsplit(x, " "))
return(paste(y[1], "を", y[3], "で層別"))
}
else if (any(grep("null probability", x))) {
return(gsub("null probability", "帰無仮説における...
}
else if (any(grep("using scores", x))) {
return(gsub("using scores", "使用したスコア", x))
}
else {
return(x)
}
}
conv6 <- function(x) ##### 検定統計量の名前 STATISTIC
{
if (x == "X-squared") {
x <- "カイ二乗値"
}
else if (x == "t") {
x <- "t値"
}
else if (x == "Z") {
x <- "Z値"
}
else if (x == "Bartlett's K-squared") {
x <- "バートレットのK二乗値"
}
else if (x == "number of successes") {
x <- "成功数"
}
else if (x == "Friedman chi-squared") {
x <- "フリードマンのカイ二乗値"
}
else if (x == "Mantel-Haenszel X-squared") {
x <- "マンテル・ヘンツェルのカイ二乗値"
}
else if (x == "Cochran-Mantel-Haenszel M^2") {
x <- "コクラン・マンテル・ヘンツェルのM二乗値"
}
else if (x == "Fligner-Killeen:med chi-squared") {
x <- "Fligner-Killeenのカイ二乗値"
}
else if (x == "Kruskal-Wallis chi-squared") {
x <- "クラスカル・ウォリスのカイ二乗値"
}
else if (x == "McNemar's chi-squared") {
x <- "マクネマーのカイ二乗値"
}
return(x)
}
# cat("\n")
#@1
writeLines(strwrap(conv1(x$method)))
cat("\n")
#@5
cat("データ: ", conv5(x$data.name), "\n")
out <- character()
if (!is.null(x$statistic))
#@6
out <- c(out, paste(conv6(names(x$statistic)), "...
4))))
if (!is.null(x$parameter))
#@2
out <- c(out, paste(conv2(names(x$parameter)), "...
3))))
if (!is.null(x$p.value)) {
fp <- format.pval(x$p.value, digits = digits)
out <- c(out, paste("P値", if (substr(fp, 1, 1) ...
"<") fp else paste("=", fp)))
}
writeLines(strwrap(paste(out, collapse = ", ")))
if (!is.null(x$alternative)) {
cat("対立仮説: ")
if (!is.null(x$null.value)) {
if (length(x$null.value) == 1) {
alt.char <- switch(x$alternative, two.si...
less = "より小さい", greater = "より大...
#@3
cat(conv3(names(x$null.value)), "は，", ...
"\n", sep="")
}
else {
cat(x$alternative, "\nnull values:\n")
print(x$null.value, ...)
}
}
else {
alt.char <- switch(x$alternative, two.sided = "...
less = "小さい", greater = "大きい")
cat(alt.char, "\n")
}
}
if (!is.null(x$conf.int)) {
cat(format(100 * attr(x$conf.int, "conf.level"))...
format(c(x$conf.int[1], x$conf.int[2])), "\n")
}
if (!is.null(x$estimate)) {
cat("標本推定値: \n")
#@4
names(x$estimate) <- conv4(x$estimate)
print(x$estimate, ...)
}
cat("\n")
invisible(x)
}
* 応用例
上の関数を有効にした後。
t.test(rnorm(10), rnorm(14, mean=0.6))
とすれば，以下のように表示される
二標本 t 検定（Welchの方法）
データ: rnorm(10) と rnorm(14, mean = 0.6)
t = -1.2367, 自由度 = 19.745, P 値 = 0.2307
対立仮説: 母平均の差は，0ではない
95パーセント信頼区間: -1.2538991 0.3209893
標本推定値:
平均値x 平均値y
-0.01261903 0.45383583
また，
t.test(rnorm(100, mean=0.5))
では
一標本 t 検定（母平均の検定）
データ: rnorm(100, mean = 0.5)
t = 5.2999, 自由度 = 99, P 値 = 7.037e-07
対立仮説: 母平均は，0ではない
95パーセント信頼区間: 0.3507048 0.7704521
標本推定値:
平均値x
0.5605784
終了行:
[[RjpWiki]]
ほとんど冗談でですが。。。
というか，どんな些細なことでも，日本語で書かれていれば。...
ということで，
以下のような関数を定義しておけば，base の統計的検定関数（...
すべての検定関数，全てのオプション（要するに出力される全...
大本の検定関数を日本語化するという方向性もあるのだけど。...
* 関数
print.htest <- function (x, digits = 4, quote = TRUE, pr...
{
conv1 <- function(x) ##### 検定の名称
{
# print(x)
if(any(grep("[0-9]+-sample test for equality of proport...
x <- "割合の一様性の検定（連続性の補正なし）"
}
else if(any(grep("[0-9]+-sample proportions test wi...
x <- "割合の検定（連続性の補正なし）"
}
else if (x == "Box-Pierce test") { # Box.test
x <- "Box-Pierce検定"
}
else if (x == "Box-Ljung test") {
x <- "Box-Ljung検定"
}
else if (x == "Phillips-Perron Unit Root Test") { #...
x <- "Phillips-Perronの単位根検定"
}
else if (x == "Ansari-Bradley test") { # ansari.test
x <- "アンサリ・ブラドレイ検定"
}
else if (x == "Bartlett test of homogeneity of var...
x <- "分散の一様性の検定（バートレット検定）"
}
else if (x == "Exact binomial test") { # binom.test
x <- "二項検定"
}
else if (x == "Pearson's Chi-squared test with Yate...
x <- "ピアソンのカイ二乗検定（イエーツの連続性補正...
}
else if (x == "Pearson's Chi-squared test") { # ch...
x <- "ピアソンのカイ二乗検定（連続性補正なし）"
}
else if (x == "Chi-squared test for given probabilit...
x <- "理論比が与えられたときのカイ二乗検定（適合度...
}
else if (x == "Pearson's product-moment correlation...
x <- "ピアソンの積率相関係数"
}
else if (x == "Spearman's rank correlation rho") { ...
x <- "スピアマンの順位相関係数"
}
else if (x == "Kendall's rank correlation tau") { #...
x <- "ケンドールの順位相関係数"
}
else if (x == "Fisher's Exact Test for Count Data")...
x <- "計数データにおけるフィッシャーの正確確率検定"
}
else if (x == "Fligner-Killeen test of homogeneity ...
x <- "分散の一様性の検定（Fligner-Killeen検定）"
}
else if (x == "Friedman rank sum test") { # fredman...
x <- "フリードマン検定"
}
else if (x == "Kruskal-Wallis rank sum test") { # k...
x <- "クラスカル・ウォリス検定"
}
else if (x == "Two-sample Kolmogorov-Smirnov test")...
x <- "二標本コルモゴロフ・スミルノフ検定"
}
else if (x == "One-sample Kolmogorov-Smirnov test")...
x <- "一標本コルモゴロフ・スミルノフ検定"
}
else if (x == "Mantel-Haenszel chi-squared test wit...
x <- "マンテル・ヘンツェルのカイ二乗検定（連続性の...
}
else if (x == "Exact conditional test of independen...
x <- "2 x 2 x k 分割表における条件付き独立性の正確...
}
else if (x == "Cochran-Mantel-Haenszel test") { # m...
x <- "コクラン・マンテル・ヘンツェル検定"
}
else if (x == "McNemar's Chi-squared test with cont...
x <- "マクネマー検定（連続性の補正）"
}
else if (x == "Mood two-sample test of scale") { # ...
x <- "尺度についての二標本ムード検定"
}
else if (x == "One-way analysis of means") { # one...
x <- "一元配置分散分析"
}
else if (x == "One-way analysis of means (not assu...
x <- "一元配置分散分析（等分散を仮定しない場合）"
}
else if (x == "Chi-squared Test for Trend in Propo...
x <- "割合の傾向についてのカイ二乗検定（傾向検定）"
}
else if (x == "Quade test") { # quade.test
x <- "Quade検定"
}
else if (x == "Shapiro-Wilk normality test") { # s...
x <- "シャピロ・ウィルクの正規性検定"
}
else if (x == "Welch Two Sample t-test") { # t.test
x <- "二標本t検定（Welchの方法）"
}
else if (x ==" Two Sample t-test") { # t.test
x <- "二標本t検定（分散が等しいと仮定できるとき）"
}
else if (x == "One Sample t-test") { # t.test
x <- "一標本t検定（母平均の検定）"
}
else if (x == "Paired t-test") { # t.test
x <- "対応のある場合のt検定"
}
else if (x == "F test to compare two variances") { ...
x <- "二群の等分散性の検定"
}
else if (x == "Wilcoxon signed rank test") { # wilc...
x <- "ウィルコクソンの符号付順位和検定"
}
else if (x == "Wilcoxon rank sum test") { # wilcox....
x <- "ウィルコクソンの順位和検定（マン・ホイットニ...
}
else if (x == "Wilcoxon rank sum test with continui...
x <- "ウィルコクソンの順位和検定（連続性の補正）"
}
else if (x == "Wilcoxon signed rank test with conti...
x <- "ウィルコクソンの符号付順位和検定（連続性の補...
}
return(x)
}
conv2 <- function(x)
{
if (length(x) == 2) {
if (x[1] == "num df") {
x[1] <- "第1自由度"
}
if (x[2] == "denom df") {
x[2] <- "第2自由度"
}
}
else if (x == "df") {
x <- "自由度"
}
else if (x == "Truncation lag parameter") {
x <- "切り捨てラグ・パラメータ"
}
else if (x == "number of trials") {
x <- "試行数"
}
return(x)
}
conv3 <- function(x) ##### 検定対象の名前
{
if (x =="difference in means") {
x <- "母平均の差"
}
else if (x == "mean") {
x <- "母平均"
}
else if (x == "mu") {
x <- "母平均"
}
else if (x == "correlation") {
x <- "母相関"
}
else if (x == "rho") {
x <- "母相関（ρ）"
}
else if (x == "tau") {
x <- "母相関（τ）"
}
else if (x == "probability of success") {
x <- "成功確率（母比率）"
}
else if (x == "ratio of scales") {
x <- "尺度の比"
}
else if (x == "common odds ratio") {
x <- "共通オッズ比"
}
else if (x == "odds ratio") {
x <- "オッズ比"
}
else if (x == "ratio of variances") {
x <- "分散比"
}
return(x)
}
conv4 <- function(x) ##### 推定値の名前
{
names(x) <- sub("mean of the differences", "差の平...
names(x) <- sub("mean of ", "平均値", names(x))
names(x) <- sub("prop ", "割合", names(x))
if (any(grep("mean in group [0-9]+", names(x)[1]))...
names(x) <- paste("グループ", 1:length(x), "の平均...
}
names(x) <- sub("cor", "相関係数", names(x))
names(x) <- sub("rho", "ρ", names(x))
names(x) <- sub("tau", "τ", names(x))
names(x) <- sub("probability of success", "成功確率...
names(x) <- sub("ratio of scales", "尺度の比", name...
names(x) <- sub("common odds ratio", "共通オッズ比"...
names(x) <- sub("odds ratio", "オッズ比", names(x))
names(x) <- sub("ratio of variances", "分散比", nam...
names(x) <- sub("difference in location", "位置母数...
}
conv5 <- function(x)
{
if (any(grep(" and ", x))) {
return(gsub(" and ", " と ", x))
}
else if (any(grep(" by ", x))) {
y <- unlist(strsplit(x, " "))
return(paste(y[1], "を", y[3], "で層別"))
}
else if (any(grep("null probability", x))) {
return(gsub("null probability", "帰無仮説における...
}
else if (any(grep("using scores", x))) {
return(gsub("using scores", "使用したスコア", x))
}
else {
return(x)
}
}
conv6 <- function(x) ##### 検定統計量の名前 STATISTIC
{
if (x == "X-squared") {
x <- "カイ二乗値"
}
else if (x == "t") {
x <- "t値"
}
else if (x == "Z") {
x <- "Z値"
}
else if (x == "Bartlett's K-squared") {
x <- "バートレットのK二乗値"
}
else if (x == "number of successes") {
x <- "成功数"
}
else if (x == "Friedman chi-squared") {
x <- "フリードマンのカイ二乗値"
}
else if (x == "Mantel-Haenszel X-squared") {
x <- "マンテル・ヘンツェルのカイ二乗値"
}
else if (x == "Cochran-Mantel-Haenszel M^2") {
x <- "コクラン・マンテル・ヘンツェルのM二乗値"
}
else if (x == "Fligner-Killeen:med chi-squared") {
x <- "Fligner-Killeenのカイ二乗値"
}
else if (x == "Kruskal-Wallis chi-squared") {
x <- "クラスカル・ウォリスのカイ二乗値"
}
else if (x == "McNemar's chi-squared") {
x <- "マクネマーのカイ二乗値"
}
return(x)
}
# cat("\n")
#@1
writeLines(strwrap(conv1(x$method)))
cat("\n")
#@5
cat("データ: ", conv5(x$data.name), "\n")
out <- character()
if (!is.null(x$statistic))
#@6
out <- c(out, paste(conv6(names(x$statistic)), "...
4))))
if (!is.null(x$parameter))
#@2
out <- c(out, paste(conv2(names(x$parameter)), "...
3))))
if (!is.null(x$p.value)) {
fp <- format.pval(x$p.value, digits = digits)
out <- c(out, paste("P値", if (substr(fp, 1, 1) ...
"<") fp else paste("=", fp)))
}
writeLines(strwrap(paste(out, collapse = ", ")))
if (!is.null(x$alternative)) {
cat("対立仮説: ")
if (!is.null(x$null.value)) {
if (length(x$null.value) == 1) {
alt.char <- switch(x$alternative, two.si...
less = "より小さい", greater = "より大...
#@3
cat(conv3(names(x$null.value)), "は，", ...
"\n", sep="")
}
else {
cat(x$alternative, "\nnull values:\n")
print(x$null.value, ...)
}
}
else {
alt.char <- switch(x$alternative, two.sided = "...
less = "小さい", greater = "大きい")
cat(alt.char, "\n")
}
}
if (!is.null(x$conf.int)) {
cat(format(100 * attr(x$conf.int, "conf.level"))...
format(c(x$conf.int[1], x$conf.int[2])), "\n")
}
if (!is.null(x$estimate)) {
cat("標本推定値: \n")
#@4
names(x$estimate) <- conv4(x$estimate)
print(x$estimate, ...)
}
cat("\n")
invisible(x)
}
* 応用例
上の関数を有効にした後。
t.test(rnorm(10), rnorm(14, mean=0.6))
とすれば，以下のように表示される
二標本 t 検定（Welchの方法）
データ: rnorm(10) と rnorm(14, mean = 0.6)
t = -1.2367, 自由度 = 19.745, P 値 = 0.2307
対立仮説: 母平均の差は，0ではない
95パーセント信頼区間: -1.2538991 0.3209893
標本推定値:
平均値x 平均値y
-0.01261903 0.45383583
また，
t.test(rnorm(100, mean=0.5))
では
一標本 t 検定（母平均の検定）
データ: rnorm(100, mean = 0.5)
t = 5.2999, 自由度 = 99, P 値 = 7.037e-07
対立仮説: 母平均は，0ではない
95パーセント信頼区間: 0.3507048 0.7704521
標本推定値:
平均値x
0.5605784
ページ名:
![[RjpWiki]](image/../Rlogo.png "[RjpWiki]")
