[R] climatological standard deviation

2008-09-24 Thread Yogesh Tiwari
Hello R users, I have a montly time series over a several year period. It's easy to compute a monthly climatology (12 values), Now, I would like to calculate the corresponding standard deviation, ie the 12 values calculated from the january values, february values, etc. What's the best way for

Re: [R] How to draw the graph of f(x,y) = x * y ?

2008-09-24 Thread Robin Hankin
Paul you might find the view() function in the 'elliptic' package useful. This function implements various methods to visualize functions over the complex plane. HTH rksh Paul Smith wrote: Dear All, The function curve() draws the graph of functions from R to R. Is there some homologous

Re: [R] Options for viewing / using results from lm

2008-09-24 Thread Michael Just
Hello, Thanks for the response. However, I have this problem with view. View(fit1) Error: could not find function View I have R 2.4.1 and even tried loading 'utils'. As for saving residuals in a spreadsheet I used something like this: res.fit1-fit1$residuals res.fit2-fit2$residuals

[R] Antwort: Re: WG: Problem during porting R-2.7.2 on HP-UP 11.11 PA-Risc

2008-09-24 Thread claus-juergen . neumann
Dear Peter, thank you so much for the information and quick help, this was my problem to get the configure run. please see the results: R is now configured for hppa2.0w-hp-hpux11.11 Source directory: . Installation directory:/opt/R C compiler:gcc -std=gnu99

Re: [R] cannot allocate memory

2008-09-24 Thread Uwe Ligges
DumpsterBear wrote: I am getting Error: cannot allocate vector of size 197 MB. I know that similar problems were discussed a lot already, but I didn't find any satisfactory answers so far! Details: *** I have XP (32bit) with 4GB ram. At the time when the problem appeared I had 1.5GB of

Re: [R] Antwort: Re: WG: Problem during porting R-2.7.2 on HP-UP 11.11 PA-Risc

2008-09-24 Thread Peter Dalgaard
[EMAIL PROTECTED] wrote: Make: Don't know how to make #. Stop. *** Error exit code 1 Sounds like you are using a make which doesn' t understand comments (??!). You might need to install GNU make and possibly set MAKE=gmake or some such to make sure that you are using the right one. At least

[R] Logistic regression

2008-09-24 Thread Annie Bissonnette
Dear all, I am currently learning to run logistic regression models with R. Would someone tell me what this exactly means: Estimated scale (compare to 1 ) 1.746724 If the value is higher or lower than 1, what should I do? The complete results of the model were as follows: Generalized linear

Re: [R] Bug in is ?

2008-09-24 Thread Keith Jewell
7 is an integer, but it's also a real. In R '?is' and '?is.integer' are clear that you're testing the class(es) of objects, not their values. I can't comment on the relationship with S Programming hth Keith J [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] This is really

[R] match.call in a function

2008-09-24 Thread Camarda, Carlo Giovanni
Dear UseRs, please find below a simple example which exemplifies my question. In words, I'd like to write a function which prints formula and outcomes only when results are not assigned or recalling the already assigned object, as e.g. glm. I'm sure a solution is rather simple and

Re: [R] cannot allocate memory

2008-09-24 Thread Bernardo Rangel Tura
Em Ter, 2008-09-23 às 21:42 -0400, DumpsterBear escreveu: I am getting Error: cannot allocate vector of size 197 MB. I know that similar problems were discussed a lot already, but I didn't find any satisfactory answers so far! Details: *** I have XP (32bit) with 4GB ram. At the time when

Re: [R] Modifying output to Google Docs

2008-09-24 Thread Jim Lemon
Ajay ohri wrote: Hi List, Graphical output to PDF's ,RTF ,CSV is known through R. Can it be modified for outputting to Google Docs (which is basically uploaded files ,published to become html pages) Is there any package on this ? Hi Ajay, Do you mean transforming R output into HTML? The

Re: [R] Modifying output to Google Docs

2008-09-24 Thread Ajay ohri
Do I mean R output to Google docs, like in Google spreadsheets which are different from a spreadsheet published as an html or Presentations in Google Docs, which are much different tahn ordinary html presentations. See examples here of Google Docs-

Re: [R] R2HTML: output from for-loops

2008-09-24 Thread Jim Lemon
Werner Wernersen wrote: Hi, I am trying to prepare a report with R2HTML using HTMLStart(outdir=./html, filename=report, echo=T, HTMLframe=F) then, for instance, I want to get the output of a loop: for (i in 1:20) print(summary(rnorm(1000))) but only the first of summaries really ends up in

[R] extract species in a phylog tree

