Dear Arun,Thanks indeed

> Date: Thu, 3 Oct 2013 10:22:38 -0700
> From: smartpink111@yahoo.com
> To: r-help@r-project.org
> Subject: Re: [R] vif
> 
> Hi Eliza,
> 
> Then, "res" needs a slight modification
> 
> 
> library(car)
>  res<- lapply(colnames(h),function(x) {x1<- h[,x];dat1<- do.call(rbind,lapply(seq_len(ncol(mat1)),function(i){ x2<- m[,mat1[,i]];GG<- lm(x1~x2[,1]+x2[,2]+x2[,3]+x2[,4]);GGsum<- summary(GG); data.frame( Models=paste(colnames(x2),collapse=","), Multiple_Rsq= GGsum$r.squared, Adjusted_Rsq = GGsum$adj.r.squared, Pval = paste(GGsum$coef[-1,4],collapse=","),Vif=paste(vif(GG),collapse=","),stringsAsFactors=FALSE)  })); dat1[rev(order(dat1[,3])),][1:10,]})
> 
> names(res)<- colnames(h)
> A.K.
> 
> 
> 
> 
> 
> ________________________________
> From: eliza botto <eliza_botto@hotmail.com>
> To: "smartpink111@yahoo.com" <smartpink111@yahoo.com> 
> Sent: Thursday, October 3, 2013 12:42 PM
> Subject: vif
> 
> 
> 
> Dear Arun,
> There is one small question however. 
> what if i also want in the table a column for "vif" values of each model. 
> vif values can be generated for any model in the following way
> 
> GG<-lm(h[,any column]~m[,any column]+m[,any other column] +m[,any other column] +m[,any other column])
> library(car)
> vif(GG)
> 
> Here will be get 4 "vif" values. I want to make a new column which could contain these values seperated by comma, very much similar to the way we did with Pr(>|t|) values.
> 
> thanks in advance
> 
> elisa
> 
> 
> > Date: Thu, 3 Oct 2013 09:01:53 -0700
> > From: smartpink111@yahoo.com
> > To: r-help@r-project.org
> > Subject: Re: [R] a simple question
> > 
> > Hi,
> > Try:
> > 
> > 
> > set.seed(494)
> >  h<- matrix(sample(1:40,4*124,replace=TRUE),ncol=4)
> > 
> >  set.seed(39)
> >  m<- matrix(sample(1:100,10*124,replace=TRUE),ncol=10)
> >  colnames(h)<- paste0("h",1:4)
> >  colnames(m)<- paste0("m",1:10)
> > mat1<-combn(colnames(m),4)
> > 
> >  res<- lapply(colnames(h),function(x) {x1<- h[,x];dat1<- do.call(rbind,lapply(seq_len(ncol(mat1)),function(i){ x2<- m[,mat1[,i]];GG<- lm(x1~x2[,1]+x2[,2]+x2[,3]+x2[,4]);GGsum<- summary(GG); data.frame( Models=paste(colnames(x2),collapse=","), Multiple_Rsq= GGsum$r.squared, Adjusted_Rsq = GGsum$adj.r.squared, Pval = paste(GGsum$coef[-1,4],collapse=","),stringsAsFactors=FALSE)  })); dat1[rev(order(dat1[,3])),][1:10,]})
> > 
> > names(res)<- colnames(h)
> > 
> > 
> > A.K.
> > 
> > 
> > 
> > 
> > 
> > ________________________________
> > From: eliza botto <eliza_botto@hotmail.com>
> > To: "smartpink111@yahoo.com" <smartpink111@yahoo.com> 
> > Sent: Thursday, October 3, 2013 11:07 AM
> > Subject: a simple question
> > 
> > 
> > 
> > 
> > Dear Arun,
> > I hope you are fine. I actually
> > wanted to discuss the following problem.
> > I have a linear model of the
> > following form. 
> > GG<-lm(h[,any column]~m[,any
> > column]+m[,any other column] +m[,any other column] +m[,any other column])
> > where,
> > h is matrix with 4 columns and
> > 124 rows
> > m is matrix with 10 columns and
> > 124 rows
> > what I want is the following
> > make a loop command to run the
> > linear model of all the possible combinations of columns of “m” with each
> > column of “h”. 
> > more precisely, if i take column
> > 1 of matrix “h”, it should be linear modeled with every combination of 10 (210
> > combinations) columns of “m”.
> > All the columns of “h” & “m”
> > have certain names (you can suppose any).  The summary(GG) will give Multiple R-squared,    Adjusted R-squared  and 4 values of Pr(>|t|). 
> > I want in the end a table in the
> > following format.
> > 
> > Models                                                                           Multiple R-squared        Adjusted R-squared       Pr(>|t|)
> > Name of columns of m separated by
> > comma      Multiple R-squared         Adjusted R-squared       Pr(>|t|) separated by comma
> > 
> > For Example 
> > 
> > Models                                                                           Multiple R-squared        Adjusted R-squared       Pr(>|t|)
> > eliza, allen, murphy, jack                                                 0.544                                     0.56                         0.000114,0.000112,0.01114,0.002114
> > 
> > where,
> > eliza, allen, murphy, jack are column names.
> > 
> > The models are to be enlisted in the order of their Adjusted R-squared values. The models with highest Adjusted R-squared value should be on the top and so on. i m only interested in top 10 models. so the remaining should be ignored. 
> > 
> > I tried to put in my question everything but if there is anything wrong plz inform me.
> > 
> > Thankyou very much in advance,
> > 
> > Eliza
> > 
> > ______________________________________________
> > R-help@r-project.org mailing list
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
> 
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
 		 	   		  
	[[alternative HTML version deleted]]

