r - Adding an column for the category of glm coeffients in broom results -


is there way add column result of broom package's tidy function can act relate term column both original names used in formula argument , columns in data argument.

for example if run following get:

library(ggplot2) library(dplyr)  mod <- glm(mpg ~ wt + qsec + as.factor(carb), data = mtcars)  tidy(mod)  #               term     estimate std.error   statistic      p.value # 1      (intercept) 21.132995090 7.5756463  2.78959633 1.017187e-02 # 2               wt -4.916303175 0.6747590 -7.28601380 1.584408e-07 # 3             qsec  0.843355538 0.3930252  2.14580532 4.221188e-02 # 4 as.factor(carb)2  0.004133826 1.5321134  0.00269812 9.978695e-01 # 5 as.factor(carb)3 -0.755346006 2.3451222 -0.32209239 7.501715e-01 # 6 as.factor(carb)4 -0.489721798 2.0628564 -0.23739985 8.143615e-01 # 7 as.factor(carb)6 -0.886846134 3.4443957 -0.25747510 7.990068e-01 # 8 as.factor(carb)8 -0.894783610 3.7496630 -0.23863041 8.134180e-01

what looking this:

#               term     estimate std.error   statistic      p.value   term_base # 1      (intercept) 21.132995090 7.5756463  2.78959633 1.017187e-02  # 2               wt -4.916303175 0.6747590 -7.28601380 1.584408e-07          wt # 3             qsec  0.843355538 0.3930252  2.14580532 4.221188e-02        qsec # 4 as.factor(carb)2  0.004133826 1.5321134  0.00269812 9.978695e-01        carb # 5 as.factor(carb)3 -0.755346006 2.3451222 -0.32209239 7.501715e-01        carb # 6 as.factor(carb)4 -0.489721798 2.0628564 -0.23739985 8.143615e-01        carb # 7 as.factor(carb)6 -0.886846134 3.4443957 -0.25747510 7.990068e-01        carb # 8 as.factor(carb)8 -0.894783610 3.7496630 -0.23863041 8.134180e-01        carb

not bothered if first row in new column empty, intercept or 1. need can match term column original variable names passed formula?

edit

would if didn't depend on using as.factor in formula, e.g. work on:

mod <- glm(mpg ~ wt + qsec + carb, data = mtcars %>% mutate(carb = factor(carb)))  tidy(mod)  #          term     estimate std.error   statistic      p.value # 1 (intercept) 21.132995090 7.5756463  2.78959633 1.017187e-02 # 2          wt -4.916303175 0.6747590 -7.28601380 1.584408e-07 # 3        qsec  0.843355538 0.3930252  2.14580532 4.221188e-02 # 4       carb2  0.004133826 1.5321134  0.00269812 9.978695e-01 # 5       carb3 -0.755346006 2.3451222 -0.32209239 7.501715e-01 # 6       carb4 -0.489721798 2.0628564 -0.23739985 8.143615e-01 # 7       carb6 -0.886846134 3.4443957 -0.25747510 7.990068e-01 # 8       carb8 -0.894783610 3.7496630 -0.23863041 8.134180e-01

we can use regex create 'term_base' column

tidy(mod) %>%         mutate(term_base = sub("intercept", "", gsub(".*\\(|\\).*", "", term))) #              term     estimate std.error   statistic      p.value term_base #1      (intercept) 21.132995090 7.5756463  2.78959633 1.017187e-02           #2               wt -4.916303175 0.6747590 -7.28601380 1.584408e-07        wt #3             qsec  0.843355538 0.3930252  2.14580532 4.221188e-02      qsec #4 as.factor(carb)2  0.004133826 1.5321134  0.00269812 9.978695e-01      carb #5 as.factor(carb)3 -0.755346006 2.3451222 -0.32209239 7.501715e-01      carb #6 as.factor(carb)4 -0.489721798 2.0628564 -0.23739985 8.143615e-01      carb #7 as.factor(carb)6 -0.886846134 3.4443957 -0.25747510 7.990068e-01      carb #8 as.factor(carb)8 -0.894783610 3.7496630 -0.23863041 8.134180e-01      carb

the as.factor can removed 'term' if mutate 'carb' factor before glm step

mtcars %>%      mutate(carb = factor(carb)) %>%       glm(formula = mpg ~wt + qsec + carb, data = .) %>%       tidy(.) %>%      mutate(term_base = sub("\\(.*\\)|\\d+", "", term)) #     term     estimate std.error   statistic      p.value term_base #1 (intercept) 21.132995090 7.5756463  2.78959633 1.017187e-02           #2          wt -4.916303175 0.6747590 -7.28601380 1.584408e-07        wt #3        qsec  0.843355538 0.3930252  2.14580532 4.221188e-02      qsec #4       carb2  0.004133826 1.5321134  0.00269812 9.978695e-01      carb #5       carb3 -0.755346006 2.3451222 -0.32209239 7.501715e-01      carb #6       carb4 -0.489721798 2.0628564 -0.23739985 8.143615e-01      carb #7       carb6 -0.886846134 3.4443957 -0.25747510 7.990068e-01      carb #8       carb8 -0.894783610 3.7496630 -0.23863041 8.134180e-01      carb

-

Comments

Post a Comment

Popular posts from this blog

amazon web services - S3 Pre-signed POST validate file type? -

i know it's possible put limit on file size content-length-range header. possible validate file type? https://docs.aws.amazon.com/amazons3/latest/dev/httppostforms.html#policyconditions i see there content-type header, if set say, audio/mp3 allow mp3 files , return error if file not mp3? i found previous question answers mention validating file size: s3 direct upload restricting file size , type you can specify content-type @ post request. you can specify @ signature, post must done content-type: { "expiration": "2007-12-01t12:00:00.000z", "conditions": [ {"acl": "public-read" }, {"bucket": "johnsmith" }, {"content-type: "audio/mp3"} ] } creating html form (using aws signature version 4) edit: @ previous question found, checking content-type.
Read more

c# - Check Keyboard Input Winforms -

i wondering if instead of doing this protected override void onkeydown(keyeventargs e) { if (e.keycode == keys.a) console.writeline("the key down."); } i set bool method , this: if(keydown(keys.a)) // whatever here i've been sitting here ages trying figure out how it. can't wrap head around it. in case wondering, plan call bool inside different method, check input. since want perform action after pressing key, using keydown event enough. but in cases suppose want check if specific key down in middle of process, can use getkeystate method way: [dllimport("user32.dll", charset = charset.auto, exactspelling = true)] public static extern short getkeystate(int keycode); public const int key_pressed = 0x8000; public static bool iskeydown(keys key) { return convert.toboolean(getkeystate((int)key) & key_pressed); } you should know, each time check key state using example iskeydown(keys.a) method returns true if ke...
Read more