r - prevent the name of model.matrix to appear in the results on a regression -

-

is possible use model matrix in regression without having name of model matrix in regression results?

i need go through such procedure have interactions have no observation. (i.e) results of interaction na.

a related question can found here.

here data illustrate point:

mydata <- read.csv("http://www.ats.ucla.edu/stat/data/binary.csv") str(mydata)  gre_ <- mydata$gre-mean(mydata$gre)  <- model.matrix(~-1+gre_:factor(rank),data=mydata)[,-c(2)]  summary(glm(admit~gpa+gre+factor(rank)+a,data=mydata, family=binomial))

results

call: glm(formula = admit ~ gpa + gre + rank + a, family = binomial,  data = mydata)  deviance residuals:      min       1q   median       3q      max   -1.6449  -0.8886  -0.6332   1.1706   2.1949    coefficients:                       estimate std. error z value pr(>|z|)     (intercept)         -3.0039781  1.4012928  -2.144   0.0321 *   gpa                  0.7634679  0.3297215   2.315   0.0206 *   gre                  0.0016098  0.0016634   0.968   0.3332     rank                -0.5584921  0.1288588  -4.334 1.46e-05 *** agre_:factor(rank)1  0.0014010  0.0028001   0.500   0.6168     agre_:factor(rank)3  0.0010074  0.0025007   0.403   0.6871     agre_:factor(rank)4  0.0009936  0.0034111   0.291   0.7708     ---

how rid of model.matrix name a in results?

if run formula syntax, r going put "a" there. can extract names of coefficients , remove first "a" if gsub() or use substr() remove first letter. depends on how want process them down line.

the other option use glm.fit , specify complete model matrix yourself. like

a <- model.matrix(~-1+gre_:factor(rank),data=mydata)[,-c(2)] b <- model.matrix(~gpa+gre+rank, data=mydata) mm<-cbind(b,a)  ff<-glm.fit(mm,mydata$admit, family=binomial()) class(ff)<-c("glm","lm") summary(ff)

will return

call: null  deviance residuals:      min       1q   median       3q      max   -1.6449  -0.8886  -0.6332   1.1706   2.1949    coefficients:                      estimate std. error z value pr(>|z|)     (intercept)        -3.0039781  1.4012928  -2.144   0.0321 *   gpa                 0.7634679  0.3297215   2.315   0.0206 *   gre                 0.0016098  0.0016634   0.968   0.3332     rank               -0.5584921  0.1288588  -4.334 1.46e-05 *** gre_:factor(rank)1  0.0014010  0.0028001   0.500   0.6168     gre_:factor(rank)3  0.0010074  0.0025007   0.403   0.6871     gre_:factor(rank)4  0.0009936  0.0034111   0.291   0.7708     --- signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1  (dispersion parameter binomial family taken 1)      null deviance: 499.98  on 399  degrees of freedom residual deviance: 459.13  on 393  degrees of freedom aic: 473.13  number of fisher scoring iterations: 4

here estimates same , variable names untouched. since not technically real glm object fibbing bit adding class information, have same properties (you can see "call" missing) , in cases should behave regular glm object functions, including summary().


Comments

Post a Comment

Popular posts from this blog

C# random value from dictionary and tuple -

-
i have following code, , need random item contains string "region_1" dictionary. public static dictionary<int, tuple<string, int, currencytype>> itemarray = new dictionary<int, tuple<string, int, currencytype>>() { { 0xb3e, new tuple<string, int, currencytype>("region_1", 1500, currencytype.fame) } }; public static int generateitemid(string shopid) { var generateditem = itemarray.elementat(new random().next(0, itemarray.count)).key; } how select this? it's not possible efficiently... first create static random @ class level... prevent non-random behaviour if run query on short period of time... (it's seeded discrete clock) static random rnd = new random(); then: var item = itemarray.values .where(t => t.item1 == shopid) .orderby(_ => rnd.next()) .firstordefault()
Read more

cgi - How do I interpret URLs without extension as files rather than missing directories in nginx? -

-
so i'm trying install zoneminder , have run under nginx, has few compiled files need run using fast-cgi. these files lack extension. if file has extension, there no issue , nginx interprets file , return 404 if can't find it. if has no extension, assume it's directory , return 404 when can't find sort of index page. here have now: # security camera domain. server { listen 888; server_name mydomain.com; root /srv/http/zm; # enable php support. include php.conf; location / { index index.html index.htm index.php; } # enable cgi support. location ~ ^/cgi-bin/(.+)$ { alias /srv/cgi-bin/zm/; fastcgi_pass unix:/run/fcgiwrap.sock; fastcgi_param script_filename $document_root/$1; include fastcgi.conf; } } the idea is, under cgi-bin directory goes through fastcgi pipe. any idea on howi can fix this? upon closer inspection, realized zoneminder has 2 cgi sc...
Read more

.htaccess - htaccess convert request to clean url and add slash at the end of the url -

-
ok , searched more 4 days , more 20 posts here in stach or in websites without kind of success first here htaccess file rewriteengine on rewriterule ^([a-za-z0-9_]+)/?$ index.php?page=$1 [nc] rewriterule ^([a-za-z0-9_]+)\/([a-za-z0-9_]+)\/?$ index.php?page=$1&action=$2 [nc] i tried many codes 1 worked me can call request , works me www.example.com/pagename => www.example.com/index.php?page=pagename here questions , please provide me information "in detail" understand code well 1- how convert request automatically www.example.com/index.php?page=pagename www.example.com/pagename 2- if tried add slash @ end shows me page without images or styling mean www.example.com/pagename => works www.example.com/pagename/ => doesn't work that's :) for first part need 301 rule this: rewritecond %{the_request} /index\.php\?page=([^\s&]+) [nc] rewriterule ^ /%1? [r=301,l] for second part use absolute path in css, js, images files r...
Read more