Re: [R] Trouble with LSD.test / agricolae package

2015-01-07 Thread Uwe Ligges
On 07.01.2015 01:45, Jeff Newmiller wrote: The R for Windows FAQ item 2.2 States that installing to a network share is not supported. I don't know how you ended up trying to use \\ccofs1\Shared\Users\Dustin Batt for your home directory, but that appears to be a problem for R. Somehow you

[R] Select best maximum likelihood tree with raxml (package ips)

2015-01-07 Thread Luigi Marongiu
Dear all, I am trying to build a maximum likelihood (ML) tree using the raxml function of the package ips. I implemented the following script: tr-raxml(seq.align,exec = /usr/bin/raxmlHPC,b = 100, N = 10) where tr is the resulting tree, seq.align is the DNA alignment created with

Re: [R] ggplot with sparse layout

2015-01-07 Thread Tom Wright
Thanks, this is pretty good. Unfortunately I made an error in generating the sample dataframe, this code better represents the situation: data-data.frame(id=factor(rep(1:6,each=3)), group=factor(rep(c('a','a','b','c','c','c'),each=3)), x=runif(18,0,10),

Re: [R-es] Resumen de R-help-es, Vol 71, Envío 7

2015-01-07 Thread Carlos Andres Perez Angarita
Muchas gracias javier Si javier cada columna representa un individuo y todos los individuos tienen curvas con comportamientos muy parecidas entre si y lo q me referia de los eventos es poder dividir cada curva en 4 secciones donde pueda hallar la posicion o el indice de cada fila los maximos o

Re: [R] #library(CHAID) - Cross validation for chaid

2015-01-07 Thread Rodica Coderie via R-help
So, if I understood correctly, the mod$finalModel should be the resulted tree after cross-validation, right? It's the same as the original one. names(mod) mod$finalModel chaidUS Thanks! Rodica - Original Message - From: Rodica Coderie rodikgeorgi...@yahoo.com To: Max Kuhn

[R] plsrglm package

2015-01-07 Thread Evan Kransdorf
Hello, Is anyone using plsrglm to do PLS regressions? I am having trouble figuring out where my factor loadings are stored. For example, if I do a PLS regression and have an object from this, what is the list element of the object with the factor loadings? Thanks very much

Re: [R] Getting column names in the titles of histograms generated automatically

2015-01-07 Thread Richard M. Heiberger
tmp - data.frame(a=rnorm(10), b=rnorm(10)) for (i in names(tmp)) hist(tmp[[i]], main=paste(histogram of, i)) On Wed, Jan 7, 2015 at 10:13 PM, John Sorkin jsor...@grecc.umaryland.edu wrote: I am trying to automatically produce a series of histograms from every column of data frame with many

Re: [R] set.cor/mat.regress

2015-01-07 Thread William Revelle
Dear Andrea. First of all, the usual advice. If you have problems with a particular function in a particular package, first try the package author. I just happened to find your query. Then, to your question. On Jan 7, 2015, at 6:19 PM, Andrea Lamont alamont...@gmail.com wrote: I am

Re: [R] Getting column names in the titles of histograms generated automatically

2015-01-07 Thread John Sorkin
Richard, WOW! A totally new way to think about loop indices, many, many, thanks!John John David Sorkin M.D., Ph.D. Professor of Medicine Chief, Biostatistics and Informatics University of Maryland School of Medicine Division of Gerontology and Geriatric Medicine Baltimore VA Medical Center 10

[R] Getting column names in the titles of histograms generated automatically

2015-01-07 Thread John Sorkin
I am trying to automatically produce a series of histograms from every column of data frame with many columns. The columns names of a shortened from of the data frame follow: colnames-names(smdata) colnames [1] X13594_pre X15568_pre I want to have X13594_pre in the title of the first

[R] Matrix element-by-element multiplication

2015-01-07 Thread Steven Yen
I like to multiple the first and second column of a 10 x 3 matrix by 100. The following did not work. I need this in an operation with a much larger scale. Any help? aa-matrix(1:30,nrow=10,ncol=3); aa bb-matrix(c(100,100,1),nrow=1,ncol=3); bb dim(aa) dim(bb) aa*bb Results:

Re: [R] Matrix element-by-element multiplication

2015-01-07 Thread John McKown
On Wed, Jan 7, 2015 at 5:05 PM, Steven Yen sye...@gmail.com wrote: I like to multiple the first and second column of a 10 x 3 matrix by 100. The following did not work. I need this in an operation with a much larger scale. Any help? aa-matrix(1:30,nrow=10,ncol=3); aa

Re: [R] Matrix element-by-element multiplication

2015-01-07 Thread Peter Langfelder
You can create a suitable matrix bb as below (note the byrow = TRUE argument) aa-matrix(1:30,nrow=10,ncol=3); aa bb-matrix(c(100,100,1),nrow=10,ncol=3, byrow = TRUE); bb dim(aa) dim(bb) aa * bb You can also use matrix multiplication, but that;s slightly more involved:

Re: [R] Matrix element-by-element multiplication

2015-01-07 Thread Peter Langfelder
On Wed, Jan 7, 2015 at 3:15 PM, Peter Langfelder peter.langfel...@gmail.com wrote: You can create a suitable matrix bb as below (note the byrow = TRUE argument) aa-matrix(1:30,nrow=10,ncol=3); aa bb-matrix(c(100,100,1),nrow=10,ncol=3, byrow = TRUE); bb dim(aa) dim(bb) aa * bb You can

Re: [R] Matrix element-by-element multiplication

2015-01-07 Thread Steven Yen
Thank you both. Both John and Peter's suggestions work great!! At 06:17 PM 1/7/2015, John McKown wrote: On Wed, Jan 7, 2015 at 5:05 PM, Steven Yen mailto:sye...@gmail.comsye...@gmail.com wrote: I like to multiple the first and second column of a 10 x 3 matrix by 100. The following did not

[R] set.cor/mat.regress

2015-01-07 Thread Andrea Lamont
I am having a problem with the function mat.regress in the psych package that I was hoping someone would be able to give guidance on. I have a correlation matrix that is [76,76] with the first column representing the Y variable, the rest are the Xs. This is what my data looks like:

Re: [R] Duplicated date values aren't duplicates

2015-01-07 Thread bhradsky
I had the same error Error in as.ltraj(xy, id, date = da) : non unique dates for a given burst, but from a different cause, that I'd like to share in case others have the same issue. My data came from GPS collars that collected dates and time in Greenwich Mean Time (GMT). However, the package

Re: [R] #library(CHAID) - Cross validation for chaid

2015-01-07 Thread Rodica Coderie via R-help
Thanks Max! You are right! I used the train function below and no model was built. Do you know what can I use instead? library(caret) mod - train(x = USvoteS[,-1], y = USvoteS$vote3, method = modelInfo, trControl = trainControl(method = cv)) Thanks! Rodica

Re: [R] R example codes for direct standardization of rates

2015-01-07 Thread Therneau, Terry M., Ph.D.
The pyears() and survexp() routines in the survival package are designed for these calculations. See the technical report #63 of the Mayo Biostat group for examples

Re: [R] R example codes for direct standardization of rates

2015-01-07 Thread Muhuri, Pradip (SAMHSA/CBHSQ)
Hello Terry, Thank you so much for sending me this reference. Pradip Pradip K. Muhuri, PhD SAMHSA/CBHSQ 1 Choke Cherry Road, Room 2-1071 Rockville, MD 20857 Tel: 240-276-1070 Fax: 240-276-1260 From: Therneau, Terry M., Ph.D. [mailto:thern...@mayo.edu] Sent: Wednesday, January 07, 2015 4:39

[R] help with looping a specific code

2015-01-07 Thread Haznadar, Majda (NIH/NCI) [E]
Hi, I would like to loop a specific code through multiple files ending in .CDF in the same directory, but I am unsure of how to do that (new to R). Can someone please help me with this-thank you. This is the code I would like to loop: xr-xcmsRaw(mh141210xg2_sccs_00101.CDF, profstep=0)

Re: [R] Trouble with LSD.test / agricolae package

2015-01-07 Thread Dustin Batt
Thanks everyone for your help! This seems to have worked! Best regards, Dustin -Original Message- From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de] Sent: Wednesday, January 07, 2015 2:35 AM To: Jeff Newmiller; Dustin Batt; r-help@r-project.org Subject: Re: [R] Trouble with

