similar to: [OT] Modeling strategies

The Design library has been fully ported to R except for Cox proportional hazards regression modeling (using Therneau's survival package) which will be available in about two weeks. It will take much longer to make all the example code executable, is it currently contains many examples for which data are not provided. Thanks to Xiao Gang Fan <xiao.gang.fan1 at libertysurf.fr> who

The Design library has been fully ported to R except for Cox proportional hazards regression modeling (using Therneau's survival package) which will be available in about two weeks. It will take much longer to make all the example code executable, is it currently contains many examples for which data are not provided. Thanks to Xiao Gang Fan <xiao.gang.fan1 at libertysurf.fr> who

Hi everybody. I'm studying the Harrell's book on regression modeling strategies, and I'd like to replicate the examples using R (instead of S+). So how's going the porting of the hmisc and design libraries? Are there available, though in beta version? TIA, Stefano Calza -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Has anyone tried to download Hmisc and used ldBands function for calculating Lan-Demets group sequential boundaries? The write-up in F.Harrell's website indicates that, besides downloading the package Hmisc, one needs to copy the progra ld98 from the University of Wisconsin website. As suggested, I did this but received another error message regarding the search path. I think I have fixed

Dear listserve members, especially Prof. Harrell: I am trying to create a factor variable that has fewer levels than the original. I have a factor: >rosa$risk1 [1] 2 2 5 1 ... [1799] 3 3 1 3 1 6 3 3 1 5 3 5 3 3 3 0 3 3 3 1 1 3 Levels: 0 1 2 3 4 5 6 8 But when I do this: rosa2 <- upData(rosa,

Not exactly an R question but.......... Does anyone have a dataset with a nice story behind it that would be suitable for regression analysis but is missing data in one or two covariates ? I am playing around with different missing data routines and need a nice dataset to play with. Thanks, Steve __________________________________________________ Yahoo! News - Today's headlines

Hi all, I'm currently working with a dataset that has quite a few missing values and after some investigation I figured that multiple imputation is probably the best solution to handle the missing data in my case. I found several references to functions in S-Plus that perform multiple imputation (NORM, CAT, MIX, PAN). Does R have corresponding functions? I searched the archives but was not

Hello, No luck in loading Frank Harrell's packages -- did anyone encountered the same problem? R : Copyright 2002, The R Development Core Team Version 1.4.1 (2002-01-30) > library(design, T) Design library by Frank E Harrell Jr, Version of Wed Apr 17 17:07:30 EDT 2002 Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library

Hi.. I'm using gam() to fit a spline model for a data set that has two predictor variables (say A and B). The results indicate that the higher order interaction terms are significant. The R^2 jumps from .5 to .9 when I change the maximum order for the interaction from 10 to 15 (i.e. (AB)^10 to (AB)^15). Is there any way of finding out which of the terms in the model are really

Hello R-user I want to compute a multiple regression but I would to include a check for collinearity of the variables. Therefore I would like to use a ridge regression. I tried lm.ridge() but I don't know yet how to get p-values (single Pr() and p of the whole model) out of this model. Can anybody tell me how to get a similar output like the summary(lm(...)) output? Or if there is

In a .First.lib I want to issue two require()s to insure that two other packages are loaded. But I want the package being loaded by .First.lib using library.dynam("mypackage",pkb,lib) to be higher in the search order than the two required packages, because I want to have a couple of functions from the two required packages overridden. What is the best way to do that? Thanks in advance

I am unclear on whether to run R CMD check pgkname as user or as root on Linux. When running as user, after all the latex, html, and man files are created, I get the following error message: Rdconv(): Couldn't open '': Permission denied Has anyone dealt with that message? When I run R CMD check a second time, all latex, html, etc. are recreated which takes quite a while. Is there

Hello I am sure this is elementary, but I can't figure it out..... I am using SAS v8.2 and R 1.4.1 on a Windows platform. I have a large (6 megabyte) file in .sd2 format. I want to import it into R. If necessary, I can create a smaller file, as I don't need all the variables that are in this file I downloaded the package foreign. I then tried to use PROC CPORT in SAS to create a

Dear lister Few months ago I posted a question about the drawing of a nomogram through R routines. Any news? thanks _________________________________________________________________ http://join.msn.com/?page=features/junkmail -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

I am using GLM to calculate logit models based on cross-sectional data. I am now down to the hard work of making the results intelligible to very average readers. Is there any way to calculate a psuedo analoque to the R^2 in standard linear regression for use as a purely descriptive statistic of goodness of fit? Most of the readers of my report will be vaguely familiar and more comfortable with

Dear all, I'm trying to solve the following problem, and hoping to get some advise here from the group I have a dataframe in which the same sample was measured more than once on the same day. I would reorganize the dataframe to get a single value (mean /median) for one day Patient Day Sample Test A 1 A 23 A 1 A 36 A 5 B 44 A 5 B 23 B 2 C 10 B 2 C 5 mean Patient Day Sample Test.mean A

Hi. I've a dataframe with about 46000 women with about 500 cases (cancers). I want to define a case-cohort sampling scheme, matching by age and hospital centre. Is there anyone who has already written a code for that? It should be something similar to the stcacoh macro in Stata. TIA, Stefano -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list

I'm interested in finding a package in R to do multiple imputations, similar to MICE in S-Plus. Does one exist? Please email directly to me. Thanks. Linda Andrews Social & Scientific Systems (301) 628-3234 -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

i make me thoughts about a "advanced tabulation" package similar to commercial software products like Quantum or Wincross. Before i'm beginning to fight with coding - is in the mailing-List anybody doing something similar in the past and have a good starting point and/or suggestions for me ? My purpose ist to define for a dataset headers (i.e. sex,age-groupes..) which should

Dear list members, I have a statistical question, that doesn't belong to this list, and I apologise for that in advance but I would appreciate your help very much. Is there some convention for selecting the a level for significance testing in scientific (e.g. chemical processes) studies? Most people use the 0.05 level but I could not find a reference to justify this. Why not 0.01 or 0.1?