作者tutonto (@@)
看板Statistics
标题[程式] R Box Cox & Shapiro test
时间Thu May 13 14:41:31 2010
------------------------------------------------------------------------
[软体程式类别]: 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
1F：→ Wush978:Box Cox 是有假设的，有时候资料就是没办法normalize 05/13 14:44
2F：→ Wush978:实务的话，我建议可以画个 QQ-plot 或是看一下 histogram 05/13 14:45
3F：→ Wush978:QQ-plot 可以看一下是 long tail 还是 short tail 05/13 14:45
4F：→ Wush978:histogram 可以看有没有左偏或右偏 05/13 14:47
5F：→ tutonto:谢谢!! 05/14 03:44