Re: [R] Robust ANCOVA

2016-03-27 Thread David Winsemius
> On Mar 27, 2016, at 9:52 AM, HAMID REZA ASHRAFI via R-help > wrote: > > HiI have a set of data with two independent variables, a pretest (covariate) > and posttest (dependent variable).Can anyone help me run robust ANCOVA?If you > wish I can send you the data file.yours > [[alternativ

[R] How to form groups for this specific problem?

2016-03-27 Thread Satish Vadlamani
Hello All: I would like to get some help with the following problem and understand how this can be done in R efficiently. The header is given in the data frame. *Component, TLA* C1, TLA1 C2, TLA1 C1, TLA2 C3, TLA2 C4, TLA3 C5, TLA3 Notice that C1 is a component of TLA1 and TLA2. I would like to

Re: [R] help with function calls

2016-03-27 Thread Bert Gunter
... and here is a maybe slightly neater approach using ?mapply (again with the method column changed to character(): f <- function(meth,i,fr) do.call(meth,list((fr[i,]))) mapply(FUN=f,meth=input.df[,4],seq_len(nrow(input.df)), MoreArgs = list(fr = input.df[,1:3]) ) Cheers, Bert Bert G

Re: [R] Robust ANCOVA

2016-03-27 Thread Jim Lemon
Hi Hamid, This looks a bit like a repeated measures analysis, but for a simple introduction to ANCOVA using R see the latter part of the following: http://www.stat.columbia.edu/~martin/W2024/R8.pdf Jim On Mon, Mar 28, 2016 at 3:52 AM, HAMID REZA ASHRAFI via R-help wrote: > HiI have a set of dat

Re: [R] 'split-lapply' vs. 'aggregate'

2016-03-27 Thread Fox, John
Dear Massimo, The difference is in the handling of NAs. Try, e.g., airquality <- na.omit(airquality) and compare again. Best, John - John Fox, Professor McMaster University Hamilton, Ontario Canada L8S 4M4 web: socserv.mcmaster.ca/jfox

[R] 'split-lapply' vs. 'aggregate'

2016-03-27 Thread Massimo Bressan
this might be a trivial question (eventually sorry for that!) but I definitely can not catch the problem here... please consider the following reproducible example: why of different results through 'split-lapply' vs. 'aggregate'? I've been also through a check against different methods (e.g. d

Re: [R] help with function calls

2016-03-27 Thread Bert Gunter
OOPS! I forgot to tell you that I first changed the "method" column, which is a factor, to character, with input.df$method <- as.character(input.df$method) Then things will work properly. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking t

Re: [R] help with function calls

2016-03-27 Thread Bert Gunter
1. return() is not needed in R functions (it's harmless, however). You might wish to go through an R function tutorial (many good ones are on the web) to learn about what slick things you can do with functions in R. 2. The following is just a brute force loop, so more elegant approaches are likely

Re: [R] loading

2016-03-27 Thread Uwe Ligges
Err, you already loaded it, there is not necessarily another message. Just try the functions you want to use. .. Best, Uwe Ligges On 27.03.2016 19:10, farzana akbari wrote: hi I install plm and pglm packages but I can not load no one of them. the massage of loading plm is library(plm)

[R] loading

2016-03-27 Thread farzana akbari
hi I install plm and pglm packages but I can not load no one of them. the massage of loading plm is > library(plm) Loading required package: Formula and for pglm is > library(pglm) Loading required package: maxLik Loading required package: miscTools Please cite the 'maxLik' package as: Hen

[R] Robust ANCOVA

2016-03-27 Thread HAMID REZA ASHRAFI via R-help
HiI have a set of data with two independent variables, a pretest (covariate) and posttest (dependent variable).Can anyone help me run robust ANCOVA?If you wish I can send you the data file.yours [[alternative HTML version deleted]] __ R-help@r-

Re: [R] Open source project that needs performance optimizations

2016-03-27 Thread Robert Sherry
I am not up on the internals of R but there does seem some run for parallelism. Are we talking about special hardware? or running this on an Intel Box? If it is the second, then I am thinking threads would be the way to go. Please consider the following R statements: for( i in 1:30 )a[i

Re: [R] Issue with var command in stats package

2016-03-27 Thread David Winsemius
> On Mar 27, 2016, at 9:14 AM, John, Larry wrote: > > Am using R version 3.2.4 in a fully updated version of Windows 7 and the most > current versions of coorplot, FactoMineR and factoextra to support multiple > correspondence analysis. However, today, a line of code that worked just fine > o

Re: [R] Issue with var command in stats package

2016-03-27 Thread Sarah Goslee
Somewhere you'd named an R object var, which is also the name of the function. There's no way to access part of a function with $ so var$contrib is throwing the closure error, which is telling you that you can't index a function. Sarah On Sun, Mar 27, 2016 at 12:14 PM, John, Larry wrote: > Am

[R] Issue with var command in stats package

2016-03-27 Thread John, Larry
Am using R version 3.2.4 in a fully updated version of Windows 7 and the most current versions of coorplot, FactoMineR and factoextra to support multiple correspondence analysis. However, today, a line of code that worked just fine on one set of data produced an error message on a different set

[R] Open source project that needs performance optimizations

2016-03-27 Thread PSATHAS NILOS-HRISTOS
Hello, i am an undergraduate student on computer engineering and im considering to do my thesis to an open source project and make performance optimizations and/or add parallelism to it where possible (or even better make use of GPU). Do you think that R-project is a good candidate? Than