Re: [R] A problem with svIDE in Tinn-R?

2007-04-30 Thread Philippe Grosjean
Hello, I konw I have some corrections to do in svIDE,... I simply don't have the time right now. Will try next week. Thank for reporting the problem. Best, Philippe Grosjean ..°})) ) ) ) ) ) ( ( ( ( (Prof. Philippe Grosjean ) ) ) ) )

Re: [R] Converting list of data frame to data frame

2007-04-30 Thread Ajit Pawar
Greetings Douglas/Jim/John/R-help, Thanks for your help so far. Answering your questions - doing an str' on the list reveals following (head) information : --- $ : Factor w/ 729 levels XX1,YY1,..: 6 9 10 12 13 14 19 22 29 30 ... $ : int [1:109] 19950201

Re: [R] Inflate/Expand/Resize an array

2007-04-30 Thread Mark Payne
Gudday, Robin's reply obviously set something in motion in my brain, because I woke up with the answer on Saturday morning (which was actually pretty frustrating because I had to wait until Monday morning to try it out!) - essentially create a dummy array of 1s, and then use sweep to fill it

[R] - ICA for less common data

2007-04-30 Thread Monica Pisica
Hi everyone, I am not sure this is the appropriate list I should put this question to, but I hope you will re-direct me to the most appropriate one if necessary. I am doing an independent component analysis on a dataset that represents different metrics for patchreefs such as depth, area,

[R] general question about use of list

2007-04-30 Thread Y G
Hi, is this list only related to R issues or it has a broader context regarding questions and discussions about statistics. Is there any other email list or forum for that? For example, I have a question regarding variance. It is defined as: variance = sum(sq(Xi-mean)) / (N-1) and I never

[R] S-plus coding

2007-04-30 Thread T . Kounouni
Hi, i would like to ask, what is the Trigg monitor and the CUSUM monitor? Thank you __ R-help@stat.math.ethz.ch 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

Re: [R] general question about use of list

2007-04-30 Thread Leeds, Mark \(IED\)
it's also not unbiased. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Y G Sent: Monday, April 30, 2007 8:15 AM To: r-help@stat.math.ethz.ch Subject: [R] general question about use of list Hi, is this list only related to R issues or it has a broader

Re: [R] weight

2007-04-30 Thread Natalie O'Toole
Hi, These are the variables in my file. I think the variable i'm having problems with is WTPP which is of the Factor type. Does anyone know how to fix this, please? Thanks, Nat data.frame': 290 obs. of 5 variables: $ PROV : num 48 48 48 48 48 48 48 48 48 48 ... $ REGION: num 4 4 4 4

Re: [R] Converting list of data frame to data frame

2007-04-30 Thread jim holtman
This does not look like a list of dataframes where each dataframe is an element of the list. This appears to just be a list with all the elements of the data frames at the same levels. A list of data frames would look like this: x - list(data.frame(a=1:5, b=1:5), data.frame(a=6:10, b=6:10))

Re: [R] general question about use of list

2007-04-30 Thread gatemaze
On 30/04/07, Y G [EMAIL PROTECTED] wrote: Hi, is this list only related to R issues or it has a broader context regarding questions and discussions about statistics. Is there any other email list or forum for that? Well, just to quickly self-reply my question. From

[R] Query about using rowSums/ColSums on table results

