Re: [R] Unsubscribe please

Hello, Don't reply only to me. 1) Filter the unwanted mails, 2) It takes few days to unsubscribe you. Regards, Pascal On 04/17/2013 02:59 PM, bert verleysen (beverconsult) wrote: I did this, but still I receive to much mails Bert Verleysen 00 32 (0)477 874 272 samen zoekend naar

Re: [R] Transformation of a variable in a dataframe

On Apr 16, 2013, at 10:33 PM, jpm miao wrote: HI, I have a dataframe with two variable A, B. I transform the two variable and name them as C, D and save it in a dataframe dfcd. However, I wonder why can't I call them by dfcd$C and dfcd$D? Because you didn't assign them to dfab$C.

[R] Merging big data.frame

Hi all, I am trying to merge 2 big data.frame. The problem is merge is memory intensive so R is going out of memory error: cannot allocate vector of size 360.1 Mb. To overcome this, I am exploring option of using data.table package. But its not helping in term of memory as merge in data.table is

Re: [R] Merging big data.frame

check out the sqldf package --- Jeff NewmillerThe . . Go Live... DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go... Live:

Re: [R] Merge

Thanks a lot:) Sent from my iPad On Apr 16, 2013, at 10:15 PM, arun smartpink...@yahoo.com wrote: Hi Farnoosh, YOu can use either ?merge() or ?join() DataA- read.table(text= ID v1 1 10 2 1 3 22 4 15 5 3 6 6 7 8 ,sep=,header=TRUE) DataB-

Re: [R] the joy of spreadsheets (off-topic)

