similar to: Simple suggestion for improvement

Hi, I'm trying to add points to 2 plots on the fly using par(mfg=vector) so switch between them. However, the appropriate scales aren't switched when changing from one plot to another, e.g. par(mfcol=c(2,1)) plot(1,1, col="blue") # blue plot plot(1.2,1.2, col="red") # red plot points(1.1,1.1) # appears to bottom left of red point par(mfg=c(1,1)) #

Just a quick point which may be easy to correct. Whilst typing the wrong thing into R 2.2.1, I noticed the following error messages, which seem to have some stray quotation marks and commas in the list of available families. Perhaps they have been corrected in the latest version (sorry, I don't want to upgrade yet, but it should be easy to check)? > glm(1 ~ 2,

Although the 'spatial' documentation doesn't mention that extractAIC works, it does seem to give an output. I may have misunderstood, but shouldn't the following give at least the same d.f.? > library(spatial) > data(topo, package="MASS") > extractAIC(surf.ls(2, topo)) [1] 46.0000 437.5059 > extractAIC(lm(z ~ x+I(x^2)+y+I(y^2)+x:y, topo)) [1]

Hi, I have an R program which takes several days to run, and sometimes hangs while running, presumably stuck in some sort of loop within a package function. I don't know where in the program code it is hanging: it is likely to be within a routine that is iterated many hundreds of thousands of times, each time with different random numbers used. I'd like to debug the R code, but

Hi, Using library(lattice), is there any way to tell xyplot to plot panels top to bottom, then left to right (i.e. panels are appended vertically, then horizontally). as.table changes the plot direction from left-to-right then top-to-bottom, to right-to-left then bottom- to-top, but that's not quite what I want to do. Thanks Yan

Hello, Is this a bug in the lmer routine? > library(lme4) > ### test case based on rats data from Crawley > a<-rnorm(36);b<-rep(1:3,each=12);c<-rep(1:2,each=6,3);d<-rep (1:3,each=2,6) > > ### mixed model works when c & d are numeric, lmer assumes they are factors > m <- lmer(a ~ b + (1|c/d)) > > ### but bails out when they are actually

Full_Name: Yan Wong Version: 2.8.1 OS: Mac OS X 10.4 Submission from: (NULL) (78.149.167.246) When printing a plot to and svg file, text intended to be in italics is rendered in bold (svg file viewed using both Safari 3.2.1 and Firefox 2.0.0.9. For example, try viewing the test.svg file generated by svg("test.svg", width=4,4,8) plot.new();

Full_Name: Yan Wong Version: 2.8.0 OS: Mac OS X 10.4 Submission from: (NULL) (78.149.183.231) When plotting multiple time series in a single plot, via plot.ts(plot.type="multiple"), the cex.lab, col.lab, and font.lab arguments are ignored > plot(ts(data.frame(a=1:10, b=1:10)), plot.type="single", cex.lab=0.5, col.lab="red") #tiny red axis labels >

Hi, just a quick question: Should weighted.mean be able to cope with the specific case where one weight is Inf? I came across this when trying to implement a simple weighted moving average algorithm for surface smoothing: these algorithms often result in a single infinite weight when predicting the surface value at known data points. e.g. > weighted.mean(c(77,88,99), c(Inf, 1, 2))

Hi, If I used drop1 in a weighted lm fit, it seems to ignore the weights in the AIC calculation of the dropped terms, see the example below. Can this be right? Yan -------------------- library(car) > unweighted.model <- lm(trSex ~ (river+length +depth)^2- length:depth, dno2) > Anova(unweighted.model) Anova Table (Type II tests) Response: trSex Sum Sq Df F value

----- Original Message ----- From: MARK LEEDS To: r-help-request@stat.math.ethz.ch Sent: Sunday, May 14, 2006 9:49 PM Subject: simple .bat question but i can't find out where to post. my apologiesfor posting a non R question I have what I think is a simple .bat question but i can't figure out where to post it ( i spent 6 hours looking through tutorials on batch commands but i

I've also tried adding the pgsql and mysql RPM's and they're not available either. On 20/03/2023 12:58, Gary Stainburn wrote: > Apologies.? This is the correct screen grab. > > The extra errors in the OP were because I had been experimenting, to > try to fix the issue. > > [root at testsvr ~]# ./ldapAuth.php gary.stainburn fake-password > PHP Fatal error:?

I'm sorry if people are getting fed up with my question, but I didn't get a response from the first post, and I didn't even see the second post appear in the summary. Hi All, this is probably a common problem and FAQ but I can't the solution. I've got two subnets, 10.1.1.0/255.255.0.0 at Leeds and 10.2.0.0/255.255.0.0 at Doncaster. At Leeds, I've got samba running as

Hi there, Does anyone know of a statistical test for heteroscedasticity for an object of class "gls"? (or alternative objective methods). Thanks in advance, Ben Gillespie, Research Postgraduate o-------------------------------------------------------------------o School of Geography, University of Leeds, Leeds, LS2 9JT o-------------------------------o http://www.geog.leeds.ac.uk/

Hi, Just a quick query - if I'm creating a function to produce a number of histograms per page of output (one per column from a matrix), how can I pass the column name of the matrix into the title (or indeed to form part of the x-axis label)? TIA, Laura Laura Quinn Institute of Atmospheric Science School of Earth and Environment University of Leeds Leeds LS2 9JT tel: +44 113 343 1596 fax:

Dear R users I have built the following model m1<-lmer(y~harn+foodn+(1|ass%in%pop%in%fam),family = "quasibinomial") where y<-cbind(alive,dead) where harn and foodn are categorical factors and the random effect is a nested term to represent experimental structure e.g. Day/Block/Replicate ass= 5 level factor, pop= 2 populations per treatment factor in each assay, 7 reps

Ok, I have reinstalled R-1.9.0 and this appears to have fixed the problems I was having with png(). However, I have a further question regarding png() Is it possible to pass a par() argument to the png() command? I am wanting to produce 4 plots per object, which I normally acheive on an X window by par(mfrow=c(1,4)). I have tried calling a new plot and setting par in this way but this has no

Hi, I was wondering if anyone could tell me why prcomp() will "Invent" modes of variation in a PCA on identical replicates of data? I would have expected 50 (or whatever number) of identical replicates to return a null score in such an analysis (or at the least, all variables would share the same PC score). This is not the case and I was wondering could someone point me in the direction

Hi, I have a table with data, as below: dput(table): structure(list(Adrian.Cole = c(0L, 0L, 0L, 0L, 0L, 0L), Alison.Wong = c(0L, 0L, 0L, 0L, 0L, 0L), Andrei.Savu = c(0L, 0L, 0L, 0L, 0L, 0L), Bruno.Dumon = c(0L, 0L, 0L, 0L, 0L, 0L), Edward.J..Yoon = c(0L, 0L, 0L, 0L, 0L, 0L), Eugene.Koontz = c(0L, 0L, 0L, 0L, 0L, 0L), Jakob.Homan = c(0L, 0L, 0L, 0L, 0L, 0L), Kelvin.Kakugawa = c(0L, 0L,

Hi everyone Can anyone tell me how I normalise a power spectral density (PSD) plot of a periodical time-series. At present I get the graphical output of spectrum VS frequency. What I want to acheive is period VS spectrum? Are these the same things but the x-axis scale needs transformed ? Any help would be greatly appreciated Tom