similar to: CI using ci.numeric

Hi all: My data has 3 variables: age(3levels : <30y=1 30-50y=2, >50y=3) gender(Male=0, Female=1) CD4 cell count(raw lab measurement) y(1:death 0:alive) I perform logistic regression to find out the factors that influence y. result<-glm(y ~ factor(age) + factor(gender) + CD4,family = binomial) >From the result,I can get OR(Odds Ratio) of gender via exp(Estimate of Female,

Dear All, I am new to ordinal logistic regression. Using ordinal regression within the R Commander GUI, I have generated an independent variable that is significant, but whose 95% confidence intervals slightly crosses "1". Is this possible? Here is the syntax and output: polr(formula = CDIcat ~ Employment, data = CDIallvariables, Hess = TRUE, method = "logistic")

Dear all, I need to do a loop in R, but I am not sure the software is generating "n" times the variables I request differently. When I ask to print the last matrix created, I just can see the loop for n=1. To be more precise, supose I need to simulate 10 times one variable and I want to fit the 10 variables simulated in a matrix. I dont really know what I am doing wrong, but I just

Dear list, I am calculating the 95th percentile of a set of values with R and with SPSS In R: > normal200<-rnorm(200,0,1) > qnorm(0.95,mean=mean(normal200),sd=sd(normal200),lower.tail =TRUE) [1] 1.84191 In SPSS, if I use the same 200 values and select Analyze -> Descriptive Statistics -> Frequencies and under "Statistics", I type in '95' under Percentiles,

Dear sir, Dear Madam I begin with R and I couldn't find in the manuels how I can solve this this problem find n such as |1/ni*sum(Wi)|<0.2 with Wi<-c(rnorm(200,0,1)) Thanks for helping Best regards ------------------------------------------ [[alternative HTML version deleted]]

Hello, We have a question about how to retrieve the nonparametric curve in gam() function. Now, we get the estimates and draw the fitted curve using the following code. You can see that our fitted curve is parallel to the output from gam() function but differs by a constant shift. However, the magnitude of the shift is not the intercept term. Also, gam's estimate appears more close to the

Dear List-Members, I have a problem with the function lqs() from package MASS. In some cases it produces different results for the same settings and needs a random seed to be set, in other cases not. I really cannot understand, why this happens. As well I do not understand what exactly you need the random seed for. Is it a starting point for iterations? Or do different results occur because of