2014-04-23 39 views
0

我想在權重做R邏輯迴歸,但我真的不知道它是如何工作的。當我應用權重時,發生了一些奇怪的事情,所有值都出現在1處,但我不明白爲什麼? (還有我該如何適應通過點的線?) 我嘗試計算觀測值與預測值的相關係數。另外,我的目標是在y軸上的「fra」(從0-1範圍),x軸上的溫度值,plot中的fra值和一條迴歸線(如下例子:http://imgur.com/FWevi36) 謝謝!加權logistic迴歸在R - 初學者水平

到目前爲止我什麼(由編碼):

#Dataframe 
temp=c(1,1,2,2,3,4,4,5,5,6,6,7,7,8,8) 
fra=c(0.0,0.0,0.0,0.0,0.0,0.0,0.5,0.2,0.2,0.3,0.1,0.3,0.4,0.0,0.5) 
bin=c(0,0,0,0,0,0,1,1,1,1,1,1,1,0,1) 
test1 <- as.data.frame(cbind(temp,bin,fra)) 

#Overview      
plot(test1$temp, test1$bin) 
plot(test1$fra) 
boxplot(test1$temp ~ test1$bin, horizontal=TRUE) 
#Logistic Regression without weight 
glmt1 <- glm(test1$bin~test1$temp, family=binomial) 
coefficients(summary(glmt1)) 
fit1 <- fitted(glmt1) 
#plot 
plot(test1$temp, fit1, ylim=range(0,1)) 
#line should go to points..??? 
lines(test1$bin, glmt1$fitted, type="l", col="red") 

#with weighted 
glmt2 <- glm(test1$bin~test1$temp, family=binomial, weights=test1$fra) 
coefficients(summary(glmt2)) 
fit2 <- fitted(glmt2) 
plot(test1$temp, fit2, ylim=range(0,1)) 
+0

的值如何,您的模型總是預測1,因爲它涉及估計技術,所以應該在交叉驗證(stats.stackexchange.com)中提出這個問題。 – ilir

回答

1

你只給予了積極的權重情況下bin == 1。這將消除響應變量中的所有變化(這次您在LHS中有fit1$bin)。這意味着無論temp1$temp