similar to: (coxph, se) Obtaining standard errors of coefficients from coxph to store

Dear R Help list, My question is regarding extracting the standard error or Z-score from a cph or coxph call. My Cox model is: - modz=cph(Surv(TSURV,STATUS)~RAGE+DAGE+REG_WTIME_M+CLD_ISCH+POLY_VS, data=kidneyT,method="breslow", x=T, y=T) I've used names(modz) but can't see anything that will let me extract the Z scores for each coefficient or the standard errors in the same

Dear R users, Below I have written 4 functions CIT1, CIT2a and CIT2b and CIT3 which recode a variable CLD_ISCH into 3 new variables(T1 T2 T3), I wish to use T1, T2 and T3 based on the values of tf1 and tf2. (NOTE:- T2a is used to create T2 in a long winded manner due to my lack of programming experience) I then attach T1 T2 and T3 to a dataset KidneyT that contains other variables i wish to use

Dear R users, I am looking for the source code of the implementation of ppr (Projection Pursuit Regression) in R. It will be great if citations of the source papers on which the implementation is based, are also provided. Thank you, Arvind Iyer, Grad student, Deptt. of Biomedical Engineering Viterbi School of Engineering University of Southern California, Los Angeles [[alternative HTML

Hi, Recently, I installed the 1.1.0 version of R (for Windows), since it includes an implementation of Projection Pursuit (I failed to write my own version of PP as a standalone C++ program). As far as I know, R offers two interfaces/sintax for the ppr() function. The first one requieres a regression formula and a data frame. The other requieres X, a matrix with the explanatory variables, and Y,

Hey, I'm curious about two of the design questions in LLVM after reading the language reference. If there's preexisting material explaining this, just point me at it (looked and couldn't find any). - Why unsigned types rather than signed-only like Java/JVM? If I understand correctly, their behavior is only distinguishable in overflow situations, and the availability of

Hi R users I try to use the lme but I can?t!!!!! My script is (some words in french, sorry!!): rm(list=ls(all=TRUE)) #Efface tous les objets en m?moire pour ?viter des erreurs library(MASS) #Chargement des Librairies library(car) library(Hmisc) library(tkWidgets) library(svDialogs) library(multtest) library(nlme) #Rep <- "C:/Documents and Settings/U3M/Bureau/steph/Scripts

Hi, I'm attempting to add a "Julian Day" column to a data frame. Here is my code and the resulting data frame: vic.data <- read.table("C:/VIC/data/vic.data.csv", header=F) names(vic.data) <- c("year", "month", "day", "precip", "evap", "runoff", "baseflow", "Tsup",

Hola a todos! perdón por molestarlos nuevamente con los contrastes, pero estoy trantando de entender que es lo que está haciendo R y de donde vienen los valores que informa pero  no lo logro. Creí haberlo entendido pero a la hora de usar mis datos los resultados no dan como deberían. Tengo dos variables explicativas que son factores con 3 niveles cada uno. Esta es la tabla de medias de la

event.nab.2 is 0/1 and I dichotomized va to get va.2 to see if I could get geeglm to work. glm has no problem with the data but geeglm chokes. Each subject (patient.id) has at most 2 observations and more than 3/4 of the subjects have 2 observations. I have even worse problems trying to use glmmPQL from MASS and worse still trying to use lmer from lme4. But I figured a marginal model would work.

Dear alltogether, two lme's, the data are available at: http://www.anicca-vijja.de/lg/hlm3_nachw.Rdata explanations of the data: nachw = post hox knowledge tests over 6 measure time points (= equally spaced) zeitn = time points (n = 6) subgr = small learning groups (n = 28) gru = 4 different groups = treatment factor levels: time (=zeitn) (n=6) within subject (n=4) within smallgroups

Hi everyone, I am fitting a simple linear regression model in R. My command is j=lm( Y ~ Sex + begsal + time + int) Call: lm(formula = Y ~ Sex + begsal + time + int) Coefficients: (Intercept) Sex begsal time int 191.916 -241.805 3.969 5.003 3.040 Now I wish to access the values of these coefficients for other purposes

Could someone identify a function that I might use to perform a components of variance analysis? In addition to the variance attributable to each factor, I would also like to obtain the SE of the variances. Thank you, John John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics Baltimore VA Medical Center GRECC and University of Maryland School of Medicine Claude Pepper OAIC University of

Dear All: Could you please tell me how I can draw figure formatted like figure 4.18 of MASS (4th) with the attached data set? Thanks Zhongming Yang --------------------------------- -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: sample.txt Url: https://stat.ethz.ch/pipermail/r-help/attachments/20050316/abfdb85e/sample.txt

List, I'm working on fitting a logistic model for a well known dataset (which is given below in case anyone wants to try to reproduce). I used both R and SAS to fit the model and have some differences in the parameter estimates. I'm wondering if R calculates the ML estimates differently. I'm making NO accusations as to which program is "right or wrong". That is not the