[R] Odp: GRR - Best methods in R

2007-09-17 Thread Petr PIKAL
Hi below you can find a function I wrote to mimic a Minitab R+R procedure. It is intended for myself only so it is not commented. If you want to know how it functions contact me off-list. Regards Petr [EMAIL PROTECTED] napsal dne 17.09.2007 05:55:00: Dear List, Which is the most

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

2007-09-17 Thread Ptit_Bleu
Hi, 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 numerical value of the integration of this graph. I looked at

[R] vector name

2007-09-17 Thread livia
I have got a list named filtered, I would like to construct alist named fdata as following: fdata - cbind(matrix(unlist(filtered),ncol=28), myregime) If I try names(filtered), it gives all the correct name for each vector, but if I try names(fdata), it appears filtered[[1]]

Re: [R] vector name

2007-09-17 Thread Stephen Tucker
Hi, Either of the following should work (I assume there are 28 elements in your list): fdata - cbind(as.matrix(as.data.frame(filtered)), myregime) fdata - cbind(colnames-(matrix(unlist(filtered),ncol=28), names(filtered)), myregime) --- livia [EMAIL PROTECTED] wrote: I have got a

Re: [R] programming question

2007-09-17 Thread Adrian Dusa
On Monday 17 September 2007, you wrote: R 2.6.0 has Reduce; myvec - c(2, 8, 24, 26, 51, 57, 58, 78, 219) Reduce(function(myvec, p) setdiff(myvec, findSubsets2(p)), myvec, myvec) Thanks Gabor, at first I jumped off my chair but... for many input variables it takes ages to reduce the vector.

[R] aggregate function oddity

2007-09-17 Thread Mihalicza Péter
__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] Matrix entries not recognised as numeric

2007-09-17 Thread snapperball
Hello, I am using a number of large matrices in a project. I read a couple of the matrices using the read.csv command. However in some instances, R does not recognize the resulting matrices as consisting of numerical values (I suspect R considers them as strings or factors) for instance when I

[R] add group to boxplot

2007-09-17 Thread Karin Lagesen
I have two sets of data in data frames (different dimensions). Now, I am able to make boxplots for one of these. I am using a formula, which gives me three boxplots which is automatically placed at at = c(1:3), since there are three groups. However, I would now like to add data from the other

Re: [R] Matrix entries not recognised as numeric