Can you resend this link please? Thanks On Tue, Apr 16, 2013 at 10:33 PM, Jim Lemon j...@bitwrit.com.au wrote: On 04/17/2013 03:25 AM, Sarah Goslee wrote: ... Ouch. (Note: I know nothing about the site, the author of the article, or the study in question. I was pointed to it by someone

Re: [R] I don't understand the 'order' function

There is a blog post about this: http://www.portfolioprobe.com/2012/07/26/r-inferno-ism-order-is-not-rank/ And proof that it is possible to confuse them even when you know the difference. Pat On 16/04/2013 19:10, Julio Sergio wrote: Julio Sergio juliosergio at gmail.com writes: I thought

Re: [R] the joy of spreadsheets (off-topic)

On Apr 17, 2013, at 10:16 , Shane Carey wrote: Can you resend this link please? Psst: https://stat.ethz.ch/pipermail/r-help/2013-April/351669.html -- Peter Dalgaard, Professor Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone:

[R] Creating a vector with repeating dates

Dear R forum I have a data.frame df = data.frame(dates = c(4/15/2013, 4/14/2013, 4/13/2013, 4/12/2013), values = c(47, 38, 56, 92)) I need to to create a vector by repeating the dates as Current_date, 4/15/2013, 4/14/2013, 4/13/2013, 4/12/2013,  Current_date, 4/15/2013, 4/14/2013,

Re: [R] Creating a vector with repeating dates

?rep On Wed, Apr 17, 2013 at 11:11 AM, Katherine Gobin katherine_go...@yahoo.com wrote: Dear R forum I have a data.frame df = data.frame(dates = c(4/15/2013, 4/14/2013, 4/13/2013, 4/12/2013), values = c(47, 38, 56, 92)) I need to to create a vector by repeating the dates as

Re: [R] Understanding why a GAM can't suppress an intercept

hi Andrew. gam does suppress the intercept, it's just that this doesn't force the smooth through the intercept in the way that you would like. Basically for the parameteric component of the model '-1' behaves exactly like it does in 'lm' (it's using the same code). The smooths are 'added on'

Re: [R] Creating a vector with repeating dates

Dear Andrija Djurovic, Thanks for the suggestion. Ia m aware of rep. However, here I need to repeat not only dates, but a string Current_date. Thus, I need to create a vector ( to be included in some other data.frame) with the name say dt which will contain dt Current_date 4/15/2013

Re: [R] Creating a vector with repeating dates

Hello, Try the following. rep(c(Current_date, as.character(df$dates)), 3) Hope this helps, Rui Barradas Em 17-04-2013 10:11, Katherine Gobin escreveu: Dear R forum I have a data.frame df = data.frame(dates = c(4/15/2013, 4/14/2013, 4/13/2013, 4/12/2013), values = c(47, 38, 56, 92)) I

Re: [R] Creating a vector with repeating dates

Hi. Here are some examples that can maybe help you: a - Current date b - Sys.Date()-1:5 a b class(a) class(b) c(a,b) mode(b) as.numeric(b) class(c(a,b)) c(a, as.character(b)) class(c(a,b)) class(c(a,as.character(b))) Hope this helps. On Wed, Apr 17, 2013 at 11:21 AM, Katherine Gobin

Re: [R] Creating a vector with repeating dates

On 04/17/2013 07:11 PM, Katherine Gobin wrote: Dear R forum I have a data.frame df = data.frame(dates = c(4/15/2013, 4/14/2013, 4/13/2013, 4/12/2013), values = c(47, 38, 56, 92)) I need to to create a vector by repeating the dates as Current_date, 4/15/2013, 4/14/2013, 4/13/2013, 4/12/2013,

Re: [R] I don't understand the 'order' function

On Apr 17, 2013, at 10:41 , Patrick Burns wrote: There is a blog post about this: http://www.portfolioprobe.com/2012/07/26/r-inferno-ism-order-is-not-rank/ And proof that it is possible to confuse them even when you know the difference. It usually helps to remember that x[order(x)] is

Re: [R] Spatial Ananlysis: zero.policy=TRUE doesn't work for no neighbour regions??

Molo kurz_m at uni-hohenheim.de writes: ... *As there are some regions without neighbours in my data I use the following code to create the Weights Matrix:* W_Matrix- nb2listw(location_nbq, style=W, zero.policy=TRUE) W_Matrix *And get this Output:* ... /(Error in

[R] Problem with DateVisit-gives wrong year?

Hi I have the following factor of dates that I want to converted to Date class so I can extract the month test.date [1] 14/05/2012 14/05/2012 14/05/2012 14/05/2012 14/05/2012 14/05/2012 [7] 14/05/2012 14/05/2012 14/05/2012 14/05/2012 201 Levels: 01/10/2012 01/11/2012 01/12/2012 02/07/2012 ...

[R] Regularized Regressions

Hi all, I would greatly appreciate if someone was so kind and share with us a package or method that uses a regularized regression approach that balances a regression model performance and model complexity. That said I would be most grateful is there is an R-package that combines Ridge (sum of

Re: [R] Problem with DateVisit-gives wrong year?

On 04/17/2013 09:18 PM, Pancho Mulongeni wrote: Hi I have the following factor of dates that I want to converted to Date class so I can extract the month test.date [1] 14/05/2012 14/05/2012 14/05/2012 14/05/2012 14/05/2012 14/05/2012 [7] 14/05/2012 14/05/2012 14/05/2012 14/05/2012 201

Re: [R] Problem with DateVisit-gives wrong year?

HI, test.date- rep(c(14/05/2012,01/10/2012,28/09/2012),each=6)  as.Date(test.date,%d/%m/%Y) # [1] 2012-05-14 2012-05-14 2012-05-14 2012-05-14 2012-05-14  #[6] 2012-05-14 2012-10-01 2012-10-01 2012-10-01 2012-10-01 #[11] 2012-10-01 2012-10-01 2012-09-28 2012-09-28 2012-09-28 #[16] 2012-09-28

Re: [R] Transformation of a variable in a dataframe

Hi, You may also try: dfab-data.frame(A,B) library(plyr)  dfcd-subset(mutate(dfab,C=A*2,D=B*3),select=-c(A,B)) #or  dfcd1-subset(within(dfab,{D-B*3;C-A*2}),select=-c(A,B)) dfcd$C #[1] 2 4 6  dfcd$D #[1] 12 18 21 A.K. - Original Message - From: jpm miao miao...@gmail.com To: r-help

Re: [R] the joy of spreadsheets (off-topic)

On Tue, Apr 16, 2013 at 4:33 PM, Jim Lemon j...@bitwrit.com.au wrote: On 04/17/2013 03:25 AM, Sarah Goslee wrote: The final point does relate to Excel and any application that hides what is going on to the casual observer. I will treasure this URL to give to anyone who chastises my moaning

[R] remove higher order interaction terms

Dear all, Consider the model below: x - lm(mpg ~ cyl * disp * hp * drat, mtcars) summary(x) Call: lm(formula = mpg ~ cyl * disp * hp * drat, data = mtcars) Residuals: Min 1Q Median 3Q Max -3.5725 -0.6603 0.0108 1.1017 2.6956 Coefficients: Estimate

Re: [R] remove higher order interaction terms

On Apr 17, 2013, at 7:23 AM, Liviu Andronic landronim...@gmail.com wrote: Dear all, Consider the model below: x - lm(mpg ~ cyl * disp * hp * drat, mtcars) summary(x) Call: lm(formula = mpg ~ cyl * disp * hp * drat, data = mtcars) Residuals: Min 1Q Median 3Q Max

Re: [R] compate tow different data frame

Hi, dat1- read.table(text= V1    V2 A  1 B  2 C  1 D  3 ,sep=,header=TRUE,stringsAsFactors=FALSE) dat2- read.table(text= V3    V2 AAA   1 BBB   2 CCC   3 ,sep=,header=TRUE,stringsAsFactors=FALSE) library(plyr) join(dat1,dat2,by=V2,type=full) #  V1 V2  V3 #1  A  1

Re: [R] remove higher order interaction terms

On Wed, Apr 17, 2013 at 2:33 PM, Marc Schwartz marc_schwa...@me.com wrote: If you only want up to say second order interactions: summary(lm(mpg ~ (cyl + disp + hp + drat) ^ 2, data = mtcars)) This is what I was looking for. Thank you so much. This is covered in ?formula Indeed. I tried to

[R] On matrix calculation

Hello again, Let say I have a matrix: Mat - matrix(1:12, 4, 3) And a vector: Vec - 5:8 Now I want to do following: Each element of row-i in 'Mat' will be divided by i-th element of Vec Is there any direct way to doing that? Thanks for your help

Re: [R] Regularized Regressions

Merhaba, Hello to you too Mehmet (Yasu ki sena) Thank you for your email and especially for sharing this package. I appreciate it. However, my feeling is that this package does not have the third component of Best Subsets (pls correct me if I am wrong). It uses only a combination of Ridge and

Re: [R] On matrix calculation

Hi, Try: sweep(Mat,1,Vec,/)  #     [,1] [,2] [,3] #[1,] 0.200    1 1.80 #[2,] 0.333    1 1.67 #[3,] 0.4285714    1 1.571429 #[4,] 0.500    1 1.50 do.call(rbind,lapply(seq_len(nrow(Mat)),function(i) Mat[i,]/Vec[i])) #  [,1] [,2] [,3] #[1,] 0.200   

Re: [R] how to change the date into an interval of date?

Hi, Try: evt_c.1- read.table(text= patient_id   responsed_at 1    2010-5 1    2010-7 1    2010-8 1    2010-9 2    2010-5 2    2010-6 2    2010-7 ,sep=,header=TRUE,stringsAsFactors=FALSE)  lst1-split(evt_c.1,evt_c.1$patient_id)  

Re: [R] On matrix calculation

On 17-04-2013, at 14:39, Christofer Bogaso bogaso.christo...@gmail.com wrote: Hello again, Let say I have a matrix: Mat - matrix(1:12, 4, 3) And a vector: Vec - 5:8 Now I want to do following: Each element of row-i in 'Mat' will be divided by i-th element of Vec Is there

Re: [R] the joy of spreadsheets (off-topic)

On Tue, Apr 16, 2013 at 1:25 PM, Sarah Goslee sarah.gos...@gmail.com wrote: Given that we occasionally run into problems with comparing Excel results to R results, and other spreadsheet-induced errors, I thought this might be of interest.

Re: [R] Problem with DateVisit-gives wrong year?

Thank you, I see this was due to me using %y instead of %Y, See ?strptime -Original Message- From: arun [mailto:smartpink...@yahoo.com] Sent: 17 April 2013 12:31 To: Pancho Mulongeni Cc: R help Subject: Re: [R] Problem with DateVisit-gives wrong year? HI, test.date-

Re: [R] how to change the date into an interval of date?

Hi, I hope this is what you are looking for: library(plyr)   mutate(evt_c.1,t=ave(as.numeric(gsub(.*\\-,,responsed_at)),patient_id,gsub(-.*,,responsed_at),FUN=function(x) c(0,cumsum(diff(x) #   patient_id responsed_at t #1   1   2010-5 0 #2   1   2010-7 2 #3  

[R] Package VIF

Hi Could you perhaps possibly help me. I would like to use the package VIF but cannot get results I attach the .csv file and my R code. What do I have to do ? Any help is greatly appreciated. library(VIF) coal - read.csv(e:/freekvif/cqa1.csv,header=TRUE) y - as.numeric(coal$AI) x -

Re: [R] Regularized Regressions

Perhaps I am wrong, but I think there are only a few packages supporting Elastic Net, and none of them also perform Best Subsets. On Wed, Apr 17, 2013 at 8:46 AM, Christos Giannoulis cgiann...@gmail.comwrote: Merhaba, Hello to you too Mehmet (Yasu ki sena) Thank you for your email and

Re: [R] odfWeave: Some questions about potential formatting options

Paul, #1: I've never tried but you might be able to escape the required tags in your text (e.g. in html you could write out the b in your text). #3: Which output? Is this in text? #2: I may be possible and maybe easy to implement. So if you want to dig into it, have at it. For me, I'm

Re: [R] how to change the date into an interval of date?

Hi, Try this: library(mondate) mutate(evt_c.1,t=ave(round(as.numeric(mondate(paste(evt_c.1[,2],-01,sep=,patient_id,FUN=function(x) c(0,cumsum(diff(x)  #  patient_id responsed_at t #1   1   2010-5 0 #2   1   2010-7 2 #3   1   2010-8 3 #4   1  

[R] Bug in VGAM z value and coefficient ?

Dear, When i multiply the y of a regression by 10, I would expect that the coefficient would be multiply by 10 and the z value to stay constant. Here some reproducible code to support the case. *Ex 1* library(mvtnorm) library(VGAM) set.seed(1) x=rmvnorm(1000,sigma=matrix(c(1,0.75,0.75,1),2,2))

[R] Error: could not find function invlogit and bayesglm

I have installed the arm package and its dependents (e.g MATRIX, etc), but cannot use the functions invlogit and bayesglm because it gives me the error message Error: could not find function invlogit or Error: could not find function invlogit. What could be the problem. Regards Carrington

Re: [R] Creating a vector with repeating dates

Dear Sir, Thanks a lot for your valuable suggestions and help. Regards Katherine --- On Wed, 17/4/13, Jim Lemon j...@bitwrit.com.au wrote: From: Jim Lemon j...@bitwrit.com.au Subject: Re: [R] Creating a vector with repeating dates To: Katherine Gobin katherine_go...@yahoo.com Cc:

[R] Anova unbalanced

Hello everybody, I have got a data set with about 400 companies. Each company has a score for its enviroment comportment between 0 and 100. These companies belong to about 15 different countries. I have e.g. 70 companies from UK and 5 from Luxembourg,- so the data set is pretty unbalanced and

[R] Pathological case for the function 'brunch' of CAPER?

Dear R-enthusiats, While using the regression fonction 'brunch' of CAPER (with R v2.15.4), in a simple case (binary variable Yes/No vs. a continuous variable) I ended with an unexplained error: Error in if (any(stRes robust)) { : missing value where TRUE/FALSE needed I simplified my code

Re: [R] Error: could not find function invlogit and bayesglm

Have you loaded it? library(arm) John Kane Kingston ON Canada -Original Message- From: masan...@uniswa.sz Sent: Wed, 17 Apr 2013 10:08:39 +0200 To: r-help@r-project.org Subject: [R] Error: could not find function invlogit and bayesglm I have installed the arm package and its

Re: [R] Error: could not find function invlogit and bayesglm

Hi Carrington, You also need the boot package (see http://stat.ethz.ch/R-manual/R-patched/library/boot/html/inv.logit.html ) As for the other function, please load the arm package, e.g., require(arm) require(boot) and then you will be able to use the functions mentioned below. HTH, Jorge.-

Re: [R] Error: could not find function invlogit and bayesglm

On 17-04-2013, at 10:08, S'dumo Masango masan...@uniswa.sz wrote: I have installed the arm package and its dependents (e.g MATRIX, etc), but cannot use the functions invlogit and bayesglm because it gives me the error message Error: could not find function invlogit or Error: could not find

Re: [R] Overlay two stat_ecdf() plots

Hi Your files have only one column, so melted data is virtualy the same. When I read them as test1 and test2 i can do plot(ecdf(test1$Down)) plot(ecdf(test2$Up), add=T, col=2) Or using previously ustated ggplot2 package test-rbind(melt(test1),melt(test2)) p-ggplot(test, aes(x=value,

[R] Mancova with R

Dear all, I'm trying to compare two sets of variables, the first set is composed exclusively of numerical variables and the second regroups factors and numerical variables. I can't use a Manova because of this inclusion of numerical variables in the second set. The solution should be to

Re: [R] plot 2 y axis

Exel is hardly the epitome of good graphing. There are a couple of ways to do what you want and Jim has shown you one , but ... What about using two panels to present the data as in opar - par(mfrow = c(2, 1)) plot(dat1$Date, dat1$Weight, col = red, xlab = , ylab = Weight) plot(dat1$Date,

[R] spatial graph and its boundary

Dear r-helpers, I have a graph created using the library 'spatgraphs'. library(spatstat) library(spatgraphs) xy - rbind(c(28.39, -16.27), c(30.62, -20.13), c(32.25, -28.7), c(22.43, -27.22), c(27.5, -21.17), c(31.22, -24.52), c(17.93, -26.92), c(18.72, -17.95), c(24.15, -17.82), c(29.23, -22.85))

[R] use of names() within lapply()

Dear all, List g has 2 elements names(g) [1] 2009-10-07 2012-02-29 and the list plot lapply(g, plot, main=names(g)) results in equal plot titles with both list names, whereas distinct titles names(g[1]) and names(g[2]) are sought. Clearly, lapply is passing 'g' in stead of consecutively

Re: [R] Mancova with R

Dear Remi, Take a look at the Anova() function in the car package. In your case, you could use Anova(lm(as.matrix(Y) ~ x + z)) or, for more detail, summary(Anova(lm(as.matrix(Y) ~ x + z))) I hope this helps, John John Fox Sen. William

Re: [R] use of names() within lapply()

On 17/04/2013 11:04 AM, Ivan Alves wrote: Dear all, List g has 2 elements names(g) [1] 2009-10-07 2012-02-29 and the list plot lapply(g, plot, main=names(g)) results in equal plot titles with both list names, whereas distinct titles names(g[1]) and names(g[2]) are sought. Clearly, lapply

Re: [R] use of names() within lapply()

Hi, Try: set.seed(25) g- list(sample(1:40,20,replace=TRUE),sample(40:60,20,replace=TRUE)) names(g)- c(2009-10-07,2012-02-29) pdf(Trialnew.pdf)  lapply(seq_along(g),function(i) plot(g[[i]],main=names(g)[i])) dev.off() A.K. - Original Message - From: Ivan Alves pap: u...@me.com To:

[R] proper way to handle obsolete function names

Dear R community, what would be the proper R way to handle obsolete function names? I have created several packages with functions and sometimes would like to change the name of a function but would like to create a mechanism that other scripts of functions using the old name still work.

Re: [R] how to change the date into an interval of date?

Hi. No problem. cc:ing to Rhelp. A.K. From: GUANGUAN LUO guanguan...@gmail.com Sent: Wednesday, April 17, 2013 10:25 AM Subject: Re: how to change the date into an interval of date? Thank you so much . That is exactly the things i want. GG Hi, Try

Re: [R] Mancova with R

Dear John, Thanks for your comments! But when I tried your suggestion, the output was as the following:  Response Dist_arbre :             Df     Sum Sq    Mean Sq F value    Pr(F)     Poids        1 0.00010398 0.00010398  6.2910 0.0364733 *   Age          1

Re: [R] use of names() within lapply()

Dear Duncan and A.K. Many thanks for your super quick help. The modified lapply did the trick, mapply died with a error Error in dots[[2L]][[1L]] : object of type 'builtin' is not subsettable. Kind regards, Ivan On 17 Apr 2013, at 17:12, Duncan Murdoch murdoch.dun...@gmail.com wrote: On

Re: [R] Mancova with R

Dear Remi, On Wed, 17 Apr 2013 08:23:07 -0700 (PDT) Rémi Lesmerises remilesmeri...@yahoo.ca wrote: Dear John, Thanks for your comments! But when I tried your suggestion, the output was as the following:  Response Dist_arbre :             Df     Sum Sq    Mean Sq F value    Pr(F)    

Re: [R] Anova unbalanced

Dear Claudia, Your question has been posed on many previous occasions. The (short) answer has always been the same: have a look at the Anova function in the car package but before doing that, get a copy of John Fox's Applied Regression Analysis and Generalized Linear Models book. Best, José

Re: [R] Regularized Regressions

Hi Levi, Thank you very much for your reply and concern. There was a package till 2011...or even 2012...entitled Generalized Path Seeking Regression in R which was using a combination of all of these methods...as a way to regularize regression. However, I presently can't find it. Thus, I was

Re: [R] odfWeave: Some questions about potential formatting options

Le mardi 16 avril 2013 à 10:15 -0700, Paul Miller a écrit : Hi Milan and Max, Thanks to each of you for your reply to my post. Thus far, I've managed to find answers to some of the questions I asked initially. I am now able to control the justification of the leftmost column in my tables,

[R] Best way to calculate averages of Blocks in an matrix?

Folks, I recently was given a simulated data set like the following subset: sim_sub-structure(list(V11 = c(0.01, 0, 0, 0.01, 0, 0.01, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), V12 = c(0, 0, 0, 0.01, 0.03, 0, 0, 0, 0, 0, 0, 0.01, 0, 0.01, 0, 0, 0, 0, 0, 0.04), V13 = c(0, 0, 0, 0.01, 0, 0, 0,

Re: [R] use of names() within lapply()

Dear Ivan, No problem. If you want it in a single plot: matplot(do.call(cbind,g),ylab=value,pch=1:2,main=Some plot,col=c(red,orange),type=o)   legend(topleft,inset=.01,lty=c(1,1),title=Plot,col=c(red,orange),names(g),horiz=TRUE) A.K. From: Ivan Alves

Re: [R] Best way to calculate averages of Blocks in an matrix?

do.call(rbind,lapply(split(sim_sub,((seq_len(nrow(sim_sub))-1)%/% 5)+1),colMeans))   #  V11   V12   V13   V14   V15  V16   V17   V18   V19   V20 #1 0.004 0.008 0.002 0.008 0.004 0.01 0.004 0.014 0.022 0.010 #2 0.002 0.000 0.002 0.006 0.006 0.01 0.018 0.006 0.010 0.014 #3 0.000 0.004 0.002 0.012

Re: [R] Merge

Hi, I have a quick question here. Lets say he has three data frames and he needs to combine those three data frame using merge. Can we simply use merge to join three data frames ? I remember I had some problem using merge for more than two dataframes. Thanks. On Wed, Apr 17, 2013 at 1:05 AM,

Re: [R] Best way to calculate averages of Blocks in an matrix?

Also, do.call(rbind,lapply(split(sim_sub,rep(1:(1+nrow(sim_sub)/5),each=5)[seq_len(nrow(sim_sub))]),colMeans)) #    V11   V12   V13   V14   V15  V16   V17   V18   V19   V20 #1 0.004 0.008 0.002 0.008 0.004 0.01 0.004 0.014 0.022 0.010 #2 0.002 0.000 0.002 0.006 0.006 0.01 0.018 0.006 0.010 0.014

Re: [R] Best way to calculate averages of Blocks in an matrix?

Hello, Try the following. blocks - rep(1:(1 + nrow(sim_sub) %/% 5), each = 5)[seq_len(nrow(sim_sub))] aggregate(sim_sub, list(blocks), FUN = mean) Hope this helps, Rui Barradas Em 17-04-2013 18:04, arun escreveu: do.call(rbind,lapply(split(sim_sub,((seq_len(nrow(sim_sub))-1)%/%

Re: [R] Mancova with R

On Apr 17, 2013, at 16:47 , Rémi Lesmerises wrote: Dear all, I'm trying to compare two sets of variables, the first set is composed exclusively of numerical variables and the second regroups factors and numerical variables. I can't use a Manova because of this inclusion of numerical

Re: [R] Merge

HI Janesh, YOu can use: library(plyr) ?join_all() #From the help page:  dfs - list(    a = data.frame(x = 1:10, a = runif(10)),    b = data.frame(x = 1:10, b = runif(10)),    c = data.frame(x = 1:10, c = runif(10)) ) join_all(dfs) join_all(dfs, x)  join_all(dfs,

[R] mgcv: how select significant predictor vars when using gam(...select=TRUE) using automatic optimization

I have 11 possible predictor variables and use them to model quite a few target variables. In search for a consistent manner and possibly non-manual manner to identify the significant predictor vars out of the eleven I thought the option select=T might do. Example: (here only 4 pedictors) first

Re: [R] Merge

Hello Arun, Thank you so much for the prompt reply. I have one simple question here. DOes three dots (...) in the reduce function means we are applying for three dataframes here ? So, if we were to combine four would that dots be four dots ? Thanks. On Wed, Apr 17, 2013 at 12:16 PM, arun

Re: [R] Merge

No, you don't have to use four dots. Please check these links for further details: http://stackoverflow.com/questions/5890576/usage-of-three-dots-or-dot-dot-dot-in-functions http://cran.r-project.org/doc/manuals/R-lang.pdf A.K. From: Janesh Devkota

Re: [R] proper way to handle obsolete function names

On Apr 17, 2013, at 10:17 AM, Jannis bt_jan...@yahoo.de wrote: Dear R community, what would be the proper R way to handle obsolete function names? I have created several packages with functions and sometimes would like to change the name of a function but would like to create a

Re: [R] mgcv: how select significant predictor vars when using gam(...select=TRUE) using automatic optimization

Jan, What mgcv version are you using, please? (Older versions have a poor p-value approximation when select=TRUE, but of course it's possible that you've managed to break the newer approximation as well) The 'select=TRUE' option adds a penalty to each smooth, to allow it to be penalized out

Re: [R] proper way to handle obsolete function names

On Wed, Apr 17, 2013 at 10:36 AM, R. Michael Weylandt michael.weyla...@gmail.com It sounds like you want .Deprecate ?.Deprecate Perhaps you meant Deprecated? ?Deprecated Best, Peter __ R-help@r-project.org mailing list

[R] simulation\bootstrap of list factors

Dear R experts, I am trying to simulate a list containing data matrices. Unfortunately, I don't manage to get it to work. A small example: n=5 nbootstrap=2 subsets-list() for (i in 1:n){ subsets[[i]] - rnorm(5, mean=80, sd=1) for (j in 1:nbootstrap){ test-list()

Re: [R] Unsubscribe please

On Apr 16, 2013, at 11:02 PM, Pascal Oettli wrote: Hello, Don't reply only to me. 1) Filter the unwanted mails, There is also an option to get Rhelp postings in digest-format. -- David. 2) It takes few days to unsubscribe you. Regards, Pascal On 04/17/2013 02:59 PM, bert

Re: [R] use of names() within lapply()

On 17/04/2013 11:33 AM, Ivan Alves wrote: Dear Duncan and A.K. Many thanks for your super quick help. The modified lapply did the trick, mapply died with a error Error in dots[[2L]][[1L]] : object of type 'builtin' is not subsettable. That's due to a typo: I should have said mapply(plot,

Re: [R] Best way to calculate averages of Blocks in an matrix?

On Apr 17, 2013, at 9:54 AM, Keith S Weintraub wrote: Folks, I recently was given a simulated data set like the following subset: sim_sub-structure(list(V11 = c(0.01, 0, 0, 0.01, 0, 0.01, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), V12 = c(0, 0, 0, 0.01, 0.03, 0, 0, 0, 0, 0, 0, 0.01, 0,

Re: [R] Best way to calculate averages of Blocks in an matrix?

 tapply(t(data.matrix(sim_sub)),rep( rep(1:4, each=5), each=10),mean)    #  1  2  3  4 #0.0086 0.0074 0.0082 0.0108 unlist(lapply(split(sim_sub,((seq_len(nrow(sim_sub))-1)%/%5)+1),function(x) mean(unlist(x #    1  2  3  4 #0.0086 0.0074 0.0082 0.0108 A.K.

[R] Full Information Maximum Likelihood estimation method for multivariate sample selection problem

Dear R experts/ users Full Information Maximum Likelihood (FIML) estimation approach is considered robust over Seemingly Unrelated Regression (SUR) approach for analysing data of multivariate sample selection problem. The zero cases in my dependent variables are resulted from three sources:

[R] Multi-core processing in glmulti

Dear list, I am trying to do an automated model selection of a glmm (function glmer; package: lme4) containing a large number of predictors. As far as i understand, glmulti is able to devide the process into chuncks and proceed by parallel processing on on multiple cores. Unfortunately this does

[R] t-statistic for independent samples

Hi, Typical things you read when new to stats are cautions about using a t-statistic when comparing independent samples. You are steered toward a pooled test or welch's approximation of the degrees of freedom in order to make the distribution a t-distribution. However, most texts give no

Re: [R] t-statistic for independent samples

On 04/17/2013 06:24 PM, David Arnold wrote: Hi, Typical things you read when new to stats are cautions about using a t-statistic when comparing independent samples. You are steered toward a pooled test or welch's approximation of the degrees of freedom in order to make the distribution a

Re: [R] t-statistic for independent samples

Dear David, On Wed, Apr 17, 2013 at 6:24 PM, David Arnold dwarnol...@suddenlink.net wrote: Hi, [snip] D. Before posting to StackExchange, check out the Wikipedia entry for Behrens-Fisher problem. Cheers, Jay -- G. Jay Kerns, Ph.D. Youngstown State University

Re: [R] t-statistic for independent samples

OK,although the variance ratio was already 2.25 to 1, tried sigma1=10, sigma2=25, which makes the ratios of the variances 6.25 to 1. Still no change. See: http://msemac.redwoods.edu/~darnold/math15/R/chapter11/DistributionForTwoIndependentSamplesPartII.html D. -- View this message in

Re: [R] Understanding why a GAM can't suppress an intercept

Simon, Many thanks as always for your help. I see and appreciate the example that you cited, but I'm having a hard time generalizing it to a multivariate case. A bit about my context -- my dataset is response ratios; the log of a treatment over a control. One of my explanatory variables is

[R] Memory usage reported by gc() differs from 'top'

In help(gc) I read, ...the primary purpose of calling 'gc' is for the report on memory usage. What memory usage does gc() report? And more importantly, which memory uses does it NOT report? Because I see one answer from gc(): used (Mb) gc trigger (Mb) max used (Mb) Ncells

Re: [R] t-statistic for independent samples

I just looked more carefully at your code. You are computing the unequal-variance (Welch) version of the t-test, so that's why there isn't a problem. Compare it with the equal-variance t-test, using the pooled variance estimate, which does have a problem, as below -thomas tstat4 -

[R] Using different function (parameters) with apply

Hi All, I have the following problem (read the commented bit below): a-matrix(1:9,nrow=3) a [,1] [,2] [,3] [1,]147 [2,]258 [3,]369 div-1:3 apply(a,2,function(x)x/div) ##want to divide each column by div- instead each row is divided## [,1]