[R] Converting ncdf to geoTIFF

2015-01-07 Thread John Wasige
Dear all, This is to request for help on how to convert ncdf data to geoTIFF in R with my script below; my script is: library(raster) library(ncdf) setwd('D:/Lund2015/CRU1982_2011') prec1 - open.ncdf('D:/Lund2015/CRU1982_2011/ b8_cruts3.20_pre_19820101_20111231_1month_mean.nc') print(prec1)

Re: [R] glmmADMB model fails when fixed effect interaction term is removed

2015-01-07 Thread Peter Claussen
Cara, Perhaps I a bit ignorant about mixed models, but I don’t understand why you have the dependent variable included as a random effect. The variable names you are using suggest an analysis of combined trials, so (1 | Site/Block/Strip.ID) seems more appropriate. Is this a factorial or a

Re: [R] ggplot with sparse layout

2015-01-07 Thread ONKELINX, Thierry
Dear Tom, You can’t do that with facet_grid() or facet_wrap(). You could try combine several plots with grid.arrange(). p1 - ggplot(subset(data, group == a),aes(x=x,y=y))+geom_point(aes(color=group))+facet_wrap(~id, nrow = 1) + scale_colour_discrete(drop = FALSE) p2 - ggplot(subset(data,

Re: [R] Converting ncdf to geoTIFF

2015-01-07 Thread Michael Sumner
raster() wants the filename, not a NetCDF connection. (raster is one of *very* few high-level interfaces to NetCDF, so it's somewhat unfamiliar to the norm) Cheers, Mike [[alternative HTML version deleted]] __ R-help@r-project.org mailing

Re: [R] help with looping a specific code

2015-01-07 Thread PIKAL Petr
See ?list.files You can use it to populate object with files in working directory. Then you can use tis object in loop. Cheers Petr -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Haznadar, Majda (NIH/NCI) [E] Sent: Wednesday, January 07, 2015

Re: [R] help with looping a specific code

2015-01-07 Thread Franklin Bretschneider
On 2015-01-07 , at 16:29, Haznadar, Majda (NIH/NCI) [E] majda.hazna...@nih.gov wrote: Hi, I would like to loop a specific code through multiple files ending in .CDF in the same directory, but I am unsure of how to do that (new to R). Can someone please help me with this-thank you.

Re: [R-es] Eventos en una serie de tiempo

2015-01-07 Thread Marcuzzi, Javier Rubén
Estimado Carlos Andres Perez Angarita Entiendo lo que busca, o creo entenderlo, yo uso algo de esa tecnología, podría ser adaptable a sus requerimientos. En primer lugar es más simple trabajar con los datos ordenados de otra forma, lo que usted dice columnas deberían ser filas, no es algo