Displaying 1 result from an estimated 1 matches for "control_indices".
2011 Jul 11
2
problem with 'predict'
Hi,
I would like to tabulate the likelihood for an affection. For this, I retrieve indices of affected people and controls for my data set and proceed as follows:
flags <- c(rep(1, length(patient_indices)), rep(0, length(control_indices)))
# dataset is a data.frame and param the parameter to be analysed:
data1 <- dataset[,param][c(patient_indices, control_indices)]
fit1 <- glm(flags ~ data1, family = binomial)
new.data <- seq(0, 300, 10)
new.p <- predict(fit1, data.frame(newdata = new.data), type = "respons...