r - Logical Subscript too long -

i realized question has been asked before, having looked @ answers, specific questions, , not find answer unique situation.

i entered following r, , worked first example, not second, , cannot understand why.

setting data glm :

setwd("p:/stat319") ucb2<-read.table('berkeley.poissontwo.txt',header=true) attach(ucb2) 

ucb2 following :

count   admit department    gender      313 false         female      512 true          female      19  false         male        89  true          male        207 false     b     female      353 true      b     female      8   false     b     male        17  true      b     male        205 false     c     female      120 true      c     female      391 false     c     male        202 true      c     male        279 false     d     female      138 true      d     female      244 false     d     male        131 true      d     male        138 false     e     female      53  true      e     female      299 false     e     male        94  true      e     male        351 false   f       female      22  true      f     female      317 false     f     male        24  true      f     male 

using factor variable, true , false admit , notadmit :

admit<-c(0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1,0,1) fadmit<-factor(admit) radmit<-factor(admit,labels=c("false","true")) glm2<-glm(count~admit+department+gender,family=poisson) glm2 

preparing way leave 1 out cross validation

library(car) vif(glm2) # gvif df gvif^(1/(2*df)) # admit         1  1               1 # department    1  5               1 # gender        1  1               1 step(glm2) # start:  aic=2272.73 # count ~ admit + department + gender #  # df deviance    aic # <none>            2097.7 2272.7 # - department  5   2257.2 2422.2 # - gender      1   2260.6 2433.6 # - admit       1   2327.7 2500.8 #  # call:  glm(formula = count ~ admit + department + gender, family = poisson) #  # coefficients: #   (intercept)        admit  departmentb  departmentc   # 5.82785     -0.45674     -0.46679     -0.01621   # departmentd  departmente  departmentf   gendermale   # -0.16384     -0.46850     -0.26752     -0.38287    # degrees of freedom: 23 total (i.e. null);  16 residual # null deviance:        2650  # residual deviance: 2098   aic: 2273  library(ipred) errorest(count~admit+department+gender,data=ucb2,model=glm,est.para=control.errorest(k=24))  # call: #   errorest.data.frame(formula = count ~ admit + department + gender,  #                       data = ucb2, model = glm, est.para = control.errorest(k = # 24)) #  # 24-fold cross-validation estimator of root mean squared error #  # root mean squared error:  180.5741  

so first 1 worked data shown. same study, had rearrange data, , perform logistic regression :

ucb1<-read.table('monday.late.txt',header=true) attach(ucb1) # following object masked _by_ .globalenv: #    #   admit  # following objects masked ucb2: #    #   admit, department, gender  y<-cbind(ucb1[,1],ucb1[,2]) glm1<-glm(y~gender+department,family=binomial) 

the data follows :

admit   notadmit    gender  department      512 313 female        353 207 female  b       120 205 female  c       138 279 female  d       53  138 female  e       22  351 female  f       89  19  male          17  8   male    b       202 391 male    c       131 244 male    d       94  299 male    e       24  317 male    f     

setting new data leave 1 out :

vif(glm1) # gvif df gvif^(1/(2*df)) # gender     1.384903  1        1.176819 # department 1.384903  5        1.033099 step(glm1) # start:  aic=103.14 # y ~ gender + department  # df deviance    aic # - gender      1    21.74 102.68 # <none>             20.20 103.14 # - department  5   783.61 856.55 #  # step:  aic=102.68 # y ~ department #  # df deviance    aic # <none>             21.74 102.68 # - department  5   877.06 948.00 #  # call:  glm(formula = y ~ department, family = binomial) #  # coefficients: #   (intercept)  departmentb  departmentc  departmentd   # 0.59346     -0.05059     -1.20915     -1.25833   # departmente  departmentf   # -1.68296     -3.26911   #  # degrees of freedom: 11 total (i.e. null);  6 residual # null deviance:        877.1  # residual deviance: 21.74  aic: 102.7 

so far, good, problem arises :

errorest(y~gender+department,data=ucb1,model=glm,est.para=control.errorest(k=12)) error in xj[i, , drop = false] : (subscript) logical subscript long 

so why happen ? tried other values k, not sure value k # meant take - assume meant of number of rows

i try same data, arranged different way :

ucb1a<-read.table('berkeley.rearranged.txt',header=true) attach(ucb1a) ucb1a 

this rearrangement of data before

admitted not_admit depart genders 1       512       313       female 2        89        19         male 3       353       207      b  female 4        17         8      b    male 5       120       205      c  female 6       202       391      c    male 7       138       279      d  female 8       131       244      d    male 9        53       138      e  female 10       94       299      e    male 11       22       351      f  female 12       24       317      f    male 

and

y<-cbind(ucb1[,1],ucb1[,2]) glm1a<-glm(y~genders+depart,family=binomial) vif(glm1a) # gvif df gvif^(1/(2*df)) # gender     1.384903  1        1.176819 # department 1.384903  5        1.033099  step(glm1a) # start:  aic=103.14 # y ~ gender + department #  # df deviance    aic # - gender      1    21.74 102.68 # <none>             20.20 103.14 # - department  5   783.61 856.55 #  # step:  aic=102.68 # y ~ department #  # df deviance    aic # <none>             21.74 102.68 # - department  5   877.06 948.00 #  # call:  glm(formula = y ~ department, family = binomial) #  # coefficients: #   (intercept)  departmentb  departmentc  departmentd   # 0.59346     -0.05059     -1.20915     -1.25833   # departmente  departmentf   # -1.68296     -3.26911   #  # degrees of freedom: 11 total (i.e. null);  6 residual # null deviance:        877.1  # residual deviance: 21.74  aic: 102.7 

again, far good, once more, occurs :

errorest(y~gender+department,data=ucb1a,model=glm,est.para=control.errorest(k=12)) error in xj[i, , drop = false] : (subscript) logical subscript long 

and believe me, tried other numbers again k, , cannot understand why 1 going wrong. if has ideas, specific example of (subscript) logical subscript being long, please reply this.

this issue arise when objects of different size. think problem comes attach() i'm not certain.. try code without it, or try with(). should check why have use attach() first before using nicola pointed out. also, i'm not trying achieve it.

you can see in section of function following : practice

attach has side effect of altering search path , can lead wrong object of particular name being found. people forget detach databases.

in interactive use, preferable use of attach/detach, unless save()-produced file in case attach() (safety) wrapper load().