Stats

Popular Posts

Followers

Mathematica 教學 統計檢定表格

戴忠淵 於 2014年1月27日星期一 下午6:17 發表


表格多的時候,這倒是個好方法~~

myTTest1[data_]:=Flatten@{ToString[NumberForm[Mean@data[[All,2]],{4,3}]]<>
"\[PlusMinus]"<>ToString[NumberForm[StandardDeviation@data[[All,2]],{4,3}]],
ToString[NumberForm[Mean@data[[All,3]],{4,3}]]<>"\[PlusMinus]"<>
ToString[NumberForm[StandardDeviation@data[[All,3]],{4,3}]],
If[Abs[#]>0.00002,NumberForm[#,{4,3}],"0.000"]&/@TTest[{data[[All,2]],
data[[All,3]]},0,{"TestStatistic","PValue"},VerifyTestAssumptions->"EqualVariance"]};

myTTest2[data_]:=Flatten@{ToString[NumberForm[Mean@data[[1,All,2]],{4,3}]]
<>"\[PlusMinus]"<>ToString[NumberForm[StandardDeviation@data[[1,All,2]],{4,3}]],
ToString[NumberForm[Mean@data[[2,All,2]],{4,3}]]<>"\[PlusMinus]"
<>ToString[NumberForm[StandardDeviation@data[[2,All,2]],{4,3}]],
If[Abs[#]>0.00002,NumberForm[#,{4,3}],"0.000"]&/@TTest[{data[[1,All,2]],
data[[2,All,2]]},0,{"TestStatistic","PValue"},VerifyTestAssumptions->"EqualVariance"]}


mydata1[data_,cond1_,cond2_]:=Block[{height,weight,bmi,x1,x2,x3,x4},
height=Select[data,#[[3]]==cond1&&#[[1]]==cond2&][[All,{1,5,8}]];
weight=Select[data,#[[3]]==cond1&&#[[1]]==cond2&][[All,{1,6,9}]];
bmi=Select[data,#[[3]]==cond1&&#[[1]]==cond2&][[All,{1,7,10}]];
x1=Heigt=Select[data,#[[3]]==cond1&&#[[1]]==cond2&][[All,{1,11,15}]];
x2=Heigt=Select[data,#[[3]]==cond1&&#[[1]]==cond2&][[All,{1,12,16}]];
x3=Heigt=Select[data,#[[3]]==cond1&&#[[1]]==cond2&][[All,{1,13,17}]];
x4=Heigt=Select[data,#[[3]]==cond1&&#[[1]]==cond2&][[All,{1,14,18}]];
{height,weight,bmi,x1,x2,x3,x4}
]

mydata2[data_,cond2_]:=Block[{height,weight,bmi,x1,x2,x3,x4},
height=Select[data,#[[1]]==cond2&][[All,{3,8}]];
weight=Select[data,#[[1]]==cond2&][[All,{3,9}]];
bmi=Select[data,#[[1]]==cond2&][[All,{3,10}]];
x1=Select[data,#[[1]]==cond2&][[All,{3,15}]];
x2=Select[data,#[[1]]==cond2&][[All,{3,16}]];
x3=Select[data,#[[1]]==cond2&][[All,{3,17}]];
x4=Select[data,#[[1]]==cond2&][[All,{3,18}]];
SplitBy[#,(#[[1]])&]&/@{height,weight,bmi,x1,x2,x3,x4}
]

mylabel1={{"身高","體重","BMI","姿體前彎","立定跳遠","仰臥起坐","心肺適能"},
{"前測","前測","T 值","P-Value"}};
mylabel2={{"身高","體重","BMI","姿體前彎","立定跳遠","仰臥起坐","心肺適能"},
{"男生","女生","T 值","P-Value"}};

TableForm[myTTest1/@mydata1[data,"男","實驗組"],
TableHeadings->mylabel1]

TableForm[myTTest1/@mydata1[data,"女","實驗組"],
TableHeadings->mylabel1]

TableForm[myTTest1/@mydata1[data,"男","對照組"],
TableHeadings->mylabel1]

TableForm[myTTest1/@mydata1[data,"女","對照組"],
TableHeadings->mylabel1]

TableForm[myTTest2/@mydata2[data,"實驗組"],TableHeadings->mylabel2]

TableForm[myTTest2/@mydata2[data,"對照組"],TableHeadings->mylabel2]


Tags: , ,

讀者回應 ( 0 意見 )

發佈留言

Please leave your name and tell me what you thought about this site. Comments, suggestions and views are welcomed.

如果這篇文章對你有幫助,那請留個訊息給我~