2007-09-17 Thread Paul Smith
On 9/17/07, snapperball [EMAIL PROTECTED] wrote: I am using a number of large matrices in a project. I read a couple of the matrices using the read.csv command. However in some instances, R does not recognize the resulting matrices as consisting of numerical values (I suspect R considers them

[R] graphs with gradients of colors

2007-09-17 Thread Van Dongen Stefan
Hi All, I would like to fill the area under a curve with a gradient of colors. Are there any packages or trick I could use Thanks Stefan Stefan Van Dongen Antwerp [[alternative HTML version deleted]] __ R-help@r-project.org mailing

Re: [R] Creating Hmisc labels inside a function

2007-09-17 Thread Steve Powell
Thanks, Frank - it doesn't work though. Your suggestion was: Test=function(obj,labe) { label(obj)=labe #at this point add the line: obj } #...but just returning the object does not permanently change the label: obj That is why I wanted to use assign. But assign will not work with attributes,

Re: [R] Matrix entries not recognised as numeric

2007-09-17 Thread Ted Harding
On 17-Sep-07 13:18:10, snapperball wrote: Hello, I am using a number of large matrices in a project. I read a couple of the matrices using the read.csv command. However in some instances, R does not recognize the resulting matrices as consisting of numerical values (I suspect R considers

Re: [R] help not working

2007-09-17 Thread Uwe Ligges
Bill Szkotnicki wrote: Hi, When we install windows R 2.5.1 on a network drive ( the P: ' drive, actually a samba public share ) R runs well but the windows help is unavailable. It complains about Navigation to the webpage was cancelled Does anyone know how to fix that? This has been

[R] Converting zoo dates to ts dates

2007-09-17 Thread Edna Bell
Hi R Gurus I'm using get.hist.quote to obtain monthly historical data. This results in a zoo series. I would like to convert the zoo series to a time series within a function. How do I get the starting date for the ts from the zoo series, please? Thanks, Edna Bell

Re: [R] system.time behavior using source

2007-09-17 Thread Gabor Grothendieck
Try this: source(myfile.R, echo = TRUE) On 9/17/07, Leeds, Mark (IED) [EMAIL PROTECTED] wrote: If I type the line below at my Rprompt, it works fine. But, if I put the same line in a file called systest.R and then source that file at the prompt, nothing happens. IF I use R CMD BATCH, it also

Re: [R] Class probabilities in rpart

2007-09-17 Thread Terry Therneau
The predict.rpart() function from the rpart library allows for calculating the class probabilities for a given test case instead of a discrete class label. How are these class probabilities derived? Is it simply the proportion of the majority class to all cases in a leaf node? That

Re: [R] side bars on dendrograms with latticeExtra

2007-09-17 Thread deepayan . sarkar
On 9/17/07, audrey [EMAIL PROTECTED] wrote: Dear all, I am using the heatmap representations of latticeExtra package and I would be interested to draw color side bars representing the groups of a factor of interest. From my understanding of : help(dendrogramGrob) ... The 'add'

Re: [R] help for high-quality plot (wmf) in files

2007-09-17 Thread Bert Gunter
May I ask a related question to tack on to this thread. One can also output a wmf plot to a file through the win.metafile() call. But as the help file says, due to Windows limitations only one plot per file is allowed and one must use parameterized filenames to produce multiple plots, which I

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

Re: [R] removing a specific number of digist from a character string

2007-09-17 Thread Jeffrey Robert Spies
test - c(060907_17_3_5_1_1_2909.tif, 060907_17_3_5_2_1_2910.tif, 060907_17_3_5_3_1_2911.tif) sub('[[:digit:]][[:digit:]][[:digit:]][[:digit:]]\.tif', '', test) or test - c(060907_17_3_5_1_1_2909.tif, 060907_17_3_5_2_1_2910.tif, 060907_17_3_5_3_1_2911.tif) sub('[[:digit:]]{4}\.tif', '', test)

[R] data frame

2007-09-17 Thread Alfredo Alessandrini
Hi everybody, If I've a data frame like this: dataframe a X0 X2 X4 X6 X8 X10 X12 X14 X16 1957 0 0 0 0 0 0 0 0 0 1958 0 0 0 0 0 0 0 0 0 1959 831 0 0 0 0 0 0 0 0 1960 544 282 0 0 0 0 0 0 0 1961 446 365 0 0 0 0 0

Re: [R] data frame

2007-09-17 Thread Stefano Calza
Ciao, do you mean the number of records? something like nrow(a) I would also suggest ncol(a) instad of length, which would work for a matrix as well. [Italian] E' quello che cercavi? [/Italian] Ciao Stefano On Mon, Sep 17, 2007 at 06:05:58PM +0200, Alfredo Alessandrini wrote: AlfredoHi

Re: [R] data frame

2007-09-17 Thread Jeffrey Robert Spies
I believe you're looking for: dim(a) dim(a)[1] # Number of observations, in your example, 12 dim(a)[2] # Number of variables per observation, in your example, 9 --Jeff. On Sep 17, 2007, at 12:05 PM, Alfredo Alessandrini wrote: Hi everybody, If I've a data frame like this: dataframe a

Re: [R] removing a specific number of digist from a character string

2007-09-17 Thread Jeffrey Robert Spies
For the sake of absolute correctness: sub('[[:digit:]]{4}\.tif', '', test) should be sub('[[:digit:]]{4}\\.tif', '', test) -- Jeff. On Sep 17, 2007, at 11:59 AM, Jeffrey Robert Spies wrote: test - c(060907_17_3_5_1_1_2909.tif, 060907_17_3_5_2_1_2910.tif, 060907_17_3_5_3_1_2911.tif)

Re: [R] What's the corresponding function in R for lo() function in S-PLUS?

2007-09-17 Thread Gavin Simpson
On Tue, 2007-09-18 at 01:19 +0800, zhijie zhang wrote: Dear friends, In S-PLUS, we can use the following argument, but not in R. mode12 - gam(score1 ~ lo(latitude) + lo(longitude)) I searched the help in S-PLUS, it says lo() Allows the user to specify a Loess fit in a GAM formula, but i

Re: [R] How to produce a macron symbol (overline / overbar) on a plot

2007-09-17 Thread Vladimir Eremeev
xlab=expression(bar(temp)) or xlab=expression(bar( temp)) The overbar looks like it begins above e in the first case. See demo(plotmath). Jason Horn wrote: Can anyone tell me how to produce a macron (or overbar or overline) symbol on an R plot. I have an axis that is labeled temp, but

Re: [R] How to produce a macron symbol (overline / overbar) on a plot

2007-09-17 Thread Marc Schwartz
On Mon, 2007-09-17 at 14:00 -0400, Jason Horn wrote: Can anyone tell me how to produce a macron (or overbar or overline) symbol on an R plot. I have an axis that is labeled temp, but I need the overbar symbol over the entire word. See ?plotmath and the examples therein. Example:

[R] Problems with compiling rimage on macOS

2007-09-17 Thread Thomas Heine
hello, i want to install the rimage package (0.5-7) into my R (2.5.1). i installed the i686-apple-darwin8-g++-4.0.1 (GCC) 4.0.1 (Apple Computer, Inc. build 5367) from the macOS cd and created symbolic links to gcc and g++. well if i try now to install the package in R the console says: g++

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 the

[R] removing a specific number of digist from a character string

2007-09-17 Thread Kim Milferstedt
Hello, I would like to remove the last 8 digists of character strings in a vector. Below I added a couple of elements from that vector. I have a problem defining a pattern to replace the digits using for example sub. Removing the .tif part works fine using sub('.tif',,x), but how do I get rid

[R] What's the corresponding function in R for lo() function in S-PLUS?

2007-09-17 Thread zhijie zhang
Dear friends, In S-PLUS, we can use the following argument, but not in R. mode12 - gam(score1 ~ lo(latitude) + lo(longitude)) I searched the help in S-PLUS, it says lo() Allows the user to specify a Loess fit in a GAM formula, but i didn't find the correponding function in R. Anybody knows how

Re: [R] MAD

2007-09-17 Thread Deepayan Sarkar
On 9/17/07, Nair, Murlidharan T [EMAIL PROTECTED] wrote: I am calculating the median absolute deviation using mad function, and it tends to ignore the parameter constant=1, when I am calculating it for x=seq(1:5). Am I missing something here? x-seq(1:5) mad(x)# gives [1] 1.4826 mad(x,

Re: [R] applying math/stat functions to rows in data frame

2007-09-17 Thread Gerard Smits
Hi all, Thanks for the suggestions. I have not yet tried the apply () approach, but have tried to get the indexed version working, so far with limited success. I realize that a transpose, as suggested, would work, but want to avoid that for something simpler. To repeat, the task is to

[R] help for high-quality plot again

2007-09-17 Thread Zheng Lu
Dear all: Thank you very much for your help. Actually, I save the plot(mfrow=c(2,3))as metafile, however, when I insert picture from the file in word to import this plot picture, everything seems to be shrinked. Is there anyone having good experience with import graph from R to Word document?

[R] map issues

2007-09-17 Thread Alexander Nervedi
Hi ! I have a shapefile that I can easily read into R using library(maptools). My problem stems from some warning messages that come even though everything seems to work fine. library(maptools) districts - read.shape(filen = a_ds, dbf.data = TRUE) length(districts$Shapes) so far so good.