Re: [R] AIC in MuMIn

2010-08-19 Thread Gavin Simpson
On Thu, 2010-08-19 at 04:42 +0800, elaine kuo wrote: Yes, I tried the example in the ?dredge and agreed that something else caused the mistake. Aside from the cause which takes time to clarify (16 explanatory variables in the model), I would like to ask another question. Please

Re: [R] AIC in MuMIn

2010-08-19 Thread elaine kuo
Please kindly advise if it is possible to show the singular model with only one certain variable using the command subset. (or maybe others) I tried the command subset=X3 but it returned multiple models including X3. The above demand might look unnecessary when visual inspection

Re: [R] AIC in MuMIn

2010-08-19 Thread Gavin Simpson
On Fri, 2010-08-20 at 06:31 +0800, elaine kuo wrote: snip / Oh, actually, I suppose you could automate this, so it will return all models with single variable: dd - dredge(lm1) parms - !is.na(dd[, -c(1, (ncol(dd) -

Re: [R] AIC in MuMIn

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 05:42 +0800, elaine kuo wrote: Thank you. Most of the answers solved the puzzles. Q2 I tried to display sub-model with only temp_ran using the code below but failed. Please kindly suggest the potential failure cause.

Re: [R] AIC in MuMIn

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 08:51 +0100, Gavin Simpson wrote: On Wed, 2010-08-18 at 05:42 +0800, elaine kuo wrote: Thank you. Most of the answers solved the puzzles. snip / = Please explain why fitted lm is better for GLM. Seriously? A GLM specified as glm(, family = gaussian) is the

Re: [R] AIC in MuMIn

2010-08-18 Thread elaine kuo
Hmmm. Thinking some more, I might not have answered your (unstated) question. What do your mean by GLM? = I meant generalized linear model as well. Thanks for the references. The first one was mentioned first in my life time after keeping asking the same question. I mean the Generalized

Re: [R] AIC in MuMIn

2010-08-18 Thread elaine kuo
Please suggest how to define subset in my case How would I know? I still haven't seen your data. You seem to be mistaken on what is and is not included in your model and you fitted it. What hope do we have...? However, given the model 'mig.stds' from above in this email: mig.stds

Re: [R] AIC in MuMIn

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 17:37 +0800, elaine kuo wrote: snip / = I posted it for help, after following the manual with the command dredge but receiving an error message two days ago. command target-dredge(mig.stds, subset = temp_max) error in eval(expr,

Re: [R] AIC in MuMIn

2010-08-18 Thread elaine kuo
A cause other than data based on standardized regression was identified. It is that the manual command added with target - at the left hand side. C1 did not work but C2 did. C1 target-dredge(mig.stds, subset = temp_max) C2 dredge(mig.stds, subset = temp_max) Elaine On Wed, Aug 18, 2010

Re: [R] AIC in MuMIn

2010-08-18 Thread Gavin Simpson
On Wed, 2010-08-18 at 21:11 +0800, elaine kuo wrote: A cause other than data based on standardized regression was identified. It is that the manual command added with target - at the left hand side. C1 did not work but C2 did. C1 target-dredge(mig.stds, subset = temp_max) C2

Re: [R] AIC in MuMIn

2010-08-18 Thread elaine kuo
Yes, I tried the example in the ?dredge and agreed that something else caused the mistake. Aside from the cause which takes time to clarify (16 explanatory variables in the model), I would like to ask another question. Please kindly advise if it is possible to show the singular model with only

Re: [R] AIC in MuMIn

2010-08-17 Thread Gavin Simpson
On Tue, 2010-08-17 at 16:05 +0800, elaine kuo wrote: Hello, Why did you decide to post the exact same message from two different email addresses?? I am using package MuMIn to calculate AIC for a full model with 10 explanatory variables. Thanks in advance in sharing your experience. Q1 In

Re: [R] AIC in MuMIn

2010-08-17 Thread elaine kuo
Thank you. Most of the answers solved the puzzles. Q2 I tried to display sub-model with only temp_ran using the code below but failed. Please kindly suggest the potential failure cause. code library(MuMIn) datam -read.csv(c:/migration/Mig_ratio_20100817.csv,header=T,