similar to: some related problems

Greetings, Quick Summary: How does R function non-interactively? I can't print graphics to an X11 window. I set up what appears to be called a "noninteractive" R session using the Perl code at the end of this message. Everthing seems to work fine as written. I get a nice postscript file with my graphics. However, I would like to display X11 R-graphics for the Perl-script

I noticed some odd behavior when I transpose a data frame containing NA's. It seems to cast all the elements as "character" including the NA's. Bug? > t(data.frame(x=1:10,y=1:10,z=rep(NA,10))) 1 2 3 4 5 6 7 8 9 10 x " 1" " 2" " 3" " 4" " 5" " 6"

Greetings, I'm having trouble finding a simple way to calculate a weighted histogram where there may be zero raw counts in a given interval. Given equal-length vectors of data 'data' and weights 'w', and breaks (intervals) for the histogram, I calculate a weighted histogram as follows (see MASS's 'truehist' for an unweighted histogram): bin <- cut(data,

Dear all, I'm a beginner in spatial regression, therefore, excuse me in advance for the stupid question. I'm using the function Tps of field library. As I read in help, the assumed model for this regression is additive Y = f(X) + e where f(X) is a d dimensional surface. I have understood that d in not an argument to specify in Tps(x, Y, ...). How can I obtain d ? Thanks a lot.

Hi! Like to write quite a lot (ca 100) of objects from my envirovment with save. The names of the objects are in a list nam. nam<-dir() nam<-grep("E",nam,value=T) length(nam) 20 for(x in nam) { #reads the objects and assigns the names. assign(x,simFromEmboss(Simmatrix(),x)) nnam<-paste(x,".rda",sep="") print(nnam) save(x,file=nnam) } I knew that

Dear developers, I wonder would you consider making a save to a graphics file format (as opposed to ps)? What prompts this is that we have just finished my wife's thesis using R heavily for stats and graphs. The combination of latex, bibtex and R generated .eps worked a treat and we were very pleased with both the final outcome and the efficiency of gernerating it, especially when all the

Hi, I have finally installed TeachingDemos trough zip file and installed on lbrary. One new question: I want to draw a line (like abline) on diagonal or the result of running a ols. How can I do it. On the other hand, I want to mark the length on each sub plot y axis. It is possible? My function is sca<-function(){ z2<-read.table("Ase.txt",header=T)

I know this topic has had plenty of discussion in the last couple of days, but.... I've been trying to compare the effects of different fitted methods for systems of equations (OLS, SUR, 2SLS, 3SLS ) and would like to compare the residual plots (easy) and the qqnorm/qqplot of the fits for the different fitted methdos. For example, qqnorm( residuals( lm( q ~ p + f + a ) ) ) par( new = TRUE )

Greetings, Summary: Has anyone else started a package (or have code) with the capability to write Vis5D-readable netCDF files containing irregularly-spaced data? I've recently put together what might be the beginnings of a (GPL'ed) R interface package to the Vis5D volumetric visualization program (http://vis5d.sourceforge.net). Vis5D is commonly used in the atmospheric science and

I want to qqnorm every column in a table. When I try the first column using qqnorm(table$column1), it worked. But when I use qqnorm(table[1]), it tells me "Error in stripchart(x1, ...) : invalid plotting method". What happen? How can I make a function that qqnorms every column? thanks a lot. -jiong The email message (and any attachments) is for the sole use of the intended

Hi all, I have just started using R to produce qqnorm plots. I am trying to switch the x and y axes so that the theoretical values are plotted on the y axis and my data on the x axis. Can anyone help me with this? Thanks -- View this message in context: http://www.nabble.com/Flipping-axes-of-qqnorm-tp23248007p23248007.html Sent from the R help mailing list archive at Nabble.com.

# I can not draw a 45 degree line on a qqnorm() plot, jj <- sample(c(1:100), 10) qqnorm(jj) abline() don't work. Thank you.

Does anyone know how to extract residuals in lmer? Here's the error I get: > crop.lme=lmer(response~variety*irrigation*pesticide+(1|rep)+(1|rep: pesticide)+(1|rep:pesticide:irrigation), crop.data) > qqnorm(crop.lme) Error in qqnorm.default(crop.lme) : y is empty or has only NAs > resid(crop.lme) NULL Thanks! --Jake

Stock returns and other financial data have often found to be heavy-tailed. Even Cauchy distributions (without even a first absolute moment) have been entertained as models. Your qq function subtracts numbers on the scale of a normal (0,1) distribution from the input data. When the input data are scaled so that they are insignificant compared to 1, say, then you get essentially the

Referring to "Using R for Data Analysis and Graphics" by J H Maindonald, and available from the R site, I found the example on p.30 non-working: > stem(qqnorm(possum$hdlngth)) Error in stem(qqnorm(possum$hdlngth)) : 'x' must be numeric Since qqnorm(possum$hdlngth) plots, and > class(possum$hdlngth) [1] "numeric" , the problem must be here: >

Hello R friends, I'm wondering why I get funny qqnorm() results. It seems that they should all be reflected in the normal qqline(). For instance: if I qqnorm() bimodal or uniform data I get a sigmoidal in which the qqnorm() points lie above the qqline() at -ve theoretical quantiles, and the qqnorm() points lie below the qqline() at +ve theoretical quantiles. Yet I expect such platykurtic

Dear list-reader, by running the following script: library(zoo) sessionInfo() search() packageDescription("zoo") data(EuStockMarkets) dax <- as.zoo(EuStockMarkets[1:10, "DAX"]) daxr <- diff(log(dax)) identical(as.vector(qnorm(daxr)), qnorm(coredata(daxr))) qqnorm(coredata(daxr)) qqnorm(daxr) qqnorm() produces an error: > qqnorm(daxr) Fehler in if (xi == xj) 0L

Hello, I am trying to do some plotting to check random effect assumptions for a model I fit using lme. I want to use qqnorm and pairs (similarly to examples given in Pinheiro & Bates p. 188), but it's not working. Here's some relevant code and the error message: library(nlme) data(Machines) m1 <- lme(fixed=score~Machine,random=~1|Worker/Machine, data=Machines) qqnorm(m1,

First my excuses if I keep bugging everyone in this list, but I am a newbie, and tend to find some behaviour that looks unexpected to me; and I would really appreciate to be pointed to some location that allows me to understand more about this software. Here is my next question: > par(mfrow=c(1,1)) > qqnorm(rnorm(20)) > qqmath(rnorm(20)) > par(mfrow=c(3,4)) > for(i in

Hi all, I'm modeling using lme in the nlme package. qqnorm makes plots just find, but when I try to add a line with qqline, I get the following error: "Error in sort.list(x, partial = unique(c(lo, hi))) : 'x' must be atomic for 'sort.list' Have you called 'sort' on a list?" For example > modelincr10<-lme(lruin~can.pos.num, data=rwushi30, random=(~1|