similar to: catch error/warning message in a loop

Hoy, I have the following problem. I have to import the dataset (survey) from SPSS into R-programm. By using the "read.spss" the all of variables have been readed as factor. I need the variables "income" and "age" as numeric for the regreesion analysis. Could you plese help me. Thanks a lot, best regards, Silika

Dear All, I am learning R so it's a very simple problem but I do not understand while I am not able to generate a graph from two vectors. when I type this code, it generates a very nice graph. pdf("mygraph.pdf") > attach(mtcars) > plot(wt,mpg) > abline(lm(mpg~wt)) > title("Regreesion of mpg") > detach(mtcars) > dev.off() But I am trying to create a

  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 =

I hope one and all will allow a stats question: When running a cox proportional hazards model ,there are two ways to deal with age, including age as a covariate, or to include age as part of the follow-up time, viz, Age as a covariate: tetest1 <- list(time= c(4, 3,1,1,2,2,3), status=c(1,NA,1,0,1,1,0), age= c(0, 2,1,1,1,0,0),

Hi all, I'd like fit this function: G(t)=A*1/(1+t/B)*1/sqrt(1+t/C) where A, B and C are fitting constants that one like search. I have got a fcs<-(t,G(t)) value table in a data frame. Any one can help me? Tahnks Horacio.

Hello, I would like to make my R script more interactive. Well, I want the script to work like this: I run the script from R Console and it offers me some options (press 1 for something, or press 2 for something else). After pressing 1 or 2, do corresponding procedure. I am not able to find out how to do this. I guess the functions menu() or switch() are the keys, but I didn't manage to use

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

Hello Rene, Maybe have you already upgrade you openssl/openssh packages if you using debian (or *buntu). http://article.gmane.org/gmane.linux.debian.security.announce/1614 We have found some weak keys : - nightlycomp@rio (sure) - comp@rio (maybe) Are you up to date ? Can Fabrice sent you new ssh public ? Tahnks Best Regards -- Gabriel CORRE gac@4js.com - Four J's Development Tools -

Hi all, I'm wanting to be able to find and store the z-score of coxph below: - modz=coxph(Surv(TSURV,STATUS)~RAGE+DAGE+REG_WTIME_M+CLD_ISCH+POLY_VS, data=kidneyT,method="breslow") I know summary(modz) will give me this, but how do i extract the standard error or z-score values in a similar way to obtaining the coefficients by coef(modz) ? I think it must be something to do with

Hi All, I have a very simple problem. I have several files in a same directory. I would like to send for an object only the files that finish in ".sens.". I execute the command below, files <- dir(pattern="*.sens") but it includes all of the files that have "sens", independent of they be in the end or in the middle of the name of the file. How could I solve

Hey guys. Thanks for the help on the Pattern matching, I got that working pretty nicely. the next problem I have is that I'm using an include file, but its not really working... In my extensions.conf: [incoming] exten => _NXXNXXXXXX,1,SetCallerID("Unknown Called Number") #include "numbers.conf" exten => _NXXNXXXXXX,3,Macro(Number,1000,${EXTEN}) [macro-Brand]

Hi, I used "princomp" for PCA analysis based on correlation matrix (cor=T). I would like to reproduce the scores for each observation by first standardizing the data matrix (mean=0, std err=1), and then multiplied by the loadings of each variable for each principle components. I get very close numbers, but not exactly the same. anything I missed here? tahnks

Hi all, I am trying to construct a pairwise coocurrence matrix for certain terms appearing in a number of documents. For example I have the following table with binary values showing the presence or absence of a certain term in a document: term1 term2 term3 term4 term5 doc1 1 1 0 0 1 doc2 1 1 0 1 1 doc3 1 0 0 0 1 And I want to have a matrix with the number of the pairwise coocurrencies. So,

Dear All: I have a question on using coxph for multiple genes: I have written code to loop through all 22283 genes in the Hgu-133A and apply coxph on survival data. However, I don't know how to work with the result for each gene: survtest<-coxph(Surv(pcc.primary.stg.3.cox[,'fup_interval'],pcc.primary.stg.

Hello, I am hoping for some advice regarding warning/error messages I received when running a Cox regression # message 1 - obtained while creating a plot of residuals > plot (NV.zph, main = "groupNUSM - UNFIT", var= 'groupNUSM') Warning messages: 1: In approx(xx, xtime, seq(min(xx), max(xx), length.out = 17)[2 * : collapsing to unique 'x' values 2: In

hallo all, i would like to know, as i would suspect, nativ bridiging should work also, if only one iax party is behind an nat router and the other has a public ip. when i connect to iax clients, which have both pubic ip's nativ bridging is working. if one of the clients is behind an nat, the iax2 channels always get routed through the asterisk server (latest stable version from cvs) ?? i

Hi everyone, I am fitting a cox proportional hazard model with a continuous variable "x" as the covariate: fit<-coxph(Surv(time, status)~x) Now I wanted to make a plot of survival probability vs. the covariate, and the 95% confidence interval for the survival probability. It's just like a Kaplan-Meier Survival curve, except now the x axis represents the value of covariate, not

Hi, I am attempting to evaluate the prediction error of a coxph model that was built after feature selection with glmnet. In the preprocessing stage I used na.omit (dataset) to remove NAs. I reconstructed all my factor variables into binary variables with dummies (using model.matrix) I then used glmnet lasso to fit a cox model and select the best performing features. Then I fit a coxph model

Hi, How can I specify a Cox proportional hazards model with a covariate which i believe its strength on survival changes/diminishes with time? The value of the covariate was only recorded once at the beginning of the study for each individual (e.g. at the diagnosis of the disease), so I do not have the time course data of the covariate for any given individual. For example, I want to state at the

Hi, I am new to time-dependent Cox model to estimate time dependent hazard ratios. Let me use aml dataset from survival package: > aml3<-survSplit(aml2,cut=c(5,10,20),end="time",start="start", event="status",episode="i") If I want to esimate hazard ratio for each of the time intervals 0-5, 5-10, 10-20 and >=20, would the following calculate