[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525258/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525259/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525263/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525269/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525261/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525260/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525268/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

Re: [R] How to display an image on RGL plot?

2009-11-26 Thread Peter Dalgaard
Vladimir Eremeev wrote: Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. The r-project would prefer to see one and only one copy of your post, though. I count five copies till now!!! -- O__ Peter Dalgaard Øster Farima

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525264/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

[R] Arrhenius Plot 2 with lattice

2009-11-26 Thread Markus Häge
Hello, sorry for incomplete code... with this I read the file and calculate my stuff. I have a plenty of them 80-300K every 5 Kelvin. I start with 79K, 80K...300K test<-read.table("T300_both.txt",header=FALSE,sep="") RH2<-c(RH2,2.5e7*.32e-4/100e-6/5100*(test$V3[c(2)]-test$V3[c(1)]+test $V3[c(4)]

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525266/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

Re: [R] How to display an image on RGL plot?

2009-11-26 Thread Barry Rowlingson
On Thu, Nov 26, 2009 at 7:14 AM, Vladimir Eremeev wrote: > The underlying picture is a JPEG image, loaded with the rimage package and > coerced to the matrix. > Spheres denote control points, collected from this picture and must be > situated over the certain points of the image. > I display the

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525271/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

Re: [R] arg min ???

2009-11-26 Thread Peterko
Petr that is what i was looking for. Thanks p.s. "ak si nahodou z Prahy mas u mna Pivo" Petr Pikal wrote: > > Hi > > r-help-boun...@r-project.org napsal dne 25.11.2009 16:50:49: > >> >> This is nice, but i have to define vector of possible theta, this is not > what >> i want to do. >> I hav

[R] How to display an image on RGL plot?

2009-11-26 Thread Vladimir Eremeev
Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525270/rgl-device.png rgl-device.png The underlying picture is a JPEG image, loaded with the rimage package and co

Re: [R] Concave hull

2009-11-26 Thread Corrado Topi
Dear David and other concave-hull-ists, yes, I meant concave hulls indeed. I know about the algorithm mentioed (www.concavehull.com) but it is not open source, so you cannot integrate it in R, and it is apparently patented, so even if you find the description you cannot apply it to implement a

[R] Sorry for repeated posts.

2009-11-26 Thread Vladimir Eremeev
I was sending my messages through Nabble, while it has went down for maintenance. -- View this message in context: http://old.nabble.com/Sorry-for-repeated-posts.-tp26525299p26525299.html Sent from the R help mailing list archive at Nabble.com. __ R-h

Re: [R] order of panels in xyplots

2009-11-26 Thread Titus Malsburg
Peter, thanks for your response! The problem is not how indexing works, but rather the question what is being indexed here. If I understand the description correctly then it is wrong. In the special and common case where all possible levels do actually occur in the data frame it coincidentally h

Re: [R] difference of two rows

2009-11-26 Thread clion
Yes, thanks, that works perfectly! great command b. jholtman wrote: > > Try this: > >> x <- read.table(textConnection("ID YEAR > + 13 2007 > + 15 2003 > + 15 2006 > + 15 2008 > + 21 2006 > + 21 2007"), header=TRUE) >> x$diff <- ave(x$YEAR, x$ID, FUN=function(a) c(diff(a), NA)) >> >> x >

Re: [R] Arrhenius Plot 2 with lattice

2009-11-26 Thread joris meys
Please provide minimal and self-contained code. Without being able to read in the dataset, it's impossible to check the code. Cheers Joris On Thu, Nov 26, 2009 at 9:13 AM, Markus Häge wrote: > Hello, > > sorry for incomplete code... > > with this I read the file and calculate my stuff. I have a

Re: [R] R help with princomp and pam clustering

2009-11-26 Thread joris meys
On Thu, Nov 26, 2009 at 1:04 AM, Tyler82 wrote: > > Hi all! > I am working with R package cluster and I have a little problem: > let's say I have two datasets...first one ("A") is divided into 4 clusters > by means of Pam algorythm. > Let's say I want to project the second database ("B") onto the

[R] analyse tab delimited textfile microarray data(help)

2009-11-26 Thread kumar kk
Hello I have microarray data in tab delimited text file.The main headers and some fields are as: following are the header: FEATURES FeatureNum Row Col SubTypeMask ControlType ProbeName SystematicName PositionX PositionY LogRatio LogRatioError PValueLogRatio gProcessedSignal rProcessedSignal gProce

Re: [R] multi variate plot with string data

2009-11-26 Thread frenchcr
is there something wrong with the way i asked the question frenchcr wrote: > > > Ive got two columns in data_set that are strings > > the first column is called "character" and has levels: > good, bad, ugly > > the second column is called "abusive" and has levels: > aggressive, moderately ag

[R] JMP <-> R ?

2009-11-26 Thread Robert Kinley
A wild thought ... ... anyone out there know if it's possible to call R code from JMP ? Bob Kinley - Eli Lilly & co, UK [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

[R] Error in do.call("expand.grid", dimnames(x)) :

2009-11-26 Thread deepak m r
Hi Friends, I am getting this error when I am running my R scripts. Can u please tell me how this error will happened. My script is given below. Errot: "Error in do.call("expand.grid", dimnames(x)) :" Script:

Re: [R] Concave hull

2009-11-26 Thread Duncan Murdoch
On 26/11/2009 4:02 AM, Corrado Topi wrote: Dear David and other concave-hull-ists, yes, I meant concave hulls indeed. I know about the algorithm mentioed (www.concavehull.com) but it is not open source, so you cannot integrate it in R, and it is apparently patented, so even if you find the des

[R] homals (level numerical, scoremat and catscores...)

2009-11-26 Thread Krusty the Klown
Just run this code: library(homals) data(sleeping) sleeping res <- homals(sleeping[,2:4], level = c(rep("numerical",2),"nominal"), ndim=1) plot(sleeping[,2],res$scoremat[,,1][,1],axes=F) axis(1,sleeping[,2]) abline(h=res$scoremat[,,1][,1]) box() It looks strange, innit? There is something which

[R] R: RE: R: Re: R: Re: chol( neg.def.matrix ) WAS: Re: Choleski and Choleski with pivoting of matrix fails

2009-11-26 Thread simona.racio...@libero.it
Thanks for your message! Actually it works quite well for me too. If I then take the trace of the final result below, I end up with a number made up of both a real and an imaginary part. This does not probably mean much if the trace of the matrix below givens me info about the degrees of freedom

Re: [R] How to display an image on RGL plot?

2009-11-26 Thread Duncan Murdoch
On 26/11/2009 2:14 AM, Vladimir Eremeev wrote: Hi all. It's been a long time since I wrote to this list. Glad to see the R project well and working. I am working with a 3D plot similar to this: http://old.nabble.com/file/p26525177/rgl-device.png rgl-device.png The underlying picture is a JP

Re: [R] How to display an image on RGL plot?

2009-11-26 Thread Martin Maechler
> Barry Rowlingson > on Thu, 26 Nov 2009 08:18:51 + writes: > On Thu, Nov 26, 2009 at 7:14 AM, Vladimir Eremeev wrote: >> The underlying picture is a JPEG image, loaded with the rimage package and >> coerced to the matrix. >> Spheres denote control points, colle

[R] Testing for strength of fit using R

2009-11-26 Thread Steve Murray
Dear all, I am trying to validate a model by comparing simulated output values against observed values. I have produced a simple X-y scatter plot with a 1:1 line, so that the closer the points fall to this line, the better the 'fit' between the modelled data and the observation data. I am now

[R] Vista truncated Sys.getenv("PATH") when and only when running as administrator

2009-11-26 Thread Keith Ponting
This is rather strange and I suspect is a Windows issue rather than an R issue, but it does not seem to be mentioned on the lists so I am logging it in case it gives anyone else problems. Running Vista Business Version 6.0.6002 Service Pack 2 Build 6002, with the latest R 2.10.0 patched. A

Re: [R] Kerning issues with CairoPDF

2009-11-26 Thread Karl Ove Hufthammer
Karl Ove Hufthammer wrote: > I'm using CairoPDF to generate PDF (because of its font embedding and > support for transparent colours). However, at least on my (Windows) > system, the text it outputs seems to have completely wrong kerning. > Here's an example: > > CairoPDF("test.pdf") > plot(rno

Re: [R] order of panels in xyplots

2009-11-26 Thread Peter Ehlers
You're right, Titus. I misunderstood. It looks like index.cond has to be in 1:(number of panels being plotted for factor f). While this can arguably be covered by the phrase "valid indexing vector", I agree that this could be made more explicit. -Peter Ehlers Titus Malsburg wrote: Peter, thank

Re: [R] ow to have R automatically print traceback upon errors

2009-11-26 Thread Hao Cen
Hi Jim, Thanks for your suggestion. It is working now. The earlier message " recover called non-interactively; > frames dumped, use debugger() to view" happened within a R-IDE but doesn't occur in a R console. Jeff -Original Message- From: jim holtman [mailto:jholt...@gmail.com] Sent: M

Re: [R] makefile for sweave

2009-11-26 Thread Wolfgang Raffelsberger
Dear all, I can't get texi2dvi working right. Basically I'd like to convert a .lex to .pdf without having to fiddle with the issue Sweave.sty not being in my current directory (as this was sugested in other posts on this list). When I'm in the R-Gui I can get the help via ?texi2dvi (So I conc

Re: [R] makefile for sweave

2009-11-26 Thread Sundar Dorai-Raj
Is texi2dvi in your PATH? What happens if you open a CMD window and type texi2dvi at the prompt? --sundar On Thu, Nov 26, 2009 at 6:14 AM, Wolfgang Raffelsberger wrote: > Dear all, > > I can't get texi2dvi working right. Basically I'd like to convert a .lex to > .pdf without having to fiddle wit

Re: [R] questions on the ff package

2009-11-26 Thread Hao Cen
Hi Jens, Thanks for your prompt and informative answers. ff is a fabulous package and your suggestions helped me solve my problems at hands. As I need to incrementally increase each of several large matrices (about 1000 rows *1 columns, 1000 matrices) by a row every day. I wonder how efficien

[R] Moving quantile()?

2009-11-26 Thread Rob Steele
Is there a faster way to get moving quantiles from a time series than to run quantile() at each step in the series? Thanks, Rob __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www

[R] Best way to preallocate numeric NA array?

2009-11-26 Thread Rob Steele
These are the ways that occur to me. ## This produces a logical vector, which will get converted to a numeric ## vector the first time a number is assigned to it. That seems ## wasteful. x <- rep(NA, n) ## This does the conversion ahead of time but it's still creating a ## logical vector first,

Re: [R] analyse tab delimited textfile microarray data(help)

2009-11-26 Thread Martin Morgan
Kumar -- ask on the Bioconductor mailing list http://bioconductor.org/docs/mailList.html see basic workflows at http://bioconductor.org/docs/workflows/oligoarrays/ You cannot analyze the data without understanding its source or questions you are asking, so determine where this data is comi

Re: [R] Best way to preallocate numeric NA array?

2009-11-26 Thread Douglas Bates
On Thu, Nov 26, 2009 at 10:03 AM, Rob Steele wrote: > These are the ways that occur to me. > > ## This produces a logical vector, which will get converted to a numeric > ## vector the first time a number is assigned to it.  That seems > ## wasteful. > x <- rep(NA, n) > > ## This does the conversio

[R] Missing and number of days

2009-11-26 Thread Val
Hi all, Assume that I have a data set (“xcv”) with several variables and some of the variables have a missing observation represented by -9 as shown below. I want to exclude these observations from the analysis ( as a NA). Is there a command that I can do it for the entire data set rather

Re: [R] Best way to preallocate numeric NA array?

2009-11-26 Thread Henrique Dallazuanna
You can try this also: rep(NA_integer_, 10) On Thu, Nov 26, 2009 at 2:03 PM, Rob Steele wrote: > These are the ways that occur to me. > > ## This produces a logical vector, which will get converted to a numeric > ## vector the first time a number is assigned to it.  That seems > ## wasteful. > x

Re: [R] Testing for strength of fit using R

2009-11-26 Thread David Winsemius
On Nov 26, 2009, at 9:48 AM, Steve Murray wrote: Dear all, I am trying to validate a model by comparing simulated output values against observed values. I have produced a simple X-y scatter plot with a 1:1 line, so that the closer the points fall to this line, the better the 'fit' betwe

Re: [R] Best way to preallocate numeric NA array?

2009-11-26 Thread Henrique Dallazuanna
Or best: rep(NA_real_, 10) On Thu, Nov 26, 2009 at 2:31 PM, Henrique Dallazuanna wrote: > You can try this also: > > rep(NA_integer_, 10) > > On Thu, Nov 26, 2009 at 2:03 PM, Rob Steele > wrote: >> These are the ways that occur to me. >> >> ## This produces a logical vector, which will get conv

[R] mcp.area vs. kernel.area Output... ?

2009-11-26 Thread T.D.Rudolph
I am trying to estimate home range size using 2 different methods in the adehabitat package, but I am slightly confounded by the results. ## Attached is an R object file containing animal relocations with a field for "id", and "x" & "y" coordinates ## (in metres) load("temp") require(adehabitat)

Re: [R] Missing and number of days

2009-11-26 Thread R Help
For your first question, I believe that dat[dat==-9]=NA should do the trick For your second, are you dates in the Date() format? If not, try ?Date Once you get them in Date format, then you can simply subtract them d1 = as.Date('20090604',format="%Y%m%d") d2 = as.Date('20080604',format="%Y%m%

[R] different fits for geese and geeglm in geepack?

2009-11-26 Thread Achaz von Hardenberg
Dear all, I am trying to fit a GEE model on eagle productivity (number of hatched offspring per nest) using the geeglm function in the library geepack and I found an odd result. My understanding is that the function geese and geeglm should give the same fits, as actually geeglm uses geese to

Re: [R] Best way to preallocate numeric NA array?

2009-11-26 Thread Rob Steele
Douglas Bates wrote: > On Thu, Nov 26, 2009 at 10:03 AM, Rob Steele > wrote: >> These are the ways that occur to me. >> >> ## This produces a logical vector, which will get converted to a numeric >> ## vector the first time a number is assigned to it. That seems >> ## wasteful. >> x <- rep(NA, n)

[R] Does nargin and nargout work with R functions?

2009-11-26 Thread Jason Rupert
I am porting some MATLAB functions over to R and hopefully into a package, so I am curious if nargin and nargout work with R functions. Here is kind of an example of where I need to head in order to port "control-1.0.11" from Octave over to R. The Octave "control-1.0.11" package has the capabi

Re: [R] Best way to preallocate numeric NA array?

2009-11-26 Thread Rob Steele
Douglas Bates wrote: > On Thu, Nov 26, 2009 at 10:03 AM, Rob Steele > wrote: >> These are the ways that occur to me. >> >> ## This produces a logical vector, which will get converted to a numeric >> ## vector the first time a number is assigned to it. That seems >> ## wasteful. >> x <- rep(NA, n)

Re: [R] Testing for strength of fit using R

2009-11-26 Thread Peter Ehlers
Steve Murray wrote: Dear all, I am trying to validate a model by comparing simulated output values against observed values. I have produced a simple X-y scatter plot with a 1:1 line, so that the closer the points fall to this line, the better the 'fit' between the modelled data and the obser

[R] Save workspace with ff objects

2009-11-26 Thread Yue Sheng
Hi All, My script generates a mixture of normal and ff objects. I need to run the script for different parameter settings. Very often I would like to save the workspace for each parameter setting, so that I can get back to it later on. Is there an easy way to do this, instead of needing to save in

Re: [R] Does nargin and nargout work with R functions?

2009-11-26 Thread baptiste auguie
Hi, I think you can use match.call() to retrieve the number of arguments passed to a function (see below), but I don't think nargout makes sense in R like it does in Matlab. foo <- function(...){ print(match.call()) nargin <- length(as.list(match.call())) -1 print(nargin) } foo(a=1, b=2)

[R] How shall one present LRT test statistic in a scientific journal ?

2009-11-26 Thread JVezilier
Hello !! I'm recently having a debate with my PhD supervisor regarding how to write the result of a likelihood ratio test in an article I'm about to submit. I analysed my data using "lme" mixed modelling. To get some p-values for my fixed effect I used model simplification and the typical outpu

Re: [R] Grouped Boxplot

2009-11-26 Thread Gary
Thanks Petr! On Thu, Nov 26, 2009 at 12:31 AM, Petr PIKAL wrote: > Hi > > r-help-boun...@r-project.org napsal dne 24.11.2009 18:44:43: > > > Back in March Soren Vogel asked exactly the same thing: > > > > Here is the solution that was offered then. (He offered a dataset as > > requested in the P

Re: [R] How shall one present LRT test statistic in a scientific journal ?

2009-11-26 Thread David Winsemius
On Nov 26, 2009, at 12:14 PM, JVezilier wrote: Hello !! I'm recently having a debate with my PhD supervisor regarding how to write the result of a likelihood ratio test in an article I'm about to submit. I analysed my data using "lme" mixed modelling. To get some p-values for my fixed

Re: [R] Help: Beanplots calculating wrong average

2009-11-26 Thread Michael Hopgood
Hi Tom, Thank you for the friendly and informative answer. It does explain a lot of things, actually. As with any good answer, it inevitably leads to other questions. In the first place, I need the arithmetic mean. It's what we base our calculations on... My code is currently this: Metall<-c

[R] Export kde object as shapefile

2009-11-26 Thread T.D.Rudolph
I am trying to estimate home range size using the plug-in method with kernel density estimation in the kernel smoothing (ks) package. Unless there is another way I am not familiar with, in order to calculate spatial area under the space I need to convert my kde () object into a spatial object som

Re: [R] How shall one present LRT test statistic in a scientific journal ?

2009-11-26 Thread Peter Dalgaard
David Winsemius wrote: > > On Nov 26, 2009, at 12:14 PM, JVezilier wrote: > >> >> Hello !! >> >> I'm recently having a debate with my PhD supervisor regarding how to >> write >> the result of a likelihood ratio test in an article I'm about to submit. >> >> I analysed my data using "lme" mixed mod

Re: [R] makefile for sweave

2009-11-26 Thread Duncan Murdoch
On 26/11/2009 9:14 AM, Wolfgang Raffelsberger wrote: Dear all, I can't get texi2dvi working right. Basically I'd like to convert a .lex to .pdf without having to fiddle with the issue Sweave.sty not being in my current directory (as this was sugested in other posts on this list). When I'm in

Re: [R] How shall one present LRT test statistic in a scientific journal ?

2009-11-26 Thread David Winsemius
On Nov 26, 2009, at 12:46 PM, Peter Dalgaard wrote: David Winsemius wrote: On Nov 26, 2009, at 12:14 PM, JVezilier wrote: Hello !! I'm recently having a debate with my PhD supervisor regarding how to write the result of a likelihood ratio test in an article I'm about to submit. I anal

[R] barchart() {Lattice} help.

2009-11-26 Thread Peng Cai
Hi R Users, I'm trying to plot a stacked barplot. Here is data: Sample Col1 Col2 Col3 Row1 -2 4 -1 Row2 3 -2 4 Row3 3 5 -2 Row4 4 1 -1 I'm using following R code: library(lattice) dta<-read.table("data.txt", header=TRUE, row.names="Sample") barchart(data.matrix(dta), horizontal=FALSE,

Re: [R] Export kde object as shapefile

2009-11-26 Thread David Winsemius
On Nov 26, 2009, at 12:40 PM, T.D.Rudolph wrote: I am trying to estimate home range size using the plug-in method with kernel density estimation in the kernel smoothing (ks) package. Unless there is another way I am not familiar with, in order to calculate spatial area under the space I

Re: [R] Export kde object as shapefile

2009-11-26 Thread T.D.Rudolph
http://old.nabble.com/file/p26533942/temp temp David Winsemius wrote: > > > On Nov 26, 2009, at 12:40 PM, T.D.Rudolph wrote: > >> >> I am trying to estimate home range size using the plug-in method >> with kernel >> density estimation in the kernel smoothing (ks) package. Unless >> ther

Re: [R] Insert elements into a vector in a defined positions

2009-11-26 Thread Manuel Ramon
That's good . Your solution works for me. Than you Rolf. Rolf Turner-3 wrote: > > > On 26/11/2009, at 10:46 AM, Manuel Ramon wrote: > >> >> Dear R users, >> I have a vector of length n and I want to insert some elements (in >> my case >> the NA string) into a defined positions. For example

Re: [R] Feature request for as.Date() function

2009-11-26 Thread nabble . 30 . miller_2555
On Thu, Nov 26, 2009 at 12:08 AM, jim holtman - jholt...@gmail.com <+nabble+miller_2555+9dc9649aca.jholtman#gmail@spamgourmet.com> wrote: > An easy way is just to write your own function that will accept "NA", > convert it to NA and then call as.Date. > I have written such a function, which ha

Re: [R] rattle(): unable to load shared library

2009-11-26 Thread Graham Williams
Sometimes I've seen this under MS/Windows and the solution has been to reinstall the gtk libraries independent of R and to make sure this is the one installed: http://downloads.sourceforge.net/gladewin32/gtk-2.12.9-win32-2.exe. It looks like glade was not in the gtk libraries you installed? Hope t

[R] lattice --- different properties of lines corresponding to type=c("l", "a") respectively

2009-11-26 Thread Kjetil Halvorsen
I think the subject says it all. I want to make a simple lattice plot, using xyplot with the argument type=c("l","a"). The problem then is that in the resulting plot it is difficult/impossible to see which plot corresponds to the average and which to the individual profiles. I triedthings like

Re: [R] Has anyone had success with RGTK2/rattle and windows 7 64-bit?

2009-11-26 Thread Graham Williams
I don't have a Windows 7 to test this on yet - works on Vista and XP. Did you install the GTK libraries (separately to R)? Regards, Graham 2009/11/21 Tetrick, Scott > I have been unable to get rattle to run in my new Windows7-64 bit > configuration. For wither Rgtk2 or rattle, I get an error:

Re: [R] different fits for geese and geeglm in geepack?

2009-11-26 Thread Søren Højsgaard
It looks strange. If you send me data for a small reproducible example then I'll look into the issue. Regards Søren Fra: r-help-boun...@r-project.org [r-help-boun...@r-project.org] På vegne af Achaz von Hardenberg [achaz.hardenb...@gmail.com] Sendt: 26.

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread Peter Ehlers
Peng Cai wrote: Hi R Users, I'm trying to plot a stacked barplot. Here is data: Sample Col1 Col2 Col3 Row1 -2 4 -1 Row2 3 -2 4 Row3 3 5 -2 Row4 4 1 -1 I'm using following R code: library(lattice) dta<-read.table("data.txt", header=TRUE, row.names="Sample") barchart(data.matrix(dta),

Re: [R] Moving quantile()?

2009-11-26 Thread Charles C. Berry
On Thu, 26 Nov 2009, Rob Steele wrote: Is there a faster way to get moving quantiles from a time series than to run quantile() at each step in the series? Yes. Run help.request() Since you have already done the first 4 items below (right?) you will answer 'y', but when you get to

[R] r code to combine random forest and glm

2009-11-26 Thread Sharma, Dhruv
Hi, is it possible to have r code which runs random forests and takes the top 7 variables based on variable importance and builds a model of trying interactions terms for each variable. Then builds 7 models and keeps the one with the lowest aic. this would be like a new form of step wise regr

Re: [R] Concave hull

2009-11-26 Thread Ted Harding
Raising a rather general question here. This is a tantalising discussion, but the notion of "concave hull" strikes me as extremely ill-defined! I'd like to see statement of what it is (generically) supposed to be. The examples discussed so far seem to rely on a person's inner feelings of what it

[R] {ggplot2} Adding Mean to (grouped) Boxplot.

2009-11-26 Thread Peng Cai
Hi R Users, I am using following R code to plot a "grouped boxplot". I'm hoping if I can add MEAN to these boxplots. Data is copied below and attached as text file. install.packages("ggplot2") library(ggplot2) dta<-read.table("Sample.txt",header=T) attach(dta) p <- ggplot(dta, aes(factor(month),

Re: [R] How shall one present LRT test statistic in a scientific journal ?

2009-11-26 Thread Frank E Harrell Jr
David Winsemius wrote: On Nov 26, 2009, at 12:46 PM, Peter Dalgaard wrote: David Winsemius wrote: On Nov 26, 2009, at 12:14 PM, JVezilier wrote: Hello !! I'm recently having a debate with my PhD supervisor regarding how to write the result of a likelihood ratio test in an article I'm abo

Re: [R] Concave hull

2009-11-26 Thread baptiste auguie
2009/11/26 Ted Harding : > Raising a rather general question here. > > This is a tantalising discussion, but the notion of "concave hull" > strikes me as extremely ill-defined! > > I'd like to see statement of what it is (generically) supposed to be. I'm curious too, but I can imagine the followin

[R] Barplot with confidence intervals

2009-11-26 Thread Jonhnny Weslley
I have a file in the following format: Scenario1 Scenario1CIL Scenario1CIU Scenario2 Scenario2CIL Scenario2CIU 60 57 62 45 48 50 110 101 111 51 50 52 120 117 122 64 62 66 192 190 194 79 75 79 where: First column = Scenario1 mean value Second column = Scenario1 Low Confidence Interval Third c

Re: [R] {ggplot2} Adding Mean to (grouped) Boxplot.

2009-11-26 Thread Peng Cai
Hi Dennis, Thanks for your reply. Answers are below: On Thu, Nov 26, 2009 at 4:14 PM, Dennis Murphy wrote: > Two questions: > (1) Do you want monthly means (averaged over groups), group means > (averaged over > months), overall average, other??? > I need mean (of nail) for each tre

Re: [R] Concave hull

2009-11-26 Thread Ted Harding
On 26-Nov-09 21:11:02, baptiste auguie wrote: > 2009/11/26 Ted Harding : >> Raising a rather general question here. >> >> This is a tantalising discussion, but the notion of "concave hull" >> strikes me as extremely ill-defined! >> >> I'd like to see statement of what it is (generically) supposed t

Re: [R] Concave hull

2009-11-26 Thread Barry Rowlingson
On Thu, Nov 26, 2009 at 9:45 PM, Ted Harding wrote: > So it is still an undefined solution. As is yours -- since you might > want to use different radii of spheres from different directions. I think the formal and rigorous definition is "a nice polygon that goes round my points". Barry __

[R] How to suppress errors generated by readHTMLTable?

2009-11-26 Thread Peng Yu
library(XML) download.file('http://polya.umdnj.edu/polya_db2/gene.php?llid=109079&unigene=&submit=Submit','index.html') tables=readHTMLTable("index.html",error=function(...){}) tables readHTMLTable gives me the following errors. Could somebody let me know how to suppress them? Opening and endi

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread Peng Cai
Thanks Peter and Dennis for your help! On Thu, Nov 26, 2009 at 3:50 PM, Dennis Murphy wrote: > Hi: > > In addition to Peter's suggestions, I would suggest that you get rid of the > row names as the > first column of your data frame, as it confuses barchart (and the > interpretation of the plot).

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread Peng Cai
Hi Peter, I'm not sure but it seems "scales" command works only with integer values. If the y-axis values are very small (such as -0.03, -0.02, -0.01, 0, 0.01,..., 0.08). My current plot has values 0, 0.05, and 0.10 only. But I need it to extend it to negative numbers and reduce the scale width (

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread Peter Ehlers
Peng Cai wrote: Hi Peter, I'm not sure but it seems "scales" command works only with integer values. If the y-axis values are very small (such as -0.03, -0.02, -0.01, 0, 0.01,..., 0.08). My current plot has values 0, 0.05, and 0.10 only. But I need it to extend it to negative numbers and reduce

Re: [R] Concave hull

2009-11-26 Thread Kjetil Halvorsen
There is a package on CRAN implementing such an idea: alphahull, phull is other package, kjetil On Thu, Nov 26, 2009 at 6:11 PM, baptiste auguie wrote: > 2009/11/26 Ted Harding : >> Raising a rather general question here. >> >> This is a tantalising discussion, but the notion of "concave hull" >

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread Peng Cai
Thanks a lot Peter! One more help, is there a similar function abline() for barchart(). I'm trying to add a (light gray colored) horizontal lines, one for each y-value. Peng On Thu, Nov 26, 2009 at 5:59 PM, Peter Ehlers wrote: > Peng Cai wrote: > >> Hi Peter, >> >> I'm not sure but it seems "s

Re: [R] How to suppress errors generated by readHTMLTable?

2009-11-26 Thread Duncan Temple Lang
Just this morning, I made suppressing these parser messages the default behavior for htmlParse() and that will apply to readHTMLTable() also. Until I release that (along with another potentially non-backward compatible change regarding character encoding), you can use readHTMLTable(htmlParse("i

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread David Winsemius
On Nov 26, 2009, at 6:12 PM, Peng Cai wrote: Thanks a lot Peter! One more help, is there a similar function abline() for barchart(). ?panel.abline I'm trying to add a (light gray colored) horizontal lines, one for each y-value. Peng On Thu, Nov 26, 2009 at 5:59 PM, Peter Ehlers w

Re: [R] Help on Html format R 2.10

2009-11-26 Thread Ruben Cabrera
Hello: I´m trying to change the plain text help to html help but I don´t know if this is possible with R 2.10.0 (Windows XP). The installation process always set to plain text. I found this: http://old.nabble.com/change-in-default-behavior-of--functionname-in-R-2.10--td26115436.html Thanks in ad

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread Peng Cai
Thanks David, I tried panel.abline(h=somevalue) -- both inside and outside of barchart() function but its not working. Any suggestions? Peng On Thu, Nov 26, 2009 at 6:42 PM, David Winsemius wrote: > > On Nov 26, 2009, at 6:12 PM, Peng Cai wrote: > > Thanks a lot Peter! One more help, is there a

Re: [R] Export kde object as shapefile

2009-11-26 Thread David Winsemius
On Nov 26, 2009, at 2:25 PM, T.D.Rudolph wrote: methods can be used to estimate home range size. The function kernel.area in the adehabitat package, for example, produces a spatial area I interpret as being equivalent to estimated home range size using one of two possible smoothing para

Re: [R] Help on Html format R 2.10

2009-11-26 Thread Ruben Cabrera
Hello Again: Let me answer myself because I already solved my problem: I found the file "...R\R-2.10.0\Rprofile.site" and used options(help_type="html"). Thanks. Ruben Cabrera On Thu, Nov 26, 2009 at 7:41 PM, Ruben Cabrera wrote: > Hello: > > I´m trying to change the plain text help to html he

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread Peter Ehlers
Peng Cai wrote: Thanks David, I tried panel.abline(h=somevalue) -- both inside and outside of barchart() function but its not working. Any suggestions? Peng Here's some code related to the data you posted earlier. barchart(data.matrix(dta), horizontal = FALSE, stack = TRUE, par.se

Re: [R] barchart() {Lattice} help.

2009-11-26 Thread David Winsemius
Cannot reproduce. Seems to work as documented. Maybe your code (whatever it might be) is not correct? -- David. On Nov 26, 2009, at 6:53 PM, Peng Cai wrote: Thanks David, I tried panel.abline(h=somevalue) -- both inside and outside of barchart() function but its not working. Any suggestions

Re: [R] Moving quantile()?

2009-11-26 Thread Rob Steele
Charles C. Berry wrote: > On Thu, 26 Nov 2009, Rob Steele wrote: > >> Is there a faster way to get moving quantiles from a time series than to >> run quantile() at each step in the series? > > > Yes. > > Run > > help.request() > > Since you have already done the first 4 items below (right

Re: [R] multi variate plot with string data

2009-11-26 Thread frenchcr
Thanks!! this works perfectly dat = data.frame(Aggressive = data_set$var1 , Behaviour=data_set$var2) plot(table(dat)) ...its better than a stacked boxplot as it gives also different width of columns on the diagram...an extra layer of information. frenchcr wrote: > > > Ive got

Re: [R] Moving quantile()?

2009-11-26 Thread Rob Steele
Charles C. Berry wrote: > On Thu, 26 Nov 2009, Rob Steele wrote: > >> Is there a faster way to get moving quantiles from a time series than to >> run quantile() at each step in the series? > > > Yes. > > Run > > help.request() > > Since you have already done the first 4 items below (right

  1   2   >