[R] question re: "summarry.lm" and NA values

Petr Pikal petr.pikal at precheza.cz
Tue Aug 15 17:32:16 CEST 2006


Hi

just as a quick workaround you probably can use aliased value from 
summary

fff<-rep(summary(reg)$aliased,4)
dim(fff)<-c(no.of.your.variables,4)
fff[which(fff)]<-NA
fff[which(!fff)]<-coef(summary(reg))

to get coef matrix with NA values

HTH
Petr


On 15 Aug 2006 at 17:15, r user <ruser2006 at yahoo.com>, wrote:

From:           	Petr Pikal <petr.pikal at precheza.cz>
To:             	r user <ruser2006 at yahoo.com>, rhelp <r-help at stat.math.ethz.ch>
Subject:        	Re: [R] question re: "summarry.lm" and NA values
Date sent:      	Tue, 15 Aug 2006 17:15:01 +0200

> Hi
> 
> On 15 Aug 2006 at 7:01, r user wrote:
> 
> Date sent:      	Tue, 15 Aug 2006 07:01:13 -0700 (PDT)
> From:           	r user <ruser2006 at yahoo.com>
> To:             	rhelp <r-help at stat.math.ethz.ch>
> Subject:        	[R] question re: "summarry.lm" and NA values
> 
> > Is there a way to get the following code to include
> > NA values where the coefficients are  NA ?
> > 
> > ((summary(reg))$coefficients)
> 
> better
> coef(reg)
> 
> > 
> > explanation:
> > 
> > Using a loop, I am running regressions on several
> >  subsets  of  data1 .
> > 
> >  reg <- ( lm(lm(data1[,1] ~., data1[,2:l])) ) 
> > 
> > My regression has 10 independent variables, and I
> > therefore expect 11 coefficients.
> > After each regression, I wish to save the coefficients
> > and standard errors of the coefficients in a table
> > with 22 columns.
> > 
> > I successfully extract the coefficients using the
> > following code:
> >  reg$coefficients 
> > 
> > I attempt to extract the standard errors using :
> > 
> > aperm((summary(reg))$coefficients)[2,]
> > 
> > ((summary(reg))$coefficients)
> > 
> > My problem:
> > For some of my subsets, I am missing data for one or
> > more of the independent variables.  This of course
> > causes the coefficients and standard erros for this
> > variable to be  NA .
> 
> ??%^&*^??
> 
> What version? My lm behaves in accordance with na.action and it 
> throws an error in case na.fail, computes a value in case of na.omit
> or na.exclude and again throws an error if the variable consist
> exclusively from NA values. 
> 
> The only way how to get NA in coeficient is when a variable is either
> constant or linear combination of other variable(s). Then coef(reg)
> will give you correctly NA in the variable which appears constant and
> in this case you could use it for setting standard error also as NA
> let say by using ifelse statement and matching of names.
> 
> HTH
> Petr
> 
> > 
> > Is there a way to include the NA standard errors, so
> > that I have the same number of standard erros and
> > coefficients for each regression, and can then store
> > the coefficients and standard erros in my table of 22
> > columns?
> > 
> > ______________________________________________
> > R-help at stat.math.ethz.ch 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.
> 
> 

Petr Pikal
petr.pikal at precheza.cz



More information about the R-help mailing list