Displaying 1 result from an estimated 1 matches for "log_int".
Did you mean:
log_init
2009 Aug 25
1
Clogit or LRM?
...two variables:
model <- lrm(label ~ val1 + val2, data = traindata )
This gives me the following results:
Coef S.E. Wald Z P
Intercept 1.8065 0.05137 35.16 0
val1 0.8105 0.02567 31.57 0
val2 0.5218 0.04308 12.11 0
2) I then calculate a conditional logit:
testdata$log_int <- exp( model$coefficients[2] * model$val1 +
model$coefficients[3] * model$val2)
for(race in testdata$races){
testlogdata$c_prob[testdata$code== race] <-
testdata$log_int[testdata$race== race] /
sum(testdata$log_int[testlogdata$race == race])
}
---------------------------------------...