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>