similar to: xyplot() - can you control how the plots are ordered?

I have a need to call and pass arguments to nlme() from within another function. I use R version 1.8. I have found an apparent way to make this work, but I would appreciate some comments on whether this fix is really appropriate, or there is another way to do it that does not involve changing the source code. I don't have enough experience to start changing the sorurce code of a library

I apologize for posting on this question again, but unfortunately, I don't have and can't get access to MASS for at least three weeks. I have found some code on the web however which implements the prediction error algorithm in cv.glm. http://www.bioconductor.org/workshops/NGFN03/modelsel-exercise.pdf Now I've tried to adapt it to my purposes, but since I'm not deeply familiar

Hi, I'm using metafor to perform a mixed-effects meta-analysis. I'd like to test whether the effect is different for animals and plants/whether "group" (animal/plant) influences the effect size, but am having trouble interpreting the results I get. I've read previous posts about QM in metafor, but I'm still a bit confused. I've dummy-coded the factors:

Hi, I'm trying to fit a multinomial logistic regression to my data which consists of 5 discrete variables (scales 1:10) and 1000 observations. I get the following error: Error in `row.names<-.data.frame`(`*tmp*`, value = c("NA.NA", "NA.NA", : duplicate 'row.names' are not allowed In addition: Warning message: non-unique value when setting

*I have the following function:* /plot_mi_time = function(mdata, miname) { mdata2 = mdata[row.names(hakat) == miname, ] print(mdata2) xcoords <- c(1,1,2,2,3,3,4,4,5,5,6,6) plot(c(xcoords), mdata2, xaxt="n", ylab="Expression", xlab="Time(h)", , main=miname) axis(1, at=xcoords,

Hi, I have the following function that i use to plot graphs. plot_mi_time = function(mdata, miname) { mdata2 = mdata[row.names(mir_test3) == miname, ] # print(mdata2) xcoords <- c(1,1,2,2,3,3) plot(xcoords, mdata2, xaxt="n", ylab="Expression", xlab="Time(h)", , main=miname) axis(1, at=xcoords,

Dear R users, I would like one of my groups in xyplot to appear to lie 'behind' the other groups. I have searched for help and find many, many topics about panel order (e.g., "as.table"), but that is not what I need. What is odd is that the group that I want to appear 'behind' the others (PROFILE) is the last in alphabetical order but xyplot is plotting it first. Any

Apologetic prologue: I've looked through the mailing list for an answer to this (since I'm sure it's trivial) but I have not been able to find a fix. So the problem is that I want each group to have a different type of plot. "Probes" should be points and "Segments" should be lines (preferably using the segment plot command, but I've just been trying --

Hi all, I need to replace missing values in a matrix by 10 % of the lowest available value in the matrix. I've got a function I've used earlier to replace negative values by the lowest value, in a data frame, but I'm not sure how to modify it... nonNeg = as.data.frame(apply(orig.df, 2, function(col) # Change negative values to a small value, close to zero { min.val =

This post is NOT a question, but an answer. For readers please disregard all earlier posts by myself about this question. I'm posting for two reasons. First to say thanks, especially to Dimitris, for suggesting the use of errorest in the ipred library. Second, so that the solution to this problem is in the archives in case it gets asked again. If one wants to run a k-fold cross-validation

Hi I was trying to import a dataset from access database into R My code looks like below >library(RODBC) >setwd("//STATESFPSOP1/docs/CRMPO/DEPT/Retail PD") >channel<-odbcConnectAccess("MortgagePD.mdb") >Mdata<-sqlFetch(channel, "MortgagePD") But I got error messages saying Error: cannot allocate vector of size 27.9 Mb In addition: Warning

Dear R users, I encounter a problem regarding space control in xyplot. Basically, I want to control spaces between label, tick and axis. I remember there is a function called mgp in general plot. Is there a similar function for xyplot or xYplot? Below is my basic code: myplotkid<-xyplot(expected_offspringnumber~afr|decade,groups=SES,data1,

Page 140 of MASS uses the function xyplot. But I don't find it in R. Is there a package that I should load to use xyplot. Or there is a function with a different name in R that does the same thing as xyplot in S. xyplot(Gas ~ Temp | Insul, whiteside, panel = function(x, y, ...) { panel.xyplot(x, y, ...) panel.lmline(x, y, ...) }, xlab = "Average external temperature (deg. C)", ylab

I'm trying to add a set of reference lines to a multipanel xyplot xyplot(y ~ x | Visit, panel = function(x, y, ...){ panel.xyplot(x, y, ...) abline(v = c(0.5, 1)) }) However, the reference lines are different for different visits. For example, for the first 2 visits, I'd like vertical lines at x = 0.5 and 1. For visits 3 and four, I'd like vertical lines at x = 1 and 1.5. I can

Hello all, I am trying to use xyplot to give a six panel plot and to highlight only points (in any panel) that meet a certain criterion. With the plot command I would do something like: plot.default(filein$Site,filein$circ.conc) points(filein$Site,filein$circ.conc,type="p", pch=ifelse(filein$p.value<5e-02,19,21)) I had thought I could just stick in the pch line from above into

Hello R-Users, I am trying to label points on an xYplot. I can label the points, but then I lose the error bars. I suspect that I should be using panel.xYplot rather than panel.xyplot, but that didn't work either... Can anybody help me out? Example code is below. First xYplot has error bars, second tries to label the points and the bars disappear. Thanks, John a=c(1:4)

Hello >From the following script I received a grafic output when I called: - xyplot.test( 'green3' ) - call.xyplot.test( 'blue3' ) I did NOT receive a grafic output when I called: - loop.xyplot.test( 'red3' ) What's the Problem? NB: I am using R 1.4.1 on Linux. --------- START OF SCRIPT ---------------- n <- 1000 x <- seq( 1, n ) y <- rnorm( n )

In xyplot, I would like to get the "data into the panel function" in the following sense: Consider xyplot(scc~time|cowidp, data=cow.s,type=c("l"), panel=function(x,y,subscripts,...){ panel.xyplot(x,y,...) vvv<-cow.s[which(!is.na(cow.s[subscripts,"mastreat"])),"time"] panel.abline(v=vvv,col="red",lwd=2) } ) If I want to use a

Hallo, I'd like to add grid lines to a lattice graph having 2 series of Y data. See these 2 examples: data(iris) [1] xyplot(Sepal.Length + Sepal.Width ~ Petal.Length , data = iris, allow.multiple = TRUE, scales = "same",type="l", ) [2] xyplot(Sepal.Length + Sepal.Width ~ Petal.Length , data = iris, allow.multiple = TRUE, scales =

Dear R list members, I have a question regarding xyplot. I managed to make a xyplot of all the IDs by using the syntax below: xyplot(PA ~ CRPC + CRPT | ID, data = redinteract) Now, I'd like to make a graph of a specific ID number (e.g., only ID number 301). I thought I could use "subset", but it seems to be not working. Could anyone let me know how I can make a graph of a