Re: [R] Modifying a built-in R function

2009-03-02 Thread Uwe Ligges
japal wrote: Hello, Something incredible (at least for me) has happen. Yesterday night I downloaded biplot.R to edit this function and add new features I wished. Namely I wanted to plot points belonging to different groups using different colors and symbols. I identified which part of the

Re: [R] nls does not accept start values

2009-03-02 Thread Uwe Ligges
Petr PIKAL wrote: Hi to all OK as I did not get any response and I really need some insight I try again with different subject line I have troubles with correct evaluating/structure of nls input Here is an example # data x -1:10 y -1/(.5-x)+rnorm(10)/100 # formula list form -

[R] Loop over several Garch-models

2009-03-02 Thread thomas.schwander
Dear all, I'm using R 2.8.1 on Windows XP. I want to loop over several garch-models to evaluate the best model-fit. After loading the package fGarch I execute the following loop: For(i in 1:5){ for(j in 1:5){ garchFit(~ arma(0,0) + garch(i,j), stetige_renditen[,7],

Re: [R] nls does not accept start values

2009-03-02 Thread Prof Brian Ripley
It is not very clear what you are trying to do here, and form - structure(list(a = list(quote(y ~ 1/(a - x)), list(a=mean(y, .Names = a) is using a historic anomaly (see the help page). I am gussing you want to give nls an object containing a formula and an expression for the starting

[R] ave and grouping

2009-03-02 Thread Patrick Hausmann
Dear list, # I have a DF like this: sleep$b - c(rep(8,10), rep(9,10)) sleep$me - with(sleep, ave(extra, group, FUN = mean)) sleep # I would like to create a new variable # holding the b-th value of group 1 and 2. # This is not what I want, it takes always the '8' from group '1' # and not

Re: [R] Does The Randvar package contain a virus(Malware) ?

2009-03-02 Thread Bernardo Rangel Tura
On Mon, 2009-03-02 at 08:41 +0100, Matthias Kohl wrote: there were only minor changes in the latest version of RandVar (from 0.6.6 to 0.6.7). Might this be a mirror problem? Best Matthias Timthy Chang wrote: Today,I update the packages in R. but AntiVir Guard dectects the Randvar

Re: [R] rounding problem

2009-03-02 Thread Bernardo Rangel Tura
On Sun, 2009-03-01 at 17:51 -0800, tedzzx wrote: Yes, round(1.5)=2. but round(2.5)=2. I want round(2.5)=3 just like the what the excel do. Can we change the setting or do some trick so that the computer will work like what we usually do with respect to rounding. My system is R 2.8.1, winXP,

[R] handle graph size in eps

2009-03-02 Thread Benoit Boulinguiez
Hi all, I've got a density graph made with the following commands: win.graph(width=13,height=6) par ( fin=c(13,3) ,mai=c(1,1,0.5,0.5) ,mfrow=c(1,2) ,cex.axis=1.5 ,cex.lab=1.5) dens-density(DATA1.y[2,]-mean(DATA1.y[2,]),kernel=gaussian) xlimit-range(dens$x) ylimit-range(dens$y) hist(

Re: [R] probleme with savePlot (to eps)

2009-03-02 Thread Christophe Genolini
I am using windows XP pack3 and R 2.8.1 so I guess this is savePlot for Windows. As Gabor say, savePlot offer to export graph also in eps. If you want to copy a screen device to postscript portably, use dev.copy2eps(). I am exporting graph in a package, so I don't want to copy screen device

Re: [R] nls does not accept start values

2009-03-02 Thread Petr PIKAL
Thank you It was simplified version of my problem. I want to elaborate a function which can take predefined list of formulas, some data and evaluate which formulas can fit the data. I was inspired by some article in Chemical engineering in which some guy used excel solver for such task. I was

Re: [R] Need help extracting info from XML file using XML package

2009-03-02 Thread Wacek Kusnierczyk
Don MacQueen wrote: I have an XML file that has within it the coordinates of some polygons that I would like to extract and use in R. The polygons are nested rather deeply. For example, I found by trial and error that I can extract the coordinates of one of them using functions from the XML

Re: [R] nls does not accept start values

2009-03-02 Thread Uwe Ligges
Petr PIKAL wrote: Thank you It was simplified version of my problem. I want to elaborate a function which can take predefined list of formulas, some data and evaluate which formulas can fit the data. I was inspired by some article in Chemical engineering in which some guy used excel solver

[R] Large data set in R

2009-03-02 Thread Hardi
Hello, I'm trying to use R statistical packages to do ANOVA analysis using aov() and lm(). I'm having a problem when I have a large data set for input data from Full Factorial Design Experiment with replications. R seems to store everything in the memory and it fails when memory is not enough

Re: [R] issue with varSel.svm.rfe in package MCRestimate

2009-03-02 Thread Elizabeth McClellan
Unfortunately no one ever answered my question and I was never able to implement the varSel.svm.rfe function in R. I found a 2005 book that mentions the RFE package: Analyzing Microarray Gene Expression Data, Chapter 6.13.4http://books.google.com/books?id=gt8JNQfpnMICpg=PA206lpg=PA206dq= but the

Re: [R] Modifying a built-in R function

2009-03-02 Thread japal
japal wrote: Hello, Something incredible (at least for me) has happen. Yesterday night I downloaded biplot.R to edit this function and add new features I wished. Namely I wanted to plot points belonging to different groups using different colors and symbols. I identified which part of

[R] xlsReadWrite package repository for Ubuntu

2009-03-02 Thread reverend33
Hi, I'm trying to install R on Ubuntu. I succeeded at installing the r-recommended package that is present in the synaptics, but i can't find the xlsReadWrite package in the repositories included in my synaptics manager. Does anybody know a liable repository in which this package is present.

Re: [R] nls does not accept start values

2009-03-02 Thread Petr PIKAL
Hi Prof Brian Ripley rip...@stats.ox.ac.uk napsal dne 02.03.2009 10:24:52: It is not very clear what you are trying to do here, and form - structure(list(a = list(quote(y ~ 1/(a - x)), list(a=mean(y, .Names = a) is using a historic anomaly (see the help page). I am gussing you

Re: [R] Inefficiency of SAS Programming

2009-03-02 Thread Gerard M. Keogh
Yes Greg, but if you're buying SAS they'll throw in IML pretty cheaply - SAS think it's only for a few nerds out there who wan to do funny stuff. G Greg Snow

Re: [R] handle graph size in eps

2009-03-02 Thread Prof Brian Ripley
On Mon, 2 Mar 2009, Benoit Boulinguiez wrote: Hi all, I've got a density graph made with the following commands: win.graph(width=13,height=6) The preferred name is windows(). par ( fin=c(13,3) ,mai=c(1,1,0.5,0.5) ,mfrow=c(1,2) ,cex.axis=1.5 ,cex.lab=1.5)

Re: [R] Adjusting confidence intervals for paired t-tests of multiple endpoints

2009-03-02 Thread Erich Studerus
Thanks, this comes close to what I need, but the problem is that I have repeated measures in each comparison and SimComp does not seem to handle this. In my experiments, each subject received both placebo and the drug in a randomized order. So, I want to do multiple paired t-test corrected for

[R] R-code help for filtering with for loop

2009-03-02 Thread John Antonydas Gaspar
Dear Sir / Madam, I am new for R coding. Kindly help me out in sorting out the following problem. There are 50 rows with six coloumns(you could see in the attached .txt file). I wish to go for filtering this 50 rows for any one of the six coloumns satisfying the value = 64. I need to have a

Re: [R] xlsReadWrite package repository for Ubuntu

2009-03-02 Thread Prof Brian Ripley
On Mon, 2 Mar 2009, reverend33 wrote: Hi, I'm trying to install R on Ubuntu. I succeeded at installing the r-recommended package that is present in the synaptics, but i can't find the xlsReadWrite package in the repositories included in my synaptics manager. Does anybody know a liable

[R] Can no longer use resize graphics X windows

2009-03-02 Thread Daniel Brewer
Hi, I use the ubuntu deb packages for R. In previous versions of R I could resize the graphics windows by dragging the bottom corner, but now I can't. Anyone got any idea why this has happened and whether it is possible to change it back? Dan --

[R] Distance between clusters

2009-03-02 Thread Corrado
Dear friends I reformulate the question. I think I did not formulate it properly. I have some data on some sites. I can define a dissimilarity between each pair of sites. Using this dissimilarity, I have clustered the sites using the hclust algorithm, with method ward. I then obtain 48

Re: [R] Modifying a built-in R function

2009-03-02 Thread Prof Brian Ripley
I am not sure what you intended by biplotes - function(x, ...) UseMethod(biplot) That does nothing different from biplot(). You need to call your modified functions 'biplotes.default' and 'biplotes.princomp' and call those via biplotes - function(x, ...) UseMethod(biplotes) Anything

Re: [R] xlsReadWrite package repository for Ubuntu

2009-03-02 Thread Uwe Ligges
reverend33 wrote: Hi, I'm trying to install R on Ubuntu. I succeeded at installing the r-recommended package that is present in the synaptics, but i can't find the xlsReadWrite package in the repositories included in my synaptics manager. Does anybody know a liable repository in which this

Re: [R] Using very large matrix

2009-03-02 Thread Corrado
Thanks a lot! Unfortunately, the R package I have to sue for my research was only released on 32 bit R on 32 bit MS Windows and only closed source I normally use 64 bit R on 64 bit Linux :) I tried to use the bigmemory in cran with 32 bit windows, but I had some serious problems.

Re: [R] R-code help for filtering with for loop

2009-03-02 Thread ONKELINX, Thierry
Dear John, It looks like you are stuck in both the second and the third circle of the R inferno (http://www.burns-stat.com/pages/Tutor/R_inferno.pdf) You problem is easy to vectorise. #could the number of columns = 64 in each row NumCols - rowSums(datax = 64) #select rows with at least one

Re: [R] Inefficiency of SAS Programming

2009-03-02 Thread Thomas Levine
R depends on all of those things to run, but you only have to use those programs through R. The software depends on these other tools, but the human doesn't have to switch interfaces. Tom! On Fri, Feb 27, 2009 at 9:22 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Fri, Feb 27, 2009

Re: [R] Need help extracting info from XML file using XML package

2009-03-02 Thread Romain Francois
Hi, You also might want to check R4X: # install.packages(R4X, repos=http://R-Forge.R-project.org;) require( R4X ) x - xml(http://code.google.com/apis/kml/documentation/KML_Samples.kml;) coords - x[Polygon///coordinates/# ] data - sapply( strsplit( coords, (,|\\s+) ), as.numeric ) Romain

[R] RES: object .trPaths not found

2009-03-02 Thread Leandro Marino
Did you do the configuration of Tinn-R after the installation? Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 R (21) 2103-9600 R.:236 ( (21) 8777-7907 ( lean...@cesgranrio.org.br Aquele que

[R] density 1?

2009-03-02 Thread Johannes Elias
Dear R-Gurus, I wonder why 'density' values as shown in hist or plot(density(x)) are sometimes over 1. How can that be? Example hist(rnorm(1000,sd=.5),freq=FALSE) The resulting plot shows density values below 1 on the y-axis. However, hist(rnorm(1000,sd=.1),freq=FALSE) shows density values

Re: [R] Need help extracting info from XML file using XML package

2009-03-02 Thread Romain Francois
Romain Francois wrote: Hi, You also might want to check R4X: # install.packages(R4X, repos=http://R-Forge.R-project.org;) require( R4X ) x - xml(http://code.google.com/apis/kml/documentation/KML_Samples.kml;) coords - x[Polygon///coordinates/# ] data - sapply( strsplit( coords, (,|\\s+) ),

Re: [R] density 1?

2009-03-02 Thread Bill.Venables
Because densities are not probabilities. It is the area under the density curve that represents probability. Example: the chi-squared density with 1 degree of freedom has a singularity at the zero and is unbounded. The area under the curve, however, is still 1. (This is a distressingly

Re: [R] density 1?

2009-03-02 Thread Eik Vettorazzi
Hi Johannes, ist more a statistical issue. In short: densities are not probabilities! With a continuous random variable probability statements are typically over intervals not over points. A density is bound to have an integral of 1 (and to be non-negative), nothing else. Consider the uniform

Re: [R] density 1?

2009-03-02 Thread Peter Dalgaard
Johannes Elias wrote: Dear R-Gurus, I wonder why 'density' values as shown in hist or plot(density(x)) are sometimes over 1. How can that be? Example hist(rnorm(1000,sd=.5),freq=FALSE) The resulting plot shows density values below 1 on the y-axis. However,

Re: [R] ave and grouping

2009-03-02 Thread Gabor Grothendieck
Try performing ave over the indexes rather than over extra itself: sleep$newcol - with(sleep, ave(1:nrow(sleep), group, FUN = function(ix) extra[b[ix][1]])) On Mon, Mar 2, 2009 at 4:28 AM, Patrick Hausmann patrick.hausm...@uni-bremen.de wrote: Dear list, # I have a DF like this: sleep$b  

Re: [R] ave and grouping

2009-03-02 Thread Dieter Menne
Patrick Hausmann Patrick.Hausmann at uni-bremen.de writes: # I have a DF like this: sleep$b - c(rep(8,10), rep(9,10)) sleep$me - with(sleep, ave(extra, group, FUN = mean)) sleep # I would like to create a new variable # holding the b-th value of group 1 and 2. # This is not what I

Re: [R] How to normalize to a set of internal references

2009-03-02 Thread Frank E Harrell Jr
Waverley wrote: Thanks for the advice. My question is more on how to do this? Let me use a biology gene analysis example to illustrate: In biology, there are always some house keeping genes which differ little even at pathological conditions. We know that at different batches, there are

Re: [R] ave and grouping

2009-03-02 Thread Gabor Grothendieck
There was an error in the line below. It should have been the following (i.e. the ix in extra[ix] was missing): sleep$newcol - with(sleep, ave(1:nrow(sleep), group, FUN = function(ix) extra[ix][b[ix][1]])) On Mon, Mar 2, 2009 at 8:01 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote:

Re: [R] Need help extracting info from XML file using XML package

2009-03-02 Thread Duncan Temple Lang
Wacek Kusnierczyk wrote: Don MacQueen wrote: I have an XML file that has within it the coordinates of some polygons that I would like to extract and use in R. The polygons are nested rather deeply. For example, I found by trial and error that I can extract the coordinates of one of them using

Re: [R] Inefficiency of SAS Programming

2009-03-02 Thread Duncan Murdoch
On 3/2/2009 6:57 AM, Thomas Levine wrote: R depends on all of those things to run, but you only have to use those programs through R. The software depends on these other tools, but the human doesn't have to switch interfaces. In fact, it doesn't even depend on them to run. Most Windows users

Re: [R] Inefficiency of SAS Programming

2009-03-02 Thread Gabor Grothendieck
If you want to write Sweave reports you have to learn latex and R does not hide that from you. This situation is somewhat better for tcltk, especially if you use one of the higher level wrapper packages that use it, but for serious work directly with it you need tcl/tk materials. On Mon, Mar 2,

Re: [R] Need help extracting info from XML file using XML package

2009-03-02 Thread Wacek Kusnierczyk
Duncan Temple Lang wrote: Wacek Kusnierczyk wrote: Don MacQueen wrote: I have an XML file that has within it the coordinates of some polygons that I would like to extract and use in R. The polygons are nested rather deeply. For example, I found by trial and error that I can extract the

[R] Newton-Raphson method and Quasi Newton

2009-03-02 Thread lelebecks
Hello, I need to apply the Newton-Raphson method and Quasi Newton method to a maximum log-likelihood function with three parameters. Which are the functions to use in R? And how to use it? Thanks, GB. __ R-help@r-project.org mailing list

Re: [R] Large data set in R

2009-03-02 Thread Kjetil Halvorsen
install.packages(biglm, dep=TRUE) library(help=biglm) kjetil On Mon, Mar 2, 2009 at 7:06 AM, Hardi sky_dr...@yahoo.com wrote: Hello, I'm trying to use R statistical packages to do ANOVA analysis using aov() and lm(). I'm having a problem when I have a large data set for input data from

[R] Bold Face in Plot

2009-03-02 Thread Rau, Roland
Dear all, I am trying to plot some text in bold face which works fine: plot(0:1,0:1,type=n) text(x=0.5, y=0.5, labels=expression(bold(the-actual-string))) Now when I try to do the following, the displayed text reads thestring: thestring - the-actual-string plot(0:1,0:1,type=n) text(x=0.5,

Re: [R] Using very large matrix

2009-03-02 Thread Steve_Friedman
I'm very interested in the bigmemory package for windows 32-bit environments. Who do I need to contact to request the Beta version? Thanks Steve Steve Friedman Ph. D. Spatial Statistical Analyst Everglades and Dry Tortugas National Park 950 N Krome Ave (3rd Floor) Homestead, Florida 33034

Re: [R] gamm (mgvc) and time-varying coefficient model

2009-03-02 Thread Simon Wood
Looks correct to me. You get a smooth of time for each level of x, so the smooths describe the way in which each coefficient of x varies in time. Note that the time varying coefficient for some level of a factor is actually given by the sum of the smooth for that factor level, and the estimate

Re: [R] smoothing a matrix (interpolate in plane)

2009-03-02 Thread Simon Wood
the thin plate spline, or tensor product smooths built into `mgcv' might be useful here (by default mgcv does automatic bandwidth selection for these). On Sunday 01 March 2009 09:51, Žroutík wrote: Hi R-users, I'd like to smooth a matrix to dismiss spikes and to interpolate in plane

Re: [R] Using very large matrix

2009-03-02 Thread Jay Emerson
Steve et.al., The old version is still on CRAN, but I strongly encourage anyone interested to email me directly and I'll make the new version available. In fact, I wouldn't mind just pulling the old version off of CRAN, but of course that's not a great idea. !-) Jay On Mon, Mar 2, 2009 at

Re: [R] R-code help for filtering with for loop

2009-03-02 Thread David Winsemius
The apply function which can work on either a row-wise of column-wise basis can be used with max and can return a logical vector that will let you separate the rows into those with and without a maximum greater than 60. datax - matrix(rnorm(300)*30,nrow=50) datax - as.data.frame(datax)

Re: [R] Bold Face in Plot

2009-03-02 Thread Prof Brian Ripley
On Mon, 2 Mar 2009, Rau, Roland wrote: Dear all, I am trying to plot some text in bold face which works fine: plot(0:1,0:1,type=n) text(x=0.5, y=0.5, labels=expression(bold(the-actual-string))) Now when I try to do the following, the displayed text reads thestring: thestring -

Re: [R] rounding problem

2009-03-02 Thread Wacek Kusnierczyk
Prof Brian Ripley wrote: I think your subject line should read 'Excel bug'. From the R help for round() Note that for rounding off a 5, the IEC 60559 standard is expected to be used, '_go to the even digit_'. In case you did not recognize it, IEC 60559 is an international

Re: [R] Bold Face in Plot

2009-03-02 Thread Rau, Roland
Dear Prof. Ripley, -Original Message- From: Prof Brian Ripley [mailto:rip...@stats.ox.ac.uk] Sent: Monday, March 02, 2009 3:05 PM To: Rau, Roland Cc: r-help@r-project.org Subject: Re: [R] Bold Face in Plot thestring - the-actual-string plot(0:1,0:1,type=n) text(x=0.5, y=0.5,

[R] Finding Lambda in Poisson distribution

2009-03-02 Thread Saeed Ahmadi
Hi, I have a dataset. First of all, I know that my dataset shall follow the Poission distribution. Now I have two questions: 1) How can I check that my data follow the Poisson distribution? 2) How can I calculate Lambda of my data? Regards Saeed -- View this message in context:

Re: [R] Bold Face in Plot

2009-03-02 Thread Dieter Menne
Rau, Roland Rau at demogr.mpg.de writes: I am trying to plot some text in bold face which works fine: plot(0:1,0:1,type=n) text(x=0.5, y=0.5, labels=expression(bold(the-actual-string))) Now when I try to do the following, the displayed text reads thestring: thestring -

[R] Bash script that uses an R command

2009-03-02 Thread stephen sefick
I have wriiten (with the help of the internet) a bash scirpt for my debian ppc 5.0 laptop to display battery remaining in a panel on the xfce desktop environment. echo `hal-device | grep battery.remaining_time | awk '{print $3/3600}' `-battery and this nicely spits out 4.95 hours in the panel.

Re: [R] Fwd: Converting R to Sweave (Rnw)

2009-03-02 Thread Neil Shephard
Rainer M Krug-6 wrote: Hi I am thinking about using Sweave more frequently, especially for documenting code. But the syntax is slightly awkward for me (name= ... @), and I was thinking if there would be a way of importing the type of code extracted from an Rnw file back into an Rnw

Re: [R] Bash script that uses an R command

2009-03-02 Thread Gabor Grothendieck
Make suitable changes for Linux (this was done in Windows Vista) but you should be able to do without R. In particular the double echo 4.95 | gawk {print int($1) : 60*($1-int($1))} On Mon, Mar 2, 2009 at 9:48 AM, stephen sefick ssef...@gmail.com wrote: I have wriiten (with the help of the

Re: [R] Fwd: Converting R to Sweave (Rnw)

2009-03-02 Thread Duncan Murdoch
On 3/2/2009 9:06 AM, Rainer M Krug wrote: Hi I am thinking about using Sweave more frequently, especially for documenting code. But the syntax is slightly awkward for me (name= ... @), and I was thinking if there would be a way of importing the type of code extracted from an Rnw file back into

Re: [R] RES: object .trPaths not found

2009-03-02 Thread rkevinburton
Yes. But I found out I the files permissions were set so that the configuration had a null effect and the configuration silently ignored the fact that the file could not be written to. Thank you. Kevin Leandro Marino lean...@cesgranrio.org.br wrote: Did you do the configuration of

Re: [R] RES: object .trPaths not found

2009-03-02 Thread rkevinburton
Yes. But I found out I the files permissions were set so that the configuration had a null effect and the configuration silently ignored the fact that the file could not be written to. Thank you. Kevin Leandro Marino lean...@cesgranrio.org.br wrote: Did you do the configuration of

Re: [R] Gradient function for optim.

2009-03-02 Thread rkevinburton
Thank you. I saw the source. But I am not sure how to get from .Internal(optim(...)) to fmingr. Kevin Katharine Mullen k...@few.vu.nl wrote: see the fmingr function in src/main/optim.c (https://svn.r-project.org/R/trunk/src/main/optim.c) On Wed, 25 Feb 2009 rkevinbur...@charter.net

Re: [R] Bash script that uses an R command

2009-03-02 Thread Romain Francois
stephen sefick wrote: I have wriiten (with the help of the internet) a bash scirpt for my debian ppc 5.0 laptop to display battery remaining in a panel on the xfce desktop environment. echo `hal-device | grep battery.remaining_time | awk '{print $3/3600}' `-battery and this nicely spits out

[R] Optim parscale?

2009-03-02 Thread rkevinburton
I am not clear on what is happening with parscale in optim It seems that scaling the parameters will produce unpredictable results in a non-linear function (which is the purpose of optim right?) The documentation states: parscale A vector of scaling values for the parameters. Optimization is

Re: [R] Need help extracting info from XML file using XML package

2009-03-02 Thread Duncan Temple Lang
Wacek Kusnierczyk wrote: Duncan Temple Lang wrote: Wacek Kusnierczyk wrote: Don MacQueen wrote: I have an XML file that has within it the coordinates of some polygons that I would like to extract and use in R. The polygons are nested rather deeply. For example, I found by trial and error

[R] Help with Wilcoxon Test

2009-03-02 Thread Amit Patel
Hi I have 2 sets of data that I want to do a Wilcoxon test on. They are of the same dimension. One has 4 zero values and the other has 5. dim(SampA) [1] 1 10 dim(SampV) [1] 1 10 I get the folowing error Error in wilcox.test.default(SampA, SampV, na.rm = TRUE, paired = FALSE, : 'x'

Re: [R] Bash script that uses an R command

2009-03-02 Thread stephen sefick
echo 'hal-device | grep battery.remaining_time | awk '{print$3/3600}' | awk '{print int($1)}:int(60*($1-int($1)))'` here is the final shell script is anyone is interested - this is written and working in debian linux Stephen Sefick On Mon, Mar 2, 2009 at 10:13 AM, Romain Francois

[R] portable R editor

2009-03-02 Thread Werner Wernersen
Hi, I have been dreaming about a complete R environment on my USB stick for a long time. Now I finally want to realize it but what I am missing is a good, portable editor for R which has tabs and syntax highlighting, can execute code, has bookmarks and a little project file management

Re: [R] handle graph size in eps

2009-03-02 Thread Benoit Boulinguiez
Sorry for the lack of information. I'm indeed under Windows. I indeed used the menu save as in the graph window. The matter with the eps obtained was the width of the graph which is lower than what I had on the screen or what I got when I saved it as a JPEG file. I tried the postscript command

Re: [R] handle graph size in eps

2009-03-02 Thread Clint Bowman
Try adding paper = special to your postscript arguments. Clint BowmanINTERNET: cl...@ecy.wa.gov Air Dispersion Modeler INTERNET: cl...@math.utah.edu Air Quality Program VOICE: (360) 407-6815 Department of Ecology FAX:

Re: [R] combining identify() and locator()

2009-03-02 Thread Brian Bolt
I am new to R and programming in general, so although your framework wasn't exactly what I needed, it was a tremendous help in getting what I needed. Here is what I ended up doing: idloc - function(xy,tol=0.10,hitcol=red,mispch=19,miscol=red){ if (type==flag) {col=red}else if

[R] comment on this book A Handbook of Statistical Analyses Using R by Brian S. Everitt (Author), Torsten Hothorn (Author)

2009-03-02 Thread choonhong ang
Is this book a good reference to learn R for statistical analysis ? A Handbook of Statistical Analyses Using R by Brian S. Everitthttp://www.amazon.com/exec/obidos/search-handle-url/ref=ntt_athr_dp_sr_1?%5Fencoding=UTF8search-type=ssindex=booksfield-author=Brian%20S.%20Everitt(Author), Torsten

Re: [R] handle graph size in eps

2009-03-02 Thread Prof Brian Ripley
See the help for argument 'paper' in ?postscript On Mon, 2 Mar 2009, Benoit Boulinguiez wrote: Sorry for the lack of information. I'm indeed under Windows. I indeed used the menu save as in the graph window. The matter with the eps obtained was the width of the graph which is lower than what

Re: [R] portable R editor

2009-03-02 Thread Wensui Liu
I feel emacs is portable enough for me. On 3/2/09, Werner Wernersen pensterfuz...@yahoo.de wrote: Hi, I have been dreaming about a complete R environment on my USB stick for a long time. Now I finally want to realize it but what I am missing is a good, portable editor for R which has tabs

Re: [R] comment on this book A Handbook of Statistical Analyses Using R by Brian S. Everitt (Author), Torsten Hothorn (Author)

2009-03-02 Thread Mitchell Maltenfort
Yes I think it is. I used it. On Mon, Mar 2, 2009 at 11:46 AM, choonhong ang angie.bear...@gmail.com wrote: Is this book a good reference to learn R for statistical analysis ? A Handbook of Statistical Analyses Using R by Brian S.

Re: [R] Help with Wilcoxon Test

2009-03-02 Thread David Winsemius
What do you get with: str(SampA) str(SampV) The error message suggests you are not giving it numeric vectors. -- David Winsemius On Mar 2, 2009, at 10:59 AM, Amit Patel wrote: Hi I have 2 sets of data that I want to do a Wilcoxon test on. They are of the same dimension. One has 4 zero

[R] Trouble - Installing R on RedHat el5

2009-03-02 Thread Steve_Friedman
I am lucky and am now working with a new Redhat Linux 64-bit OS But I am getting the following error can anyone provide some graciously needed assistance: [r...@bluebird system-files]# rpm -i R-2.8.1-1.rh5.x86_64.rpm warning: R-2.8.1-1.rh5.x86_64.rpm: Header V3 DSA signature: NOKEY, key ID

Re: [R] Need help extracting info from XML file using XML package

2009-03-02 Thread Wacek Kusnierczyk
Duncan Temple Lang wrote: not really: the xpath pattern '//coordinates' does say 'find all coordinates nodes searching from the root', but the root here is not the original root of the whole document, but each polygon node in turn. try: root = xmlInternalTreeParse(' root

Re: [R] Trouble - Installing R on RedHat el5

2009-03-02 Thread Marc Schwartz
on 03/02/2009 11:07 AM steve_fried...@nps.gov wrote: I am lucky and am now working with a new Redhat Linux 64-bit OS But I am getting the following error can anyone provide some graciously needed assistance: [r...@bluebird system-files]# rpm -i R-2.8.1-1.rh5.x86_64.rpm warning:

Re: [R] Bash script that uses an R command

2009-03-02 Thread Wacek Kusnierczyk
stephen sefick wrote: echo 'hal-device | grep battery.remaining_time | awk '{print$3/3600}' | awk '{print int($1)}:int(60*($1-int($1)))'` here is the final shell script is anyone is interested - this is written and working in debian linux it can't work -- you have unmatched quotes here.

[R] Goldbach partitions code

2009-03-02 Thread Murali.MENON
Folks, I put up a brief note describing my naive attempts to compute Goldbach partitions, starting with a brute-force approach and refining progressively. http://jostamon.blogspot.com/2009/02/goldbachs-comet.html I'd welcome your suggestions on improvements, alternatives, other

Re: [R] Newton-Raphson method and Quasi Newton

2009-03-02 Thread Ravi Varadhan
Check out: ?optim ?nlminb Ravi. Ravi Varadhan, Ph.D. Assistant Professor, Division of Geriatric Medicine and Gerontology School of Medicine Johns Hopkins University Ph. (410) 502-2619 email: rvarad...@jhmi.edu -

[R] TinnR Philips Webcam

2009-03-02 Thread waterhouse
Hi All, I have a Philips Webcam, model SPC110NC. I had to download a diver in order for the webcam to work with Skype from the philips support page (www.philips.com/support). Now when I start up TinnR to use with R, it opens the webcam. If I manually start up R, Tinn R will no longer synch

[R] QQplot

2009-03-02 Thread kayj
Hi All, this might be an easy question but I do not know how to find R-squared for a QQplot when doing linear regression? -- View this message in context: http://www.nabble.com/QQplot-tp22288218p22288218.html Sent from the R help mailing list archive at Nabble.com.

[R] xyplot color question

2009-03-02 Thread Paul Heinrich Dietrich
Hi, I am plotting scatterplots of horsepower by torque, conditional on brand (I'm just making up the variables for this example), and the goal is to see both the scatterplot points as well as the smoothed line. When I do the following code, I get the same color for the points and line, and would

[R] FW: partial residual plots

2009-03-02 Thread Culp, Dwayne
I was advised someone might be able to help me with this. Very truly yours, Dwayne E. Culp, P.E., CFM Effective February 8, 2009, my contact information becomes: Dwayne E. Culp, P.E., CFM | JACOBS ENGINEERING GROUP : North American Infrastructure | Manager: Hydrology Hydraulics Section

[R] barplot with specific order of x axis labels

2009-03-02 Thread R User R User
Hi all, I'd appreciate your help with this problem. I need to plot a barplot with the categories in a specific order. My data might be: hours Freq AN 10 MO 14 LU 30 I need the categories to be in the order: MO LU AN Is there some way to pass sort the dataset into this order so

Re: [R] logistic regression model validation through bootstrapping

2009-03-02 Thread Vivienne_O . Ozohili
Hi, I was wondering whether this query was addressed on how to perform validation through boostrapping. I am currently trying to implement a boostrapping approach to validation but don't know where to start. Help please. Thank you and Regards, Vivienne Ozohili Risk Model Validation Manager

[R] (no subject)

2009-03-02 Thread Brajkovic J.
Greetings, I am using fGarch package to estimate and simulate GARCH models. What I would like to do is to perform Monte Carlo simulation. Unfortunately I cannot figure how to modify the code to achieve this. I use the following code to run a single simulation: spec=garchSpec(model=list(ar=

[R] RWeka dataset running

2009-03-02 Thread ahmet kocyigit
Hi to all list members, I'm a newbie to R and will work on it from now. I'm preparing my thesis based on rpart and RWeka. I've managed to run already existing test data but couldn't import a data set from out of R. (By the way I'm an ubuntu user.) Thanks for your help Ahmet KOCYIGIT

[R] initial gradient and vmmin not finite

2009-03-02 Thread June Wong
Dear Rhelpers I have the problem with initial values, could you please tell me how to solve it? Thank you June p = summary(maxLik(fr,start=c(0,0,0,1,0,-25,-0.2))) Error in maxRoutine(fn = logLik, grad = grad, hess = hess, start = start, : NA in the initial gradient p =

Re: [R] TinnR Philips Webcam

2009-03-02 Thread Stefan Grosse
On Mon, 2 Mar 2009 12:48:16 -0500 (EST) waterho...@vims.edu wrote: WE Now when I start up TinnR to use with R, it opens the webcam. If I WE manually start up R, Tinn R will no longer synch with it. More detail would help to answer those questions properly. Which Tinn-R are you using, and how

Re: [R] Fwd: Converting R to Sweave (Rnw)

2009-03-02 Thread Duncan Murdoch
On 3/2/2009 10:12 AM, Duncan Murdoch wrote: On 3/2/2009 9:06 AM, Rainer M Krug wrote: Hi I am thinking about using Sweave more frequently, especially for documenting code. But the syntax is slightly awkward for me (name= ... @), and I was thinking if there would be a way of importing the type

[R] using par funtions in bargraph.CI()

2009-03-02 Thread maxa0006
I'm trying to create a bargraph of means with standard error bars using the function bargraph.CI (in the sciplot package). Like this: bargraph.CI(x.factor, response,data,xlab, ylab, par(family=serif),font=11) However, an error message comes up when I try to use the par funtion. Does the

[R] ways to put multiple graphs on single page (using ggplot2)

2009-03-02 Thread Juliet Hannah
Hi, Here are three plots: library(ggplot2) data(diamonds) randind - sample(nrow(diamonds),1000,replace=FALSE) dsmall - diamonds[randind,] qplot(carat, data=dsmall, geom=histogram,binwidth=1) qplot(carat, data=dsmall, geom=histogram,binwidth=.1) qplot(carat, data=dsmall,

Re: [R] Finding Lambda in Poisson distribution

2009-03-02 Thread Rau, Roland
Hi, -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Saeed Ahmadi Sent: Monday, March 02, 2009 3:16 PM To: r-help@r-project.org Subject: [R] Finding Lambda in Poisson distribution Hi, I have a dataset. First of all, I

Re: [R] FW: partial residual plots

2009-03-02 Thread Mark Difford
Dwayne Culp wrote: Is it possible to produce partial residual plots in R? See ?cr.plots in package car. This does Component+Residual (Partial Residual) Plots. Regards, Mark. Culp, Dwayne wrote: I was advised someone might be able to help me with this. Very truly yours, Dwayne E.

Re: [R] FW: partial residual plots

2009-03-02 Thread Mark Difford
Dwayne, Sorry, I forgot about this: also look at ?ceres.plots, also in package car. From Prof. Fox's documentation: Ceres plots are a generalization of component+residual (partial residual) plots that are less prone to leakage of nonlinearity among the predictors. Regards, Mark. Culp, Dwayne

  1   2   >