similar to: Reporting odds ratios or risk ratios from GLM

The code library(Design) f <- lrm(y~x1+x2+x1*x2, data=data) plot(f) produces a plot of log odds vs x2 with 0.95 confidence intervals. How do I get a plot of odds ratios vs x2 instead? Thanks -- View this message in context: http://r.789695.n4.nabble.com/Odds-ratios-from-lrm-plot-tp4033340p4033340.html Sent from the R help mailing list archive at Nabble.com.

Hello, I have an example of a 2 paneled plot, with two different aspect ratios displayed on one page. An example would help n=20 x1 <- cumsum(runif(n)) x2 <- cumsum(runif(n)) d <- data.frame(val=c(x1,x2),id=c(1:n,1:n), nt=c(rep("A",n),rep("B",n))) u1 <- xyplot(val~id | nt, data=d,aspect=1,layout=c(1,2)) u2 <- xyplot(val~id|nt, data=d,aspect=0.5,layout=c(1,2))

I am trying to calculate Poisson regression based adjusted rate ratios in R, but R's default in glm does not code the intercept as the global rate. In SAS I use "cell means" coding so that the intercept is the global rate, but I do not know how to do this in R. If anyone knows a way to make glm use "cell means", or an how to find adjusted rate ratios I would be grateful.

hi, I'm looking for a package to compare two hazard ratios (and assign statistical significance) obtained from two different predictive models. I know of the hr.comp2 function from the survcomp package, but was wondering if there's any other packages out there. thanks! [[alternative HTML version deleted]]

I partly answered my question since independence_test() function in coin package apparently do Cochran-Armitage trend test just like Eric Lecoutre's function tabletrend() - slightly modified here: > independence_test(pheno ~ geno, data = dat2, teststat = "quad", scores = list(geno = c(0, 1, 2))) Asymptotic General Independence Test data: pheno by groups 1 < 2

I am trying to do logistic regression with a categorical predictor variable with the glm() function, family=binomial. Using glm() I would like to be able to calculate the confidence intervals of all three possible odds ratios for a factor (the factor has three categories). Three categories imply two columns of 0's and 1's in the design matrix, and two parameter estimates with their

Dear R users: Analysis of the impact of a time-dependent covariate (GVHD or use of steroid after bone marrow transplantation) on two competing endpoints (invasive fungal infection and death) is frequently encountered in the setting of BMT data. Coxph package can be used as the following: for the analysis of GVHD: > gvhd -> coxph(Surv(start,stop,status = =1) ~ GVHD, data=bmt.data)

My apologies if this is obvious: Is there a simple way (other than simulation or bootstrapping) to obtain a (approximate)confidence interval for the ratio of 2 variance components in a fitted lme model? -- In particular, if there are only 2 components (1 grouping factor). I'm using nlme but lme4 would be fine, too. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA

Hi, Is there a R test available that tests whether 2 hazard ratios obtained from Cox regressions on the same patient sample by 2 different classifiers are significantly different? Thanks Steve

Hi folks, I have assumed that ratios of variance components (Fst and Qst in population genetics) could be estimated using the output of mcmcsamp (the series on mcmc sample estimates of variance components). What I have started to do is to use the matrix output that included the log(variances), exponentiate, calculate the relevant ratio, and apply either quantile or or HPDinterval to get

Greetings, I am new to R and have been comparing CPH survival analysis hazard ratios between R and SAS PhReg. The binary covariates' HRs are the same, however the continuous variables, for example age, have quite different HRs although in the same direction. SAS PhReg produces HRs which are the change in risk for every one increment change in the independent variable. How do I

Hi, does anyone know of a function to calculate odds ratios in multiway tables (stratified) (+ the other usual statistics involved) i mean: say we have a table r*c*d, For every d (depth) we have a r*c table, and in this table the odds ratio's are calculated for every 2*2 subtable in it. logically this function would look like): ORs(multiwaytable) or ORs(data$var1r,data$var2c,data$var3d)

Hello there, I am trying to run an ANOVA model using a non-Standard F ratio. Imagine that the treatments (treatments 1 & 2) are applied to the row not to individual samples. Thus the row is the experimental unit. Therefore my error term in my ANOVA table should be the error associated with with row. The question is how do I check the assumptions of an ANOVA model when I have a non-standard F

Dear R users, Here is the coxme output I obtain on my survival dataset having 3 strains and 2 infection status (i: infected, ni: non infected) coxme(Surv(lay) ~ infection*strain, data=datalay, random= ~1 |block) Cox mixed-effects model fit by maximum likelihood Data: datalay n= 1194 Iterations= 3 77 NULL Integrated Penalized Log-likelihood -7270.028 -7223.859 -7218.175

so here is the code I need help to fix the line for 'st' when it reaches maximum, that's where i am getting stuck. thanks Sharad set.seed(100) d=data.frame(x=rnorm(20)+5, x1=rnorm(20)+5, x2=rnorm(20)+5, x3=rnorm(20)+5, x4=rnorm(20)+5, x5=rnorm(20)+5, x6=rnorm(20)+5, x7=rnorm(20)+5, x8=rnorm(20)+5) r2=c() cinit=0 for (col in 1:9) { cinit=cinit+1 st=cinit+1

Hello all, I am having some trouble finding a good csv file or even any kind of table with the financial ratios (i.e. P/E ratio, ROI, ROA etc.) from the Internet. I was able to find a good source with poor information. That wouold be Finviz.com. Another option would be creating a XML function loop to go through YahooFinance website, get the information and store it in a variable. This, however,

Hello All, I am writing to ask your opinion on how to interpret this case. I have two vectors "a" and "b" that I am trying to compare. The wilcoxon test is giving me a pvalue of 5.139217e-303 of a over b with the alternative "greater". Now if I make a summary on each of them I have the following > summary(a) Min. 1st Qu. Median Mean 3rd Qu.

i need o know what are the ratios used in language 'R'. i just get the sharp, Treynor, information ratios, is ther any other ratios used in it ? and also please give any reference web page link too. Thanks . -- View this message in context: http://r.789695.n4.nabble.com/ratios-tp4648826.html Sent from the R help mailing list archive at Nabble.com.

Dear useRs, I have two different length vectors: one column (1...m) and one row vector (1...n): 20 40 20 60 5 4 2 Now I have to calculate ratios between column vector elements and each row vector elements: 4 5 10 8 10 20 4 5 20 15 12 30 Thank's in advance for any suggestions, Andrej

Does anyone know of another version of the Mcnemar test that provides: 1. Odds Ratios 2. 95% Confidence intervals of the Odds Ratios 3. Sample probability 4. 95% Confidence intervals of the sample probability Obviously the Odds Ratios and Sample probabilities are easy to calculate from the contingency table, but I would appreciate any help on how to calculate the confidence