Re: [R] Zeta-squared transformation use R?

2014-07-12 Thread Jeff Newmiller
Have you tried RSiteSearch(zeta squared) ? Someone may recognize this, but it never hurts to communicate where you have already looked. --- Jeff NewmillerThe . . Go Live...

Re: [R] Merge rows

2014-07-12 Thread Jim Lemon
On Fri, 11 Jul 2014 12:19:39 PM Ryan de Vera wrote: Hello all, I have a data frame filled with senders and recipients. Some of the senders have multiple rows with different recipients and I want to merge those rows. For example I have a...@email.com b...@email.com a...@email.com

Re: [R] Zeta-squared transformation use R?

2014-07-12 Thread Miles Yang
Hi Jeff, Yes, I searched that but it comes out zeta-squared coefficient instead of transformation method. I hope if there is someone have the experience in applying zeta-squared transformation in R? Appreciate any help on this. Regards, miles On Sat, Jul 12, 2014 at 4:01 PM, Jeff Newmiller

[R] lapply returns NULL ?

2014-07-12 Thread ce
Dear all, I have a list of arrays : foo-list(A = c(1,3), B =c(1, 2), C = c(3, 1)) foo $A [1] 1 3 $B [1] 1 2 $C [1] 3 1 if( foo$C[1] == 1 ) foo$C[1] lapply(foo, function(x) if(x[1] == 1 ) x ) $A [1] 1 3 $B [1] 1 2 $C NULL I don't want to list $C NULL in the output. How I can do

Re: [R] lapply returns NULL ?

2014-07-12 Thread Uwe Ligges
On 12.07.2014 15:25, ce wrote: Dear all, I have a list of arrays : foo-list(A = c(1,3), B =c(1, 2), C = c(3, 1)) foo $A [1] 1 3 $B [1] 1 2 $C [1] 3 1 if( foo$C[1] == 1 ) foo$C[1] lapply(foo, function(x) if(x[1] == 1 ) x ) $A [1] 1 3 $B [1] 1 2 $C NULL I don't want to

Re: [R] Merge rows

2014-07-12 Thread David Stevens
This is a (very) slightly modified version of Jim's reply that takes the sender's email our of the list element and uses it as the name so it can be accessed as newdat$'senders email' or newdat[['senders email']] newdat-list() for(sndr in unique(rdvdf$sender)) { newvec-

Re: [R] lapply returns NULL ?

2014-07-12 Thread Rui Barradas
Hello, Try the following. res - lapply(foo, function(x) if(x[1] == 1 ) x ) res[!sapply(res, is.null)] Hope this helps, Rui Barradas Em 12-07-2014 14:25, ce escreveu: Dear all, I have a list of arrays : foo-list(A = c(1,3), B =c(1, 2), C = c(3, 1)) foo $A [1] 1 3 $B [1] 1 2 $C [1]

Re: [R] lapply returns NULL ?

2014-07-12 Thread Jeff Newmiller
I think that removing them is something the OP doesn't understand how to do. The lapply function ALWAYS produces an output element for every input element. If this is not what you want then you need to choose a looping structure that is not so tightly linked to the input, such as a for loop

Re: [R] lapply returns NULL ?

2014-07-12 Thread ce
Thanks Jeff et. all, This is exactly what I needed. -Original Message- From: Jeff Newmiller [jdnew...@dcn.davis.ca.us] Date: 07/12/2014 10:38 AM To: Uwe Ligges lig...@statistik.tu-dortmund.de, ce zadi...@excite.com, r-help@r-project.org Subject: Re: [R] lapply returns NULL ? I think

Re: [R] lapply returns NULL ?

2014-07-12 Thread luke-tierney
Another option is Filter(function(x) x[1] == 1, foo) Best, luke On Sat, 12 Jul 2014, ce wrote: Thanks Jeff et. all, This is exactly what I needed. -Original Message- From: Jeff Newmiller [jdnew...@dcn.davis.ca.us] Date: 07/12/2014 10:38 AM To: Uwe Ligges

[R] lm and 450k data

2014-07-12 Thread Jessica Timms
Hi, I'm quite new to R and currently trying to use lm to fit linear models but I am currently stuck my code is as follows: Model1 = function(meth_matrix,exposure, X1, X2, X3, batch) { mod = lm(meth_matrix[, methcol]~exposure+X1+X2+X3+batch) res = summary(mod)$coef[2,]

[R] multiplicative error

2014-07-12 Thread Kevin Kunzmann
Hi, I am currently trying to build a regression model for calibration of HPLC outputs. I decided to use a multiplicative error model: Y_i = (a*X_i + b)*eps_i where the eps_i ~ iid N(0, s^2). Now I am having a hard time estimating my parameters ;) So the idea was to apply log() to both

[R] Metropolis-within Gibbs sampling

2014-07-12 Thread Baba Bukar
Dear R group, I am trying to update Negative Binomial distribution parameters using MCMC (Metropolis-within Gibbs sampling) method with Gamma prior. Can someone direct me to where i can get some R codes for this? Kind regards Zakir [[alternative HTML version deleted]]

[R] Fwd: lm and 450k data

2014-07-12 Thread Raghuraman Ramachandran
Try this: mod = lm(methcol ~ exposure+X1+X2+X3+batch, data = meth_matrix) res=coef(summary(mod))[2,] It will good if you give a subset of data you are working on by dput. Raghu On Sat, Jul 12, 2014 at 3:11 PM, Jessica Timms j.a.ti...@newcastle.ac.uk wrote: Hi, I'm quite new to R and

[R] R GUI for undergraduate lab class?

2014-07-12 Thread Louise Stevenson
Hi, I'm working on a new set of simple, ecological modeling exercises for our campus' undergraduate Introductory Biology lab series. The students work with simple population models by looking at graphs and seeing how changing parameter values and initial population sizes changes how the

Re: [R] R GUI for undergraduate lab class?

2014-07-12 Thread Roy Mendelssohn
There are probably several solutions to what you want to do, but you might look at Shiny as one possibility: http://shiny.rstudio.com -Roy M. On Jul 12, 2014, at 1:47 PM, Louise Stevenson louise.steven...@lifesci.ucsb.edu wrote: Hi, I'm working on a new set of simple, ecological modeling

Re: [R] R GUI for undergraduate lab class?

2014-07-12 Thread Richard M. Heiberger
I recommend RExcel. RExcel is an add-in for Windows Excel that gives complete access to the entirety of R from Windows Excel. It is free for educational use. The program by Erich Neuwirth is at http://rcom.univie.ac.at Our book, designed as a supplement to any text, is at

Re: [R] multiplicative error

2014-07-12 Thread David Winsemius
On Jul 12, 2014, at 4:25 AM, Kevin Kunzmann wrote: Hi, I am currently trying to build a regression model for calibration of HPLC outputs. I decided to use a multiplicative error model: Y_i = (a*X_i + b)*eps_i where the eps_i ~ iid N(0, s^2). Now I am having a hard time estimating my

[R] canonical correlation

2014-07-12 Thread Monaly Mistry
Hi, I was wondering if it's possible in R to do a canonical correlation with only one dependent variable and several independent variables. I've tried using cc(X,Y) but I got an error message. In this case I had 1 dependent variable and 10 independent variables. Error in cor(X, use = pairwise)

Re: [R] R GUI for undergraduate lab class?

2014-07-12 Thread Ranjan Maitra
Hi Louise, I tried using Deducer (graphical frontend to R) in my introductory class (Stat 105 at Iowa State University) for civil and construction engineers and this was a roaring success this year. The class itself has a wee bit of software experience, which has previously been done using JMP.