[R] Use glm coefficients for other datasets
maxsilva
mmsilva3 at uc.cl
Wed Feb 9 04:16:06 CET 2011
Thank you for your kindness, but ive done what you've said and the problem
remains. What im doing is pretty straightforward,
>data
response pred1 pred2
1 1 0 1
2 0 0 0
3 1 0 0
4 1 1 1
5 1 0 1
6 0 1 1
7 1 1 0
8 1 0 1
9 0 1 1
10 0 0 1
>sdata <- data[sample(nrow(data), 5), ]
> sdata
response pred1 pred2
8 1 0 1
2 0 0 0
9 0 1 1
10 0 0 1
6 0 1 1
>model<-glm(data$response~data$pred1+datos$pred2,
family=binomial(link="logit"))
> summary(model)
###the model ran correctly
But when I ask for the predictions,
> pred<- predict(model, newdata=sdata,type="response")
Mensajes de aviso perdidos
'newdata' had 5 rows but variable(s) found have 10 rows
> length(pred)
[1] 10
And those 10 values are the fitted values corresponding for the model over
its origin dataset.
I really can't get what's the problem...
--
View this message in context: http://r.789695.n4.nabble.com/Use-glm-coefficients-for-other-datasets-tp3276626p3296561.html
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list