similar to: using weights in lrm

I needed to compute a complicated cross tabulation to show weighted means and standard deviations and the only method I could get that worked uses a series of nested for next loops. I know that there must be a better way to do so, but could use some assistance pointing the way. Here is my working, but inefficient script: library(Hmisc) rm(list=ls()) load('NHTS.Rdata') day.wt <-

  Hi there,   I am beginner in R and I have some basic question. Suppose I run a common procedure such as a t test or cox model like below:   out<-coxph( Surv(tstart,tstop, death1) ~ x1+x1:log(tstop+1) , test1,method=c("breslow"))    Which yields the following result:   Call: coxph(formula = Surv(tstart, tstop, death1) ~ x1 + x1:log(tstop +     1), data = test1, method =

Dear R People: Here is a little section of a data frame: > zzz[1:10,] DATE GENDER AGE Co DEATH1 3945 2009-04-16 M 24 Botulinic 23 3851 2009-04-16 M 35 Constitutional 23 8495 2009-04-16 F 49 Constitutional 27 10967 2009-04-16 F 47 Constitutional 28 11164 2009-04-16 F 26 Constitutional 28 15220

Hi there, I have a rather large data set and perform the following cox model: test1 <- list(tstart,tstop,death1,chemo1,radio1,horm1) out1<-coxph( Surv(tstart,tstop, death1) ~ chemo1+chemo1:log(tstop+1)+horm1+horm1:log(tstop+1)+age1+grade1+grade1:log(tstop+1)+positive1+positive1:log(tstop+1)+size1+size1:log(tstop+1), test1) out1 Up to here everything works fine (with each covariate

I am trying to use apply and a list to supply names to a set of tables I want to generate. Below is an example that I hope mimics the larger original problem. EXAMPLE aa <- c( 2,2,1,1,2) bb <- c(5,6,6,7,4) aan <- c("yes", "no") bbn <- c("a", "b", "c", "d") mynames <- c("abby", "billy") mylist <-

Hello, Does anyone know of a way I can speed this up? Basically I'm attempting to get the data item on the same row as the report date for each report date available. In reality, I have over 11k of columns, not just A, B, C, D and I have to do that over 100 times. My solution is slow, but it works. The loop is slow because of merge. # create sample data z.dates =

Here is my solution: tapply(zzz$DEATH1,zzz$GENDER,sum, na.rm=TRUE) Dear R People: Here is a little section of a data frame: > zzz[1:10,] DATE GENDER AGE Co DEATH1 3945 2009-04-16 M 24 Botulinic 23 3851 2009-04-16 M 35 Constitutional 23 8495 2009-04-16 F 49 Constitutional 27 10967 2009-04-16 F 47 Constitutional

Dear all, I'm facing a problem in estimation of glm model with weibull distribution. I run this : eqn0<-formula(fdh~cup1+cup2+cup3+cup4+fin1+vd1+cm2+cm4+milieu+cpro1+cpro2+cpro3a+cpro3b+schef+log(y)) regWeib0<-vglm(eqn0,family=weibull,subset(br, fdh<1)) I have en estimation but there is a message saying that regularity conditions are violated : Message d'avis : In eval(expr,

Respected sir/madam can you please suggest what is an unexpected symbol in the below code for running a multinomial logistic regression model <- multinom(adoption ~ age + education + HH size + landholding + Farmincome + nonfarmincome + creditaccesibility + LHI, data=newdata) [[alternative HTML version deleted]]

Dear R listers, I have a silly problem. I am trying to load a dta (Stata) file in R. The dta is about 650 MB and contains the integrated World Values Survey/ European Value Study data-set. My problem is that I don't manage to load the file. After almost 1 hour I issued the following command: data <- read.dta("http://www.stata-press.com/data/kkd/data1.dta", convert.dates = TRUE,

I have been doing some research on echoes and echo cancellation. Here is a summary of what I've learned: 1) When solving echo problems it is extremely important to locate the source of the echo, and if possible, eliminate it at the source. The most common cause of echo is an unbalanced hybrid at the PSTN to TDM interface. This is where the 2-wire POTS line is divided into a 4-wire

Hi all, After working on my first rails app and having handed over some very sophisticated coding from a lot of you guys. My app does what I want it to do, but I''m no where near the end yet, and it seems that the KISS and DRY objectives already went down the drain. It could be lack of knowledge, only doing Ruby on (and) Rails for a month, but I''m stuck on the following: