代码是老师给的，验证本福特定律，但是R一直提示错误，我不知道我少装什么程序包了


library(xtable)
> options(encoding="latin1")
>  
> pob < - read.csv("http://www.grserrano.es/datos/pobmun09.csv")
错误: 找不到对象'pob'
> head(pob)
错误于head(pob) : 找不到对象'pob'
>  
> totprov <- pob[ is.na(pob[, "CMUN"]), ]   ## Totales de provincia
错误: 找不到对象'pob'
> pob     <- pob[!is.na(pob[, "CMUN"]), ]   ## Dejo sólo los municipios
错误: 找不到对象'pob'
>  
> pob[, "digpob"]  <- as.numeric(substr(as.character(pob[, "POB09"]), 1, 1))
错误于substr(as.character(pob[, "POB09"]), 1, 1) : 找不到对象'pob'
> digitos <- c(1:9)
> (abs.freq     <- table(pob[, "digpob"]))
错误于table(pob[, "digpob"]) : 找不到对象'pob'
> (benford.prob <- log10(1+1/digitos))
[1] 0.30103000 0.17609126 0.12493874 0.09691001 0.07918125 0.06694679
[7] 0.05799195 0.05115252 0.04575749
> tabmun  <- cbind( Observada  = prop.table(abs.freq),
+                   Teorica    = benford.prob,
+                   Municipios = abs.freq )
错误于prop.table(abs.freq) : 找不到对象'abs.freq'
> xtabmun <- xtable(tabmun, digits=c(0, 4,4,0), align="r|rr|r|",
+       caption="Frecuencias observadas, esperadas y número de municipios por el primer dígito de la población de los 8112 municipios espa&ntilde;oles a 1 de enero de 2009",
+       label="tab:benfordmun")
错误于xtable(tabmun, digits = c(0, 4, 4, 0), align = "r|rr|r|", caption = "Frecuencias observadas, esperadas y número de municipios por el primer dígito de la población de los 8112 municipios espa<U+00F1>oles a 1 de enero de 2009",  :
  找不到对象'tabmun'
> print(xtabmun, type="latex", file="tabmun.tex")
错误于print(xtabmun, type = "latex", file = "tabmun.tex") :
  找不到对象'xtabmun'
> #print(xtabmun, type="HTML", file="tabmun.html")
>  
> plot(digitos, tabmun[, "Observada"],
+      main="Ley de Benford para población de municipios espa&ntilde;oles",
+      xlab="Primer dígito", ylab="Frecuencia")
错误于xy.coords(x, y, xlabel, ylabel, log) : 找不到对象'tabmun'
> lines(digitos, tabmun[, "Teorica"], lty=2, col="blue")
错误于xy.coords(x, y) : 找不到对象'tabmun'
> legend("topright", c("Observado", "Teórico"), lty=c(0,2),
+        pch=c(1,NA), col=c("black", "blue"))
> grid()
> utils:::menuInstallPkgs()