similar to: plot.ranef.lme (PR#2986)

Hi, I just wonder if R is still used for PK analysis. However, I have to use R for some purposes anyways. Here is a really strange problem I am running into and I hope to get some assistance from the group. The following code works fine in S-Plus (Population estimates and MAP estimates look reasonable ). cs.fit<-nlme ( CONC~phenoModel(CID,TIME, DOSE, lCl, lV), fixed = lCl+lV ~1,

Hi, The solution to my problem is to use 'lsegments' instead of 'segments' within lattice commands. (Although I wont forget again, a comment in the segments help file referring to 'lsegments' might help others not to make the same mistake in the future.) My thanks to Renaud Lancelot. Regards, John. John Gavin <john.gavin at ubsw.com>, Quantitative Risk Models

Hi, How might I use xyplot to plot segments where the segments are in the input data? (ie a directed acyclic forest). Here's an example in base graphics: n = data.frame(id = c(1,2,3,4), parent = c(0,1,2,2), value = c(5,5.5,7,3), date = c(1,2,3,3.5)) plot(n$date, n$value) do.call( segments, with( merge(n,n,by.x="parent", by.y="id"),

Dear R-helpers, I have a function that adds some segments to the current plot, but if I the current plot is made with any of the lattice functions (in my case, levelplot), I should use lsegments rather than segments. How can I tell if the current device was made with plot() or e.g. levelplot() or another lattice function? dev.cur() does not help me, it just tells me "windows 2".

If I have a contourplot (in the lattice package) and I want to add straight lines to it, how do I do this? I see that there are llines() and lsegement() functions for lattice plots, but they don't seem to do anything in this case: library(lattice) library(KernSmooth) x=rnorm(10000) y=x+rnorm(x,0,.5) a=bkde2D(cbind(x,y),.7) z=as.vector(a$fhat) grid=expand.grid(x=a$x1,y=a$x2) grid$z=z

Dear R-Users, how to include error bars within lattice? How should the panel = function(x,y,...){ looks like? Does panel.arrows works here as well? I appreciate any help on this. Regards, Mario AT [[alternative HTML version deleted]]

Hi I am following the model building strategy that is outlined in the Pinheiro and Bates book wrt including covariates but am having a problem with the plot. Basically I am using 4 covariates (1 of them is continuous) and 3 of them are fine but the 4th one is being shown as a scatterplot despite the fact that it is a factor. I have explicitly declared this to be a factor (pcat<-as.factor(pcat))

First, apologies for no example data but I don't think it's needed in this case, Q: Can (and if so how ) the line end style be changed for 'cloud' plots? I've tried par(lend=2), trellis.par.set(add.line = list(lend=2)) and much googling but to no avail Thanks in advance Dan P.S. the reason for this is that the round end looks bad at lwd=3 or more Daniel Alcock Malaria

Hi all, I have been previously been using scatterplot3d package to create some graphs but unfortunately it does not allow me to rotate the plot on all three axis. The cloud() function in the lattice package does allow me to do so. When I was using scatterplot3d I was using a script (Shown Below) to calculate the mean, quartiles and range limits for all three axis and I was representing that on the

Dear r-help, I need to draw xyplot() graphs with several regression lines: one line for the whole range of x (the variable on the horizontal axis) and two additional lines for subranges of x. Is it possible to make first regression line (panel.lmline(x,y,...);) to be drawn on the whole graph and regression lines of the subsets to be drawn only over their subsets? I have defined

In the example: x = rep(c(0,0,1,1),4) y = rep(c(0,1,0,1),4) z = c(1,0,1,0,0,0,1,1,0,1,0,0,1,1,1,1) f = as.factor(c(rep("a",4),rep("b",4),rep("c",4),rep("d",4))) levelplot(z~x+y|f,data.frame(x=x,y=y,z=z,f=f)) I noted that the last ("d") plot remains empty. I guess the reason for this is that the values are constant (1), but I consider it more

Dear R list members, I'd like to produce rug plots at the bottom of panels in a trellis display (using the lattice package), but par("usr") doesn't return user coordinates for panels, and consequently rug fails, as the following example (suggested to me by Georges Monette) illustrates: > x <- rnorm(50) > y <- rnorm(50) > f <-

Hi, There is a double "the" in x11.Rd: "...which is faster will depend on the the X11 connection" Thank you, Stephen -- Rochester, Minn. USA

I'm looking at Pinheiro & Bates "Mixed-Effects Models in S and S-PLUS" at the moment. Several datasets are used, one of which is called "PBIB" (a partially balanced incomplete block design). All the other datasets can be found somewhere or other in R. However, I cannot locate PBIB, and it does not seem to be mentioned in the latest edition of the R Full Reference

Hi There are new versions of the grid and lattice packages on CRAN: grid 0.2 and lattice 0.3-0 These packages are still _under development_ (they are in CRAN/src/contrib/Devel) There have been some bug fixes and, especially in lattice, there are lots of new features. Descriptions of the changes are included below. If you have any grid code, it may break under the new version because of a

Hi There are new versions of the grid and lattice packages on CRAN: grid 0.2 and lattice 0.3-0 These packages are still _under development_ (they are in CRAN/src/contrib/Devel) There have been some bug fixes and, especially in lattice, there are lots of new features. Descriptions of the changes are included below. If you have any grid code, it may break under the new version because of a

I'm looking for a way to plot lines on top of a levelplot(), where the lines are borders between cells of different values. The clines() function provides contours suitable for continuous data. I am dealing with discrete values, spatial clusters of nodes where each cluster has an integer value, and I want to plot the borderlines between these areas. So, in a levelplot having nodes with

Hi, I use panel.points to add points to a xyplot graphic. But I like to use the sunflowerplot to plot my points because this is very superimposed. It is possible to use this? I try but it dont work directly. It may be need to put this function inside a panel.??? Thanks Ronaldo -- Where there's a will, there's a relative. -- > Prof. Ronaldo Reis J?nior | .''`.

> To sum up, I can't figure out how MLWin calculates the > standardized residuals. But I understand that this is not a > question for the R list. > Nevertheless, it would help if someone could point me to some > arguments why not to use them and stick to the results > obtainable by ranef(). Hi Dirk: Well, it is interesting that mlWin and lmer generate the same exact

Hello. I'm running a generalized linear model and am interested in using the random effects that are output for further analysis. My random effect is interacting with two different fixed effects (which which are factors with two levels each). When I retrieve the random effects I get something like this: (Intercept) nutrient (Intercept) light (Intercept) Aa-0 0.59679192