[程式] R Box Cox & Shapiro test

看板Statistics作者 (@@)時間15年前 (2010/05/13 14:41), 編輯推噓0(005)
留言5則, 2人參與, 最新討論串1/1
------------------------------------------------------------------------ [軟體程式類別]: R [程式問題]: R output [軟體熟悉度]: (1年以上) [問題敘述]: 我有一筆資料 想先把他們normal後再做分析 用了Box-Cox transfromation的code box.cox.powers 然後用跑出來的power 把原來的變數變成 variable^power 接著又用shapiro.test(variable^power)做檢驗 結果p-value超小> < 這樣不就不是normal? 想請問兩個問題 1.希望我沒有誤會box.cox.powers的output的意思 ex 如下面的output 新1=舊1^(-0.4754) 新2=舊2^(1.7826) 新~Normal ? 2.ex shapiro.test(新1) p-value=0.00005 這種情況下 我是否該繼續用box-cox的結果 假設它是常態分佈? 還是有其他的建議? 不好意思 麻煩大家了 謝謝!!!! [程式範例]: > summary(box.cox.powers(cbind(newdata[,4],newdata[,6],newdata[,7]))) Box-Cox Transformations to Multinormality Est.Power Std.Err. Wald(Power=0) Wald(Power=1) [1,] -0.4754 0.9043 -0.5257 -1.6316 [2,] 1.7826 0.6311 2.8248 1.2402 [3,] 2.8328 1.0234 2.7682 1.7910 L.R. test, all powers = 0: 12.7782 df = 3 p = 0.0051 L.R. test, all powers = 1: 6.4149 df = 3 p = 0.0931 ----------------------------------------------------------------------------- -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 98.255.1.168

05/13 14:44, , 1F
Box Cox 是有假設的,有時候資料就是沒辦法normalize
05/13 14:44, 1F

05/13 14:45, , 2F
實務的話,我建議可以畫個 QQ-plot 或是看一下 histogram
05/13 14:45, 2F

05/13 14:45, , 3F
QQ-plot 可以看一下是 long tail 還是 short tail
05/13 14:45, 3F

05/13 14:47, , 4F
histogram 可以看有沒有左偏或右偏
05/13 14:47, 4F

05/14 03:44, , 5F
謝謝!!
05/14 03:44, 5F
文章代碼(AID): #1BwvyUO2 (Statistics)