[R] Errors using update.packages()

2006-02-11 Thread David Scott
When trying to update packages after the sysadmin updated R on my unix box I got errors on some packages. For example for chron: * DONE (chron) mkdir: cannot create directory `/usr/lib/R/library/00LOCK': Permission denied ERROR: failed to lock directory '/usr/lib/R/library' for modifying I

Re: [R] Transferring R results to word prosessors

2006-02-11 Thread Philippe Grosjean
The copy to clipboard feature is written (by Eric Lecoutre myself) since a couple of years. It is in the SciViews bundle, library svIO. In this library, you have: - copy() that copies an object to the clipboard in various formats (icluding HTML, by using R2HTML and LaTeX), - export() does

Re: [R] Need frequency distribution for x,y coordinates

2006-02-11 Thread Gabor Grothendieck
Check out: http://addictedtor.free.fr/graphiques/graphcode.php?graph=116 On 2/11/06, mark shanks [EMAIL PROTECTED] wrote: Hi, I have a set of data in x,y coordinates across the range of -5 to 5 in each dimension. I would like to obtain the frequency distribution of the different points, and

Re: [R] Errors using update.packages()

2006-02-11 Thread Peter Dalgaard
David Scott [EMAIL PROTECTED] writes: When trying to update packages after the sysadmin updated R on my unix box I got errors on some packages. For example for chron: * DONE (chron) mkdir: cannot create directory `/usr/lib/R/library/00LOCK': Permission denied ... From the startup:

Re: [R] Need frequency distribution for x,y coordinates

2006-02-11 Thread Roger Bivand
On Sat, 11 Feb 2006, mark shanks wrote: Hi, I have a set of data in x,y coordinates across the range of -5 to 5 in each dimension. I would like to obtain the frequency distribution of the different points, and then graph them so you can see which of the points are the most frequently

Re: [R] Need frequency distribution for x,y coordinates

2006-02-11 Thread Prof Brian Ripley
On Sat, 11 Feb 2006, Roger Bivand wrote: On Sat, 11 Feb 2006, mark shanks wrote: I have a set of data in x,y coordinates across the range of -5 to 5 in each dimension. I would like to obtain the frequency distribution of the different points, and then graph them so you can see which of the

[R] hclust(stats) merge matrix interpretation

2006-02-11 Thread Tarun Kumar Singh
Hi, We are trying to interpret the clusters generated by hclust method of R stats package. The problem here is when i get the hc$order then there is some order, while exporting to file that order is lost. Here is the example code and their results: hc - hclust(dist(USArrests), ave) plot(hc)

[R] tk mouse cursor icon widget tkwinfo tkfocus questions

