Re: [R] customizing RGui/Rconsole

2009-08-24 Thread Dieter Menne
Philip A. Viton wrote: Is it possible to change the behavior of RGui/RConsole under MS Windows so that, if you submit a bunch of commands from a script window, when they've finished, focus is given to the Console, and not (as now) to the script? That's normally the job of the

[R] Generate random sample interval

2009-08-24 Thread Tammy Ma
Hi, R users, I have a problem about how to generate random sample interval from a duration. For example, during a time duration 0-70s, I want to generate a sample which last 10s. the sample could be 0-10 or 30-40s or 25-35s etc How could I do it in R`?? Thanks a lot. Tammy

[R] Odp: When factor is better than other types, such as vector and frame?

2009-08-24 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 23.08.2009 05:00:11: Hi, It is easy to understand the types vector and frame. But I am wondering why the type factor is designed in R. What is the advantage of factor compare with other data types in R? Can somebody give an example in which

Re: [R] Convert list to data frame while controlling column types

2009-08-24 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 23.08.2009 17:29:48: On 8/23/2009 9:58 AM, David Winsemius wrote: I still have problems with this statement. As I understand R, this should be impossible. I have looked at both you postings and neither of them clarify the issues. How can you have

Re: [R] Generate random sample interval

2009-08-24 Thread Tammy Ma
I need a simple algorithm to generate random ranges. From: metal_lical...@live.com To: r-help@r-project.org Date: Mon, 24 Aug 2009 09:55:35 +0300 Subject: [R] Generate random sample interval Hi, R users, I have a problem about how to generate random sample interval from a duration.

Re: [R] Generate random sample interval

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 2:55 AM, Tammy Ma wrote: Hi, R users, I have a problem about how to generate random sample interval from a duration. For example, during a time duration 0-70s, I want to generate a sample which last 10s. the sample could be 0-10 or 30-40s or 25-35s etc How could I

[R] how to add 95% confidential interval as vertical lines to x axe in density plot

2009-08-24 Thread Mao Jianfeng
Dear R-help listers, I want to add 95% confidential interval as vertical lines to x axe in density plot. I have found the library(hdrcde) can do this work, but I do not know how to handle functions of this library when I used ggplot2 to draw the graph. Thank you in advance. The data and codes

[R] Help in building new function

2009-08-24 Thread Ana Paula Mora
Hi: I've installed the precompiled binary for Windows. I need to use an existing function, but I want to introduce some slight changes to it. 1. Is there a way for me to find the source files through windows explorer? I know I can see it using edit(object name) but I want to know if I can see it

Re: [R] Help in building new function

2009-08-24 Thread Gabor Grothendieck
1. Just enter arima at the R console to see its source code (without comments). The source tar.gz for R is found by googling for R, clicking on CRAN in left column and choosing mirror. Or to view it online or get it via svn: https://svn.r-project.org/R/ 2. You want myarima's free variables

[R] R 2.9.2 is released

2009-08-24 Thread Peter Dalgaard
I've rolled up R-2.9.2.tar.gz a short while ago. This is a maintenance release and fixes a number of mostly minor issues. See the full list of changes below. You can get it from http://cran.r-project.org/src/base/R-2/R-2.9.2.tar.gz or wait for it to be mirrored at a CRAN site nearer to you.

[R] Using R CMD Batch similarly to gnuplot pause mouse

2009-08-24 Thread W Eryk Wolski
I have a gnuplot script: plot 'gp.dat' using 1:2 with line pause mouse and would like to have R CMD BATCH script that works similarily... x-read.table(gp.dat) options(device=windows) plot(x[,1],x[,2]) what I am looking for now is an R functions which behaves like pause mouse in gnuplot, i.e.

Re: [R] how to add 95% confidential interval as vertical lines to x axein density plot

2009-08-24 Thread ONKELINX, Thierry
Show us how you extract the confidence interval from the functions in the hdrcde library and then we might be able to help you. HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research

[R] Lattice xyplot: modify line width of plot lines

2009-08-24 Thread ukoenig
# Hi all, # I want to increase the line width of the plotted lines # in a xy-lattice plot. My own attempts were all in vain. # Without the group option the line width is modified - # with the option it is funnily enough not. # Please have a look at my syntax. # # Many thanks in advance # Udo

Re: [R] Lattice xyplot: modify line width of plot lines

2009-08-24 Thread Chuck Cleland
On 8/24/2009 4:47 AM, ukoe...@med.uni-marburg.de wrote: # Hi all, # I want to increase the line width of the plotted lines # in a xy-lattice plot. My own attempts were all in vain. # Without the group option the line width is modified - # with the option it is funnily enough not. # Please

Re: [R] Lattice xyplot: modify line width of plot lines

2009-08-24 Thread ukoenig
Now it works. Many thanks, Chuck! Quoting Chuck Cleland cclel...@optonline.net: On 8/24/2009 4:47 AM, ukoe...@med.uni-marburg.de wrote: # Hi all, # I want to increase the line width of the plotted lines # in a xy-lattice plot. My own attempts were all in vain. # Without the group option the

Re: [R] compare observed and fitted GAM values

2009-08-24 Thread Simon Wood
You haven't given quite enough information to be sure, but I would guess that this is not really a problem, but rather the interesting proporty of GLMs fitted with a canonical link described in e.g. section 2.1.8 of Wood (2006) Generalized additive models: and introduction with R, or at the

[R] natural sorting a data frame /vector by row

2009-08-24 Thread Moumita Das
How to NATURAL sort a vector or data frame* by row* , in ascending order ? V1 V2V3 V4 i1 5.00e-01 1.036197e-17 4.825338e+16 0. i104.001692e-18 1.365740e-17 2.930053e-01 0.76973827 i12 -1.052843e-17

Re: [R] image() generates many border lines in pdf, not on screen (quartz) - R 2.9.1 GUI 1.28 Tiger build 32-bit (5444) - OS X 10.5.8

2009-08-24 Thread Stefan Evert
On 23 Aug 2009, at 20:26, Uwe Ligges wrote: Since it looks like nobody answered so far: Your code is not reproducible, we do not have rfc, y, zVals nor NoCols. It's much easier to reproduce: just type in the first example from the image help page x - y - seq(-4*pi, 4*pi,

[R] calculating probability

2009-08-24 Thread maram salem
Hi all, I've a trivial question. If (q) is a continous variable,actually a vector of 1000 values. how to calculate the probability that q is greater than a specific value, i.e. P(q45)?? Thanks Maram [[alternative HTML version deleted]]

[R] Import/export ENVI files

2009-08-24 Thread Lucas Sevilla García
Hi! I'm a beginner with this webpage so, I don't know if I'm sending my question to the correct site. Anyway, I'm working with R and I need to import and export ENVI files, (*.HDR files). A colleague told me that there is a package to import/export envi files but I haven't found that package,

[R] R with MPI

2009-08-24 Thread polemon
Hello, I plan to use R with my cluster with OpenMPI. I need the packaged 'snow' and 'Rmpi' for that, however, I get an error while downloading and installing them: When I do a: install.packages(Rmpi, dependencies=T) I get this error: checking for mpi.h... no Try to find libmpi.so or

[R] hdf5 package segfault when processing large data

2009-08-24 Thread Budi Mulyono
Hi there, I am currently working on something that uses hdf5 library. I think hdf5 is a great data format, I've used it somewhat extensively in python via PyTables. I was looking for something similar to that in R. The closest I can get is this library: hdf5. While it does not work the same way

[R] transforming data glm

2009-08-24 Thread Mcdonald, Grant
Dear sir, I am fitting a glm with default identity link: model-glm(timetoacceptsecs~maleage*maletub*relweight*malemobtrue*femmobtrue) the model is overdisperesed and plot model shows a low level of linearity of the residuals. The overdispersion and linearity of residulas on the normal Q-Q

[R] Copy Paste from tktext on Mac

2009-08-24 Thread Anne Skoeries
Hi there, a text Window is supposed to map the shortcuts for copying and pasting (Ctrl-C, Ctrl-V) automatically. I'm working under Mac OS X and my text window doesn't really map these functions automatically - it works fine under Windows. Is there an easy way to map copypaste functions to

[R] error in creating gantt chart.

2009-08-24 Thread rajclinasia
hi every one, i have a excel sheet like this labels starts ends 1 first task 1-Jan-04 3-Mar-04 2 second task 2-Feb-04 5-May-04 3 third task 3-Mar-04 6-Jun-04 4 fourth task 4-Apr-04 8-Aug-04 5 fifth task 5-May-04 9-Sep-04 now i converted this excel sheet into csv file and i read

[R] Text in barplot.

2009-08-24 Thread koj
Dear all, i want to display a text in a barplot. Actual I simply work with the following lines: text(0.3,0.2,substring(teammembers.paint,1,140),cex=0.45,pos=4,srt=90) text(0.48,0.2,substring(teammembers.paint,141,1000),cex=0.45,pos=4,srt=90) The problem is, that teammembers.paint could be very

Re: [R] Text in barplot.

2009-08-24 Thread jim holtman
?strwrap On Mon, Aug 24, 2009 at 3:53 AM, kojjens.k...@gmx.li wrote: Dear all, i want to display a text in a barplot. Actual I simply work with the following lines: text(0.3,0.2,substring(teammembers.paint,1,140),cex=0.45,pos=4,srt=90)

Re: [R] discriminant analysis

2009-08-24 Thread Rubén Roa Ureta
Beatriz Yannicelli wrote: Dear all: Is it possible to conduct a discriminant analysis in R with categorical and continuous variables as predictors? Beatriz Beatriz, Simply doing this in the R console: RSiteSearch(discriminant) yields many promising links. In particular, check documentation

Re: [R] image() generates many border lines in pdf, not on screen (quartz) - R 2.9.1 GUI 1.28 Tiger build 32-bit (5444) - OS X 10.5.8

2009-08-24 Thread Uwe Ligges
Stefan Evert wrote: On 23 Aug 2009, at 20:26, Uwe Ligges wrote: Since it looks like nobody answered so far: Your code is not reproducible, we do not have rfc, y, zVals nor NoCols. It's much easier to reproduce: just type in the first example from the image help page x - y -

Re: [R] Problem with passing a string to subset

2009-08-24 Thread Sebastien Bihorel
Thanks Baptiste The eval(parse()) combination is just what I need. baptiste auguie wrote: Try this, mystr -c==1 subset(foo, eval(parse(text = mystr)) ) library(fortunes) fortune(parse) # try several times # I prefer this, but there is probably a better way mycond- quote(c==1) subset(foo,

Re: [R] calculating probability

2009-08-24 Thread Uwe Ligges
maram salem wrote: Hi all, I've a trivial question. If (q) is a continous variable,actually a vector of 1000 values. how to calculate the probability that q is greater than a specific value, i.e. P(q45)?? Do you want to estimate any distribution or do you just want the empirical

Re: [R] calculating probability

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 8:11 AM, maram salem wrote: Hi all, I've a trivial question. If (q) is a continous variable,actually a vector of 1000 values. how to calculate the probability that q is greater than a specific value, i.e. P(q45)?? sum(q45)/1000 # if no NA's in vector sum(q45,

Re: [R] natural sorting a data frame /vector by row

2009-08-24 Thread Ottorino-Luca Pantani
Moumita Das ha scritto: How to NATURAL sort a vector or data frame* by row* , in ascending order ? V1 V2V3 V4 i1 5.00e-01 1.036197e-17 4.825338e+16 0. i104.001692e-18 1.365740e-17 2.930053e-01 0.76973827

Re: [R] Copy Paste from tktext on Mac

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 8:21 AM, Anne Skoeries wrote: Hi there, a text Window is supposed to map the shortcuts for copying and pasting (Ctrl-C, Ctrl-V) automatically. I'm working under Mac OS X and my text window doesn't really map these functions automatically - it works fine under Windows.

[R] how to specify between group variance in lme

2009-08-24 Thread Bond, Stephen
Hello r-help, I am using lme with two specs for the variance func varComb(varFixed(~1/n)),varPower(~Age)) this produces worse forecasts than the lm model with simple weights=n I think due to the fact that the lme spec works on variance inside the group. I need to show it that 1/n scales the

Re: [R] image() generates many border lines in pdf, not on screen (quartz) - R 2.9.1 GUI 1.28 Tiger build 32-bit (5444) - OS X 10.5.8

2009-08-24 Thread Stefan Evert
Your code is not reproducible, we do not have rfc, y, zVals nor NoCols. It's much easier to reproduce: just type in the first example from the image help page x - y - seq(-4*pi, 4*pi, len=27) r - sqrt(outer(x^2, y^2, +)) image(z = z - cos(r^2)*exp(-r/6), col=gray((0:32)/32)) then

Re: [R] Import/export ENVI files

2009-08-24 Thread Barry Rowlingson
2009/8/24 Lucas Sevilla García luckocto...@hotmail.com: Hi! I'm a beginner with this webpage so, I don't know if I'm sending my question to the correct site. Anyway, I'm working with R and I need to import and export ENVI files, (*.HDR files). A colleague told me that there is a package to

Re: [R] how to specify between group variance in lme

2009-08-24 Thread Bond, Stephen
Clarification: Lm is much better than the base forecast from lme level=0, Level=1 produces a much tighter fit than lm. I was expecting that level=0 would produce something very close to lm, but it does not. [[alternative HTML version deleted]]

Re: [R] error in creating gantt chart.

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 8:22 AM, rajclinasia wrote: hi every one, i have a excel sheet like this labels starts ends 1 first task 1-Jan-04 3-Mar-04 2 second task 2-Feb-04 5-May-04 3 third task 3-Mar-04 6-Jun-04 4 fourth task 4-Apr-04 8-Aug-04 5 fifth task 5-May-04 9-Sep-04 now i

[R] R survival package error message - bug?!

2009-08-24 Thread Damjan Krstajic
Dear all, I have encountered a weird behaviour in R survival package which seems to me to be a bug. The weird behaviour happens when I am using 100 variables in the ridge function when calling coxph with following formula Surv(time = futime, event = fustat, type = right) ~ ridge(X1, X2,

Re: [R] R with MPI

2009-08-24 Thread polemon
On Mon, Aug 24, 2009 at 2:20 PM, polemon pole...@gmail.com wrote: Hello, I plan to use R with my cluster with OpenMPI. I need the packaged 'snow' and 'Rmpi' for that, however, I get an error while downloading and installing them: When I do a: install.packages(Rmpi, dependencies=T) I

Re: [R] Convert list to data frame while controlling column types

2009-08-24 Thread Alexander Shenkin
On 8/24/2009 2:06 AM, Petr PIKAL wrote: Hi r-help-boun...@r-project.org napsal dne 23.08.2009 17:29:48: On 8/23/2009 9:58 AM, David Winsemius wrote: I still have problems with this statement. As I understand R, this should be impossible. I have looked at both you postings and neither of

Re: [R] Selecting groups with R

2009-08-24 Thread Michael A. Miller
To drop empty factor levels from a subset, I use the following: a.subset - subset(dataset, Color!='BLUE') ifac - sapply(a.subset,is.factor) a.subset[ifac] - lapply(a.subset[ifac],factor) Mike dataset Color Score 1 RED10 2 RED13 3 RED12 4 WHITE22 5 WHITE27 6 WHITE

Re: [R] error in creating gantt chart.

2009-08-24 Thread Stefan Grosse
On Mon, 24 Aug 2009 05:22:17 -0700 (PDT) rajclinasia r...@clinasia.com wrote: R my.gantt.info-read.csv(C:/Documents and R Settings/balakrishna/Desktop/one.csv). R R and for create gantt chart i used below code. R R gantt.chart(my.gantt.info). This again is why others have pointed you first to

[R] CRAN (and crantastic) updates this week

2009-08-24 Thread Hadley Wickham
CRAN (and crantastic) updates this week New packages Updated packages New reviews --- This email provided as a service for the R community by http://crantastic.org. Like it? Hate it? Please let us know: crana...@gmail.com.

Re: [R] CRAN (and crantastic) updates this week

2009-08-24 Thread hadley wickham
Sorry, we had some problems with the initial sending of our weekly digest which resulted in a rather empty email. Here is the correct version: CRAN (and crantastic) updates this week New packages * atm (0.1.0) Charlotte Maia http://crantastic.org/packages/atm An R package for

Re: [R] CRAN (and crantastic) updates this week

2009-08-24 Thread Ronggui Huang
I like it. Thanks. 2009/8/24 hadley wickham crana...@gmail.com: Sorry, we had some problems with the initial sending of our weekly digest which resulted in a rather empty email. Here is the correct version: CRAN (and crantastic) updates this week New packages * atm (0.1.0)

[R] Re current neural network (RNN) package?

2009-08-24 Thread J_Laberga
Hello, I'm trying to tackle a problem that would require the implementation of a recurrent NN. However, even though the CRAN is very big, I can’t seem to find a package for this. Does anybody here know if one exits? BR, John -- View this message in context:

Re: [R] CRAN (and crantastic) updates this week

2009-08-24 Thread John Kerpel
Great idea - thx! On Mon, Aug 24, 2009 at 9:30 AM, hadley wickham crana...@gmail.com wrote: Sorry, we had some problems with the initial sending of our weekly digest which resulted in a rather empty email. Here is the correct version: CRAN (and crantastic) updates this week New packages

[R] Two lines, two scales, one graph

2009-08-24 Thread Rick
First of all, thanks to everyone who answers these questions - it's most helpful. I'm new to R and despite searching have not found an example of what I want to do (there are some good beginner's guides and a lot of complex plots, but I haven't found this). I would like to plot two variables

[R] robust method to obtain a correlation coeff?

2009-08-24 Thread Christian Meesters
Hi, Being a R-newbie I am wondering how to calculate a correlation coefficient (preferably with an associated p-value) for data like: d[,1] [1] 25.5 25.3 25.1 NA 23.3 21.5 23.8 23.2 24.2 22.7 27.6 24.2 ... d[,2] [1] 0.0 11.1 0.0 NA 0.0 10.1 10.6 9.5 0.0 57.9 0.0 0.0 ...

Re: [R] robust method to obtain a correlation coeff?

2009-08-24 Thread Ted Harding
On 24-Aug-09 14:47:02, Christian Meesters wrote: Hi, Being a R-newbie I am wondering how to calculate a correlation coefficient (preferably with an associated p-value) for data like: d[,1] [1] 25.5 25.3 25.1 NA 23.3 21.5 23.8 23.2 24.2 22.7 27.6 24.2 ... d[,2] [1] 0.0 11.1 0.0 NA

Re: [R] robust method to obtain a correlation coeff?

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 11:26 AM, (Ted Harding) wrote: On 24-Aug-09 14:47:02, Christian Meesters wrote: Hi, Being a R-newbie I am wondering how to calculate a correlation coefficient (preferably with an associated p-value) for data like: d[,1] [1] 25.5 25.3 25.1 NA 23.3 21.5 23.8 23.2 24.2

Re: [R] robust method to obtain a correlation coeff?

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 11:38 AM, David Winsemius wrote: On Aug 24, 2009, at 11:26 AM, (Ted Harding) wrote: On 24-Aug-09 14:47:02, Christian Meesters wrote: Hi, Being a R-newbie I am wondering how to calculate a correlation coefficient (preferably with an associated p-value) for data like:

[R] Multiply List by a Numeric

2009-08-24 Thread Brigid Mooney
I apologize for what seems like it should be a straighforward query. I am trying to multiply a list by a numeric and thought there would be a straightforward way to do this, but the best solution I found so far has a for loop. Everything else I try seems to throw an error non-numeric argument to

Re: [R] transforming data glm

2009-08-24 Thread Ben Bolker
Mcdonald, Grant wrote: Dear sir, I am fitting a glm with default identity link: model-glm(timetoacceptsecs~maleage*maletub*relweight*malemobtrue*femmobtrue) the model is overdisperesed and plot model shows a low level of linearity of the residuals. I don't see how the

[R] table function

2009-08-24 Thread Inchallah Yarab
hi, i want to use the function table to build a table not of frequence (number of time the vareable is repeated in a list or a data frame!!) but in function of classes I don t find a clear explnation in  examples of  ?table !!! example x      y    z 1    0   100 5    1   1500 6   

Re: [R] robust method to obtain a correlation coeff?

2009-08-24 Thread Bert Gunter
Inline below. Bert Gunter Genentech Nonclinical Biostatisics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of David Winsemius Sent: Monday, August 24, 2009 8:53 AM To: David Winsemius Cc: r-help@r-project.org Help;

Re: [R] Multiply List by a Numeric

2009-08-24 Thread Marc Schwartz
On Aug 24, 2009, at 10:58 AM, Brigid Mooney wrote: I apologize for what seems like it should be a straighforward query. I am trying to multiply a list by a numeric and thought there would be a straightforward way to do this, but the best solution I found so far has a for loop. Everything

Re: [R] Multiply List by a Numeric

2009-08-24 Thread Peter Ehlers
Try lapply(abc, function(x) x*3) Peter Ehlers Brigid Mooney wrote: I apologize for what seems like it should be a straighforward query. I am trying to multiply a list by a numeric and thought there would be a straightforward way to do this, but the best solution I found so far has a for

[R] Combining matrices

2009-08-24 Thread Daniel Nordlund
If I have two matrices like x - matrix(rep(c(1,2,3),3),3) y - matrix(rep(c(4,5,6),3),3) How can I combine them to get ? 1 1 1 4 4 4 1 1 1 5 5 5 1 1 1 6 6 6 2 2 2 4 4 4 2 2 2 5 5 5 2 2 2 6 6 6 3 3 3 4 4 4 3 3 3 5 5 5 3 3 3 6 6 6 The number of rows and the actual numbers above are unimportant,

[R] problem with BRugs

2009-08-24 Thread R Heberto Ghezzo, Dr
Hello, I am sorry, I have this problem before and Uwe send me the answer but I misplaced it and can not find it. writing a model for BRugs library(BRugs) Loading required package: coda Loading required package: lattice Welcome to BRugs running on OpenBUGS version 3.0.3 setwd(c:/tmp) Error in

Re: [R] table function

2009-08-24 Thread Erik Iverson
You need to create a factor that indicates which group the values in 'z' belong to. The easiest way to do that based on your situation is to use the 'cut' function to construct the factor, and then call 'table' using the result created by 'cut'. See ?cut and ?factor -Original

Re: [R] table function

2009-08-24 Thread Marc Schwartz
On Aug 24, 2009, at 10:59 AM, Inchallah Yarab wrote: hi, i want to use the function table to build a table not of frequence (number of time the vareable is repeated in a list or a data frame!!) but in function of classes I don t find a clear explnation in examples of ?table !!! example

Re: [R] table function

2009-08-24 Thread Glen Sargeant
Inchallah Yarab wrote: i want to do a table summerizing the number of variable where z is in [0-1000],],[1000-3000], [ 3000] You can use cut to create a new vector of labels and tabulate the result. Options control closed/open endpoints (see ?cut): z -

Re: [R] [Rd] Formulas in gam function of mgcv package

2009-08-24 Thread Gavin Simpson
[Note R-Devel is the wrong list for such questions. R-Help is where this should have been directed - redirected there now] On Mon, 2009-08-24 at 17:02 +0100, Corrado wrote: Dear R-experts, I have a question on the formulas used in the gam function of the mgcv package. I am trying to

Re: [R] Combining matrices

2009-08-24 Thread Henrique Dallazuanna
Try this; do.call(rbind, lapply(split(x, seq(nrow(x))), cbind, y)) On Mon, Aug 24, 2009 at 1:16 PM, Daniel Nordlund djnordl...@verizon.netwrote: If I have two matrices like x - matrix(rep(c(1,2,3),3),3) y - matrix(rep(c(4,5,6),3),3) How can I combine them to get ? 1 1 1 4 4 4 1 1 1 5

Re: [R] table function

2009-08-24 Thread Daniel Malter
Your question is a little vague. Do you just want to know how often z falls in one the three classes? If so, you could either code an indicator variable (e.g. z.cat) that expresses the three categories and then do table(z.cat). Alternatively, you could just do sum(z0z1000) sum(z1z3000)

Re: [R] image plot

2009-08-24 Thread Paul Hiemstra
ogbos okike schreef: Hi, I am trying to use the image function to do a color plot. My matrix columns are labeled y and x. I tried image(y, x) but I had error message (Error in image.default(y, x) : increasing 'x' and 'y' values expected). Could anybody please tell me how to add these increasing

Re: [R] Is there a fast way to do several hundred thousand ANOVA tests?

2009-08-24 Thread Benilton Carvalho
have you tried: fits - lm(a~b) fstat - sapply(summary(fits), function(x) x[[fstatistic]][[value]]) it takes 3secs for 100K columns on my machine (running on batt) b On Aug 23, 2009, at 9:55 PM, big permie wrote: Dear R users, I have a matrix a and a classification vector b such that

Re: [R] Combining matrices

2009-08-24 Thread Marc Schwartz
On Aug 24, 2009, at 11:16 AM, Daniel Nordlund wrote: If I have two matrices like x - matrix(rep(c(1,2,3),3),3) y - matrix(rep(c(4,5,6),3),3) How can I combine them to get ? 1 1 1 4 4 4 1 1 1 5 5 5 1 1 1 6 6 6 2 2 2 4 4 4 2 2 2 5 5 5 2 2 2 6 6 6 3 3 3 4 4 4 3 3 3 5 5 5 3 3 3 6 6 6 The

Re: [R] Selecting groups with R

2009-08-24 Thread Glen Sargeant
jlwoodard wrote: Each of the above lines successfully excludes the BLUE subjects, but the BLUE category is still present in my data set; that is, if I try table(Color) I get RED WHITE BLUE 82 151 0 How can I eliminate the BLUE category completely so I can do a t-test

[R] help with recalling data points in a specific region of the plot

2009-08-24 Thread Edward Chen
Hi all, Is there a quick way to display or recall data points from a specific region on the plot? For example I want the points from x5 and y5? Thank you very much! -- Edward Chen [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] Combining matrices

2009-08-24 Thread Marc Schwartz
On Aug 24, 2009, at 11:46 AM, Marc Schwartz wrote: On Aug 24, 2009, at 11:16 AM, Daniel Nordlund wrote: If I have two matrices like x - matrix(rep(c(1,2,3),3),3) y - matrix(rep(c(4,5,6),3),3) How can I combine them to get ? 1 1 1 4 4 4 1 1 1 5 5 5 1 1 1 6 6 6 2 2 2 4 4 4 2 2 2 5 5 5 2 2

Re: [R] robust method to obtain a correlation coeff?

2009-08-24 Thread John Kane
I may be misunderstanding the question but would cor(d1, use='complete.obs') or some other variant of use help? --- On Mon, 8/24/09, Christian Meesters meest...@imbie.uni-bonn.de wrote: From: Christian Meesters meest...@imbie.uni-bonn.de Subject: [R] robust method to obtain a correlation

Re: [R] Combining matrices

2009-08-24 Thread Gabor Grothendieck
Try this: kronecker(cbind(x, y), rep(1, 3)) On Mon, Aug 24, 2009 at 12:16 PM, Daniel Nordlunddjnordl...@verizon.net wrote: If I have two matrices like x - matrix(rep(c(1,2,3),3),3) y - matrix(rep(c(4,5,6),3),3) How can I combine  them to get ? 1 1 1 4 4 4 1 1 1 5 5 5 1 1 1 6 6 6 2 2

Re: [R] Combining matrices

2009-08-24 Thread Daniel Nordlund
-Original Message- From: Marc Schwartz [mailto:marc_schwa...@me.com] Sent: Monday, August 24, 2009 9:57 AM To: Daniel Nordlund Cc: r help Subject: Re: [R] Combining matrices On Aug 24, 2009, at 11:46 AM, Marc Schwartz wrote: On Aug 24, 2009, at 11:16 AM, Daniel Nordlund

Re: [R] help with recalling data points in a specific region of the plot

2009-08-24 Thread John Kane
I am assuming that you want to desplay all the data and highlight the subset. Set up a vector to indicate the breakdown of the data and you can do it fairly easily in ggplot2 if you treat the vector as a factor. library(ggplot) mydata - data.frame(x=1:21, y= -10:10) z - ifelse(mydata[,1]5

Re: [R] RMySQL - overwrite record, not table

2009-08-24 Thread Adrian Dusa
whizvast wrote: Hi, Adrian- If you use overwrite=T parameter, you will overwrite the entire table, not each record. this is the essence of my problem and i still haven't found out right solution. i am thinking of writing my own MySQLwriteTable function... Thank you for your answer

Re: [R] hdf5 package segfault when processing large data

2009-08-24 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Budi Mulyono Sent: Monday, August 24, 2009 3:38 AM To: r-help@r-project.org Subject: [R] hdf5 package segfault when processing large data Hi there, I am currently working

Re: [R] trouble with Vista reading files

2009-08-24 Thread spencerg
Dear Mike: I don't know. 1. What specific error message do you get? 2. Your example is too long for me to parse, especially with the color being stripped before I saw it. 3. Have you tried using debug(cro.etest.grab), then walking through

Re: [R] help with recalling data points in a specific region of the plot

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 12:55 PM, Edward Chen wrote: Hi all, Is there a quick way to display or recall data points from a specific region on the plot? For example I want the points from x5 and y5? Thank you very much! Your question is pretty light on specifics but assuming that you have

[R] lme, lmer, gls, and spatial autocorrelation

2009-08-24 Thread Timothy_Handley
Hello folks, I have some data where spatial autocorrelation seems to be a serious problem, and I'm unclear on how to deal with it in R. I've tried to do my homework - read through 'The R Book,' use the online help in R, search the internet, etc. - and I still have some unanswered questions. I'd

[R] Unique command not deleting all duplicate rows

2009-08-24 Thread Mehdi Khan
Hello everyone, when I run the unique command on my data frame, it deletes the majority of duplicate rows, but not all of them. Here is a sample of my data. How do I get it to delete all the rows? 6 -115.38 32.894 195 162.94 D 8419 D 7 -115.432 32.864 115 208.91 D 8419 D 8 -115.447 32.773

Re: [R] Unique command not deleting all duplicate rows

2009-08-24 Thread Erik Iverson
I really don't think this is the issue. I think the issue is that some columns of the data.frame, specifically V1, V2, and V4 should be checked versus R FAQ 7.31. -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Don McKenzie

Re: [R] lme, lmer, gls, and spatial autocorrelation

2009-08-24 Thread Bert Gunter
Have you looked at the Spatial task view on CRAN? That would seem to me the logical first place to go. Bert Gunter Genentech Nonclinical Biostatisics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of timothy_hand...@nps.gov Sent:

Re: [R] Unique command not deleting all duplicate rows

2009-08-24 Thread Mehdi Khan
Duplicated did not work, I agree with Erik. Is there any way I can specify a tolerance limit and then delete? On Mon, Aug 24, 2009 at 11:41 AM, Erik Iverson eiver...@nmdp.org wrote: I really don't think this is the issue. I think the issue is that some columns of the data.frame, specifically

Re: [R] Between-group variance from ANOVA

2009-08-24 Thread Mark Difford
Hi Emma, R gives you the tools to work this out. ## Example set.seed(7) TDat - data.frame(response = c(rnorm(100, 5, 2), rnorm(100, 20, 2))) TDat$group - gl(2, 100, labels=c(A,B)) with(TDat, boxplot(split(response, group))) summary(aov(response ~ group, data=TDat)) Regards, Mark. emj83

[R] create list entry from variable

2009-08-24 Thread rami jiossy
Hi; assume i-10 how can i create a list having key=10 and value=11 list(i=11) generates a list with 'i' [1] 11 and not 10 [1] 11 any help? Thanks _ Facebook. :ON:WL:en-US:SI_SB_facebook:082009 [[alternative

Re: [R] Unique command not deleting all duplicate rows

2009-08-24 Thread Bert Gunter
?round Bert Gunter Genentech Nonclinical Biostatisics -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Mehdi Khan Sent: Monday, August 24, 2009 11:52 AM To: Erik Iverson Cc: r-help@r-project.org Subject: Re: [R] Unique command not

Re: [R] Unique command not deleting all duplicate rows

2009-08-24 Thread Don McKenzie
duplicated() test.df V1 V2 V3 V4 V5 V6 V7 1 -115.380 32.894 195 162.940 D 8419 D 2 -115.432 32.864 115 208.910 D 8419 D 3 -115.447 32.773 1170 264.570 D 8419 D 4 -115.447 32.773 1170 264.570 D 8419 D 5 -115.447 32.773 1170 264.570 D 8419 D 6 -115.447 32.773 1170

Re: [R] lme, lmer, gls, and spatial autocorrelation

2009-08-24 Thread Timothy_Handley
Bert - I took a look at that page just now, and I'd classify my problem as spatial regression. Unfortunately, I don't think the spdep library fits my needs. Or at least, I can't figure out how to use it for this problem. The examples I have seen all use spdep with networks. They build a graph,

Re: [R] create list entry from variable

2009-08-24 Thread Henrique Dallazuanna
Try this: l - list(i + 1) names(l) - i On Mon, Aug 24, 2009 at 4:01 PM, rami jiossy sra...@hotmail.com wrote: Hi; assume i-10 how can i create a list having key=10 and value=11 list(i=11) generates a list with 'i' [1] 11 and not 10 [1] 11 any help? Thanks

Re: [R] create list entry from variable

2009-08-24 Thread rami jiossy
Yep; great thanks :) Date: Mon, 24 Aug 2009 16:11:20 -0300 Subject: Re: [R] create list entry from variable From: www...@gmail.com To: sra...@hotmail.com CC: r-help@r-project.org Try this: l - list(i + 1) names(l) - i On Mon, Aug 24, 2009 at 4:01 PM, rami jiossy sra...@hotmail.com wrote:

Re: [R] create list entry from variable

2009-08-24 Thread William Dunlap
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of rami jiossy Sent: Monday, August 24, 2009 12:01 PM To: R-Help Subject: [R] create list entry from variable Hi; assume i-10 how can i create a list having key=10 and

Re: [R] Two lines, two scales, one graph

2009-08-24 Thread Patrick Connolly
On Mon, 24-Aug-2009 at 08:00AM -0700, Rick wrote: First of all, thanks to everyone who answers these questions - it's most helpful. I'm new to R and despite searching have not found an example of what I want to do (there are some good beginner's guides and a lot of complex plots, but I

Re: [R] Help on comparing two matrices

2009-08-24 Thread Michael Kogan
David: Well, e.g. the first row has 2 ones in your output while there were no rows with 2 ones in the original matrix. Since the row and column sums can't be changed by sorting them, the output matrix can't be equivalent to the original one. But that means nothing, maybe it's intended and just

[R] plotting a grid with grid() ?

2009-08-24 Thread John Kane
I am trying to come up with a way of shading-in a grid for a simple pattern So far I can draw a square where I want but I cannot seem to draw a complete grid. I am just drawing them along the diagonal!! Clearly I am missing something simple but what? Any suggestions gratefully accepted.

Re: [R] Two lines, two scales, one graph

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 11:00 AM, Rick wrote: First of all, thanks to everyone who answers these questions - it's most helpful. I'm new to R and despite searching have not found an example of what I want to do (there are some good beginner's guides and a lot of complex plots, but I haven't

Re: [R] Help on comparing two matrices

2009-08-24 Thread David Winsemius
On Aug 24, 2009, at 4:01 PM, Michael Kogan wrote: David: Well, e.g. the first row has 2 ones in your output while there were no rows with 2 ones in the original matrix. Since the row and column sums can't be changed by sorting them, the output matrix can't be equivalent to the original

  1   2   >