similar to: p value for joint probability

Dear all, I was intrigued by the ctree command and wanted to check it out. I first ran the demo with example(ctree) and did get the survival graphs in the end. Upon doing this with my own data and yielding a "Invalid operation on a survival time" I tried to rerun example(ctree) and now I also get "Invalid operation on a survival time" after the example runs plot(GBSG2ct)...

Dear R?ers I?m trying to get a summary table using latex and summary in the rms package to no avail. I?m running R 2.10.1, Mac OS X snow leopard and I have the mactex 2009 distribution installed. Any obvious things I?m missing? //M options(digits=3) set.seed(173) sex <- factor(sample(c("m","f"), 500, rep=TRUE)) age <- rnorm(500, 50, 5) treatment <-

Dear R´ers.. In this mock dataset how can I generate a logical variable based on whether just tes or tes3 are NA in each row?? test<-sample(c("A",NA,"B"),100,replace=T) test2<-sample(c("A",NA,"B"),100,replace=T) test3<-sample(c("A",NA,"B"),100,replace=T) tes<-cbind(test,test2,test3)

Dear R?ers. I?ve got a dataset with age and year of diagnosis. In order to age-standardize the incidence I need to transform the data into a matrix with age-groups (divided in 5 or 10 years) along one axis and year divided into 5 years along the other axis. Each cell should contain the number of cases for that age group and for that period. I.e. My data format now is ID-age (to one

glm(A~B+C+D+E+F,family = binomial(link = "logit"),data=tre,na.action=na.omit) Error in `contrasts<-`(`*tmp*`, value = "contr.treatment") : contrasts can be applied only to factors with 2 or more levels however, glm(A~B+C+D+E,family = binomial(link = "logit"),data=tre,na.action=na.omit) runs fine glm(A~B+C+D+F,family = binomial(link =

Dear all, Using the SetInternet2(TRUE) option works wonders with R in my sealed down work-environment. However, I'd like to use RCurl and apparently the proxy settings are not carried over. Is it possible to figure out the proxy-IP and port number from R after invoking SetInternet2? //M

Dear all, I?m having trouble getting a list of regression variables back into a dataframe. mydf <- data.frame(x1=rnorm(100), x2=rnorm(100), x3=rnorm(100)) mydf$fac<-factor(sample((0:2),replace=T,100)) mydf$y<- mydf$x1+0.01+mydf$x2*3-mydf$x3*19+rnorm(100) dlply(mydf,.(fac),function(df) lm(y~x1+x2+x3,data=df))->dl here I?d like to use ldply(dl,coef(summary)) or something

Dear R´ers, How can I create one single factor variable from two variables incorporating all possible combinations of the values?? test<-sample(c("A",NA,"B"),100,replace=T) test2<-sample(c("E",F,"A"),100,replace=T) tes<-cbind(test,test2) pseduocode: r<-function(test,test2) r AE AF AA NAE NAF NAA BE BF BA //M [[alternative HTML

Dear All, I?ve been trying to find an option to scatter3D from rcmdr to remove the individual points from the plots but to no help so far. Removing the residuals is easy, but I cannot find a similar point option. Is there such an option that can be set to FALSE? Best, //M

Dear all, With the latest update of Hmisc I no longer have any problems with latex. However using the ctable option produces latex code that at least on both the miktex distribution at work and mactex distribution at home refuses to run due to an extra blank line inserted between the multicolumn lines in the latex code... It runs fine if the line is deleted or if the ctable option is left out.

Dear R family, I want to calculate the joint probability (distribution) of two random continuous variables X and Y. Could to please tell me how to do it?Thanks in advance.. elisa [[alternative HTML version deleted]]

Hello, I want to fit a tree parameter distribution to given data. I tried it with sample data using the "fitdistr" function. Here my workflow that didn't had any result: I started with the generalized gamma distr, which is: r*dgamma(x^r,shape,rate) The R-function is: ggamma = function (x,r,shape,rate) r*dgamma(x^r,shape,rate=rate) For the first step I assumed r = 1 and I

I?m trying to run a logistic joinpoint regression utilising the ljr package. I?ve been using the forward selection technique to get the number of knots for the analysis, but I?m uncertain as to my results and the interpretation. The documentation is rather brief ( in the package and the stats in medicine article is quite technical) and without any good examples. At the moment I?m thinking 1)find

Hi Roscoe, I couldn't make it though that excel worksheet without my eyes glazing over... but do you want something like data(mtcars) tab <- xtabs(~cyl + gear, data=mtcars) tab/sum(tab) ? Best, Isat On Mon, May 3, 2010 at 10:08 PM, Roscoe <Roscoe at exemail.com.au> wrote: > Dear R-Help, > > First post, be gentle... > > I am new to R: I can perform basic functions

Dear all, Is there an implementation of Simon & Makuch method of plotting the survival function with time-dependent variables. I?m only able to find event.chart in Hmisc for the purpose and I would prefer the Simon and Makuch method. I believe stata has it implemented for this purpose, but I cannot find it on CRAN. Simon R, Makuch RW. A non-parametric graphical representation of the

Dear all, I just received a file from a colleague in spss. The read.spss could not finish the file due to an error (Unrecognized record type 7, subtype 18 encountered in system file) so instead I converted the file using stat-transfer. Looking at my data I see that most labels are in the attributes and I?d love to access them and assign the pertinent variables to factors without doing the whole

Dear all, I did post this more or less identical mail in a follow up to another question I posted, but under another heading. I try again, but now under the correct header. upon running this code (from the Hmisc library-latex function) I believe the call to summary.formula is allright and produces wonderful tables, but the latex command results in a correct formatted table but where all the

Hello, I am using sm.density() to find the bivariate joint PDFof events: For eg, x<-cbind(rnorm(30),rnorm(30)) den<-sm.density(x) Then I get the joint pdf from den$estimate in order to constructthe joint cumulative PDF. However, summing up all the values from den$estimateisnot equal to 1(have multipliedby the grid size). Anyone could help? Thanks. mc [[alternative HTML version

Hello list. I need help (e.g., a reference, code, package, etc.) in calculating the joint entropy of many variables (some sure highly mutually-informative and some not). Is there anyone here who knows a computationally-efficient solution (such as an R package)? I appreciate you help ... Best, Reza [[alternative HTML version deleted]]

Strange. Both systems are Linux. I compiled the LDAP slave Samba with CFLAGS="-O2 -march=i686" and SuSE's GCC 3.3.1. The fine working master got the same CFLAGS but SuSE's GCC 2.95.3. Without CFLAGS the GCC 3.3.1 works fine too. Didn't know that optimization (and GCC>3) is that dangerous.... Thanks anyway Daniel ------ Hi! We have a replicated LDAP Samba network. I