similar to: coxph and variables

Hi, I used Dr. Harrell's rms package to make a nomogram. Below is my code for nomogram and calculate total points and probability *in original data set* used for building the nomogram. *My question is how I get the formula for calculating the survival probability for this nomogram. Then I can use this formula to do validation by using other data set. * f1 <- cph(Surv(retime,dfs) ~

Hello Everyone, I have what I think is a complex text parsing task. I've provided some sample data below. There's a relatively simple version of the coding that needs to be done and a more complex version. If someone could help me out with either version, I'd greatly appreciate it. Here are my sample data. haveData <- structure(list(profile_key = structure(c(1L, 1L, 2L, 2L, 2L,

Hello, I am using Dr. Harrell's rms package to make a nomogram. I was able to make a beautiful one. However, I want to change 5-year survival probability to 5-year failure probability. I couldn?t get hazard rate from Hazard(f1) because I used cph for the model. Here is my code: library(rms) f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11,

Hello, I am using Dr. Harrell's design package to make a nomogram. I was able to make a beautiful one. However, I want to change 5-year survival probability to 5-year failure probability. I couldn?t get hazard rate from Hazard(f1) because I used cph for the model. Here is my code: f1 <- cph(Surv(retime,dfs) ~ age+her2+t_stage+n_stage+er+grade+cytcyt+Cyt_PCDK2 , data=data11, surv=T,

COMPUTATIONAL ANALYSIS OF NEXT GENERATION SEQUENCE DATA (http://bioinformatics.nki.nl/vacancies.php) THREE BIOINFORMATICS POSITIONS PROJECT OUTLINE Genomic alterations are major determinant of responses to (targeted) therapies in cancer. In fact, the best positive and negative predictors of responses to targeted therapies are alterations in kinases or their direct downstream effectors. To

Hello, I am using the clogit (conditional logistic regression) on a simple data set which is not related to survivorship, which I understand to be fine. I have trouble understanding the output. I would like to find parameter estimates of the logistic models I am constructing with R, and am unclear as to what the coefficients represent that are provided. The below is a copy of my output.

Dear list, I have the following matrix. How can I make the following plot? 1. The x-axis has index 1:7, and the first column is plotted against index 1, second against 2, and so on. 2. I want the points from the left upper conner including the antidiagonal to be plotted with col=2, and the rest with col=3 [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 0.589 0.857 0.923 0.944 0.954 0.963

Hello, i'm a french student of medical oncology and i'm working on breast cancer. I have a variable with the histologic type of tumor wich is between 1 and 5. I use as.factor function to make some variable with level between 1 and 5. When i put it in the cox model i have only the level between 2 and 5. The level 1 doesn't appear. I think i have to change the number of level but i

Hello, I have tried reading the documentation and googling for the answer but reviewing the online matches I end up more confused than before. My problem is apparently simple. I fit a glm model (2^k experiment), and then I would like to predict the response variable (Throughput) for unseen factor levels. When I try to predict I get the following error: > throughput.pred <-

I have a situation with a large dataset (3000+ observations), where I'm doing lags as regressors, where I get: Call: lm(formula = rj ~ rM + rM.1 + rM.2 + rM.3 + rM.4) Residuals: 1990-06-04 1994-11-14 1998-08-21 2002-03-13 2005-09-15 -5.64672 -0.59596 -0.04143 0.55412 8.18229 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) -0.003297 0.017603

Hola que tal alguien de ustedes me podría auxiliar, lo que pasa es que estoy tratando de hacer un programa en R utilizando la siguiente tabla: x<1 km x>1km Estabilidad a c d f c d f A 213 440.8 1.94 9.27 459.70 2.094 -9.6 B 156 106.6 1.15 3.30 108.20 1.098 2.0 C 104 61.0 0.91 0.00 61.00 0.911 0.0 D 68 33.2 0.73 -1.70 44.50 0.516 -13.0 E 51 22.8 0.68 -1.30 55.40 0.305

Hi list members, I saw a couple lavaan posts here so I think I?m sending this to the correct list. I am trying to run a multigroup analysis with lavaan in order to compare behavioural correlations across two populations. I?m following the method suggested in the paper by Dingemanse et al. (2010) in Behavioural Ecology. In one of the groups, lavaan returns negative variance for one path and I?m

Hi, I would like to fit parametric survival models to time-to-event data that are left truncated. I have checked the help page for survreg and looked in the R-help archive, and it appears that the R function survreg from the survival library (version 2.16) should allow me to take account of left truncation. However, when I try the command

Dear R fans, I would like to create a scatterplot showing the relationship between 4 continuous variables. I thought of using the package "scatterplot 3d" to have a 3-dimensional plot and then using the size of the symbols to represent the 4th variable. Does anybody know how to do this? I already tried to create this graph using the colour of the symbols, but I was unable to generate

Tobias, Thanks for taking the time to summarize all this. It's a great writeup. I'm moving the thread to llvm-dev. My responses below. On May 21, 2012, at 5:06 AM, Tobias Grosser <tobias at grosser.es> wrote: > First of all some information about the RegionInfo pass: > > ======================================================================= > The very first paper I

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Hi, I've data from an incomplete fatorial design. One level of a factor doesn't has the levels of the other. When I use lm(), the summary() return NA for that non estimable parameters. Ok, I understant it. But I use contrast::contrast(), gmodels::estimable(), multcomp::glht() and all these fail when model has NA estimates. This is becouse vcov() and coef() has different dimensions. Is

... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal. On May 15, 2018 7:04:31 AM PDT, Bert Gunter

Unless there is good reason not to, always cc the list -- there are lots of smarter folks than I on it who can help. I may or may not have time to look at this. Hopefully someone else will. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip