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=