similar to: can't get summary() in loop

I use AfterStep on Linux. I made an R.xpm file (attached) and put this icon in the "wharf". To invoke R I click on the R icon. I did not find an offical R.xpm, and so I am submitting this one for your collective consideration.... The idea behind the icon is that it looks like a scatterplot, which is "statistical". My main beefs with my R.xpm are that it doesn't have a

Full_Name: Version: 1.0.0 OS: linux Submission from: (NULL) (193.62.250.209) ?par lists and describes xlog and ylog, though these are obsolete. It does not list log. I found out about log from ?plot.defaults--so that help is OK. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

I was wondering if anyone out there has created a series of images (e.g. jpegs) using R, then animated them on a web page. I have looked around a bit and it seems fairly complicated. One approach that seems sensible is to use a java script that displays the list of images. I have only used html on web pages before so I don't know how to use java... Ideally the final result is easily ported to

I have not really used lm before and I was hoping for some help on a simple problem. Here is a toy version of the problem I want to solve. y x grp -.9 1 a -.8 2 a -.7 3 a -.7 1.5 b -.5 2.5 b -.3 3.5 b -.19 2.7 c -.11 3.7 c -.41 4.7 c I want to fit a model that has one y-intercept and three slopes, one for

Full_Name: Bill Simpson Version: 65.1 OS: linux Submission from: (NULL) (193.62.250.209) Try this code as separate lines entered interactively, then try doit() doit<-function() { x<-seq(1,10) y1<-x+rnorm(10,1,1.5) y2<-x+rnorm(10,1,1) t.test(x,y1,paired=TRUE) t.test(x,y2,paired=TRUE) } doit() apparently executes only the last of a series of t.test()s. Maybe this is no bug,

There were some questions about hist() a couple of days ago which triggered this post. My question/suggestion is about the y-axis in hist. There are reasons to prefer making the y-axis density=relative frequency/bin width. One reason is that the height of the plot does not depend on the bin width; another is that if your histogram is in density then you can easily superimpose a smooth theoretical

Full_Name: Bill Simpson Version: 0.64.1 OS: linux Submission from: (NULL) (193.62.250.209) Two problems with dev.print help. 1. Under Examples: dev.print(width=6, height=6, horizontal=F) ^ FALSE 2. dev.print help does not describe any of its parameters (e.g. height or width). Not described on the other dev.xxx pages either

I just installed 99.0a, and I can't set mar and mgp: > par(mar(2,2,.1,.1)) Error in par(mar(2, 2, .1, .1)) : couldn't find function "mar" I had the same happen with mgp. Are these bugs or a problem in my installation (everything else seems OK)? Is R core going to kill mar and mgp? If so, how should I fix the huge margins and huge space between axis title and axis? Try this

Full_Name: Bill Simpson Version: 65.1 OS: Linux Redhat 6.1 Submission from: (NULL) (193.62.250.209) Try this: #default plot symbols and lettering are too small, need to scale up par(cex=2,mex=.7) par(mar=c(5,5,1,1)) spdiff<-c(1,2,3,4,5) dpdet<-c(1,2,3,4,5)/10 dpsp<-c(2,3,4,5,6)/10 dpdir<-c(2,4,6,8,10)/10 plot(spdiff,dpdet,pch=15,ylim=c(0,max(dpdet,dpsp,dpdir)),xlab="Speed

Full_Name: Bill Simpson Version: 1.2.3 OS: linux Submission from: (NULL) (193.62.250.209) I think I have found an accuracy problem with sum() > a<-1:1000 > b<-1:1000 > > crossprod(a,a)*crossprod(b,b) [,1] [1,] 1.114448e+17 > drop(.Last.value) [1] 1.114448e+17 > #this should give same answer as above > sum(a*a)*sum(b*b) [1] -652010736 #seems to be due to

I want to display a series of plots: display plot 1 wait for user input (keypress or mouse click) display plot 2 etc How to do it? Thanks for any help. (I have been using locator().) Bill Simpson -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

I have a data file like this: std cf hit miss fa cr 920 980 40 15 14 31 950 1010 24 23 23 30 1190 1250 26 21 27 26 1010 1070 33 10 28 29 1040 1100 35 10 11 44 I use read.table to read it in. My problem is that read.table makes std and cf into factors. I want them just

If I source in the function (see below) calib(), I get: > source("papers/helle/threshold.r") > calib() Error in eval(expr, envir, enclos) : Object "energy" not found But if I cut and paste the code for calib() one line at a time into the R window it works fine. calib<-function() { contrast<-c(.01,.02,.0325,.055,.0775,.1,.125,.15,.175,.2)

I fit the model fit<-lm(thresh~cond*Ne) where thresh is the reponse cond is a factor with levels a, b, and c Ne is a continuous indep var I think of this full model as having three lines: thresh as a function of Ne for each condition. Thus we have slopea, slopeb, slopec, inta, intb, intc. lm output my params ------------------------- (Intercept) inta condb intb - inta condc

A while ago there was a thread on graphics. I complained that the plots produced by dev.print had symbols and fonts that were too small. I said I had been using par(cex=2), but the advice was to use pointsize. I have tried ps, is that what was meant? par(ps=40) x<-rnorm(100) y<-rnorm(100) dev.print(file="~/junk.ps", horizontal=FALSE) No matter how big I make ps (40,60...), the

This is probably a trivial question, but I did not see it in the FAQs. Perhaps I simply need to be pointed to some URLs. In any case, I have some sampled data that has peaks, and I would like to extract peak height, location, and width, according to almost any statistical model. Gaussian would seem to make the most sense. There can also be "steps" increasing the amplitude over a

I was just playing around with the R-streams package, reading in some data from an old experiment we did here, and was again amazed by the flexibility of the R environment. The eaze with wich I was able to generate simple plots of my EEG data made me wonder if someone tried to implement an ERP package. For the non-cognitive scientists ou there: EEG is electrical activity measured from the scalp

I have some Bernoulli data something like this: x<-sort(runif(100,1,20)) p<-pnorm(x,10,3) y<-as.numeric(runif(x)<p) plot(x,y) lines(x,p) This plot is not very satisfactory because the ogive does not visually fit the (0,1) points very well, and also because the points tend to fall on top of one another. The second problem can be eliminated by adding vertical jitter. However I was

This example, x<-1:10 y1<-x y2<-x+2 plot(x,y1,pch=21,ylim=range(c(y1,y2)),type="o",lty="solid",bg="white") points(x,y2,pch=22,type="o",lty="dashed",bg="white") legend(7,3,legend=c("y1","y2"),pch=c(21,22),lty=c("solid","dashed"),bg="white") works if I use pch = one of 21-25. But

I wanted to produce a layout like this: |----------| | | | 1 | | | |-----------------------| | | | | 2 | 3 | | | | |-----------------------| Three equal sized plots with the top one centred about the other two. So far I am using nf<-layout(matrix(c(1,1,2,3), 2, 2,