2007-04-30 Thread lalitha viswanath
Hi I have data of the form class age A 0.5 B 0.4 A 0.5 C 0.785 D 0.535 A 0.005 C 0.015 D 0.205 A 0.605 etc etc... I tabulated the above as tab -table(data$class, cut(data$age, seq(0,0.6,0.02)) I wish to view the results in individual bins as a

[R] Definition of Hinges in Boxplot

2007-04-30 Thread fredrik odegaard
Hi All, Is the definition and/or algorithm for the 'Hinges' (upper and lower) in Boxplot available somewhere online? I would like to know a little more than what is described in boxplot.stats. In particular I am trying to understand exactly how it differs from the default (Type 7) of the quantile

[R] Definition of Hinges in Boxplot

2007-04-30 Thread fredrik odegaard
Hi all, sorry if this posting appears twice, but I seemed to have some trouble on my first attempt. I am wondering if the definition and/or algorithm for the 'Hinges' (upper and lower) in Boxplot available somewhere online? I would like to know a little more than what is described in

Re: [R] general question about use of list

2007-04-30 Thread John Kane
You might find the usenet groups sci.stat.math sci.stat.consult to be what you're looking for. --- Y G [EMAIL PROTECTED] wrote: Hi, is this list only related to R issues or it has a broader context regarding questions and discussions about statistics. Is there any other email list or

[R] Definition of Hinges in Boxplot

2007-04-30 Thread fredrik odegaard
Hi all, sorry if this posting appears twice, but I seemed to have some trouble on my first attempt. I am wondering if the definition and/or algorithm for the 'Hinges' (upper and lower) in Boxplot available somewhere online? I would like to know a little more than what is described in

[R] Definition of Hinges in Boxplot

2007-04-30 Thread fredrik odegaard
Hi all, sorry if this posting appears twice, but I seemed to have some trouble on my first attempt. I am wondering if the definition and/or algorithm for the 'Hinges' (upper and lower) in Boxplot available somewhere online? I would like to know a little more than what is described in

[R] thousand separator (was RE: weight)

2007-04-30 Thread Liaw, Andy
I've run into this occasionally. My current solution is simply to read it into Excel, re-format the offending column(s) by unchecking the thousand separator box, and write it back out. Not exactly ideal to say the least. If anyone can provide a better solution in R, I'm all ears... Andy

Re: [R] weight

2007-04-30 Thread Natalie O'Toole
__ Hi, These are the variables in my file. I think the variable i'm having problems with is WTPP which is of the Factor type. Does anyone know how to fix this, please? Thanks, Nat data.frame': 290 obs. of 5 variables: $ PROV : num 48 48 48 48 48 48 48 48 48 48

Re: [R] Query about using rowSums/ColSums on table results

2007-04-30 Thread Gabor Grothendieck
Try prop.table(tab, 1) prop.table(tab, 2) prop.table(tab) for the three ways of taking fractions or see ?CrossTable in the gmodels package. On 4/30/07, lalitha viswanath [EMAIL PROTECTED] wrote: Hi I have data of the form class age A 0.5 B 0.4 A 0.5 C 0.785 D

[R] R GUI in Ubuntu Feisty?

2007-04-30 Thread Andy Weller
Dear all, After an update from Ubuntu Edgy to Feisty, I seem to have lost package JGR()!? I have updated my sources.list to point to the Feisty repos at http://cran.ch.r-project.org/ and re-installed JGR() via: $ sudo rm -rf /usr/local/lib/R/site-packages/* $ sudo R CMD javareconf $ sudo R

Re: [R] thousand separator (was RE: weight)

2007-04-30 Thread Marc Schwartz
One possibility would be to use something like the following post-import: WTPP [1] 1,106.8250 1,336.5138 str(WTPP) Factor w/ 2 levels 1,106.8250,1,336.5138: 1 2 as.numeric(gsub(,, , WTPP)) [1] 1106.825 1336.514 Essentially strip the ',' characters from the factors and then coerce the

Re: [R] thousand separator (was RE: weight)

2007-04-30 Thread Bert Gunter
Except this doesn't work for 1,123,456.789 Marc. I hesitate to suggest it, but gregexpr() will do it, as it captures the position of **every** match to ,. This could be then used to process the vector via some sort of loop/apply statement. But I think there **must** be a more elegant way using

Re: [R] thousand separator (was RE: weight)

2007-04-30 Thread Marc Schwartz
Bert, What am I missing? print(as.numeric(gsub(,, , 1,123,456.789)), 10) [1] 1123456.789 FWIW, this is using: R version 2.5.0 Patched (2007-04-27 r41355) Marc On Mon, 2007-04-30 at 10:13 -0700, Bert Gunter wrote: Except this doesn't work for 1,123,456.789 Marc. I hesitate to suggest

[R] ROracle issues

2007-04-30 Thread Matt Anthony
Hi all - I am trying to install ROracle for linux machines ... I have read the INSTALL documentation and followed the directions for setting the paths as follows: export PATH=/app/oracle/product/10.2.0/db_1/bin:$PATH [EMAIL PROTECTED] ~]$ export

