sorry, I had a typo there, it's datadist(b) for the
analysis of data frame "b".
--- Robert Balshaw <Robert.Balshaw@syreon.com> wrote:
> Not sure if this will help, but did you mean to use
> datadist(a) for
> the analysis of B?
>
> Rob
>
> > -----Original Message-----
> > From: r-help-bounces@stat.math.ethz.ch
> > [mailto:r-help-bounces@stat.math.ethz.ch]On Behalf
> Of array chip
> > Sent: Tuesday, December 14, 2004 4:28 PM
> > To: s-news@lists.biostat.wustl.edu;
> r-help@stat.math.ethz.ch
> > Subject: [R] Re: [S] using Hmisc and Design
> library
> >
> >
> > actually data frame "b" has 177 rows, the script
> ran
> > ok on the first 166 rows as a subset, but started
> to
> > break down if subset of the first 177 rows were
> used
> > as the dataset, or the first 166 rows plus 168th
> row,
> > ....
> >
> > the data in those rows are:
> > time status scores
> > 165 172 0 -123.3
> > 166 105 0 -138.4
> > 167 166 0 -128.8
> > 168 140 0 -114.2
> > 169 163 0 -117.0
> > 170 141 0 -115.8
> >
> > please advise!
> >
> >
> > --- array chip <arrayprofile@yahoo.com> wrote:
> >
> > > Hi, I encountered a weird problem when using the
> > > Design and Hmisc problem. I have 2 data frame
> called
> > > "a" and "b", both have 3 columns: "time",
> "status"
> > > and
> > > "scores", a sample of the data frame is like:
> > >
> > > data frame "a":
> > > time status scores
> > > 1 21 1 99.61
> > > 2 38 0 101.11
> > > 3 51 0 100.62
> > > 4 48 0 87.52
> > > 5 78 0 97.18
> > >
> > > data frame "b":
> > > time status scores
> > > 1 27 0 -126.7
> > > 2 24 0 -135.6
> > > 3 30 0 -139.5
> > > 4 49 0 -137.6
> > > 5 27 0 -136.9
> > >
> > > when I ran the following using data frame "a",
> > > everything is ok.
> > >
> > > > library(Hmisc,T);library(Design,T)
> > > > dd <- datadist(a)
> > > > options(datadist='dd')
> > > > fit <- cph(Surv(time,status) ~ scores,
> > > data=a,surv=T, x=T, y=T)
> > > > fit
> > > Cox Proportional Hazards Model
> > >
> > > cph(formula = Surv(time, status) ~ scores, data
> =
> > > data.combined.testing.set.scores, x = T, y = T,
> > > surv = T)
> > >
> > > Obs Events Model L.R. d.f. P Score Score P
> > > R2
> > >
> > > 85 31 7.47 1 0.0063 7.21 0.0072
> > > 0.088
> > >
> > >
> > > coef se(coef) z p
> > > scores 0.0766 0.0287 2.67 0.00752
> > >
> > >
> > > but when I ran the same script using data frame
> "b",
> > > even though no error message was produced, no
> fit
> > > object was generated:
> > >
> > > > rm(fit)
> > > > library(Hmisc,T);library(Design,T)
> > > > dd <- datadist(a)
> > > > options(datadist='dd')
> > > > fit <- cph(Surv(time,status) ~ scores,
> > > data=a,surv=T, x=T, y=T)
> > > > fit
> > > Problem: Object "fit" not found, while calling
> > > subroutine S_agsurv2
> > > Use traceback() to see the call stack
> > >
> > > can anyone tell me why? no missing value is
> present
> > > in
> > > either data frame.
> > >
> > > Thanks
> > >
> > >
> > >
> > >
> > >
> > > __________________________________
> > > Do you Yahoo!?
> > > The all-new My Yahoo! - What will yours do?
> > > http://my.yahoo.com
> > >
> >
>
----------------------------------------------------------
> > ----------
> > > This message was distributed by
> > > s-news@lists.biostat.wustl.edu. To
> > > unsubscribe send e-mail to
> > > s-news-request@lists.biostat.wustl.edu with
> >
> > >
> >
> >
> >
> >
> > __________________________________
> >
> > Dress up your holiday email, Hollywood style.
> Learn more.
> >
> > ______________________________________________
> > R-help@stat.math.ethz.ch mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
>
>
__________________________________
Do you Yahoo!?
Jazz up your holiday email with celebrity designs. Learn more.
http://celebrity.mail.yahoo.com
|