Displaying 3 results from an estimated 3 matches for "repcnd".
2010 Jul 09
1
Appropriate tests for logistic regression with a continuous predictor variable and Bernoulli response variable
I have a data with binary response variable, repcnd (pregnant or not) and one predictor continuous variable, svl (body size) as shown below. I did Hosmer-Lemeshow test as a goodness of fit (as suggested by a kind “R-helper” previously). To test whether the predictor (svl, or body size) has significant effect on predicting whether or not a female sna...
2010 Jul 07
1
Different goodness of fit tests leads to contradictory conclusions
...ced p-values well above 0.05. Would anyone please teach me why deviance test produce p-values such a small value (0.001687886) that suggest lack of fit, while other tests suggest good fit? Did I do something wrong?
Thank you for your time and help!
Kiyoshi
mod.fit <- glm(formula = no.NA$repcnd ~ no.NA$svl, family = binomial(link = logit), data = no.NA, na.action = na.exclude, control = list(epsilon = 0.0001, maxit = 50, trace = F))
> # Option 1: Hosmer-Lemeshow test
> mod.fit <- glm(formula = no.NA$repcnd ~ no.NA$svl, family = binomial(link = logit), data = no.NA, na.actio...
2010 Jul 15
0
Histogram with two groups on the same graph (not on separate panels)
...9.2 9 8.5 9.5 8.8 ...
$ bm : num 142 128 148.3 192.3 70.5 ...
$ defensiveness: int NA NA NA NA NA 1 1 1 1 2 ...
$ ftr : int NA NA NA NA NA 5 10 60 18 30 ...
$ latency : num NA NA NA NA NA NA NA NA NA NA ...
$ dvc : int NA NA NA NA NA 25 25 38 28 28 ...
$ RepCnd : Factor w/ 4 levels "Male","Nonpregnant",..: 4 4 4 2 1 2 2 2 2
2 ...
$ repcnd : int 1 1 1 0 2 0 0 0 0 0 ...
$ repstatus : int 1 1 1 1 NA 0 0 0 0 0 ...
$ LS.w_egg : int NA NA NA NA NA NA NA NA NA NA ...
$ estimatedLS : int NA NA NA NA NA NA NA NA NA...