[問題] 卡方檢定時如何加權

看板R_Language作者 (Wen)時間8年前 (2017/09/20 21:38), 編輯推噓3(309)
留言12則, 2人參與, 最新討論串1/1
[問題類型]:程式諮詢 [軟體熟悉度]:新手 [問題敘述]: 簡單的做兩組類別資料的卡方檢定 但因為抽樣問題,需進行加權 有用sas算過,程式碼舉例如下 proc freq data=a; tables survey*sex/chisq; weight wt; run; 但是現在要把他轉成r我不知道要怎麼加權@@ chisq.test() 好像不能直接使用加權 麻煩大家幫我解答一下了QQ,謝謝~ -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 101.139.144.45 ※ 文章網址: https://www.ptt.cc/bbs/R_Language/M.1505914710.A.E34.html

09/20 21:45, , 1F
weight的用法可能不是如你想像的。是否可能將(假)資
09/20 21:45, 1F

09/20 21:45, , 2F
料及其SAS的結果附上?
09/20 21:45, 2F

09/20 22:14, , 3F
如果我沒記錯的話,此時sas裡的weight就只是頻率而已,
09/20 22:14, 3F

09/20 22:15, , 4F
而在R中並沒有「加權」這樣的想法。
09/20 22:15, 4F

09/20 22:16, , 5F
你看看 https://wp.me/pMm6L-zTx 是否已經可以解決問題
09/20 22:16, 5F

09/20 23:14, , 6F
Sas weight後結果p值=0.2195
09/20 23:14, 6F

09/20 23:14, , 7F
未weight,p值=0.2551 (r與sas結果相同)
09/20 23:14, 7F

09/20 23:14, , 8F
剛剛試著用weights這個package裡的wtd.chi.sq(survey
09/20 23:14, 8F

09/20 23:14, , 9F
,sex,weight=wt),但是與sas的結果不同><
09/20 23:14, 9F

09/21 00:16, , 10F
建議你確定到底要什麼檢驗方法。
09/21 00:16, 10F

09/21 07:45, , 11F
好~謝謝你!
09/21 07:45, 11F

09/21 19:30, , 12F
我仍然懷疑你對sas裡weight的理解。歡迎再討論。
09/21 19:30, 12F
文章代碼(AID): #1PmczMuq (R_Language)
文章代碼(AID): #1PmczMuq (R_Language)