2008-09-24 Thread Christine N. Meynard
Hi, I am working with a phylog tree and I would like to extract a subset of the tree based on the species names (conserving the evolutionary distance and relationships between the pairs of species I am interested in). I see there is an option to select a subset of the tree using node names

Re: [R] read.table readLines behaviour?

2008-09-24 Thread J . delasHeras
Quoting Peter Dalgaard [EMAIL PROTECTED]: [EMAIL PROTECTED] wrote: Hi, I have been using 'read.table' regularly to read tab-delimited text files with data. No problem, until now. Now I have a file that appeared to have read fine, and the data inside looks correct (structure etc), except I

Re: [R] read.table readLines behaviour?

2008-09-24 Thread J . delasHeras
The result was 11, 24001 times, as I expected originally. hmmm... JOse Quoting Gabor Grothendieck [EMAIL PROTECTED]: Try looking at the result of count.fields to diagnose it. On Tue, Sep 23, 2008 at 5:19 AM, [EMAIL PROTECTED] wrote: Hi, I have been using 'read.table' regularly to read

[R] lattice: Color in Barchart legend

2008-09-24 Thread Patrick Hausmann
Dear list, with the code below I produce the right graph, but the colours of the legend are different from the colours of the graph. The colours of the graph are the desired colours. Thanks for any help. Patrick library(lattice) pal1 - rgb(196, 255, 255, max = 255) pal2 - rgb( 0, 35,

Re: [R] How to convert backslash to slash?

2008-09-24 Thread Duncan Murdoch
Shengqiao Li wrote: On Tue, 23 Sep 2008, Duncan Murdoch wrote: On 23/09/2008 4:00 PM, Shengqiao Li wrote: How to use sub, gsub, etc. to replace \ in a string to /? For example,convert C:\foo\bar to C:/foo/bar. If those are R strings, there are no backslashes in the first one.

Re: [R] lattice: Color in Barchart legend

2008-09-24 Thread baptiste auguie
Hi, you can use par.settings to get a consistent color scheme, pal1 - rgb(196, 255, 255, max = 255) pal2 - rgb( 0, 35, 196, max = 255) df - data.frame( Gruppe = c(A, B, A, B), Kat = c(x1, x1, w1, w1), value= c(1,2, 4, 5)) barchart(value ~ Kat, group= Gruppe,

[R] hexbin object to vector (shapefile)

2008-09-24 Thread Paulo Cardoso
Is it possible to convert an hexbin object like hexbin(rep(1:100,each=100),rep(1:100,100),xbins=10) and convert it to a vectorial (shp) format? Paulo - This electronic mail

[R] Changing a plot

2008-09-24 Thread R Help
Hello list, I've been working on this problem for a while and I haven't been able to come up with a solution. I have a couple of functions that plot a bunch of data, then a single point on top of it. What I want is to be able to change the plot of the point without replotting all the data.

Re: [R] Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness?

2008-09-24 Thread Gabor Grothendieck
Check out JohnsonFit in SuppDists. On Tue, Sep 23, 2008 at 10:59 AM, zhijie zhang [EMAIL PROTECTED] wrote: Dear R users, I hope to explain the concepts of skewness and kurtosis by generating series of distributions with same skewness and different kurtosis or with same kurtosis and

Re: [R] Bug in is ?

2008-09-24 Thread Martin Maechler
KJ == Keith Jewell [EMAIL PROTECTED] on Wed, 24 Sep 2008 09:46:08 +0100 writes: KJ 7 is an integer, but it's also a real. KJ In R '?is' and '?is.integer' are clear that you're testing the class(es) of KJ objects, not their values. KJ I can't comment on the relationship

Re: [R] How to convert backslash to slash?

2008-09-24 Thread Shengqiao Li
On Wed, 24 Sep 2008, Duncan Murdoch wrote: Shengqiao Li wrote: On Tue, 23 Sep 2008, Duncan Murdoch wrote: On 23/09/2008 4:00 PM, Shengqiao Li wrote: How to use sub, gsub, etc. to replace \ in a string to /? For example,convert C:\foo\bar to C:/foo/bar. If those are R strings, there are

Re: [R] How to draw the graph of f(x,y) = x * y ?

2008-09-24 Thread Paul Smith
On Wed, Sep 24, 2008 at 4:05 AM, Ben Bolker [EMAIL PROTECTED] wrote: The function curve() draws the graph of functions from R to R. Is there some homologous function to curve() to draw functions from R^2 to R? There is a curve3d function in the emdbook package on CRAN. Thanks, Ben and

Re: [R] How to convert backslash to slash?

2008-09-24 Thread jim holtman
It probably depends on how you are prompting the user for input. Here is what happens with 'readline' and you will see that the string after input does have the backslashes escaped. x - readline(Input File: ) Input File: c:\dir1\dir2\file.name x [1] c:\\dir1\\dir2\\file.name gsub(, /, x)

[R] rgl: ellipse3d with axes

2008-09-24 Thread Michael Friendly
Last week I asked about data ellipses with rgl:::ellipse3d() with lines showing the principal axes. (The goal is a visual demonstration of PCA as a rotation of variable space to component space.) I was trying, unsuccessfully, to use princomp() to generate the PCA axes and plot them using

Re: [R] How to convert backslash to slash?

2008-09-24 Thread Duncan Murdoch
Shengqiao Li wrote: On Wed, 24 Sep 2008, Duncan Murdoch wrote: Shengqiao Li wrote: On Tue, 23 Sep 2008, Duncan Murdoch wrote: On 23/09/2008 4:00 PM, Shengqiao Li wrote: How to use sub, gsub, etc. to replace \ in a string to /? For example,convert C:\foo\bar to

[R] Bug in persp?

2008-09-24 Thread Paul Smith
Dear All, I am not sure whether the following can be considered a bug: x - seq(-1,1,length=20) y - x z - matrix(1,20,20) persp(x,y,z) Error in persp.default(x, y, z) : invalid 'z' limits It works with: persp(x,y,z,zlim=c(0,1.5)) Should not persp set the z limits by default when the matrix

Re: [R] bandwidth selection for locpoly

2008-09-24 Thread Liaw, Andy
Don't think it's in kernSmooth, but locfit has such facility. See Loader's book for details. Best, Andy -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Chinthaka Kuruwita Sent: Tuesday, September 23, 2008 12:18 PM To: r-help@r-project.org

[R] qplot, stacked area, own colourscheme

2008-09-24 Thread Albin Blaschka
Dear list, I would like to apply my own colours to a stacked area plot, done with qplot, but I have not succeeded... What do I have so far (I am dealing with the development of cover of specific groups of plants): library(ggplot2) library(RODBC) channel - odbcConnect(myusername,

Re: [R] How to convert backslash to slash?

2008-09-24 Thread Shengqiao Li
On Wed, 24 Sep 2008, jim holtman wrote: It probably depends on how you are prompting the user for input. Here is what happens with 'readline' and you will see that the string after input does have the backslashes escaped. x - readline(Input File: ) Input File: c:\dir1\dir2\file.name x

[R] time calculation for simulation

2008-09-24 Thread Ralph Scherer
Hello list, I want to compare time needed for a simulation between R and Ra (time optimized loops). Does anybody know, how I can count time needed for a simulation in R? Bests, Ralph __ R-help@r-project.org mailing list

Re: [R] qplot, stacked area, own colourscheme

2008-09-24 Thread hadley wickham
On Wed, Sep 24, 2008 at 8:39 AM, Albin Blaschka [EMAIL PROTECTED] wrote: Dear list, I would like to apply my own colours to a stacked area plot, done with qplot, but I have not succeeded... What do I have so far (I am dealing with the development of cover of specific groups of plants):

Re: [R] time calculation for simulation

2008-09-24 Thread jim holtman
?system.time On Wed, Sep 24, 2008 at 9:46 AM, Ralph Scherer [EMAIL PROTECTED] wrote: Hello list, I want to compare time needed for a simulation between R and Ra (time optimized loops). Does anybody know, how I can count time needed for a simulation in R? Bests, Ralph

Re: [R] climatological standard deviation

2008-09-24 Thread milton ruser
Dear Yogesh, Send us a minimal reproducive code, with a sample of your data structure. It will help us to help you :-) Best wishes, miltinho astronauta brazil On Wed, Sep 24, 2008 at 2:00 AM, Yogesh Tiwari [EMAIL PROTECTED]wrote: Hello R users, I have a montly time series over a several

Re: [R] rgl: ellipse3d with axes

2008-09-24 Thread Duncan Murdoch
Michael Friendly wrote: Last week I asked about data ellipses with rgl:::ellipse3d() with lines showing the principal axes. (The goal is a visual demonstration of PCA as a rotation of variable space to component space.) I was trying, unsuccessfully, to use princomp() to generate the PCA axes

[R] rowSums()

2008-09-24 Thread Doran, Harold
Say I have the following data: testDat - data.frame(A = c(1,NA,3), B = c(NA, NA, 3)) testDat A B 1 1 NA 2 NA NA 3 3 3 rowsums() with na.rm=TRUE generates the following, which is not desired: rowSums(testDat[, c('A', 'B')], na.rm=T) [1] 1 0 6 rowsums() with na.rm=F generates the

Re: [R] How to convert backslash to slash?

2008-09-24 Thread jim holtman
Your users will have to learn that if they are inputting quoted strings into R, then, by convention, a backslash is used to 'escape' certain character sequences and if you want a backslash, you have to escape it ('\\'). You can also have your users upgrade to a system that does not use

[R] Error results from MS Access via RODBC

2008-09-24 Thread stefan . petersson
I have an MS Access database with one table and one column holding rep(1:10) I use: library(RODBC) channel - odbcConnect(test) sqlQuery(channel, paste(SELECT col FROM tblTest), believeNRows=FALSE) and get: 49 50 51 52 53 54 55 56 57 12337 What? The above should, of course, read: 1

[R] lower / upper case letters in a plot

2008-09-24 Thread Jörg Groß
Hi, How can I generate lower case letters for my axis-titles? Thanks, Jörg __ 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

Re: [R] cannot allocate memory

2008-09-24 Thread DumpsterBear
I am getting Error: cannot allocate vector of size 197 MB. I know that similar problems were discussed a lot already, but I didn't find any satisfactory answers so far! Details: *** I have XP (32bit) with 4GB ram. At the time when the problem appeared I had 1.5GB of available physical

Re: [R] Changing a plot

2008-09-24 Thread Ben Bolker
R Help rhelp.stats at gmail.com writes: Hello list, I've been working on this problem for a while and I haven't been able to come up with a solution. I have a couple of functions that plot a bunch of data, then a single point on top of it. What I want is to be able to change the plot

Re: [R] rgl: ellipse3d with axes

2008-09-24 Thread Michael Friendly
Duncan Murdoch wrote: The normals component contains the surface normals. It is used to help in rendering the surface, but isn't much use for your purposes. Unfortunately, I'm not familiar enough with the internals of princomp to tell you how to get the axes you want. Duncan Murdoch OK,

Re: [R] qplot, stacked area, own colourscheme

2008-09-24 Thread Albin Blaschka
hadley wickham schrieb: On Wed, Sep 24, 2008 at 8:39 AM, Albin Blaschka [EMAIL PROTECTED] wrote: [snip] As I have those categories (cover grasses, cover herbs...) I would like to assign special colours to each (for example cover grasses = dark green, cover open soil = brown ...), but I did

Re: [R] lower / upper case letters in a plot

2008-09-24 Thread Dimitris Rizopoulos
check: ?tolower() Best, Dimitris Jörg Groß wrote: Hi, How can I generate lower case letters for my axis-titles? Thanks, Jörg __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] rowSums()

2008-09-24 Thread Adaikalavan Ramasamy
I guess this would be the fastest way would be: rs - rowSums( testDat, na.rm=T) rs[ which( rowMeans(is.na(testDat)) == 1 ) ] - NA since both rowSums and rowMeans are internally coded in C. Regards, Adai Doran, Harold wrote: Say I have the following data: testDat - data.frame(A =

Re: [R] Changing a plot

2008-09-24 Thread Adaikalavan Ramasamy
One way is to keep a copy of the original and then return to it when you need it. x - rnorm(100,1,0.5) y - rnorm(100,1,0.5) plot(x,y,pch=16) original - recordPlot() for( i in 1:10 ){ points( x[i], y[i], pch=19, col=yellow, cex=3) points( x[i], y[i], pch=16) Sys.sleep(1)

Re: [R] rowSums()

2008-09-24 Thread Chuck Cleland
On 9/24/2008 10:06 AM, Doran, Harold wrote: Say I have the following data: testDat - data.frame(A = c(1,NA,3), B = c(NA, NA, 3)) testDat A B 1 1 NA 2 NA NA 3 3 3 rowsums() with na.rm=TRUE generates the following, which is not desired: rowSums(testDat[, c('A', 'B')],

Re: [R] lower / upper case letters in a plot

2008-09-24 Thread Jörg Groß
Oh, ok ... perhaps it was the wrong description; I want lowered cases... In LaTeX I get them with x_1 for example. Or x^2 for upper cases. But how can I make this in R? Am 24.09.2008 um 16:17 schrieb Dimitris Rizopoulos: check: ?tolower() Best, Dimitris Jörg Groß wrote: Hi, How can I

Re: [R] How to convert backslash to slash?

2008-09-24 Thread Gabor Grothendieck
On Wed, Sep 24, 2008 at 9:59 AM, Shengqiao Li [EMAIL PROTECTED] wrote: On Wed, 24 Sep 2008, jim holtman wrote: It probably depends on how you are prompting the user for input. Here is what happens with 'readline' and you will see that the string after input does have the backslashes

Re: [R] lower / upper case letters in a plot

2008-09-24 Thread Dimitris Rizopoulos
so you mean subscripts and superscripts, then check: ?plotmath Best, Dimitris Jörg Groß wrote: Oh, ok ... perhaps it was the wrong description; I want lowered cases... In LaTeX I get them with x_1 for example. Or x^2 for upper cases. But how can I make this in R? Am 24.09.2008 um 16:17

Re: [R] rowSums()

2008-09-24 Thread Dimitris Rizopoulos
try the following: testDat - data.frame(A = c(1,NA,3), B = c(NA, NA, 3)) ind - rowSums(is.na(testDat)) == length(testDat) out - rowSums(testDat, na.rm = TRUE) out[ind] - NA out I hope it helps. Best, Dimitris Doran, Harold wrote: Say I have the following data: testDat - data.frame(A =

Re: [R] lower / upper case letters in a plot

2008-09-24 Thread Adaikalavan Ramasamy
An example would help. You generally control the titles using arguments like main, xlab, ylab, sub in the plotting functions or afterwards using title() function. You can get the upper/lower case using toupper()/tolower() functions. See help(par), help(title), help(tolower). Here is an

Re: [R] rowSums()

2008-09-24 Thread Marc Schwartz
on 09/24/2008 09:06 AM Doran, Harold wrote: Say I have the following data: testDat - data.frame(A = c(1,NA,3), B = c(NA, NA, 3)) testDat A B 1 1 NA 2 NA NA 3 3 3 rowsums() with na.rm=TRUE generates the following, which is not desired: rowSums(testDat[, c('A', 'B')],

Re: [R] rowSums()

2008-09-24 Thread Chuck Cleland
On 9/24/2008 10:38 AM, Marc Schwartz wrote: on 09/24/2008 09:06 AM Doran, Harold wrote: Say I have the following data: testDat - data.frame(A = c(1,NA,3), B = c(NA, NA, 3)) testDat A B 1 1 NA 2 NA NA 3 3 3 rowsums() with na.rm=TRUE generates the following, which is not desired:

Re: [R] *solved* lower / upper case letters in a plot

2008-09-24 Thread Jörg Groß
Great thanks and sorry for my bad english! Am 24.09.2008 um 16:32 schrieb Dimitris Rizopoulos: so you mean subscripts and superscripts, then check: ?plotmath Best, Dimitris Jörg Groß wrote: Oh, ok ... perhaps it was the wrong description; I want lowered cases... In LaTeX I get them

Re: [R] How to convert backslash to slash?

2008-09-24 Thread hadley wickham
I agree that that some sort of facility would be convenient. Creating latex output is another situation where not having to escape backslashes would be convenient. And regular expressions, of course. Hadley -- http://had.co.nz/ __

Re: [R] Counting character occurrences in data frame

2008-09-24 Thread Hutchinson,David [PYR]
Thanks Charles, ftable() works perfectly. -Original Message- From: Charles C. Berry [mailto:[EMAIL PROTECTED] Sent: Tuesday, September 23, 2008 5:06 PM To: Hutchinson,David [PYR] Cc: r-help@r-project.org Subject: Re: [R] Counting character occurrences in data frame See

[R] How do I Convert 1 to the number 1?

2008-09-24 Thread Tom La Bone
Is there an elegant way in R to change a number reported as a less-than number in text format, 1 for example, to the numeric equivalent 1? I have been trying to use as.numeric, but have not come up with anything clever yet. Tom -- View this message in context:

Re: [R] rgl: ellipse3d with axes

2008-09-24 Thread Duncan Murdoch
On 24/09/2008 10:12 AM, Michael Friendly wrote: Duncan Murdoch wrote: The normals component contains the surface normals. It is used to help in rendering the surface, but isn't much use for your purposes. Unfortunately, I'm not familiar enough with the internals of princomp to tell you how to

Re: [R] How to draw the graph of f(x,y) = x * y ?

2008-09-24 Thread Farley, Robert
You could always try F(x,y) = f(x) + 0*y That is zero out the degenerate dimensions. Of course you'll be plotting what is essentially a two dimensional object as if it were three dimensional. The degeneracy in y means a 2-D curve will be extruded along the Y dimension. Robert Farley Metro

[R] looping through variables

2008-09-24 Thread K. Fleischer
Hello everyone, I have the following problem: My analysis includes many predictor variables (50) in the form of raster maps (asc), but I am trying to avoid having to type all their names over and over again in the analysis (e.g. for vectorisation, for deletion of NA's, etc.) So ideally I

Re: [R] How to convert backslash to slash?

2008-09-24 Thread Duncan Murdoch
On 24/09/2008 10:53 AM, hadley wickham wrote: I agree that that some sort of facility would be convenient. Creating latex output is another situation where not having to escape backslashes would be convenient. And regular expressions, of course. There were proposals to do this last year (I

[R] keep the row indexes/names when do aggregate

2008-09-24 Thread zhihuali
Hi, R-users, If I have a data frame like this: x-data.frame(g=c(g1,g2,g1,g1,g2),v=c(1,7,3,2,8)) g v 1 g1 1 2 g2 7 3 g1 3 4 g1 2 5 g2 8 It contains two groups, g1 and g2. Now for each group I want the max v: aggregate(x$v,list(g=x$g),max) g x 1 g1 3 2 g2 8 Beautiful. But what if I want

Re: [R] rgl: ellipse3d with axes

2008-09-24 Thread Duncan Murdoch
On 24/09/2008 11:11 AM, Duncan Murdoch wrote: On 24/09/2008 10:12 AM, Michael Friendly wrote: Duncan Murdoch wrote: The normals component contains the surface normals. It is used to help in rendering the surface, but isn't much use for your purposes. Unfortunately, I'm not familiar enough

Re: [R] Weights for polr

2008-09-24 Thread Gregory Wawro
Thanks for directing my attention to the survey package and svyolr(). Best, Greg . Gregory Wawro [EMAIL PROTECTED] Associate Professor phone: 212-854-8540 Dept. of Political Science fax:212-222-0598 741

Re: [R] keep the row indexes/names when do aggregate

2008-09-24 Thread Gabor Grothendieck
Rather than aggregate, use order and duplicated as in this post: https://stat.ethz.ch/pipermail/r-help/2008-September/173139.html On Wed, Sep 24, 2008 at 11:21 AM, zhihuali [EMAIL PROTECTED] wrote: Hi, R-users, If I have a data frame like this:

[R] (nlme) Repeated measures with continuous covariate in lme

2008-09-24 Thread Robert Buitenwerf
Dear readers, I have a basic question about how to use lme for my design. I haven't been able to find an example in r-help that made it clear to me how to tackle this problem and unfortunately I also cannot get hold of Pinheiro Bates 2000. I hope someone can help. Data for the response

[R] Graph question

2008-09-24 Thread Georgina Sarah Humphreys
If I have a set of data comprising a list of numbers of eggs on mosquito guts that range from 1 to 157. How can I get R to draw a barchart of the distribution of the data (i.e. x axis= number of eggs on a gut, y axis=number of mosquitoes found with that number of eggs)? Data is listed below.

[R] climatological standard deviation- (question re-posted)

2008-09-24 Thread Yogesh Tiwari
Sorry for re-posting the question, I did not get any reply. Kindly reply please if any one can. ### Hello R users, I have a montly time series over a several year period. It's easy to compute a monthly climatology (12 values), Now, I would like to calculate the corresponding standard deviation,

[R] Simulations / repetitions help!

2008-09-24 Thread Lauren_G
Dear all, My question concerns using repetitions and simulations (loops?) in R. I am very new R user, so any help that can be offered would be greatly appreciated! I am using fitdistr() to determine the distribution of empirical univariate datasets, and ks.test to assess the goodness of fit.

[R] t tests/ANOVA

2008-09-24 Thread Georgina Sarah Humphreys
I have a set of data that comprises genome numbers in single eggs from three different parasite clones - 3D7, HB3, and MIX. I can draw a boxplot of the genome numbers for each clonefed but how do I carry out a t test or ANOVA to compare if the means are signifcantly different? (Data is listed

Re: [R] How do I Convert 1 to the number 1?

2008-09-24 Thread John Fox
Dear Tom, I don't know whether it qualifies as elegant or clever, but you should be able to use sub(), as in as.numeric(sub(, , c(1, -2, 1))) [1] 1 -2 1 I hope this helps, John On Wed, 24 Sep 2008 08:08:21 -0700 (PDT) Tom La Bone [EMAIL PROTECTED] wrote: Is there an elegant way in R

Re: [R] Bug in is ?

2008-09-24 Thread ctu
Thank you for all of you. Intuitively, 7 is an integer for people who live in this planet. It is just very difficult for me to believe that R does not think 7 is an integer but 7L is. is.integer(7) # R 2.7.2 [1] FALSE Thus, based on Martin's comments, I try it again on the S-PLUS 8.0 and

Re: [R] climatological standard deviation- (question re-posted)

2008-09-24 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Yogesh Tiwari Sent: Wednesday, September 24, 2008 5:43 AM To: [EMAIL PROTECTED] Subject: [R] climatological standard deviation- (question re-posted) Sorry for re-posting the question, I did not get

Re: [R] climatological standard deviation- (question re-posted)

2008-09-24 Thread Sarah Goslee
On Wed, Sep 24, 2008 at 8:43 AM, Yogesh Tiwari [EMAIL PROTECTED] wrote: Sorry for re-posting the question, I did not get any reply. Kindly reply please if any one can. Most likely that's because you didn't clearly explain the problem. Please read the posting guide (see link at bottom of each

Re: [R] Why horizontal bars instead of a line

2008-09-24 Thread Sarah Goslee
Since x$GrSe isn't a number, what do you expect R to do? On Wed, Sep 24, 2008 at 11:28 AM, Schreiber, Stefan [EMAIL PROTECTED] wrote: Hey list! It looks simple, though it's not possible for me to plot the following properly: (some made-up data) GrSeClone1 Clone2 Clone3 Clone4

Re: [R] climatological standard deviation- (question re-posted)

2008-09-24 Thread stephen sefick
package zoo rollapply if I had to guess at what you want to do, but heed the advice given above. On Wed, Sep 24, 2008 at 12:11 PM, Sarah Goslee [EMAIL PROTECTED] wrote: On Wed, Sep 24, 2008 at 8:43 AM, Yogesh Tiwari [EMAIL PROTECTED] wrote: Sorry for re-posting the question, I did not get any

Re: [R] Graph question

2008-09-24 Thread stephen sefick
look at histogram and densityplot in the lattice package On Wed, Sep 24, 2008 at 9:26 AM, Georgina Sarah Humphreys [EMAIL PROTECTED] wrote: If I have a set of data comprising a list of numbers of eggs on mosquito guts that range from 1 to 157. How can I get R to draw a barchart of the

[R] weights option in lmer

2008-09-24 Thread maria jose juan jorda
Hi all, I am trying to run a linear mixed effect models in lmer() from the lme4 package using the weights option. I am using the R version 2.7.2 (2008-08-25) and lmer version in lme4_0.999375-26, which I think it is the latest version! I am getting and error message when I add the

Re: [R] rgl: ellipse3d with axes

2008-09-24 Thread Michael Friendly
Thanks Duncan ( others) Here is a function that does what I want in this case, and tries to do it to work generally with ellipse3d. (Note that I reverse the order of centre and scale 'cause I was bitten by trying ellipse3d.axes(cov, mu)) # draw axes in the data ellipse computed by ellipse3d

Re: [R] Bug in is ?

2008-09-24 Thread Keith Jewell
Have you tried is.integer(7.0) in S-Plus? (I have) Do you think 7.0 is integer? IMHO in R there is nothing to be fixed (in this regard) except your understanding. This is a computer language, not English; intuition isn't reliable, so we have help pages. is.integer(x) is not intended to

[R] logistic regression model validation through bootstrapping

2008-09-24 Thread Van Dyke, Alison
I'm new to R and to validation of logistic regression models via bootstrapping. I see that there are various approaches to bootstrapping, and I am wondering (1) which approach would be best for internal validation of predictive discrimination, (2) are there any good resources that could be

Re: [R] logistic regression model validation through bootstrapping

2008-09-24 Thread Frank E Harrell Jr
Van Dyke, Alison wrote: I'm new to R and to validation of logistic regression models via bootstrapping. I see that there are various approaches to bootstrapping, and I am wondering (1) which approach would be best for internal validation of predictive discrimination, (2) are there any good

Re: [R] rgl: ellipse3d with axes

2008-09-24 Thread Duncan Murdoch
On 24/09/2008 12:32 PM, Michael Friendly wrote: Thanks Duncan ( others) Here is a function that does what I want in this case, and tries to do it to work generally with ellipse3d. (Note that I reverse the order of centre and scale 'cause I was bitten by trying ellipse3d.axes(cov, mu)) #

Re: [R] How to draw the graph of f(x,y) = x * y ?

2008-09-24 Thread Paul Smith
On Wed, Sep 24, 2008 at 4:14 PM, Farley, Robert [EMAIL PROTECTED] wrote: You could always try F(x,y) = f(x) + 0*y That is zero out the degenerate dimensions. Of course you'll be plotting what is essentially a two dimensional object as if it were three dimensional. The degeneracy in y

[R] Fortunes Candidates from Keith Jewell?

2008-09-24 Thread Bert Gunter
IMHO Keith Jewell has a couple of jewels here (sorry Keith -- couldn't resist): nice possible fortunes package candidates: IMHO in R there is nothing to be fixed (in this regard) except your understanding. This is a computer language, not English; intuition isn't reliable, so we have help

Re: [R] Bug in is ?

2008-09-24 Thread ctu
Hi Keith, No doubt, 7.0 is integer in math. But if people can write 7 why people need to write 7.0 (I do not see any reason to do this). My point is that R maybe can do something like S-plus. No point to argue. don't you think so? Thanks Chunhao Quoting Keith Jewell [EMAIL PROTECTED]:

Re: [R] Options for viewing / using results from lm

2008-09-24 Thread Greg Snow
I don't remember which version View was added in, but it was more recent than 2.4.1. I would suggest upgrading (2.4.1 is ancient in R terms). You can get something similar with: invisible(edit(mydata)) But you will have to close the window to continue with R (View lets you look at the data

[R] lattice xyplot symbols instead of colors and legend matching plot symbols or colors

2008-09-24 Thread stephen sefick
I would like to use the data below where the plots are close to what I want. Instead of color I would like to use different symbols, and have the symbols in the legend match the graphs. I am also going to add a regression line to these I know about the type=r (which is fine for these particular

[R] possible interesting R projects for undergrads

2008-09-24 Thread Erin Hodgess
Dear R People: I finally (Yay!) got R installed in a classroom! Anyhow, I have a respectful request, please: could anyone recommend some nice undergrad projects in R, please? This is in a statistical computation class; first time being run. Thanks, Erin -- Erin Hodgess Associate Professor

[R] Add title in sink output and possibility of plot-like output for text

2008-09-24 Thread Michael Just
Hello, I have been using sink to create text file outputs. sink(summ_model1and2.txt) summary(model1) summary(model2) sink() Q1: Is there a way I could add a line of the text above the summary to act like a title? Also, I have been using the following to save plots from the lm function:

Re: [R] Generating series of distributions with the same skewness and different kurtosis or with same kurtosis and different skewness?

2008-09-24 Thread Greg Snow
And here is one way of using JohnsonFit to display how a distribution changes with changes in the moments: library(TeachingDemos) library(SuppDists) tmpfun - function(m1=0,m2=1,m3=-0.2,m4=3.7, xmin=-3, xmax=3){ parms - JohnsonFit( c(m1,m2,m3,m4), 'use') plot(function(x)

Re: [R] Proper power computation for one-sided binomial tests.

2008-09-24 Thread Collin Lynch
Thank you Peter. That is incredibly helpfyul, and much much smaller! Best, Collin. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Modifying output to Google Docs

2008-09-24 Thread Duncan Temple Lang
The command line app is fine, although there are some potential security issues or annoyances with a smooth work flow. But it would be more flexible and powerful to do it directly in R. And we don't need to know or program in C. The RCurl (https, getForm() and postForm()) and XML packages

Re: [R] Modifying output to Google Docs

2008-09-24 Thread Ajay ohri
Wow. Dont you get to name it on me or something (like the RPackageGoogleAjayDoc ?) Just kidding This is fabulous and the reason this community is so innovative Now if only someone can create a package for me to use R from Amazon EC2 and be free of the damn RAM problem forever

Re: [R] Proper power computation for one-sided binomial tests.

2008-09-24 Thread Johannes Hüsing
Am 23.09.2008 um 23:57 schrieb Peter Dalgaard: For this kind of problem I'd go directly for the binomial distribution. If the actual probability is 0, this is essentially deterministic and you can look at binom.test(0,99,p=.03, alt=less) This means that you don't sample from the p=.03

Re: [R] Why horizontal bars instead of a line

2008-09-24 Thread Greg Snow
Remember that computers are not as smart as you. Some smart people have written instructions for the computer on what to do in certain cases, but they can't anticipate everything, so when you tell the computer to do something that was different from what is anticipated, it either gives an

Re: [R] lattice xyplot symbols instead of colors and legend matching plot symbols or colors

2008-09-24 Thread stephen sefick
xyplot(GPP~(1/Iron)+(1/TSS)+TIN.TP, groups=RiverMile,data=f, scales=list(relation=free), par.settings = list( superpose.symbol = list( pch = c(1:7), fill = c(red, blue), col = black)), auto.key=TRUE, panel=function(...){

  1   2   >