getwd()
[1] "/cloud/project"
regression1<-read.csv("incidents.csv",header=T,sep = ",")
str(regression1)
'data.frame':   16 obs. of  4 variables:
 $ area      : chr  "Boulder" "California-lexington" "Huntsville" "Seattle" ...
 $ zone      : chr  "west" "east" "east" "west" ...
 $ population: chr  "107,353" "326,534" "444,752" "750,000" ...
 $ incidents : int  605 103 161 1703 1003 527 721 704 105 403 ...
summary(regression1)
     area               zone            population          incidents     
 Length:16          Length:16          Length:16          Min.   : 103.0  
 Class :character   Class :character   Class :character   1st Qu.: 277.8  
 Mode  :character   Mode  :character   Mode  :character   Median : 654.0  
                                                          Mean   : 695.2  
                                                          3rd Qu.: 853.0  
                                                          Max.   :2072.0  
# make sure the packages for this chapter
# are installed, install if necessary
pkg <- c("ggplot2", "scales", "maptools",
              "sp", "maps", "grid", "car" )
new.pkg <- pkg[!(pkg %in% installed.packages())]
if (length(new.pkg)) {
 
}
NULL
regression1$population <- as.numeric(gsub(",","",regression1$population))
regression1$population
 [1]  107353  326534  444752  750000   64403 2744878 1600000 2333000 1572816  712091 6900000 2700000
[13] 4900000 4200000 5200000 7100000
regression2<-regression1[,-1]
G3;  File "<string>", line 1
    regression2<-regression1[,-1]
                             ^
SyntaxError: invalid syntax
g
head(regression2)
G3;Traceback (most recent call last):
  File "<string>", line 1, in <module>
NameError: name 'head' is not defined
g
#reg.fit1<-lm(regression2$incidents~regression2$population)
summary(reg.fit1)
Error: object 'reg.fit1' not found
regression2$zone<-ifelse(regression2$zone=="west",1,0)
Error: object 'regression2' not found
reg.fit4<-lm(regression2$incidents~interaction)
Error in eval(predvars, data, env) : object 'regression2' not found
LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKCgpgYGB7cn0KZ2V0d2QoKQpgYGAKYGBge3J9CnJlZ3Jlc3Npb24xPC1yZWFkLmNzdigiaW5jaWRlbnRzLmNzdiIsaGVhZGVyPVQsc2VwID0gIiwiKQpzdHIocmVncmVzc2lvbjEpCmBgYApgYGB7cn0Kc3VtbWFyeShyZWdyZXNzaW9uMSkKYGBgCmBgYHtyfQojIG1ha2Ugc3VyZSB0aGUgcGFja2FnZXMgZm9yIHRoaXMgY2hhcHRlcgojIGFyZSBpbnN0YWxsZWQsIGluc3RhbGwgaWYgbmVjZXNzYXJ5CnBrZyA8LSBjKCJnZ3Bsb3QyIiwgInNjYWxlcyIsICJtYXB0b29scyIsCiAgICAgICAgICAgICAgInNwIiwgIm1hcHMiLCAiZ3JpZCIsICJjYXIiICkKbmV3LnBrZyA8LSBwa2dbIShwa2cgJWluJSBpbnN0YWxsZWQucGFja2FnZXMoKSldCmlmIChsZW5ndGgobmV3LnBrZykpIHsKIAp9CmBgYAoKCmBgYHtyfQpyZWdyZXNzaW9uMSRwb3B1bGF0aW9uIDwtIGFzLm51bWVyaWMoZ3N1YigiLCIsIiIscmVncmVzc2lvbjEkcG9wdWxhdGlvbikpCnJlZ3Jlc3Npb24xJHBvcHVsYXRpb24KYGBgCmBgYHtweXRob259CnJlZ3Jlc3Npb24yPC1yZWdyZXNzaW9uMVssLTFdCmhlYWQocmVncmVzc2lvbjIpCmBgYApgYGB7cn0KI3JlZy5maXQxPC1sbShyZWdyZXNzaW9uMiRpbmNpZGVudHN+cmVncmVzc2lvbjIkcG9wdWxhdGlvbikKc3VtbWFyeShyZWcuZml0MSkKYGBgCmBgYHtyfQpyZWdyZXNzaW9uMiR6b25lPC1pZmVsc2UocmVncmVzc2lvbjIkem9uZT09Indlc3QiLDEsMCkKaW50ZXJhY3Rpb248LXJlZ3Jlc3Npb24yJHpvbmUqcmVncmVzc2lvbjIkcG9wdWxhdGlvbgpyZWcuZml0MzwtbG0ocmVncmVzc2lvbjIkaW5jaWRlbnRzfmludGVyYWN0aW9uK3JlZ3Jlc3Npb24yJHBvcHVsYXRpb24rcmVncmVzc2lvbjIkem9uZSkKc3VtbWFyeShyZWcuZml0MykKYGBgCmBgYHtyfQpyZWcuZml0NDwtbG0ocmVncmVzc2lvbjIkaW5jaWRlbnRzfmludGVyYWN0aW9uKQpzdW1tYXJ5KHJlZy5maXQ0KQpgYGAKCgo=