[程式] R: 關於probit/logit迴歸的問題已刪文
[軟體程式類別]:
R
[程式問題]:
probit/logit迴歸後的敘述統計問題
[軟體熟悉度]:
新手(不到1個月)
[問題敘述]:
我有一筆資料X (4個column y, x_1, x_2, x_3 有含header)
先用attach()把header當係數
然後用glm()做完了probit和logit的迴歸
y = b_0 +b_1 x_1 +b_2 x_2 +b_3 x_3 (x_2的資料也是binary)
，得出了迴歸的係數b_0, b_1, b_2, b_3 係數名稱是我用attach()得來的
我也用summary()得到x_1, x_2, x_3的敘述統計包含mean值
[1]我現在想要把x_1=mean, x_2=1 or 0, x_3=mean 代入上述已經估計出
b_0, b_1, b_2, b_3的迴歸式中，
[2]然後再把y的value代入probit/logit的cdf裡，
算出probability<-最終目標
請問我可以怎麼用R寫?
[1],[2]是需要分開寫嗎? 還是可以合在一起寫?
還是其實我需要自己寫programing?
或是有甚麼package可以解決我的問題?
----------------------------------以上問題都解決了-----------------------------------
新問題:
[3]如果今天加了一個binary的資料x_4，
也做了新的迴歸式後，我一樣想要按照[1], [2]的做法算probability
可是x_4我想要代入mean(不是代入0或1)
但是我在迴歸式那邊已經加了as.factor()了
請問需要甚麼fuction才能讓as.factor()接受非0或1的數值輸入?
或是有沒有甚麼其他作法可以讓我這樣做?
先謝謝大家了~
[程式範例]:
>X<-read.table("X.txt", header=T, sep="")
>attach(X)
>name(X)
"y", "x_1", "x_2", "x_3", "x_4"
>summary(X)
#這裡得出y, x_1, x_2, x_3, x_4共五個mean
>reg1<-glm(y ~ x_1 + as.factor(x_2) + x_3 +as.factor(x_4),
+family=binomial(link="logit"),data=X)
>summary(reg1)
#會估出截距, x_1, x_2, x_3, x_4 的係數b_0, b_1, b_2, b_3, b_4
>predict.glm(reg1, data.frame(x_1=mean(X$x_1), x_2=1, x_3=mean(X$x_3),
+x_4=?mean(X$x_4)),type="response")
麻煩大家幫忙提示了~~
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 59.121.148.160
→
02/25 07:58, , 1F
02/25 07:58, 1F
→
02/25 09:59, , 2F
02/25 09:59, 2F
→
02/25 11:13, , 3F
02/25 11:13, 3F
→
02/25 11:14, , 4F
02/25 11:14, 4F
※ 編輯: ieuser 來自: 140.112.4.191 (02/25 15:31)
→
02/25 15:32, , 5F
02/25 15:32, 5F
→
02/25 15:32, , 6F
02/25 15:32, 6F