Ошибка R не будет связана с PDF в результате "ypos" - PullRequest
Новая
       11

Ошибка R не будет связана с PDF в результате "ypos"

0 голосов
/ 26 марта 2020

Я действительно застрял, пытаясь связать свой проект в PDF. В настоящее время он не будет вязать, как говорится ypos is not a graphical function

Вот код. Пожалуйста, дайте мне знать любые советы / хитрости / указатели, которые помогут мне улучшить. 

##r code to reproduce results for part (b) 
##read data file 
dat = read.csv("parking.csv", comment.char = "#", as.is=TRUE)
##remove rows with missing data 
dat = dat[dat$prepost!=""&!is.na(dat$viol_pc), ];
##only keep variables to use 
dat = dat[, c("wbcode", "country", "corruption", "prepost", "viol_pc")]
##check data 
#print(str(dat))

##define high/low corrupt countries 
dat$hilo = ifelse(dat$corruption>0, "hi" , "lo");
##check distribution of hi/lo 
#print (table(dat$hilo))

##compute alternative measures of parking violation 
dat$viol_log = log(1+dat$viol_pc + sqrt(dat$viol_pc^2 + 1));
##check correlations between three measures 
##print(cor(dat[,c(viol_pc" , "viol_log" , "viol_ihs")]))

##generate pre/post factor where pre comes before post 
dat$fprepost = factor(dat$prepost, c("pre", "post"));

##informal comparison pf parking violation hi/lo
#boxplot(viol_pc~ fprepost, dat, xlab="",
#         ylab="parking violation/diplomat",
#         ylim=c(0, 50));

##informal comparison of parking violations hi/lo
#boxplot(viol_ihs ~ hilo, dat, xlab="", 
#       ylab=ihs("Parking violation/diplomat)")

##test difference in means pre/post for(viol_pc)
pp_log ~ t.test(dat$viol_log[dat$prepost=="pre"], 
                dat$viol_log[dat$prepost=="pos"])
#print(str(pp_log))

##plot 95% confidence interval for difference in means 
x = pp_log$conf.int;
y = c(0,0);
par(mar=c(3, 3, 7, 2), bty="n", yaxt="n")
plot(x,y, type="l", xlim=range(0, x),
     lwd=17, col=gray(1), lend=3, 
    main="95% confidence interval 
          for the difference in ihs transform of parking violations 
          before/after change in law");
##add vertical line at zero 
abline(v=0, lty="dashed");
##add pre/post means 
points(pp_log$estimate[1], 0);
text(pp_log$estimate[1], 0, "mean(pre)" , pos=3, xpd=TRUE);
points(pp_log$estimate[2],0);
text(pp_log$estimate[2], 0, "mean(post)", pos=3, xpd=TRUE);

Как вы можете видеть, ypos не упоминается, поэтому я не уверен, где искать. Может кто-нибудь указать, где мне следует искать и как исправить эту ошибку. Большое спасибо:)

...