similar to: Help on extract paramters from fitted models

Hi, thank you all for those who helped me on prediction of newdata for linear model, it is my new question on the prediction of coxph for newdata, for example, i have the model: coxout<-coxph( Surv(time, status) ~ x predict(coxout) will give the fitted values I have tried predict(coxout, newdata), it still gave me the fitted values only. Can anyone hlep me on how to do prediction for

Hi, I need to fit models and use AIC method to campare the best fitted model manually. When i extract AIC by using extractAIC, it gave me the df and AIC values. Now the problem is, how can I compare the AIC values from two models? is there anyway to extract AIC with no df so that I can compare directly? Thank you! > extractAIC(coxout) [1] 1.000 1723.038 [[alternative HTML version

Dear R-community, I have a short question: How do I interpret the result of a likelihood ratio test correctly? I am fitting a parametric survival model (with aftreg {eha}) and the output tells me the overall p-value of my model is < 0.001. My simple question is: Does the result mean my model fits the data well OR does it mean my model DOES NOT fit the data well? Some side information how the

Hi, I have a question about predicting new values from a fitted model in R. For example, if i fit a linear model fit<-lm(y~x) predict(fit) will give the fitted values for each x. I am wondering if there is a way to do some prediction of the fitted for some new values of x that are not in the data points? for example, x=0.5 is not in the original data (not a datapoint), is there a way to

Hi All: I have searched high and low with no success. Does anyone know of a way to pass the values of more than one form field through the :with parameter of observe_field? Basically I have two fields on the form with number values. When the user changes the value of one I would like to send the values of both fields to the controller for calculations and return the result to the user on

Is there a setting in gconftool-2 to that sets up all the sound and volume settings? If so I havenot found it. Can someone point me in the right direction. Or is this just done directly with alsacntl? THanks, Jerry

Hi, I would like to call a R script sometimes with: R --slave --vanilla --args Debug=1 N=25 < test.R and sometimes with R --slave --vanilla Debug=1 < test.R or similar. Inside the R script I get the arguments with: args=(commandArgs(TRUE)) if(length(args)==0){ print("=== No arguments supplied. =============") }else{ for(i in 1:length(args)){

Dear all, is this a bug or do a miss something (this concerns version 0.63.1 on Linux)? > boxplot (c (1, 2, 3, 4, 5), ylim= c (-5, 5)) Warning: parameter "ylim" couldn't be set in high-level plot() function > boxplot (c (1, 2, 3, 4, 5), pars= list (ylim= c (-5, 5))) Warning: parameter "ylim" couldn't be set in high-level plot() function > and actually in

Dear List, How do I extract the approximate Wald test for the frailty (in the following example 17.89 value)? What about the P-values, other Chisq, DF, se(coef) and se2? How can they be extracted? ######################################################> kfitm1 Call: coxph(formula = Surv(time, status) ~ age + sex + disease + frailty(id, dist = "gauss"), data = kidney)

I am calculating cox propotional hazards models with the coxph function from the survival package. My data relates to failure of various types of endovascular interventions. I can successfully obtain the LR, Wald, and Score test p-values from the coxph.object, as well as the hazard ratio as follows: formula.obj = Surv(days, status) ~ type coxph.model = coxph(formula.obj, df) fit =

I am trying to use t.test on the following data: date type INTERVAL nCASES MTF SDF MTO SDO nFST MF nOBS MO MB BIASCV BIASEV ME MAE RMSE CRCF 2001-06-15 avn GE1.00 4385 0.246 0.300 1.502 0.556 1367 1.373 4385 1.502 1.471 0.285 0.164 -1.256 1.266 1.399 0.056 2001-06-15 avn

We''ve all seen the productivity benefits of RoR. It would be worthwhile to figure out where they come from. For me, I''d say the following: 1. The number one benefit isn''t from Rails per se, but from Ruby, or any other dynamic OO language (like Python). They avoid the weaknesses of loosely typed, procedural languages like Perl & PHP, but have the swiftness

Maybe this is common knowledge, I don''t know. I figured I''d post it, hopefully it can help some people. If anyone has a better way, please share. http://jrhicks.net/120 (not my blog, I''m Mike in the comments I couldn''t get it to literally run the test in the identical way as rails does when you type rake. However, I could get it to run an individual

Dear R users I have a big matrix like 6021 1188 790 290 1174 1015 1990 6613 6288 100714 6021 1 0.658 0.688 0.474 0.262 0.163 0.137 0.32 0.252 0.206 1188 0.658 1 0.917 0.245 0.331 0.122 0.148 0.194 0.168 0.171 790 0.688 0.917 1 0.243 0.31 0.122 0.15 0.19 0.171 0.174 290 0.474

Hi, I am having a question. With multiple observations for each ID, I would like to keep only two or 3 observations for each ID. Does R have any command to do this? thank you very much! A simple example is below, How can I only keep the first 2 observations for each j? j x 1 1 0.795373270 2 1 0.326845207 3 1 0.049116967 4 1 0.673830996 5 2 0.411789618 6 2 0.628034020 7

Richard Smith <richard at metafoo.co.uk> writes: > On Fri, Mar 28, 2014 at 4:50 AM, Renato Golin <renato.golin at linaro.org>wrote: >> On 28 March 2014 11:16, Chandler Carruth <chandlerc at google.com> wrote: >> > Just the reserved part. >> >> Ok, in that case, I share you concerns. >> >> We could easily only implement the reserved ones

Hello: Please be patient, its hard to me to explain this even in my native language. I think it may be useful to add an option like lograte's prerotate and postrotate, that triggers server side scripts before and after rsync. It would help the administrator work. Escenary 1: ----------- client -> rsync server -> 2nd stage rsync server Want to rsync from rsync server to 2nd stage

Hello everybody, I am sorry if my questions are too simple or not easily understandable. I’m not a native English speaker and this is my first analysis using this function. I have a problem with a cross correlation function and I would like to understand how I have to perform it in R. I have yearly data of an independent variable (x) from 1982 to 2010, and I also have yearly data of a variable

Der r helpers, I have a following set of data. Though some other software, I know it follows 3 parameter Gamma distrinution. I need to write a R code for estimating the parameters of 3 Parameter Gamma distribution. c(0,0,18561.9,0,0,0,34400,0,0,0,0,2190,0,0,0,0,60000,0,0,0, ? 19583,0,0,0,109872.87,0,0,0,0,0,0,1244,0,0,25150,0,500,0,0,0, ?

Since the WEBrick server has been so beefed up in Rails 0.9, I think it''s a darn shame that its lingering outside the support of custom pretty URLs. Also, it seems that the connection between mod_rewrite and url_for and friends is pretty fuzzy for most people. Hence, as we''ve talked about before, the rewriting of the url should be pulled in-house. There has been a few