Thetazero Pubs
MMoore25named names factor iap yes int chr attr dont know list gss

I tried doing basic regressions with the variables ?colhomo? and ?colrac,? but was unsuccesful in the final step using Zelig. I got one of many errors that carried over throughout the assignment.


library(foreign)
GSS2014 <- read.dta ("C:/Users/clyve_000/Documents/Advanced Analytics/GSS2014OLD.dta")
library(Zelig)
library(DescTools)
library(dplyr)
library(knitr)
library(stargazer)
str(GSS2014)
## 'data.frame':    2538 obs. of  866 variables:
##  $ abany        : Factor w/ 5 levels "iap","yes","no",..: 2 3 3 NA 3 2 NA NA 3 2 ...
##  $ abdefect     : Factor w/ 5 levels "iap","yes","no",..: 2 3 3 NA 2 2 NA NA 3 2 ...
##  $ abhlth       : Factor w/ 5 levels "iap","yes","no",..: 2 2 2 NA 2 2 NA NA 2 2 ...
##  $ abnomore     : Factor w/ 5 levels "iap","yes","no",..: 2 3 3 NA 3 2 NA NA 3 2 ...
##  $ abpoor       : Factor w/ 5 levels "iap","yes","no",..: 2 3 3 NA 3 2 NA NA 3 2 ...
##  $ abrape       : Factor w/ 5 levels "iap","yes","no",..: 2 3 3 NA 2 2 NA NA 2 2 ...
##  $ absingle     : Factor w/ 5 levels "iap","yes","no",..: 2 3 3 NA 3 2 NA NA 3 2 ...
##  $ accptoth     : Factor w/ 9 levels "iap","many times a day",..: 6 4 NA NA 3 NA NA NA NA NA ...
##  $ acqntsex     : Factor w/ 5 levels "iap","sex w acquaintance",..: 2 NA NA NA NA NA NA NA NA NA ...
##  $ actassoc     : int  NA NA 5 4 NA 4 7 5 1 5 ...
##  $ actlaw       : Factor w/ 7 levels "iap","very likely",..: NA NA 4 4 NA 3 4 5 5 4 ...
##  $ adults       : int  1 2 1 3 2 2 2 3 1 2 ...
##  $ advfront     : Factor w/ 7 levels "iap","strongly agree",..: NA NA NA 2 NA NA 2 2 4 3 ...
##  $ affctlaw     : Factor w/ 7 levels "iap","very likely",..: NA NA 4 3 NA 4 5 5 5 4 ...
##  $ affrmact     : Factor w/ 7 levels "iap","strongly support pref",..: NA NA 5 4 NA 5 5 5 NA NA ...
##  $ age          : int  53 26 59 56 74 56 63 34 37 30 ...
##  $ aged         : Factor w/ 6 levels "iap","a good idea",..: 2 2 NA 4 2 NA 3 4 4 2 ...
##  $ agekdbrn     : int  NA NA 27 38 24 23 23 24 18 15 ...
##  $ aidscndm     : Factor w/ 5 levels "iap","yes","no",..: NA 2 NA 2 2 NA NA NA 3 2 ...
##  $ aidslook     : Factor w/ 5 levels "iap","yes","no",..: NA 2 NA 2 2 NA NA NA 2 2 ...
##  $ amancstr     : Factor w/ 7 levels "iap","very important",..: 5 4 NA NA 3 NA NA NA NA NA ...
##  $ ambetter     : Factor w/ 8 levels "iap","strongly agree",..: 2 4 NA NA 2 NA NA NA NA NA ...
##  $ ambornin     : Factor w/ 7 levels "iap","very important",..: 2 3 NA NA 2 NA NA NA NA NA ...
##  $ amchrstn     : Factor w/ 7 levels "iap","very important",..: 5 4 NA NA 3 NA NA NA NA NA ...
##  $ amcit        : Factor w/ 7 levels "iap","very important",..: 3 2 NA NA 2 NA NA NA NA NA ...
##  $ amcitizn     : Factor w/ 8 levels "iap","strongly agree",..: 3 2 NA NA 2 NA NA NA NA NA ...
##  $ amcult       : Factor w/ 8 levels "iap","agree strongly",..: 3 6 NA NA 6 NA NA NA NA NA ...
##  $ amenglsh     : Factor w/ 7 levels "iap","very important",..: 2 2 NA NA 2 NA NA NA NA NA ...
##  $ amfeel       : Factor w/ 7 levels "iap","very important",..: 2 2 NA NA 2 NA NA NA NA NA ...
##  $ amgovt       : Factor w/ 7 levels "iap","very important",..: 2 3 NA NA 2 NA NA NA NA NA ...
##  $ amlived      : Factor w/ 7 levels "iap","very important",..: 3 2 NA NA 2 NA NA NA NA NA ...
##  $ amownway     : Factor w/ 8 levels "iap","agree strongly",..: 3 5 NA NA 3 NA NA NA NA NA ...
##  $ amproud1     : Factor w/ 8 levels "iap","very proud",..: 2 2 NA NA 2 NA NA NA NA NA ...
##  $ amshamed     : Factor w/ 8 levels "iap","strongly agree",..: 2 6 NA NA 6 NA NA NA NA NA ...
##  $ amsports     : Factor w/ 8 levels "iap","strongly agree",..: 2 2 NA NA 2 NA NA NA NA NA ...
##  $ amtv         : Factor w/ 8 levels "iap","agree strongly",..: 3 5 NA NA 6 NA NA NA NA NA ...
##  $ arthrtis     : Factor w/ 5 levels "iap","yes","no",..: 2 3 NA 3 NA NA NA NA NA NA ...
##  $ astrolgy     : Factor w/ 5 levels "iap","yes","no",..: NA NA NA 3 NA NA 3 3 2 3 ...
##  $ astrosci     : Factor w/ 6 levels "iap","very scientific",..: NA NA NA 4 NA NA 4 4 3 3 ...
##  $ attend       : Factor w/ 10 levels "never","lt once a year",..: 1 8 4 4 5 2 4 4 1 4 ...
##  $ attrally     : Factor w/ 7 levels "iap","have done it in the past yr",..: NA NA 4 5 NA 3 5 5 5 5 ...
##  $ avoidbuy     : Factor w/ 7 levels "iap","have done it in the past yr",..: NA NA 3 5 NA 2 4 2 NA 2 ...
##  $ away1        : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away10       : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away11       : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away12       : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away13       : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away14       : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away2        : Factor w/ 4 levels "iap","yes","dk",..: NA 2 NA NA NA NA NA NA NA NA ...
##  $ away3        : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away4        : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away5        : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away6        : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away7        : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away8        : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ away9        : Factor w/ 4 levels "iap","yes","dk",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ babies       : int  0 0 0 0 0 0 0 0 0 0 ...
##  $ backpain     : Factor w/ 5 levels "iap","yes","no",..: 3 3 NA 3 NA NA NA NA NA NA ...
##  $ ballot       : Factor w/ 5 levels "iap","ballot a",..: 4 4 2 3 4 2 3 3 4 4 ...
##  $ balneg       : Factor w/ 5 levels "iap","strongly in favor",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ balpos       : Factor w/ 5 levels "iap","strongly in favor",..: NA NA NA 2 NA NA 2 2 2 3 ...
##  $ barate       : Factor w/ 5 levels "iap","0-39.9%",..: NA NA NA NA NA NA NA NA NA NA ...
##  $ befair       : Factor w/ 7 levels "iap","advntge all time",..: NA NA 3 4 NA 4 2 3 2 3 ...
##  $ belikeus     : Factor w/ 8 levels "iap","strongly agree",..: 5 3 NA NA 2 NA NA NA NA NA ...
##  $ bettrlfe     : Factor w/ 7 levels "iap","strongly agree",..: NA NA NA 3 NA NA 2 3 4 3 ...
##  $ betrlang     : Factor w/ 5 levels "iap","language 1",..: NA NA NA 3 NA 3 NA 2 NA 2 ...
##  $ bible        : Factor w/ 7 levels "iap","word of god",..: 4 3 3 3 3 3 2 3 3 3 ...
##  $ bigbang      : Factor w/ 5 levels "iap","true","false",..: NA NA NA 2 NA NA 2 NA 3 2 ...
##  $ born         : Factor w/ 5 levels "iap","yes","no",..: 2 2 2 2 2 2 2 2 2 2 ...
##  $ boyorgrl     : Factor w/ 5 levels "iap","true","false",..: NA NA NA 3 NA NA 3 2 3 2 ...
##  $ buypol       : int  NA NA 5 5 NA 6 2 7 5 4 ...
##  $ buyvalue     : int  NA NA NA NA NA NA 9 NA NA NA ...
##  $ cantrust     : Factor w/ 7 levels "iap","always trusted",..: NA NA 4 3 NA 3 5 4 5 4 ...
##  $ cappun       : Factor w/ 5 levels "iap","favor",..: 2 3 2 3 3 2 2 2 2 3 ...
##  $ careself     : Factor w/ 8 levels "iap","strongly agree",..: 4 6 NA NA 5 NA NA NA NA NA ...
##  $ carried      : Factor w/ 9 levels "iap","more than once a week",..: 7 7 NA NA 6 NA NA NA NA NA ...
##  $ childs       : int  0 0 1 2 3 1 2 2 4 3 ...
##  $ chldidel     : int  NA NA 2 3 NA 2 8 8 NA NA ...
##  $ chngeoth     : Factor w/ 7 levels "iap","often",..: NA NA 3 5 NA 5 2 5 5 4 ...
##  $ chngtme      : Factor w/ 7 levels "iap","often",..: 2 3 NA 2 NA NA NA NA NA NA ...
##  $ choices      : Factor w/ 8 levels "iap","strongly agree",..: NA NA 3 4 NA 5 3 2 5 3 ...
##  $ citizen      : Factor w/ 5 levels "iap","yes","no",..: 2 2 NA NA 2 NA NA NA NA NA ...
##  $ citworld     : Factor w/ 8 levels "iap","agree strongly",..: 5 6 NA NA 6 NA NA NA NA NA ...
##  $ class        : Factor w/ 8 levels "iap","lower class",..: 4 4 3 4 4 4 NA 3 2 3 ...
##  $ closeblk     : int  5 4 5 NA 7 5 NA NA 9 9 ...
##  $ closewht     : int  5 9 1 NA 9 9 NA NA 9 9 ...
##  $ clsenoam     : Factor w/ 7 levels "iap","very close",..: 3 2 NA NA 3 NA NA NA NA NA ...
##  $ clsestat     : Factor w/ 7 levels "iap","very close",..: 5 4 NA NA 4 NA NA NA NA NA ...
##  $ clsetown     : Factor w/ 7 levels "iap","very close",..: 5 5 NA NA 3 NA NA NA NA NA ...
##  $ clseusa      : Factor w/ 7 levels "iap","very close",..: 3 2 NA NA 2 NA NA NA NA NA ...
##  $ cntctgov     : Factor w/ 7 levels "iap","have done it in the past yr",..: NA NA 4 5 NA 3 3 5 5 5 ...
##  $ cohort       : int  1961 1988 1955 1958 1940 1958 1951 1980 1977 1984 ...
##  $ colath       : Factor w/ 5 levels "iap","allowed",..: 2 2 2 NA 3 2 NA NA 2 2 ...
##  $ colcom       : Factor w/ 5 levels "iap","fired",..: 2 2 3 NA 3 2 NA NA NA 3 ...
##  $ coldeg1      : Factor w/ 11 levels "iap","associate's",..: NA NA NA 3 NA NA NA 4 NA 2 ...
##  $ colhomo      : Factor w/ 5 levels "iap","allowed",..: 2 2 3 NA 2 2 NA NA 3 2 ...
##  $ colmil       : Factor w/ 5 levels "iap","allowed",..: 2 3 2 NA 2 2 NA NA NA 2 ...
##  $ colmslm      : Factor w/ 5 levels "iap","yes, allowed",..: 2 3 2 NA 2 3 NA NA 3 NA ...
##  $ colrac       : Factor w/ 5 levels "iap","allowed",..: 2 3 2 NA 3 3 NA NA 2 2 ...
##   [list output truncated]
##  - attr(*, "datalabel")= chr ""
##  - attr(*, "time.stamp")= chr " 8 Mar 2015 17:06"
##  - attr(*, "formats")= chr  "%8.0g" "%8.0g" "%8.0g" "%8.0g" ...
##  - attr(*, "types")= int  251 251 251 251 251 251 251 251 251 251 ...
##  - attr(*, "val.labels")= chr  "LABA" "LABA" "LABA" "LABA" ...
##  - attr(*, "var.labels")= chr  "abortion if woman wants for any reason" "strong chance of serious defect" "womans health seriously endangered" "married--wants no more children" ...
##  - attr(*, "expansion.fields")=List of 4
##   ..$ : chr  "gradtounder" "note0" "1"
##   ..$ : chr  "gradtounder" "note1" "ratio of graduate students to undergraduate students in college respondent attended"
##   ..$ : chr  "cshutyp10" "note0" "1"
##   ..$ : chr  "cshutyp10" "note1" "during april through june <u>2010</u>, were you a resident in a household ??? such as living in a house, apartment or si"
##  - attr(*, "version")= int 12
##  - attr(*, "label.table")=List of 384
##   ..$ ZODIAC  : Named int  2147483630 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "iap" "aries" "taurus" "gemini" ...
##   ..$ YEARSJOB: Named int  2147483630 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "dont know" "no answer"
##   ..$ XNORCSIZ: Named int  2147483635 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "not assigned" "city gt 250000" "city,50-250000" "suburb, lrg city" ...
##   ..$ LABCJ   : Named int 2147483630
##   .. ..- attr(*, "names")= chr "iap"
##   ..$ WRKWAYUP: Named int  2147483630 1 2 3 4 5 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "agree strongly" "agree somewhat" "neither agree nor disagree" ...
##   ..$ WRKTYPE : Named int  2147483630 1 2 3 4 5 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "independent contractor/consultant/freelance worker" "on-call, work only when called to work" "paid by a temporary agency" ...
##   ..$ WRKSCHED: Named int  2147483630 1 2 3 4 5 6 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "day shift" "afternoon shift" "night shift" ...
##   ..$ WRKHOME : Named int  2147483630 1 2 3 4 5 6 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "never" "a few times a year" "about once a month" ...
##   ..$ WRKGOVT : Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "government" "private" "dk" ...
##   ..$ WORKFOR1: Named int  1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "for-profit company" "non-profit or not-for-profit organization" "government or government agency" "don't know" ...
##   ..$ LABCI   : Named int  2147483630 1 7 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "hardworking" "lazy" "dont know" ...
##   ..$ WORDSUM : Named int  2147483630 2147483625
##   .. ..- attr(*, "names")= chr  "iap" "did not try"
##   ..$ LABCH   : Named int  2147483630 0 1 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "incorrect" "correct" "na"
##   ..$ LABCG   : Named int  2147483630 1 7 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "rich" "poor" "dont know" ...
##   ..$ LABCF   : Named int  2147483630 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "no" "dk" ...
##   ..$ WKPRAISE: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "maybe" "no" ...
##   ..$ WKCONTCT: Named int  2147483630 1 2 3 4 5 6 7 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "never" "less than once a month" "once or twice a month" ...
##   ..$ WHYWKHME: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "worker wants to work at home" "worker has to work at home to keep up with job" "other combinations and other reasons" ...
##   ..$ LABCE   : Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "no" "dont know" ...
##   ..$ LABCD   : Named int  2147483630 1 2 3 4 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "another household" "traveling" "institution" ...
##   ..$ LABCC   : Named int  2147483630 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "dk" "na"
##   ..$ WEIGHT  : Named int  2147483630 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "don't know" "no answer"
##   ..$ WEBMOB  : Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "no" "don't know" ...
##   ..$ WEALTH  : Named int  2147483630 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "iap" "less than $5,000" "$5,000 to $20,000" "$20,000 to $40,000" ...
##   ..$ WAYPAID : Named int  2147483630 1 2 3 4 5 6 7 8 2147483625 ...
##   .. ..- attr(*, "names")= chr  "iap" "salaried" "paid by the hour" "other" ...
##   ..$ LABCB   : Named int 2147483630
##   .. ..- attr(*, "names")= chr "iap"
##   ..$ VOTE12  : Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "voted" "did not vote" "ineligible" ...
##   ..$ VOTE08  : Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "voted" "did not vote" "ineligible" ...
##   ..$ VOLMONTH: Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "no" "don't know" ...
##   ..$ VOEDNME2: Named int  2147483630 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "iap" "jail/prison" "high school" "company/business" ...
##   ..$ VOEDNME1: Named int  2147483630 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "iap" "jail/prison" "high school" "company/business" ...
##   ..$ LABCA   : Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "no" "dk" ...
##   ..$ VISITORS: Named int  2147483630 0 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "no visitors" "1 visitor" "2 visitors" ...
##   ..$ VETYEARS: Named int  2147483630 0 1 2 3 4 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "none" "less than 2 yrs" "2 to 4  years" ...
##   ..$ VALGIVEN: Named int  2147483630 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "dk" "na"
##   ..$ USETECH : Named int  2147483630 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "dk" "na"
##   ..$ USESKILL: Named int  2147483630 1 2 3 4 2147483624 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "almost none" "a little" "a lot" ...
##   ..$ USEDUP  : Named int  2147483630 1 2 3 4 5 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "very often" "often" "sometimes" ...
##   ..$ USCITZN : Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "a u.s. citizen" "not a u.s. citizen" "a u.s. citizen born in puerto rico, the u.s. virgin islands, or the northern mar" ...
##   ..$ UNRELAT : Named int  2147483630 8 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "8 or more" "na"
##   ..$ UNION   : Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "r belongs" "spouse belongs" "r and spouse belong" ...
##   ..$ TRYNEWJB: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "very likely" "somewhat likely" "not at all likely" ...
##   ..$ TRUST   : Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "can trust" "cannot trust" "depends" ...
##   ..$ TRDESTCK: Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "no" "don't know" ...
##   ..$ TAX     : Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "too high" "about right" "too low" ...
##   ..$ TALKSPVS: Named int  2147483630 1 2 3 4 5 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "not at all comfortable" "a little" "somewhat" ...
##   ..$ SUPERVIS: Named int  2147483635 1 2 2147483625 2147483622
##   .. ..- attr(*, "names")= chr  "no issp" "supervises" "doesnt supervise" "dk" ...
##   ..$ STRREDPG: Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "no" "don't know" ...
##   ..$ STRESS12: Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes" "no" "don't know" ...
##   ..$ STRESS  : Named int  2147483635 1 2 3 4 5 2147483624 2147483622
##   .. ..- attr(*, "names")= chr  "no issp" "always" "often" "sometimes" ...
##   ..$ LABBZ   : Named int  2147483630 9999995 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "other" "dont know" "no answer"
##   ..$ SRCBELT : Named int  2147483635 1 2 3 4 5 6
##   .. ..- attr(*, "names")= chr  "not assigned" "12 lrgst smsa's" "smsa's 13-100" "suburb, 12 lrgst" ...
##   ..$ LABBY   : Named int  2147483630 1 2 3 4 5 6 7 8 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "working fulltime" "working parttime" "temp not working" ...
##   ..$ SPSECTOR: Named int  2147483630 0 1 2 3 4 5 6 7 8 ...
##   .. ..- attr(*, "names")= chr  "iap" "administrative unit" "public, 4-year or above" "private not-for-profit, 4-year or above" ...
##   ..$ SPRTPRSN: Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "very spiritual" "modeate spirtual" "slight spiritual" ...
##   ..$ SPREL   : Named int  2147483630 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "iap" "protestant" "catholic" "jewish" ...
##   ..$ SPOTHER : Named int  2147483630 1 2 3 5 6 7 8 9 10 ...
##   .. ..- attr(*, "names")= chr  "iap" "hungarian reformed" "evangelical congregational" "ind bible, bible, bible fellowship" ...
##   ..$ SPKMSLM : Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "yes, allowed" "not allowed" "dont know" ...
##   ..$ SPKLANG : Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "very well" "well" "not well" ...
##   ..$ LABBX   : Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "allowed" "not allowed" "dk" ...
##   ..$ SPGRADTO: Named int  2147483630 1 2 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "0-.210499" ".211+" "no answer"
##   ..$ SPEFTOTL: Named int  2147483630 1 2 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "0-467" "468+" "no answer"
##   ..$ SPDIPGED: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "high school diploma" "ged" "other" ...
##   ..$ SPBARATE: Named int  2147483630 1 2 3 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "0-39.9%" "40-57.9%" "58%+" ...
##   ..$ SPANSELF: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "have been able to do the interview easily in english" "have been able to do the interview with difficulty" "not have been able to do the interview" ...
##   ..$ SPANINT : Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "could have been interviewed in english" "would have been excluded as language problem" "dont know" ...
##   ..$ SPANENG : Named int  2147483630 1 2
##   .. ..- attr(*, "names")= chr  "iap" "english" "spanish"
##   ..$ SOLARREV: Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "one day" "one month" "one year" ...
##   ..$ LABBW   : Named int  2147483630 1 2 3 4 5 6 7 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "almost daily" "sev times a week" "sev times a mnth" ...
##   ..$ SLPPRBLM: Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "often" "sometimes" "rarely" ...
##   ..$ SIZE    : Named int 2147483635
##   .. ..- attr(*, "names")= chr "not assigned"
##   ..$ SHORTCOM: Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "agree strongly" "agree" "neither agree nor disagree" ...
##   ..$ LABBV   : Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "exclusively male" "both male and female" "exclusively female" ...
##   ..$ SEXORNT : Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "gay, lesbian, or homosexual" "bisexual" "heterosexual or straight" ...
##   ..$ SEXFREQ : Named int  2147483630 0 1 2 3 4 5 6 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "not at all" "once or twice" "once a month" ...
##   ..$ SEXEDUC : Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "favor" "oppose" "depends" ...
##   ..$ SEX     : Named int  1 2
##   .. ..- attr(*, "names")= chr  "male" "female"
##   ..$ SERVEPEO: Named int  2147483630 1 2 3 4 2147483624 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "very committed" "somewhat committed" "not very committed" ...
##   ..$ SEEKSCI : Named int  2147483630 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "iap" "newspapers" "magazines" "the internet" ...
##   ..$ SECTOR  : Named int  2147483630 0 1 2 3 4 5 6 7 8 ...
##   .. ..- attr(*, "names")= chr  "iap" "administrative unit" "public, 4-year or above" "private not-for-profit, 4-year or above" ...
##   ..$ SCITEXT : Named int  2147483630 1 2 3 4 5 6 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "formulation of theories test hypotheses" "do experiments control group" "rigorous systematic comparison" ...
##   ..$ SCISTUDY: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "clear understanding" "general sense" "little understanding" ...
##   ..$ SCINEWS3: Named int  2147483630 1 2 3 4 5 6 7 8 10 ...
##   .. ..- attr(*, "names")= chr  "iap" "online newspapers" "online magazines" "science site" ...
##   ..$ SCINEWS2: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "printed magazines" "online magazines" "other specify" ...
##   ..$ SCINEWS1: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "printed newspapers" "online newspapers" "other specify" ...
##   ..$ SCIBNFTS: Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "benefits greater" "about equal if volunteered" "harmful results greater" ...
##   ..$ SATJOB1 : Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "very satisfied" "somewhat satisfied" "not too satisfied" ...
##   ..$ SATJOB  : Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "very satisfied" "mod. satisfied" "a little dissat" ...
##   ..$ SATFIN  : Named int  2147483630 1 2 3 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "satisfied" "more or less" "not at all sat" ...
##   ..$ SAMPCODE: Named int 2147483630
##   .. ..- attr(*, "names")= chr "iap"
##   ..$ RVISITOR: Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "r. is household member" "r. is visitor" "dk" ...
##   ..$ RICHWORK: Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "continue working" "stop working" "dk" ...
##   ..$ RESPNUM : Named int  2147483630 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "iap" "1st person" "2nd person" "3rd person" ...
##   ..$ RES2012 : Named int  1 2 8 9
##   .. ..- attr(*, "names")= chr  "living in the u.s." "not living in u.s." "don't know" "refused"
##   ..$ RES2010 : Named int  2147483630 1 2 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "living in the u.s." "not living in u.s." "dont know" ...
##   ..$ RES16   : Named int  2147483630 1 2 3 4 5 6 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "country,nonfarm" "farm" "town lt 50000" ...
##   ..$ LABBU   : Named int  2147483630 1 2 3 4 5 6 7 8 9 ...
##   .. ..- attr(*, "names")= chr  "iap" "hh spouse, partner" "spouse" "partner, fiance, boyfriend, girlfriend, etc" ...
##   ..$ RELPERSN: Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "very religious" "modrte religious" "slight religious" ...
##   ..$ RELITEN : Named int  2147483630 1 2 3 4 2147483625 2147483635
##   .. ..- attr(*, "names")= chr  "iap" "strong" "not very strong" "somewhat strong" ...
##   .. [list output truncated]
GSS2014.mod<- lm(colrac ~ colhomo, data= GSS2014)
GSS2014.z <- zelig(colrac ~ colhomo, data = GSS2014, model ="normal")
## <simpleError in glm.fit(x = structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 0, 0, 1, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 1, 1, 0, 1, 1, 1, 0, 1, 1, 1, 0, 1, 0, 0, 0, 1, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 1, 0, 0, 1, 1, 0, 1, 0, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), .Dim = c(1660L, 2L), .Dimnames = list(    c("1", "2", "3", "5", "6", "9", "10", "11", "12", "13", "14",     "18", "20", "21", "23", "24", "25", "26", "27", "29", "30",     "31", "32", "33", "34", "36", "37", "39", "40", "41", "42",     "43", "45", "46", "47", "48", "49", "50", "52", "53", "55",     "56", "57", "58", "60", "61", "63", "64", "65", "66", "67",     "69", "70", "72", "75", "76", "77", "78", "79", "80", "82",     "83", "84", "85", "88", "89", "90", "92", "93", "96", "97",     "98", "99", "103", "105", "111", "113", "114", "115", "118",     "119", "120", "121", "122", "123", "125", "126", "127", "128",     "129", "130", "131", "133", "134", "136", "137", "138", "139",     "140", "142", "144", "146", "147", "148", "151", "152", "153",     "154", "157", "158", "159", "161", "163", "165", "167", "171",     "172", "173", "174", "175", "176", "177", "178", "179", "180",     "183", "184", "185", "186", "189", "190", "191", "193", "194",     "196", "198", "200", "201", "202", "204", "205", "207", "208",     "211", "212", "213", "216", "218", "220", "221", "223", "227",     "228", "229", "231", "233", "234", "235", "236", "238", "239",     "240", "241", "245", "247", "248", "249", "250", "252", "253",     "255", "262", "263", "264", "265", "271", "272", "273", "274",     "275", "276", "277", "283", "285", "287", "289", "290", "291",     "292", "293", "294", "295", "297", "298", "300", "301", "303",     "306", "307", "308", "309", "310", "311", "312", "315", "316",     "317", "318", "319", "323", "326", "327", "328", "329", "330",     "331", "332", "337", "338", "339", "340", "341", "342", "344",     "345", "347", "350", "351", "353", "355", "356", "359", "360",     "361", "362", "363", "364", "366", "368", "369", "370", "372",     "373", "374", "375", "377", "380", "381", "382", "383", "386",     "387", "389", "392", "393", "397", "399", "400", "401", "402",     "403", "405", "407", "411", "412", "413", "417", "418", "419",     "421", "424", "425", "426", "427", "429", "430", "431", "432",     "433", "434", "435", "436", "437", "438", "439", "440", "442",     "443", "444", "446", "447", "448", "451", "453", "454", "455",     "456", "457", "458", "459", "460", "461", "462", "464", "465",     "466", "468", "470", "471", "474", "475", "476", "477", "478",     "480", "481", "482", "483", "484", "485", "488", "489", "490",     "491", "492", "494", "495", "497", "498", "499", "500", "501",     "502", "503", "504", "507", "508", "510", "511", "512", "513",     "514", "516", "517", "518", "519", "522", "523", "524", "526",     "528", "530", "531", "534", "535", "536", "537", "538", "539",     "541", "542", "544", "547", "548", "549", "550", "552", "553",     "554", "555", "560", "562", "564", "567", "568", "569", "570",     "571", "573", "574", "576", "581", "582", "584", "585", "586",     "587", "588", "589", "590", "592", "593", "596", "597", "598",     "599", "603", "604", "606", "611", "612", "615", "616", "618",     "619", "621", "622", "623", "624", "626", "630", "631", "632",     "635", "636", "637", "639", "640", "643", "645", "646", "647",     "648", "649", "650", "651", "654", "655", "656", "658", "660",     "661", "662", "664", "666", "667", "669", "670", "672", "673",     "674", "676", "677", "681", "684", "685", "687", "689", "690",     "691", "694", "695", "696", "698", "701", "702", "704", "705",     "707", "709", "710", "711", "713", "714", "718", "719", "721",     "722", "723", "724", "725", "727", "728", "729", "731", "733",     "735", "736", "737", "738", "741", "742", "743", "745", "746",     "748", "749", "750", "752", "753", "754", "755", "757", "758",     "759", "760", "761", "762", "763", "764", "765", "766", "767",     "768", "769", "770", "772", "776", "777", "780", "781", "782",     "786", "788", "789", "790", "792", "793", "794", "796", "797",     "798", "799", "800", "801", "803", "804", "807", "808", "810",     "812", "813", "814", "816", "818", "819", "820", "822", "823",     "826", "827", "829", "830", "831", "832", "833", "834", "836",     "837", "839", "840", "841", "843", "844", "845", "846", "847",     "848", "850", "851", "852", "854", "856", "858", "860", "862",     "863", "870", "871", "872", "874", "876", "878", "880", "881",     "882", "885", "886", "888", "889", "891", "892", "893", "896",     "899", "902", "905", "906", "907", "908", "909", "912", "913",     "914", "915", "916", "921", "925", "926", "927", "931", "933",     "934", "935", "936", "938", "939", "940", "941", "943", "950",     "951", "953", "954", "958", "960", "961", "962", "964", "965",     "966", "967", "968", "969", "970", "971", "973", "976", "977",     "979", "980", "981", "982", "983", "984", "985", "986", "988",     "989", "990", "992", "993", "994", "996", "997", "998", "1001",     "1002", "1003", "1005", "1006", "1007", "1008", "1009", "1011",     "1012", "1014", "1018", "1019", "1022", "1025", "1026", "1029",     "1031", "1032", "1034", "1036", "1038", "1039", "1040", "1041",     "1042", "1046", "1047", "1048", "1049", "1051", "1052", "1054",     "1055", "1056", "1057", "1060", "1063", "1064", "1066", "1067",     "1068", "1069", "1071", "1072", "1073", "1075", "1076", "1077",     "1079", "1082", "1083", "1085", "1087", "1088", "1089", "1090",     "1091", "1092", "1093", "1095", "1098", "1099", "1101", "1102",     "1103", "1104", "1105", "1107", "1108", "1111", "1112", "1113",     "1114", "1115", "1116", "1118", "1119", "1121", "1124", "1125",     "1127", "1129", "1131", "1133", "1136", "1137", "1138", "1139",     "1140", "1141", "1142", "1144", "1145", "1146", "1149", "1150",     "1151", "1153", "1154", "1155", "1156", "1157", "1158", "1160",     "1161", "1165", "1167", "1168", "1169", "1171", "1172", "1173",     "1178", "1179", "1180", "1181", "1182", "1183", "1184", "1185",     "1187", "1188", "1189", "1190", "1191", "1193", "1194", "1196",     "1197", "1199", "1200", "1203", "1204", "1205", "1206", "1209",     "1210", "1212", "1214", "1219", "1220", "1221", "1224", "1225",     "1226", "1227", "1228", "1229", "1230", "1231", "1233", "1234",     "1235", "1236", "1237", "1238", "1239", "1243", "1244", "1245",     "1247", "1248", "1249", "1251", "1252", "1253", "1255", "1256",     "1258", "1259", "1260", "1261", "1265", "1266", "1267", "1271",     "1272", "1273", "1275", "1276", "1277", "1278", "1279", "1280",     "1281", "1282", "1283", "1284", "1286", "1287", "1288", "1289",     "1290", "1291", "1292", "1294", "1297", "1298", "1300", "1301",     "1304", "1307", "1308", "1310", "1312", "1313", "1315", "1316",     "1317", "1318", "1319", "1321", "1322", "1323", "1327", "1329",     "1330", "1332", "1334", "1335", "1337", "1338", "1340", "1342",     "1346", "1348", "1349", "1350", "1352", "1354", "1357", "1358",     "1362", "1364", "1365", "1367", "1368", "1369", "1370", "1371",     "1372", "1377", "1378", "1379", "1380", "1381", "1382", "1383",     "1385", "1387", "1390", "1391", "1392", "1393", "1394", "1395",     "1396", "1397", "1398", "1399", "1400", "1402", "1404", "1405",     "1406", "1412", "1414", "1415", "1416", "1417", "1418", "1420",     "1421", "1422", "1423", "1425", "1426", "1427", "1428", "1429",     "1431", "1432", "1433", "1434", "1435", "1436", "1438", "1439",     "1441", "1443", "1445", "1450", "1451", "1452", "1453", "1456",     "1458", "1461", "1465", "1466", "1469", "1470", "1471", "1472",     "1473", "1474", "1476", "1477", "1478", "1479", "1480", "1481",     "1482", "1483", "1484", "1485", "1489", "1491", "1492", "1493",     "1496", "1502", "1504", "1505", "1507", "1508", "1509", "1511",     "1512", "1514", "1515", "1516", "1517", "1519", "1520", "1524",     "1530", "1532", "1533", "1534", "1536", "1538", "1540", "1541",     "1542", "1543", "1544", "1547", "1551", "1552", "1553", "1554",     "1555", "1557", "1560", "1561", "1562", "1563", "1564", "1565",     "1567", "1568", "1569", "1572", "1573", "1574", "1575", "1577",     "1579", "1580", "1581", "1582", "1584", "1586", "1587", "1589",     "1592", "1594", "1595", "1598", "1600", "1601", "1605", "1606",     "1607", "1608", "1610", "1611", "1612", "1613", "1616", "1618",     "1619", "1620", "1623", "1624", "1625", "1627", "1628", "1629",     "1630", "1631", "1632", "1633", "1634", "1636", "1637", "1638",     "1641", "1642", "1643", "1644", "1645", "1646", "1649", "1650",     "1651", "1652", "1653", "1654", "1659", "1660", "1662", "1664",     "1665", "1668", "1670", "1671", "1672", "1673", "1674", "1677",     "1678", "1679", "1681", "1683", "1688", "1689", "1691", "1693",     "1694", "1695", "1696", "1699", "1700", "1702", "1704", "1705",     "1706", "1707", "1708", "1709", "1710", "1711", "1712", "1713",     "1715", "1717", "1719", "1721", "1722", "1723", "1724", "1725",     "1728", "1729", "1732", "1733", "1734", "1735", "1738", "1740",     "1741", "1742", "1743", "1744", "1748", "1749", "1750", "1752",     "1754", "1755", "1756", "1757", "1760", "1761", "1763", "1764",     "1765", "1767", "1768", "1769", "1770", "1771", "1773", "1774",     "1778", "1780", "1781", "1782", "1783", "1784", "1787", "1788",     "1789", "1790", "1793", "1794", "1795", "1796", "1800", "1801",     "1802", "1803", "1804", "1805", "1807", "1808", "1810", "1811",     "1812", "1814", "1815", "1816", "1817", "1818", "1819", "1820",     "1821", "1822", "1823", "1824", "1825", "1827", "1828", "1829",     "1830", "1832", "1834", "1839", "1840", "1841", "1842", "1843",     "1844", "1845", "1846", "1847", "1848", "1849", "1850", "1851",     "1853", "1854", "1855", "1856", "1857", "1858", "1859", "1860",     "1861", "1862", "1863", "1865", "1866", "1867", "1868", "1869",     "1870", "1873", "1875", "1876", "1877", "1878", "1880", "1881",     "1882", "1883", "1884", "1886", "1887", "1890", "1892", "1893",     "1895", "1898", "1899", "1900", "1901", "1902", "1903", "1907",     "1910", "1914", "1915", "1916", "1917", "1920", "1921", "1922",     "1923", "1924", "1925", "1926", "1927", "1928", "1929", "1930",     "1932", "1933", "1934", "1936", "1938", "1939", "1940", "1941",     "1945", "1946", "1949", "1951", "1952", "1953", "1954", "1955",     "1956", "1957", "1958", "1959", "1961", "1962", "1963", "1964",     "1966", "1967", "1968", "1969", "1970", "1971", "1972", "1973",     "1976", "1978", "1979", "1980", "1981", "1982", "1983", "1984",     "1985", "1986", "1987", "1988", "1989", "1991", "1994", "1996",     "1997", "1998", "1999", "2000", "2003", "2004", "2005", "2006",     "2008", "2009", "2010", "2011", "2013", "2014", "2015", "2016",     "2017", "2019", "2020", "2021", "2022", "2023", "2026", "2028",     "2031", "2033", "2036", "2038", "2039", "2040", "2042", "2044",     "2045", "2046", "2047", "2050", "2051", "2052", "2053", "2054",     "2057", "2058", "2059", "2061", "2064", "2065", "2067", "2069",     "2071", "2073", "2074", "2076", "2078", "2079", "2080", "2081",     "2082", "2083", "2084", "2086", "2088", "2089", "2091", "2092",     "2093", "2094", "2095", "2096", "2098", "2099", "2102", "2103",     "2104", "2105", "2106", "2107", "2108", "2111", "2112", "2113",     "2115", "2116", "2117", "2119", "2120", "2121", "2122", "2124",     "2125", "2126", "2132", "2133", "2134", "2136", "2137", "2139",     "2140", "2141", "2142", "2143", "2147", "2148", "2149", "2151",     "2152", "2155", "2157", "2158", "2159", "2160", "2161", "2162",     "2163", "2164", "2166", "2167", "2168", "2169", "2171", "2173",     "2174", "2175", "2176", "2177", "2178", "2179", "2180", "2181",     "2182", "2184", "2185", "2186", "2187", "2188", "2189", "2190",     "2193", "2195", "2196", "2198", "2199", "2200", "2201", "2202",     "2203", "2205", "2206", "2208", "2211", "2213", "2214", "2216",     "2218", "2219", "2222", "2223", "2224", "2225", "2227", "2228",     "2229", "2230", "2233", "2234", "2235", "2236", "2238", "2239",     "2240", "2241", "2243", "2245", "2246", "2248", "2249", "2250",     "2251", "2252", "2253", "2254", "2255", "2258", "2259", "2261",     "2262", "2264", "2267", "2268", "2270", "2271", "2272", "2275",     "2276", "2278", "2279", "2281", "2282", "2285", "2287", "2288",     "2290", "2291", "2298", "2299", "2300", "2301", "2303", "2304",     "2305", "2307", "2308", "2312", "2313", "2314", "2316", "2317",     "2318", "2319", "2320", "2321", "2322", "2323", "2325", "2327",     "2329", "2330", "2331", "2332", "2334", "2335", "2336", "2337",     "2338", "2339", "2340", "2341", "2342", "2343", "2344", "2347",     "2348", "2349", "2350", "2351", "2352", "2353", "2355", "2356",     "2358", "2359", "2360", "2361", "2365", "2366", "2368", "2369",     "2371", "2372", "2375", "2376", "2377", "2378", "2379", "2380",     "2381", "2382", "2385", "2386", "2387", "2389", "2392", "2393",     "2394", "2395", "2396", "2400", "2401", "2402", "2404", "2406",     "2408", "2409", "2410", "2411", "2412", "2413", "2414", "2415",     "2418", "2421", "2422", "2424", "2425", "2427", "2429", "2430",     "2431", "2433", "2434", "2435", "2436", "2437", "2439", "2440",     "2443", "2444", "2447", "2450", "2451", "2453", "2454", "2455",     "2457", "2458", "2459", "2460", "2462", "2464", "2465", "2466",     "2470", "2471", "2472", "2474", "2476", "2478", "2479", "2481",     "2482", "2487", "2490", "2491", "2492", "2495", "2499", "2501",     "2502", "2504", "2505", "2507", "2508", "2510", "2513", "2514",     "2516", "2517", "2518", "2520", "2521", "2522", "2524", "2529",     "2530", "2535", "2536", "2537", "2538"), c("(Intercept)",     "colhomonot allowed")), assign = 0:1, contrasts = structure(list(    colhomo = "contr.treatment"), .Names = "colhomo")), y = structure(c(1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 2L, 1L, 1L, 2L, 2L, 2L, 1L, 1L, 2L, 1L, 2L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L), .Label = c("allowed", "not allowed"), class = "factor", .Names = c("1", "2", "3", "5", "6", "9", "10", "11", "12", "13", "14", "18", "20", "21", "23", "24", "25", "26", "27", "29", "30", "31", "32", "33", "34", "36", "37", "39", "40", "41", "42", "43", "45", "46", "47", "48", "49", "50", "52", "53", "55", "56", "57", "58", "60", "61", "63", "64", "65", "66", "67", "69", "70", "72", "75", "76", "77", "78", "79", "80", "82", "83", "84", "85", "88", "89", "90", "92", "93", "96", "97", "98", "99", "103", "105", "111", "113", "114", "115", "118", "119", "120", "121", "122", "123", "125", "126", "127", "128", "129", "130", "131", "133", "134", "136", "137", "138", "139", "140", "142", "144", "146", "147", "148", "151", "152", "153", "154", "157", "158", "159", "161", "163", "165", "167", "171", "172", "173", "174", "175", "176", "177", "178", "179", "180", "183", "184", "185", "186", "189", "190", "191", "193", "194", "196", "198", "200", "201", "202", "204", "205", "207", "208", "211", "212", "213", "216", "218", "220", "221", "223", "227", "228", "229", "231", "233", "234", "235", "236", "238", "239", "240", "241", "245", "247", "248", "249", "250", "252", "253", "255", "262", "263", "264", "265", "271", "272", "273", "274", "275", "276", "277", "283", "285", "287", "289", "290", "291", "292", "293", "294", "295", "297", "298", "300", "301", "303", "306", "307", "308", "309", "310", "311", "312", "315", "316", "317", "318", "319", "323", "326", "327", "328", "329", "330", "331", "332", "337", "338", "339", "340", "341", "342", "344", "345", "347", "350", "351", "353", "355", "356", "359", "360", "361", "362", "363", "364", "366", "368", "369", "370", "372", "373", "374", "375", "377", "380", "381", "382", "383", "386", "387", "389", "392", "393", "397", "399", "400", "401", "402", "403", "405", "407", "411", "412", "413", "417", "418", "419", "421", "424", "425", "426", "427", "429", "430", "431", "432", "433", "434", "435", "436", "437", "438", "439", "440", "442", "443", "444", "446", "447", "448", "451", "453", "454", "455", "456", "457", "458", "459", "460", "461", "462", "464", "465", "466", "468", "470", "471", "474", "475", "476", "477", "478", "480", "481", "482", "483", "484", "485", "488", "489", "490", "491", "492", "494", "495", "497", "498", "499", "500", "501", "502", "503", "504", "507", "508", "510", "511", "512", "513", "514", "516", "517", "518", "519", "522", "523", "524", "526", "528", "530", "531", "534", "535", "536", "537", "538", "539", "541", "542", "544", "547", "548", "549", "550", "552", "553", "554", "555", "560", "562", "564", "567", "568", "569", "570", "571", "573", "574", "576", "581", "582", "584", "585", "586", "587", "588", "589", "590", "592", "593", "596", "597", "598", "599", "603", "604", "606", "611", "612", "615", "616", "618", "619", "621", "622", "623", "624", "626", "630", "631", "632", "635", "636", "637", "639", "640", "643", "645", "646", "647", "648", "649", "650", "651", "654", "655", "656", "658", "660", "661", "662", "664", "666", "667", "669", "670", "672", "673", "674", "676", "677", "681", "684", "685", "687", "689", "690", "691", "694", "695", "696", "698", "701", "702", "704", "705", "707", "709", "710", "711", "713", "714", "718", "719", "721", "722", "723", "724", "725", "727", "728", "729", "731", "733", "735", "736", "737", "738", "741", "742", "743", "745", "746", "748", "749", "750", "752", "753", "754", "755", "757", "758", "759", "760", "761", "762", "763", "764", "765", "766", "767", "768", "769", "770", "772", "776", "777", "780", "781", "782", "786", "788", "789", "790", "792", "793", "794", "796", "797", "798", "799", "800", "801", "803", "804", "807", "808", "810", "812", "813", "814", "816", "818", "819", "820", "822", "823", "826", "827", "829", "830", "831", "832", "833", "834", "836", "837", "839", "840", "841", "843", "844", "845", "846", "847", "848", "850", "851", "852", "854", "856", "858", "860", "862", "863", "870", "871", "872", "874", "876", "878", "880", "881", "882", "885", "886", "888", "889", "891", "892", "893", "896", "899", "902", "905", "906", "907", "908", "909", "912", "913", "914", "915", "916", "921", "925", "926", "927", "931", "933", "934", "935", "936", "938", "939", "940", "941", "943", "950", "951", "953", "954", "958", "960", "961", "962", "964", "965", "966", "967", "968", "969", "970", "971", "973", "976", "977", "979", "980", "981", "982", "983", "984", "985", "986", "988", "989", "990", "992", "993", "994", "996", "997", "998", "1001", "1002", "1003", "1005", "1006", "1007", "1008", "1009", "1011", "1012", "1014", "1018", "1019", "1022", "1025", "1026", "1029", "1031", "1032", "1034", "1036", "1038", "1039", "1040", "1041", "1042", "1046", "1047", "1048", "1049", "1051", "1052", "1054", "1055", "1056", "1057", "1060", "1063", "1064", "1066", "1067", "1068", "1069", "1071", "1072", "1073", "1075", "1076", "1077", "1079", "1082", "1083", "1085", "1087", "1088", "1089", "1090", "1091", "1092", "1093", "1095", "1098", "1099", "1101", "1102", "1103", "1104", "1105", "1107", "1108", "1111", "1112", "1113", "1114", "1115", "1116", "1118", "1119", "1121", "1124", "1125", "1127", "1129", "1131", "1133", "1136", "1137", "1138", "1139", "1140", "1141", "1142", "1144", "1145", "1146", "1149", "1150", "1151", "1153", "1154", "1155", "1156", "1157", "1158", "1160", "1161", "1165", "1167", "1168", "1169", "1171", "1172", "1173", "1178", "1179", "1180", "1181", "1182", "1183", "1184", "1185", "1187", "1188", "1189", "1190", "1191", "1193", "1194", "1196", "1197", "1199", "1200", "1203", "1204", "1205", "1206", "1209", "1210", "1212", "1214", "1219", "1220", "1221", "1224", "1225", "1226", "1227", "1228", "1229", "1230", "1231", "1233", "1234", "1235", "1236", "1237", "1238", "1239", "1243", "1244", "1245", "1247", "1248", "1249", "1251", "1252", "1253", "1255", "1256", "1258", "1259", "1260", "1261", "1265", "1266", "1267", "1271", "1272", "1273", "1275", "1276", "1277", "1278", "1279", "1280", "1281", "1282", "1283", "1284", "1286", "1287", "1288", "1289", "1290", "1291", "1292", "1294", "1297", "1298", "1300", "1301", "1304", "1307", "1308", "1310", "1312", "1313", "1315", "1316", "1317", "1318", "1319", "1321", "1322", "1323", "1327", "1329", "1330", "1332", "1334", "1335", "1337", "1338", "1340", "1342", "1346", "1348", "1349", "1350", "1352", "1354", "1357", "1358", "1362", "1364", "1365", "1367", "1368", "1369", "1370", "1371", "1372", "1377", "1378", "1379", "1380", "1381", "1382", "1383", "1385", "1387", "1390", "1391", "1392", "1393", "1394", "1395", "1396", "1397", "1398", "1399", "1400", "1402", "1404", "1405", "1406", "1412", "1414", "1415", "1416", "1417", "1418", "1420", "1421", "1422", "1423", "1425", "1426", "1427", "1428", "1429", "1431", "1432", "1433", "1434", "1435", "1436", "1438", "1439", "1441", "1443", "1445", "1450", "1451", "1452", "1453", "1456", "1458", "1461", "1465", "1466", "1469", "1470", "1471", "1472", "1473", "1474", "1476", "1477", "1478", "1479", "1480", "1481", "1482", "1483", "1484", "1485", "1489", "1491", "1492", "1493", "1496", "1502", "1504", "1505", "1507", "1508", "1509", "1511", "1512", "1514", "1515", "1516", "1517", "1519", "1520", "1524", "1530", "1532", "1533", "1534", "1536", "1538", "1540", "1541", "1542", "1543", "1544", "1547", "1551", "1552", "1553", "1554", "1555", "1557", "1560", "1561", "1562", "1563", "1564", "1565", "1567", "1568", "1569", "1572", "1573", "1574", "1575", "1577", "1579", "1580", "1581", "1582", "1584", "1586", "1587", "1589", "1592", "1594", "1595", "1598", "1600", "1601", "1605", "1606", "1607", "1608", "1610", "1611", "1612", "1613", "1616", "1618", "1619", "1620", "1623", "1624", "1625", "1627", "1628", "1629", "1630", "1631", "1632", "1633", "1634", "1636", "1637", "1638", "1641", "1642", "1643", "1644", "1645", "1646", "1649", "1650", "1651", "1652", "1653", "1654", "1659", "1660", "1662", "1664", "1665", "1668", "1670", "1671", "1672", "1673", "1674", "1677", "1678", "1679", "1681", "1683", "1688", "1689", "1691", "1693", "1694", "1695", "1696", "1699", "1700", "1702", "1704", "1705", "1706", "1707", "1708", "1709", "1710", "1711", "1712", "1713", "1715", "1717", "1719", "1721", "1722", "1723", "1724", "1725", "1728", "1729", "1732", "1733", "1734", "1735", "1738", "1740", "1741", "1742", "1743", "1744", "1748", "1749", "1750", "1752", "1754", "1755", "1756", "1757", "1760", "1761", "1763", "1764", "1765", "1767", "1768", "1769", "1770", "1771", "1773", "1774", "1778", "1780", "1781", "1782", "1783", "1784", "1787", "1788", "1789", "1790", "1793", "1794", "1795", "1796", "1800", "1801", "1802", "1803", "1804", "1805", "1807", "1808", "1810", "1811", "1812", "1814", "1815", "1816", "1817", "1818", "1819", "1820", "1821", "1822", "1823", "1824", "1825", "1827", "1828", "1829", "1830", "1832", "1834", "1839", "1840", "1841", "1842", "1843", "1844", "1845", "1846", "1847", "1848", "1849", "1850", "1851", "1853", "1854", "1855", "1856", "1857", "1858", "1859", "1860", "1861", "1862", "1863", "1865", "1866", "1867", "1868", "1869", "1870", "1873", "1875", "1876", "1877", "1878", "1880", "1881", "1882", "1883", "1884", "1886", "1887", "1890", "1892", "1893", "1895", "1898", "1899", "1900", "1901", "1902", "1903", "1907", "1910", "1914", "1915", "1916", "1917", "1920", "1921", "1922", "1923", "1924", "1925", "1926", "1927", "1928", "1929", "1930", "1932", "1933", "1934", "1936", "1938", "1939", "1940", "1941", "1945", "1946", "1949", "1951", "1952", "1953", "1954", "1955", "1956", "1957", "1958", "1959", "1961", "1962", "1963", "1964", "1966", "1967", "1968", "1969", "1970", "1971", "1972", "1973", "1976", "1978", "1979", "1980", "1981", "1982", "1983", "1984", "1985", "1986", "1987", "1988", "1989", "1991", "1994", "1996", "1997", "1998", "1999", "2000", "2003", "2004", "2005", "2006", "2008", "2009", "2010", "2011", "2013", "2014", "2015", "2016", "2017", "2019", "2020", "2021", "2022", "2023", "2026", "2028", "2031", "2033", "2036", "2038", "2039", "2040", "2042", "2044", "2045", "2046", "2047", "2050", "2051", "2052", "2053", "2054", "2057", "2058", "2059", "2061", "2064", "2065", "2067", "2069", "2071", "2073", "2074", "2076", "2078", "2079", "2080", "2081", "2082", "2083", "2084", "2086", "2088", "2089", "2091", "2092", "2093", "2094", "2095", "2096", "2098", "2099", "2102", "2103", "2104", "2105", "2106", "2107", "2108", "2111", "2112", "2113", "2115", "2116", "2117", "2119", "2120", "2121", "2122", "2124", "2125", "2126", "2132", "2133", "2134", "2136", "2137", "2139", "2140", "2141", "2142", "2143", "2147", "2148", "2149", "2151", "2152", "2155", "2157", "2158", "2159", "2160", "2161", "2162", "2163", "2164", "2166", "2167", "2168", "2169", "2171", "2173", "2174", "2175", "2176", "2177", "2178", "2179", "2180", "2181", "2182", "2184", "2185", "2186", "2187", "2188", "2189", "2190", "2193", "2195", "2196", "2198", "2199", "2200", "2201", "2202", "2203", "2205", "2206", "2208", "2211", "2213", "2214", "2216", "2218", "2219", "2222", "2223", "2224", "2225", "2227", "2228", "2229", "2230", "2233", "2234", "2235", "2236", "2238", "2239", "2240", "2241", "2243", "2245", "2246", "2248", "2249", "2250", "2251", "2252", "2253", "2254", "2255", "2258", "2259", "2261", "2262", "2264", "2267", "2268", "2270", "2271", "2272", "2275", "2276", "2278", "2279", "2281", "2282", "2285", "2287", "2288", "2290", "2291", "2298", "2299", "2300", "2301", "2303", "2304", "2305", "2307", "2308", "2312", "2313", "2314", "2316", "2317", "2318", "2319", "2320", "2321", "2322", "2323", "2325", "2327", "2329", "2330", "2331", "2332", "2334", "2335", "2336", "2337", "2338", "2339", "2340", "2341", "2342", "2343", "2344", "2347", "2348", "2349", "2350", "2351", "2352", "2353", "2355", "2356", "2358", "2359", "2360", "2361", "2365", "2366", "2368", "2369", "2371", "2372", "2375", "2376", "2377", "2378", "2379", "2380", "2381", "2382", "2385", "2386", "2387", "2389", "2392", "2393", "2394", "2395", "2396", "2400", "2401", "2402", "2404", "2406", "2408", "2409", "2410", "2411", "2412", "2413", "2414", "2415", "2418", "2421", "2422", "2424", "2425", "2427", "2429", "2430", "2431", "2433", "2434", "2435", "2436", "2437", "2439", "2440", "2443", "2444", "2447", "2450", "2451", "2453", "2454", "2455", "2457", "2458", "2459", "2460", "2462", "2464", "2465", "2466", "2470", "2471", "2472", "2474", "2476", "2478", "2479", "2481", "2482", "2487", "2490", "2491", "2492", "2495", "2499", "2501", "2502", "2504", "2505", "2507", "2508", "2510", "2513", "2514", "2516", "2517", "2518", "2520", "2521", "2522", "2524", "2529", "2530", "2535", "2536", "2537", "2538")), weights = NULL, start = NULL,     etastart = NULL, mustart = NULL, offset = NULL, family = structure(list(        family = "gaussian", link = "identity", linkfun = function (mu)         mu, linkinv = function (eta)         eta, variance = function (mu)         rep.int(1, length(mu)), dev.resids = function (y, mu,             wt)         wt * ((y - mu)^2), aic = function (y, n, mu, wt, dev)         {            nobs <- length(y)            nobs * (log(dev/nobs * 2 * pi) + 1) + 2 - sum(log(wt))        }, mu.eta = function (eta)         rep.int(1, length(eta)), initialize = expression({            n <- rep.int(1, nobs)            if (is.null(etastart) && is.null(start) && is.null(mustart) &&                 ((family$link == "inverse" && any(y == 0)) ||                   (family$link == "log" && any(y <= 0))))                 stop("cannot find valid starting values: please specify some")            mustart <- y        }), validmu = function (mu)         TRUE, valideta = function (eta)         TRUE), .Names = c("family", "link", "linkfun", "linkinv",     "variance", "dev.resids", "aic", "mu.eta", "initialize",     "validmu", "valideta"), class = "family"), control = structure(list(        epsilon = 1e-08, maxit = 25, trace = FALSE), .Names = c("epsilon",     "maxit", "trace")), intercept = TRUE): NA/NaN/Inf in 'y'>
## 
## 
##  How to cite this model in Zelig:
##   Kosuke Imai, Gary King, and Olivia Lau. 2015.
##   "normal: Normal Regression for Continuous Dependent Variables"
##   in Kosuke Imai, Gary King, and Olivia Lau, "Zelig: Everyone's Statistical Software,"
##   http://gking.harvard.edu/zelig
## 

You can also embed plots, for example:

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.

Copyright © 2016 thetazero.com All Rights Reserved. Privacy Policy