[R] thousand separator (was RE: weight)

2007-04-30 Thread Natalie O'Toole
Thank-you Andy!! That works great now Nat __ I've run into this occasionally. My current solution is simply to read it into Excel, re-format the offending column(s) by unchecking the thousand separator box, and write it back out. Not exactly ideal to say the least. If

[R] About autocorrelation

2007-04-30 Thread Atte Tenkanen
Hello R-users, I'd like to know if there is an autocorrelation function for cases when the time differences between the observations vary? I want to give both the time points (x) and values (y) for the function and get as an output an estimation of autocorrelation. Atte Tenkanen University of

Re: [R] thousand separator (was RE: weight)

2007-04-30 Thread Bert Gunter
Nothing! My mistake! gsub -- not sub -- is what you want to get 'em all. -- Bert Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Marc Schwartz Sent: Monday, April 30, 2007 10:18 AM To: Bert Gunter Cc:

Re: [R] thousand separator (was RE: weight)

2007-04-30 Thread Liaw, Andy
Still, though, it would be nice to have the data read in correctly in the first place, instead of having to do this kind of post-processing afterwards... Andy From: Bert Gunter Nothing! My mistake! gsub -- not sub -- is what you want to get 'em all. -- Bert Bert Gunter Genentech

[R] Read every second line from ASCII file

2007-04-30 Thread Roland Rau
Dear all, I have an ASCII file where records are separated by a blank. I would like to read those data; however, only the data in rows 1, 3, 5, 7, ... are important; the other lines (2,4,6,8,) contain no useful information for me. So far I used awk/gawk to do it: gawk '{if ((FNR % 2) != 0)

Re: [R] Read every second line from ASCII file

2007-04-30 Thread Ferdinand Alimadhi
Hi, You can start by reading all lines from your file: lines - readLines(pathtoyourfile) then keep only odd lines: oddlines - lines[seq (1, length(lines),2)] You have to split the line into fields, e.g res - strsplit(oddlines, split = \t) if you have tab as field seperator hth On

Re: [R] Read every second line from ASCII file

2007-04-30 Thread jim holtman
Use readLines and then just use the odd numbered lines: x.in - readLines(yourFile) x.in - x.in[seq(1, length(x.in), 2)] # every 2nd line or just to make sure, only delete blank lines: x.in - x.in[!(x.in == )] On 4/30/07, Roland Rau [EMAIL PROTECTED] wrote: Dear all, I have an ASCII file

[R] Intercept Coefficient in a Model with Orthogonal Polynomials

2007-04-30 Thread Chuck Cleland
This very likely falls in the category of an unexpected result due to user ignorance. I generated the following data: time - 0:10 set.seed(4302007) y - 268 + -9*time + .4*(time^2) + rnorm(11, 0, .1) I then fit models using both orthogonal and raw polynomials: fit1 - lm(y ~ poly(time, 2))

[R] intersect filled.contour and polygon

2007-04-30 Thread Pedro Mardones
Dear R users; Is there any way to intersect a filled contour image and a polygon? My problem is that I want to create a kriging map and the boundaries of my map are given by the coordinates of the polygon. So far I can superompose the polygon in the filled.contour image but I don't know how to

Re: [R] thousand separator (was RE: weight)

2007-04-30 Thread Gabor Grothendieck
That could be accomplished using a custom class like this: library(methods) setClass(num.with.junk) setAs(character, num.with.junk, function(from) as.numeric(gsub(,, , from))) ### test ### Input - A B 1,000 1 2,000 2 3,000 3 DF - read.table(textConnection(Input), header = TRUE,

