ˊ2020/11/18 Wed.
德國已故總理科爾有句名言，謊言有三種：謊言、善意的謊言、統計數字。
GMBA 本學期一堂必修課Statistics and Data(統計及資料分析)，本來興致勃勃地覺得重回校園就是要擁抱以前懼怕的東西，總覺得十年前我還小，社會經驗不足，不懂得思考還情有可原；然而，年過而立又三的小姐姐我，至少在江湖打滾了七八年也對商業模式有點概念的前提下，應該有機會一點就通！💪
OK，每週二晚上我都提早報到，坐在帝王方位(黑板前第一排)，整理好思緒、集中元氣專心聽講！第一堂 no problem😃、第二堂Nice🤨、第三堂OK😑、第四堂課...Uh...Fine🙄️，第五堂課....shoot.........Fxxk......🤦🏿♀️
聽不懂課，讓我非常慌張！💦
σ、 xi、 i=1~N 、μ、R squre、ANOVA、T Test、P value、昨天還來了一個F，還有很多奇奇怪怪的符號。
問了隔壁從小在美國長大的ABC同學：老師在幹嘛？他用流利的英文+閃著信心的眼神回我這些就是數學基礎，在美國讀高中就都學過了，你以前沒有學過嗎？痾...情急之下我也用流利的英文回他：喔喔我在台灣讀的高中是數乙，這些沒碰過。(自我辯解時英文突然很遛，但我是亂回答的，不確定數乙有教？不該污衊台灣的科普教育😫）
然後他又接著問我數乙是什麼？嗯....沒時間跟他解釋啊啊啊，隨口跟他說了台灣教育跟你們不一樣，這牽扯到Cultural Variation(文化差異)與Educational Gap (教育斷層)，下次再跟你講🤦🏿♀️
課後沮喪地回到家，睡前驚鴻一撇看到老公書架買了一本書叫做「用數學的語言看世界」，書皮上又有一句驚人的標：學數學，「讓你多擁有一種靈魂」。
馬上點燃我心中三把火🔥🔥🔥失去理智大吼：「朱先生！請你下次不要亂買書！」「學數學並沒有讓我多擁有一個靈魂，簡直是讓我失魂落魄，身心受創」🚑🚑🚑
平常沒什麼機會唸他，這次真是逮到機會可以讓我可以開罵，有點紓壓～～～喜悅！🤗
想起年初在美國準備考雅思的日子，當時覺得好煩題目好多，但現在比一比，還真覺得比起念數學我寧願去做雅思考題咧！(眼神死)
好吧，寫完還是希望各位強者/學霸給我一些方法愛上數學/統計，我將不勝感激。
#不諳數字的女子廢文
#投胎需要數學腦
📸 安東尼張❤️
r anova 在 Re: [問題] type III anova - 看板R_Language - 批踢踢實業坊 的推薦與評價
先解決MS不正確的問題。
因為不平衡設計而採用 Type III SS 時，
必須將類別變數轉換成正交的虛擬變數。
不過，R 內建的虛擬變數轉換方式是 simple contrast
（在 R 中稱為 contr.treatment），
而 simple contrast 並不是正交的。
R 有內建其它常用的正交轉換方式，
包括 contr.sum 和 contr.helmert。
回到例子，正確的方式是。
contrasts(A) <- contr.sum # 更變虛擬變數
contrasts(B) <- contr.sum
m <- aov(Y~A*B)
library(car)
Anova(m, type=3) # Type III SS
即可得到正確的 MS。
至於如何設立隨機因子就先不談，還蠻複雜的。
※ 引述《DrRd (就這樣吧)》之銘言：
: [問題類型]:
: 程式諮詢(我想用R 做某件事情，但是我不知道要怎麼用R 寫出來)
: [軟體熟悉度]:
: 新手(沒寫過程式，R 是我的第一次)
: [問題敘述]:
: 要進行type III的變異數分析
: 共有A、B兩個因子，其中A是固定效果，B是隨機效果
: 網路上有看到可以用car裡面的Anova(,type=3)
: 但它沒有讓我指定何者為固定何者為隨機
: 算出來的結果也不正確
: 請問除了手算之外，有其他的套件可以算嗎？
: [程式範例]:
: aov.result = aov(Y~A*B)
: summary(aov.result)
: Df Sum Sq Mean Sq F value Pr(>F)
: A 2 222.6 111.30 5.532 0.00711 **
: B 1 0.1 0.14 0.007 0.93447
: A:B 2 76.4 38.18 1.898 0.16169
: Residuals 45 905.4 20.12
: ---
: Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
: 用aov得到的結果A效果的F值是除上Residuals得到的
: 但是因為B是隨機效果，所以A的F值應該是要除以A:B
: 目前我是都用手算來得到A效果的真正的F值
: 我在網路上看到說可以用car:::Anova
: 所以我改用
: >library(car)
: >Anova(aov.result,type=3)
: Anova Table (Type III tests)
: Response: Y
: Sum Sq Df F value Pr(>F)
: (Intercept) 544.44 1 27.0598 4.683e-06 ***
: A 22.03 2 0.5475 0.5822
: B 47.45 1 2.3584 0.1316
: A:B 76.37 2 1.8979 0.1617
: Residuals 905.40 45
: ---
: Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
: 改用這個之後Sum Sq都變了
: 也許是因為不等格的關係？
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 122.117.37.172
※ 編輯: andrew43 來自: 122.117.37.172 (05/06 04:13)
... <看更多>