similar to: Modifications to an abline

In a previous post, I mentioned a loop being used to generate graphs. I have some sample code partially put together but have found one offending line of code that I cannot figure out what to do with. I have one data frame called grade4. If I do something like hist(grade4$math) I get the appropriate chart. Within the loop, however, I am doing this for multiple files and grades, so I use

Thanks for sharing the questions and responses! Is it possible to appreciate how much the coefficients matter in one or the other model? Say, using Biau's example, using coxph, as.factor(grade2 == "high")TRUE gives hazard ratio 1.27 (rounded). As clinician I can grasp this HR as 27% relative increase. I can relate with other published results. With survreg the Weibull model gives a

Dear R help list, I am modeling some survival data with coxph and survreg (dist='weibull') using package survival. I have 2 problems: 1) I do not understand how to interpret the regression coefficients in the survreg output and it is not clear, for me, from ?survreg.objects how to. Here is an example of the codes that points out my problem: - data is stc1 - the factor is dichotomous

1. The weibull is the only distribution that can be written in both a proportional hazazrds for and an accelerated failure time form. Survreg uses the latter. In an ACF model, we model the time to failure. Positive coefficients are good (longer time to death). In a PH model, we model the death rate. Positive coefficients are bad (higher death rate). You are not the first to be confused

Hi Igor, Thanks so much for troubleshooting all this while and we found out none of our configuration is the problem but the source code. Hope that the samba team will modify to a working code so that I can deploy it. Actually my dateline to deploy is coming soon and I do not know what to do now..... when do you think the code will be modified and be released? Thanks so much for your help.

Hi there, I have a rather large data set and perform the following cox model: test1 <- list(tstart,tstop,death1,chemo1,radio1,horm1) out1<-coxph( Surv(tstart,tstop, death1) ~ chemo1+chemo1:log(tstop+1)+horm1+horm1:log(tstop+1)+age1+grade1+grade1:log(tstop+1)+positive1+positive1:log(tstop+1)+size1+size1:log(tstop+1), test1) out1 Up to here everything works fine (with each covariate

I have a matrix with 3 years of data (2006, 2009, 2010). I am trying to split this matrix by year so that I have 3 separate matrices. My matrix looks like this:     Q16.1 Year Gender Grade1       3 2006      1     52       2 2006      0     53       3 2006      0     54       3 2006      0     55       3 2006      1     56       2 2006      0     57       3 2006      0     58       3 2006      1  

Hello R List, I'm trying to run a linear mixed effects analysis on fairly large datasets with lmer (from the lme4 package) on a 32-bit Windows machine running XP with 3 GB of RAM. It's not working. (details below) I've researched the ff and bigmemory packages, but it appears they won't handle the mixed mode dataset I'm analyzing. It has some character fields for the

Hi Igor, Once again, thanks for keeping up with me. I have been migrating my master ldap server to 2.1 version so to keep it the same with the PDCs version of LDAP. Now they are the same. I have rectified such that "wbinfo -u" on both sides worked now. I am made "net rpc trustdom list" worked. It was not working before. I had to put "stuadmin = root" in the

Hello fellows, I would like to create a sequence for repeated numbers in a dataset. For example: ID <- c(1:20) grade <- c(4,4,4,5,5,7,7,7,7,8,8,8,9,9,9,9,9,10,10,10) Data: ID Grade 1 4 2 4 3 4 4 5 5 5 6 7 7 7 8 7 9 7 (...) I would like to create a variable "sequence": Data: ID Grade Sequence: 1 4 1 2 4 2 3 4 3 4 5

Dear list: I have some data for which I am generating a series of barplots for percentages. One issue that I am dealing with is that I am trying to get the legend to print in a fixed location for each chart generated by the data. Because these charts are being created in a loop, with different data, my code searches the data to identify the maximum value in the data and then print the data values

Hello, I have a dataframe (var_1, var_2, ..., var_n) and I would like to export summary statistics to Latex in the form of a table. I want specific summary statistics by crossing numerous variables 2x2 AT ONCE. In each cell I would like sometimes to have the median (Q1 - Q3), or frequency and proportion, etc. CrossTable, xtab, etc... do not allow for multiple 2 by 2 crosstabulation. The table

Dear R-users, The help does not give much details on the use of the arguments 'from' and 'to' in the panel.abline function. I have looked in the archives but did not find how to implement them. My different tries failed miserably. E.g, the following code doesn't seem to work, in a sense that the line is not limited to the (0,10) range. Do these arguments really apply to

> Hi, > > I am plotting time series data, then using abline to draw a horizontal > line through the graph (for mean, etc). I would like to be able to > display the level at which the abline is drawn on either the left > (main) > Y axis or a supplemental right Y axis. > > Additionally, I would like to display 1 and 2 standard deviations of a > certain measure

Dear All, I am wondering on how to use the abline.xyplot with xyplot such that I will have different vertical lines for each panel. More sepcifically, suppose that the xyplot generates 4 panels defined by the combination of two binary variables: X_1 and X_2. i.e. xyplot(Y ~ Z | X_1*X_2, data = df) I want something like: abline(v = 5) if X_1=0 and X_2 = 0 abline(v =

Dear R-users, I am unsuccessful in trying to add an horizontal line to all graphs in the example below: library(lattice) val<-runif(15) x<-rep(seq(1:5),3) type<-c(rep("a",5), rep("b",5), rep("c",5)) xyplot(val ~ x | type, panel.abline(h=.6)) Any hints are appreciated. Best regards, ---Jean -- View this message in context:

> In a lattice plot like this: > > win.graph() > xyplot(tmx~frequ|as.factor(as.numeric(spf)),groups=as.factor(blm), > data=tmx,type="l",pch=16,xlab="frequency (N)",ylab="Area held (ha)", > auto.key=list(blm,points=F,lines=T,title="Blm factor",cex.title=0.7, > cex=0.7,corner=c(1,1)),main="Mangroves target=1573ha",layout=c(2,3))

Is it kosher to add, say, horizontal lines to a lattice xyplot using abline(h=xxx)? The lines don't appear at the right value, that is, if h=150, the line might appear at h=140? Jeff This communication, including any attachments, is intended solely for the use of the addressee and may contain information which is privileged, confidential, exempt from disclosure under applicable law or

Hi, Is there a way I can draw an abline on a heatmap? I try the abline function, but don't get the line. My sample code is: mat <- matrix(sample(1:100,40),nrow=5) heatmap(mat,col=greenred(75),trace="none", dendrogram = "column",labCol = NULL) abline(h=5,v=4) thanks! [[alternative HTML version deleted]]

A non-text attachment was scrubbed... Name: not available Type: text Size: 1947 bytes Desc: not available Url : https://stat.ethz.ch/pipermail/r-help/attachments/19970719/f75cd159/attachment.pl