search for: logt

I have data on a single variable LOGT. It has about 300,000 observations. I am trying to make a Histogram out of this data set. Following is my effort. Could anyone help me to solve this error. > hist(x) Error in hist.default(x) : `x' must be numeric > class(x) [1] "data.frame" > is.object(x) [1] TRUE > is.v...

...ls each, so about 1000 columns. The 300 factors all pass a preliminary test of association with the outcome. Solving this with cph from Design takes about 3 hours. I have created the fit with x=T, y=T to save the model data. I want to validate the PH assumption by calling survplot(fit, gender=NA, logt=TRUE, loglog=TRUE) for many of the factors (here gender is one column name). Just creating this one plot takes 40m. I'd be happy to sample from the fitted model to create these tests, or figure out another way to check assumptions in the model. Has anyone done something similar, or have othe...

...strGender, data = recidivismv) summary (cox.V) # produces a log-log plot mvfit <- survfit (Surv(recidivismv$intDaysUntilFVPO, recidivismv$Event_v) ~ MHC, data = recidivismv) neg.ll <- function(mvfit) -log(-log(mvfit)) library (Design) survplot(mvfit, fun = neg.ll, conf = "none", logt=TRUE) # hazard values mhc2 <- coxph (Surv(recidivismv$intDaysUntilFVPO, recidivismv$Event_v) ~ MHC, data = recidivismv) basehaz(mhc2)

I expect that there's something glaringly obvious that I'm overlooking, as I'm justr getting back involved in using R after a several-month hiatus (from R). So I welcome clues. When I invoke plot(), merely specifying a data.frame with 2 columns, specify the plot type ("type") of "p" ("points"), and that I want the point to be green ('col =

When I use zelig to compute parameter estimates with "model = logit.survey", I receive the following error: Nicholas Carnes. 2007. "logt.surveyWarning message: In eval(expr, envir, enclose) : non-integer #successes in a binomial glm! I believe this is because the model is not using "quasibinomial" (see R-News 2003, Analyzing Survey Data in R, by Thomas Lumley). Is there a way to change the model to use "quasibinomia...

Dear R users, I know one way to see the code for a hidden function, say function_x, is using default.function_x (e.g. summary.default). But how can I see the code for imported packages that have no namespace (in this case Design)? Many Thanks Eleni