search for: fdrop

Dear R users, Would it be difficult to change the code of stepAIC (from the MASS library) to use AICc instead of AIC? It would be great to know of someone has tried this already. Best wishes Christoph.

...;s version of glm") md <- missing(direction) direction <- match.arg(direction) backward <- direction == "both" | direction == "backward" forward <- direction == "both" | direction == "forward" if (missing(scope)) { fdrop <- numeric(0) fadd <- attr(Terms, "factors") if (md) forward <- FALSE } else { if (is.list(scope)) { fdrop <- if (!is.null(fdrop <- scope$lower)) attr(terms(update.formula(object, fdrop)), "factors...

...g code: > step(glm00, scope=formula(glm1), method="both") and another one: > step(glm00, scope=formula(glm1), method="forward") In both cases step procedure terminated after several steps with a warning: > Error in factor.scope(ffac, list(add = fadd, drop = fdrop)) : > upper scope does not include model What does the word "model" refer to? To glm00 (the null model)? To my knowledge and understanding, the initial null model should have been included in every step of the search. Am I wrong? Thanks for all comments :-) Konr...

...al+msa+racekid+racemom+sex+shotcard )^2) lower <- as.formula(haspdata~1) lst <- list( upper=upper , lower=lower ) form1a.step <- step(lg.mod1, scope=lst , trace=T , direction="both" ) I get the following error message in R: Error in factor.scope(ffac, list(add = fadd, drop = fdrop)) : upper scope does not include model Because I am trying to migrate from Splus to R, I ran that same code in Splus and do not get an error. Could someone please tell me: 1) what I'm doing wrong in R, and 2) what the distinction between R and Splus is in this case. Thank you, Ph...

...at ../../../netinet/tcp_usrreq.c:443 #11 0xc0689822 in sodisconnect (so=0x0) at ../../../kern/uipc_socket.c:576 #12 0xc0689490 in soclose (so=0xc19ec164) at ../../../kern/uipc_socket.c:457 #13 0xc0678d17 in soo_close (fp=0xc1736c60, td=0xc1730c00) at ../../../kern/sys_socket.c:317 #14 0xc062e818 in fdrop_locked (fp=0xc1736c60, td=0xc1730c00) at file.h:289 #15 0xc062e769 in fdrop (fp=0xc1736c60, td=0xc1730c00) at ../../../kern/kern_descrip.c:2112 #16 0xc062cd97 in closef (fp=0xc1736c60, td=0xc1730c00) at ../../../kern/kern_descrip.c:1932 #17 0xc062a175 in close (td=0xc1730c00, uap=0x12) at ../../../...

...914.34 225.30 Step: AIC= 204.71 y ~ x3 + x1 Df Sum of Sq RSS AIC + x1:x3 1 533.91 195.49 75.03 + x2 1 47.44 681.96 199.98 <none> 729.40 204.71 Step: AIC= 75.03 y ~ x3 + x1 + x3:x1 Error in factor.scope(ffac, list(add = fadd, drop = fdrop)) : upper scope does not include model Thanks, Min [[alternative HTML version deleted]]

Dear all: I have a shapefile in UTM coordinate system and I would like to transform it to Lat/Log coordinates (WSG84). The package PBSmapping contains function convUL to transform between the two coordinate systems when data is in the form of a data frame with attributes specifying the coordinate system. However, when shapefiles are imported using function read.shape (package maptools), a list

I can only assume I'm betraying my ignorance here, but this is not what I would expect. I'm getting the following from a stepwise selection (with both step and stepAIC): > step(lm(sqrt(Grids)~ SE + Edge + NH), scope=~ (Edge + SE + NH)^2) Start: AIC= 593.56 sqrt(Grids) ~ SE + Edge + NH Df Sum of Sq RSS AIC <none> 2147.0 593.6 + Edge:NH 1

Dear list, here is an example of stepAIC that I do not understand. The data is n=42, Lage is the only factor and there are four other variables treated as continuous. First you see the stepAIC-forward solution (fs7). The strange thing here is that apparently not all interactions are tried for inclusion, but only WQ:Lage. In particular, I think that WFL:Lage should be tried in the last two

...viance AIC + IA:Region 8 540.24 1377.73 + IM 1 586.84 1410.32 + altura:Region 8 574.21 1411.70 <none> 600.74 1422.23 Step: AIC= 1377.73 MCHNV ~ Region + IA + altura + Region:IA + offset(Offset1) Error in factor.scope(ffac, list(add = fadd, drop = fdrop)) : upper scope does not include model ------------------------------------------------------------------------------- What am I doing wrong? Thank you for your help. I put here the database I use, so you can replicate the problem. (I am using R221 on WinXP) ---------------------------...