2006-02-11 Thread Frank Samuelson
1. I want to change the mouse cursor over my window into a wait/watch icon while R computes. Can this be done directly? Some ancient tcltk mailing lists said change the cursor over every widget in the window: foreach widget [winfo children $window] {

[R] R-newbie-question, fixed effects panel model, large number of observations

2006-02-11 Thread Thomas Wilde
Hi, I'm trying to fit a fixed effect (LSDV) panelmodel with R. I have a dataset with y as dependent, x1x2 as indeps, t as time index and i as an id-variable for each individual. There are three observations for each individual (t=1, t=2, t=3). I want to try a simple regression, but with

Re: [R] lme4: Error in getResponseFormula(form) : Form must be a two sided formula

2006-02-11 Thread Douglas Bates
On 2/10/06, Spencer Graves [EMAIL PROTECTED] wrote: Hi, Doug and Chris: I just got the same error message with the lmList example in lme4: (fm1 - lmList(breaks ~ wool | tension, warpbreaks)) Call: Error in getResponseFormula(form) : Form must be a two sided formula

Re: [R] R-newbie-question, fixed effects panel model, large number of observations

2006-02-11 Thread ronggui
This is one function I wrote. http://sociology.yculblog.com/post.794856.html 2006/2/12, Thomas Wilde [EMAIL PROTECTED]: Hi, I'm trying to fit a fixed effect (LSDV) panelmodel with R. I have a dataset with y as dependent, x1x2 as indeps, t as time index and i as an id-variable for each

Re: [R] Need frequency distribution for x,y coordinates

2006-02-11 Thread nojhan
Le Sat, 11 Feb 2006 04:40:41 -0500, Gabor Grothendieck a écrit : http://addictedtor.free.fr/graphiques/graphcode.php?graph=116 Or hist2d, in the gplots package : http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=70 -- NoJhan __

[R] plots

2006-02-11 Thread statistical . model
Hi all, I have a basic question. how can i visualize two or more density curves on the same plot? ex: x1-runif(100,10,80) x2-runif(100,1,100) kernelgraf-density(x1,kernel = gaussian, width= 20) plot(kernelgraf, xlab=Probability, xlim=c(0,100), ylim=c(0,.1), col=rgb(0,1,0), main=)

Re: [R] plots

2006-02-11 Thread Romain Francois
Le 11.02.2006 19:44, [EMAIL PROTECTED] a écrit : Hi all, I have a basic question. how can i visualize two or more density curves on the same plot? ex: x1-runif(100,10,80) x2-runif(100,1,100) kernelgraf-density(x1,kernel = gaussian, width= 20) plot(kernelgraf, xlab=Probability,

Re: [R] plots

2006-02-11 Thread Charles Annis, P.E.
You could substitute lines(kernelgraf) for your last line: points(kernelgraf, xlab=Probability, xlim=c(0,100), ylim=c(0,.1), col=rgb(0,0,1), main=) See ?lines Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com

Re: [R] plots

2006-02-11 Thread Duncan Murdoch
Romain Francois wrote: Le 11.02.2006 19:44, [EMAIL PROTECTED] a écrit : Hi all, I have a basic question. how can i visualize two or more density curves on the same plot? ex: x1-runif(100,10,80) x2-runif(100,1,100) kernelgraf-density(x1,kernel = gaussian, width= 20) plot(kernelgraf,

Re: [R] plots

2006-02-11 Thread Romain Francois
Le 11.02.2006 19:59, Duncan Murdoch a écrit : Romain Francois wrote: Le 11.02.2006 19:44, [EMAIL PROTECTED] a écrit : Hi all, I have a basic question. how can i visualize two or more density curves on the same plot? ex: x1-runif(100,10,80) x2-runif(100,1,100)

Re: [R] srt --- slope text with function?

2006-02-11 Thread Patricia J. Hawkins
iw == ivo welch [EMAIL PROTECTED] writes: iw * SUGGESTION: can we please offer the ?: operator ala C in iw addition to ifelse()? This would make R code prettier. R: if (condition) xxx else yyy if (condition) xxx else if (yyy) zzz else qqq C: condition ? xxx : yyy; condition ? xxx : yyy ? zzz

[R] how do I relate tables in R?

2006-02-11 Thread Hans Gardfjell
?merge -- * Hans Gardfjell Ecology and Environmental Science Umeå University 90187 Umeå, Sweden email: [EMAIL PROTECTED] phone: +46 907865267 mobile: +46 705984464 __ R-help@stat.math.ethz.ch mailing list

[R] aggregate vs tapply; is there a middle ground?

2006-02-11 Thread Joseph LeBouton
Dear all, I'm wanting to do a series of comparisons among 4 categorical variables: a - aggregate(y, list(var1, var2, var3, var4), sum) This gets me a very nice 2-dimensional data frame with one column per variable, BUT, as help for aggregate says, empty subsets are removed. I don't see in

[R] aggregate vs tapply; is there a middle ground?

2006-02-11 Thread Hans Gardfjell
I faced a similar problem. Here's what I did tmp - data.frame(A=sample(LETTERS[1:5],10,replace=T),B=sample(letters[1:5],10,replace=T),C=rnorm(10)) tmp1 - with(tmp,aggregate(C,list(A=A,B=B),sum)) tmp2 - expand.grid(A=sort(unique(tmp$A)),B=sort(unique(tmp$B))) merge(tmp2,tmp1,all.x=T) At least

Re: [R] aggregate vs tapply; is there a middle ground?

2006-02-11 Thread hadley wickham
I faced a similar problem. Here's what I did tmp - data.frame(A=sample(LETTERS[1:5],10,replace=T),B=sample(letters[1:5],10,replace=T),C=rnorm(10)) tmp1 - with(tmp,aggregate(C,list(A=A,B=B),sum)) tmp2 - expand.grid(A=sort(unique(tmp$A)),B=sort(unique(tmp$B))) merge(tmp2,tmp1,all.x=T) At

Re: [R] aggregate vs tapply; is there a middle ground?

2006-02-11 Thread Joseph LeBouton
Thanks, Phil! I've literally spent two hours on my own trying to find something that does exactly that. Thanks for another pair of functions added to my (slowly!) growing R vocabulary. -jlb Phil Spector wrote: Joseph - I'm sure there are clearer and more efficient ways to do it, but

[R] failed installing hdf5 package--can't find zlib

2006-02-11 Thread Zepu Zhang
I use Mac, tried to install the package hdf5-1.6. Its configuration script can't find zlib, which apparently exists in my /usr directory. I used R CMD INSTALL --configure-args=--with-hdf5=/sw --configure-args=--with-zlib=/usr hdf5tar.gz It still can't find zlib. Is this a bug? Any pointer

Re: [R] srt --- slope text with function?

2006-02-11 Thread Gabor Grothendieck
On 2/11/06, Patricia J. Hawkins [EMAIL PROTECTED] wrote: iw == ivo welch [EMAIL PROTECTED] writes: iw * SUGGESTION: can we please offer the ?: operator ala C in iw addition to ifelse()? This would make R code prettier. R: if (condition) xxx else yyy if (condition) xxx else if (yyy) zzz

Re: [R] aggregate vs tapply; is there a middle ground?

2006-02-11 Thread Peter Dalgaard
hadley wickham [EMAIL PROTECTED] writes: I faced a similar problem. Here's what I did tmp - data.frame(A=sample(LETTERS[1:5],10,replace=T),B=sample(letters[1:5],10,replace=T),C=rnorm(10)) tmp1 - with(tmp,aggregate(C,list(A=A,B=B),sum)) tmp2 -

[R] Mathematical typesetting of column heads using the latex (Hmisc) function

2006-02-11 Thread Michael Kubovy
Dear r-helpers, I would very much appreciate help with the following problem: The following command (in a .Rnw file) latex(anova(e7.lmer3, e7.lmer4), file = 'e7lmer34.tex', rowname = c ('nonlinear', 'linear'), longtable = FALSE, dcolumn = T, booktabs = T, table.env = F) produces the

[R] lme, nlsList, nlsList.selfStart

2006-02-11 Thread Patrick Giraudoux
Dear listers, I am trying to fit a model using nlsList() using alternately a SSfol() selfstart function or its developped equivalent formulae. This preliminary trial works well mydata-groupedData(Conc~Tps|Organ,data=mydata) mymod1-nls(Conc~SSfol(Dose,Tps,lKe,lKa,lCl),data=mydata) as well as a