
Hi Dennis,
 
Thanks so much, I will try change the code and see what happens. 
 
Many thanks,
Candice

PhD. Student
Centre for Invasion Biology
Department of Botany and Zoology
Stellenbosch University
Cell: 083 987 2200


 



Date: Mon, 24 Oct 2011 09:18:02 -0700
From: ml-node+s789695n3933698h72@n4.nabble.com
To: candice.lyons@hotmail.com
Subject: Re: bestglm function and output in R

Hi: 

bestglmtest is your input data frame, is it not? From the names() line, 
you can see that it has no variable named BestModel that corresponds 
to a list containing a component named coefficients. 

Were you perhaps looking for 
output$BestModel$coefficients  ?? 

Dennis 

On Mon, Oct 24, 2011 at 8:35 AM, Candice <[hidden email]> wrote: 

> Hi all, 
> 
> I have been trying to run a bestglm in R for a while now and am struggling 
> to get it to run. When I thought I had succeeded, the "output" it gave me 
> was "NULL" and that's it. Below is my code: 
> 
> bestglmtest<-read.table("C:\\Documents and 
> Settings\\clyons\\Desktop\\bestglmtest.txt",header=T,row.names=1) 
>> names(bestglmtest) 
> [1] "acc"   "age"   "sex"   "ctmin" 
>> library(bestglm) 
>> output<-bestglm(bestglmtest,IC="AIC",family=gaussian(link=identity)) 
> Morgan-Tatar search since family is non-gaussian. 
>> bestglmtest$BestModel$coefficients 
> NULL 
> 
> I believ my data is Gaussian so don't get the Morgan-Tatar search 
> "error/warning" which it threw out. Trying to find info on it was also not 
> successful but I did manage to find an article with example which also had 
> the same message and figured out it was not an error message but a warning, 
> if I can say that. Asking for the best mode coefficients however, did not 
> provide me with any useful answers. 
> 
> What I want to do is predict how four factors/parameters (this example is 
> only three because I was testing the code) influence my response variable. I 
> have my response variable in the very last column of my text file with 
> predictors before it. Two of my predictors have two levels and the other two 
> have three levels. Perhaps this makes a difference to the ditribution type I 
> select? 
> 
> Any help would be greatly appreciated since I am now rightly stuck and have 
> no idea what else to try. 
> 
> Many thanks, 
> Candice 
> 
> -- 
> View this message in context: http://r.789695.n4.nabble.com/bestglm-function-and-output-in-R-tp3933547p3933547.html
> Sent from the R help mailing list archive at Nabble.com. 
> 
> ______________________________________________ 
> [hidden email] 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. 
> 
______________________________________________ 
[hidden email] 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. 






If you reply to this email, your message will be added to the discussion below:http://r.789695.n4.nabble.com/bestglm-function-and-output-in-R-tp3933547p3933698.html 
To unsubscribe from bestglm function and output in R, click here.  		 	   		  

--
View this message in context: http://r.789695.n4.nabble.com/bestglm-function-and-output-in-R-tp3933547p3935845.html
Sent from the R help mailing list archive at Nabble.com.
	[[alternative HTML version deleted]]