[R] Independent contrasts from lme with interactions

2007-04-30 Thread Ken Nussear
Hi All, I've been searching the help archives but haven't found a workable solution to this problem. I'm running an lme model with the following call: lme.fnl - lme(Max ~ S + Tr + Yr + Tr:Yr, random = ~1 |TID) anova(lme.fnl) numDF denDF F-value p-value (Intercept) 1

Re: [R] Independent contrasts from lme with interactions

2007-04-30 Thread Andrew Robinson
Ken, estimable in the gmodels package will help you. Cheers Andrew On Mon, Apr 30, 2007 at 03:11:57PM -0700, Ken Nussear wrote: Hi All, I've been searching the help archives but haven't found a workable solution to this problem. I'm running an lme model with the following call:

Re: [R] Perpendicular symbol in plotmath?

2007-04-30 Thread Paul Murrell
Hi Matthew Neilson wrote: Thanks for your response, Gabor. That works quite nicely. The documentation states that it is not possible to mix and match Hershey fonts with plotmath symbols. My *ideal* scenario would be to write the perpendicular symbol as a subscript (specifically, I

Re: [R] Perpendicular symbol in plotmath?

2007-04-30 Thread Gabor Grothendieck
On 4/30/07, Paul Murrell [EMAIL PROTECTED] wrote: though I had to switch to a single-byte encoding to get my system to pick up the symbol font Exactly what command did you use to do that? __ R-help@stat.math.ethz.ch mailing list

Re: [R] Perpendicular symbol in plotmath?

2007-04-30 Thread Paul Murrell
Hi Gabor Grothendieck wrote: On 4/30/07, Paul Murrell [EMAIL PROTECTED] wrote: though I had to switch to a single-byte encoding to get my system to pick up the symbol font Exactly what command did you use to do that? Sorry, that should have read switch to a single-byte locale, which for

[R] JGR starter problems

2007-04-30 Thread J Dougherty
I'm running R 2.4.1 on SuSE Linux 10.2. My system is an AMD based PC with 2 Gigs of ram and abundant HD space. I have always run R from the console without problems, but in the interests of broadening R's user-base, and reducing the complexity of my computing environment, I am hoping to

[R] Unsubscription Confirmation

2007-04-30 Thread Ultimate Wealth
Thank you for subscribing. You have now unsubscribed and no more messages will be sent. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] thousand separator (was RE: weight)

2007-04-30 Thread Liaw, Andy
Looks very neat, Gabor! I just cannot fathom why anyone who want to write numerics with those separators in a flat file. That's usually not for human consumption, and computers don't need those separators! Andy From: Gabor Grothendieck That could be accomplished using a custom class

[R] to draw a smooth arc

2007-04-30 Thread Paulo Barata
Dear R-list members, I would like to draw a smooth arc. I can draw an arc parametrically, but this produces an arc too coarse, even allowing for different increments in sequence t in the example below. Function symbols (graphics) does produce a smooth circle, but it cannot produce an arc.

[R] bivariate normal distribution in likelihood

2007-04-30 Thread DEEPANKAR BASU
Hi all, I am trying to do a maximum likelihood estimation. In my likelihood function, I have to evaluate a double integral of the bivariate normal density over a subset of the points of the plane. For instance, if I denote by y the y co-ordinate and by x, the x co-ordinate then the area over

[R] Concepts question: environment, frame, search path

2007-04-30 Thread Graham Wideman
Folks: I'd appreciate if someone could straighten me out on a few concepts which are described a bit ambiguously in the docs. 1. data.frame: Refan p84: 'A data frame is a list of variables of the same length with unique row names, given class data.frame.' I probably don't

Re: [R] bivariate normal distribution in likelihood

2007-04-30 Thread Robert A LaBudde
At 11:32 PM 4/30/2007, Deepankar wrote: Hi all, I am trying to do a maximum likelihood estimation. In my likelihood function, I have to evaluate a double integral of the bivariate normal density over a subset of the points of the plane. For instance, if I denote by y the y co-ordinate and by