こにしき(言葉・日本社会・教育)

関西学院大学(2016.04~)の寺沢拓敬のブログです(専門:言語社会学)。

有意確率をアステリスクに変換するRスクリプト

hoshi(x)

hoshi <- function(x,ns){

X <- NULL
for(i in 1:length(x)){
if(x[i] < .001){X[i] <- c("***")
  }else{
if(x[i] >= .001 & x[i] < .01){X[i] <- c("**")
  }else{
if(x[i] >= .01  & x[i] < .05){X[i] <- c("*")
  }else{
if(x[i] >= .05  & x[i] < .1 ){X[i] <- c("†")
  }else{X[i] <- ns}
}
}
}
}
return(X)
}

Usage

hoshi( x, ns = ... )

 x : ベクトル形式のオブジェクト
 ns: Non-significantのとき(p値が0.1以上の時)に何を代入するか

Example

>  x <- runif(10)^2   # 乱数で有意確率っぽいものを生成
>  p1 <- hoshi(x, ns = NA) # pが0.1以上のときNAを出力
>  p2 <- hoshi(x, ns = "") # pが0.1以上のとき空白を出力
>  # 表で比較
>  data.frame( round(x,4), p1, p2)   
   round.x..4.   p1  p2
1       0.0265    *   *
2       0.3300 <NA>    
3       0.9126 <NA>    
4       0.0512   †  †
5       0.7036 <NA>    
6       0.0003  *** ***
7       0.2855 <NA>    
8       0.0071   **  **
9       0.7554 <NA>    
10      0.1077 <NA>