Re: [R] spplot cuts

2010-09-29 Thread bbolker
Hmmm. Maybe a documentation typo in ?spplot. If you follow the documentation through to ?levelplot, you find that cuts: number of levels the range of ‘z’ would be divided into (no mention of actual breakpoints) but: at: numeric vector giving breakpoints along the range of ‘z’.

Re: [R] Fry Plots

2010-07-29 Thread bbolker
[cc'ing back to r-help] At this point it sounds like you may need to poke around some more on your own -- I'm guessing you're not familiar with R. Reading through the various documentation on these packages (try help(package=...) to see if there is a vignette too) and/or the "Introduction to

Re: [R] 3-dimensional graph

2007-10-01 Thread bbolker
John Sorkin wrote: > > Windows XP > R 2.3.1 > > I have a funciton > fit1<-lm(y~x+z) > Is there a function that will produce a 3-dimensional plot of y,x,z? > > I looked at the help files, but did not find a clean answer to my questio > Check out ?scatter3d in the Rcmdr package. Ben Bolker

Re: [R] resetting par() to all defaults: par(reset=TRUE) ?

2007-09-30 Thread bbolker
Michael Friendly wrote: > > In a long session, producing multiple graphs, I sometimes repeatedly > change par() settings, particularly with multi-row/col displays. > If I'm using a script, I'll do > > op <- par(newsettings) > ... plots ... > par(op) > > but sometimes I do things on the fly an

Re: [R] kurtosis

2007-09-28 Thread bbolker
S Bina wrote: > > Hi, > > I cannot find the function kurtosis. Is it sth additional I am meant to > download? I use the MacOS X version of R. > > Many thanks > Samira > > __ > R-help@r-project.org mailing list > https://stat.ethz.ch/mailman/listin

Re: [R] add points to wireframe

2007-09-26 Thread bbolker
marcg wrote: > > Hello > > with follwing code I want to produce the wireframe and add some points on > it: > > wireframe(z ~ x * y, data=dem, aspect = c(1, .5), > scales = list(arrows = FALSE), > panel.3d.wireframe = function(x, y, z,...) { > panel.3dwire(x =

Re: [R] R2Winbugs problem

2007-09-25 Thread bbolker
Sam_zhz wrote: > > Just while I use R2Winbugs, the following error will be presented. Please > tell me how to cope with. > Error in file(file, "r") : unable to open connection > In addition: Warning message: > cannot open file 'codaIndex.txt', reason 'No such file or directory' in: > file(file,

Re: [R] 3d barplot in rgl

2007-09-25 Thread bbolker
hadley wrote: > > On 9/25/07, Chris Stubben <[EMAIL PROTECTED]> wrote: >> hadley wickham gmail.com> writes: >> >> > >> > Why do you want a 3d barchart? They are generally a bad way to >> > present information as tall bars can obscure short bars, and it is >> > hard to accurately read off the

Re: [R] Root finding problem

2007-09-24 Thread bbolker
tuki <- function(u, x, a, lambda){ u <- u+0i f <- Re(x-(a*(u)^lambda-(1-(u))^lambda)/lambda) f } ## What I want to do is to find the root, but without specifying the ## interval within which to search for it. I can do it easily in MATLAB ## with fsolve() or fzero() functions. ## BB: how d

Re: [R] calculating/plotting error ellipses

2007-09-24 Thread bbolker
Jan M. Wiener wrote: > > hello, > sorry for posting what may be a simple question: > i do have a matrix of coordinates (positional judgments, see below) and > now want to calculate and plot the corresponding error ellipse. > can anyone help me with the exact steps/syntax? > Something along th

Re: [R] save 3dplot to file

2007-09-24 Thread bbolker
Birgit Lemcke wrote: > > I think it is an easy question but i haven´t found a solution since I > am r beginner. > > I did this > > plot3d(PCoA, type="p", col=rainbow(1000),size=5) > text3d(PCoA, text=Nam) > > > How can I save it including the labels to a jpg-file? > You can

Re: [R] Ambiguities in vector

2007-09-20 Thread bbolker
Birgit Lemcke wrote: > > > Perhaps you haven´t understood my question in the mail yesterday. So > I will try to describe my problem in a different way > > You see the tables. I would like to test the variables between the > tables. > I'm afraid that even before we start to deal with th

[R] R-related: problem with openvt on Ubuntu 7.04

2007-09-19 Thread bbolker
This is a configuration/OS problem, but it's affecting my use of R ... Whenever I try to open _any_ vignette (as far as I can tell) from within R, I get "Could not get a file descriptor referring to the console", coming from /usr/bin/openvt, which is pointed to by /usr/bin/open, which is cal

Re: [R] Create correlated data with skew

2007-09-18 Thread bbolker
Mike Lawrence wrote: > > Hi all, > > I understand that it is simple to create data with a specific > correlation (say, .5) using mvrnorm from the MASS library: > > > library(MASS) > > set.seed(1) > > > > a=mvrnorm( > + n=10 > + ,mu=rep(0,2) > + ,Sigma=matrix(c(1,.5,.5,1),2,2

Re: [R] Cannot get contrasts to work with aov.

2007-09-17 Thread bbolker
ThatDeadDude wrote: > > I have been trying for hours now to perform an orthogonal contrast through > an ANOVA in R. > > I have done a two-factor factorial experiment, each factor having three > levels. I converted this dataset to a dataframe with one factor with nine > treatments, as I cou

Re: [R] are hurdle logit-poisson model and posson model nested?

2007-09-17 Thread bbolker
Wensui Liu wrote: > > Dear Listers, > I have a general statistical question. Are hurdle logit-poisson model > and posson model nested? > You might have to give us a little more detail. On first glance, my impression is that the Poisson model is _not_ nested in the hurdle-Poisson, because t

Re: [R] Problem with nlm() function.

2007-09-17 Thread bbolker
Rolf Turner-3 wrote: > > > In the course of revising a paper I have had occasion to attempt to > maximize a rather > complicated log likelihood using the function nlm(). This is at the > demand of a referee > who claims that this will work better than my proposed use of a home- > grown i

Re: [R] Must be easy, but haven't found the function (numerical integration)

2007-09-17 Thread bbolker
On 9/17/07, Ptit_Bleu <[EMAIL PROTECTED]> wrote: > I have a data frame of 2 columns with the following types : > data$day char > data$value num > > And I plot my data with : > plot(strptime(donnees$day,format="%Y-%m-%d %H:%M:%S"),donnees$value, > type="l") > > And I'd just like to get the numeri