library("segmented", lib.loc="~/R/win-library/3.1")
#a: (bontott csirke ár/sertéscomb ár)*100
#b: 1 főre jutó sertés/kg/(sertés/kg+baromfi/kg)*100
a<-c(60.9, 60.8, 60.1, 70.3, 57.1, 48.6, 53.4, 59.8, 57.5, 56.6, 53.8, 63.9, 64.8, 63.4, 65.7, 64.5, 61.0, 59.6, 60.1)
b<-c(46.3, 47.4, 49.8, 46.1, 54.6, 57.6, 55.3, 54.7, 51.7, 52.7, 52.5, 51.0, 52.7, 50.7, 49.3, 49.6, 50.9, 51.9, 52.8)
husfogy<-data.frame(css_ararany=a, csfogyasztas_reszesedese=b)
attach(husfogy)
lm(csfogyasztas_reszesedese~css_ararany)
lm.out=lm(csfogyasztas_reszesedese~css_ararany)
summary(lm.out)
plot(ylim=c(44,60), xlim=c(48,73), csfogyasztas_reszesedese~css_ararany, main="A baromfi- és sertéshús fogyasztásának vizsgálata (1996-2014)",
     xlab="Bontott csirke ár/Sertéscomb ár(%)", ylab="A csirkefogyasztás részesedése(%)")
abline(lm.out, col="blue", lwd=2)
abline(v=60.148, col="green", lwd=2)
abline(v=71.605, col="red", lwd=2)
#ccs áfacsökk. előtt (2013)=60,148
#ccs áfacsökk. után (2013)=71,605
eh_lm=coefficients(lm.out)
fogyar_afacs_el_lm=eh_lm[1]+60.148*eh_lm[2]
fogyar_afacs_ut_lm=eh_lm[1]+71.605*eh_lm[2]
cor(csfogyasztas_reszesedese, css_ararany)
szak_lin_reg1 <- segmented(lm(csfogyasztas_reszesedese~css_ararany), seg.Z=~css_ararany, psi=c(60,65))#,control = seg.control(stop.if.error = FALSE))
summary(szak_lin_reg1)
eh_szlm=coefficients(szak_lin_reg1)
c<-slope(szak_lin_reg1)
slope_szak_lin_reg1<-data.frame(c)
fogyar_afacs_el_szlm = eh_szlm[1]+eh_szlm[2]*60.148
fogyar_afacs_ut_szlm = eh_szlm[1]+slope_szak_lin_reg1[1,1]*szak_lin_reg1$psi[3]+slope_szak_lin_reg1[2,1]*(szak_lin_reg1$psi[4]-szak_lin_reg1$psi[3])+slope_szak_lin_reg1[3,1]*(71.605-szak_lin_reg1$psi[4])
plot(ylim=c(44,60), xlim=c(48,73),csfogyasztas_reszesedese~css_ararany,main="A baromfihús (1996-2014)",
     xlab="Bontott csirke ár/sertéscomb ár(%)", ylab="A baromfihús-fogyasztás részesedése (%)")
plot(szak_lin_reg1, add = TRUE, col = "blue", lwd=2)
#segments(szak_lin_reg1$psi[4],eh_szlm[1]+slope_szak_lin_reg1[1,1]*szak_lin_reg1$psi[3]+slope_szak_lin_reg1[2,1]*(szak_lin_reg1$psi[4]-szak_lin_reg1$psi[3]),73, eh_szlm[1]+slope_szak_lin_reg1[1,1]*szak_lin_reg1$psi[3]+slope_szak_lin_reg1[2,1]*(szak_lin_reg1$psi[4]-szak_lin_reg1$psi[3])+slope_szak_lin_reg1[3,1]*(73-szak_lin_reg1$psi[4]),col="blue")
abline(v=60.148, col="green", lwd=2)
abline(v=71.605, col="red", lwd=2)
együtthatók=coefficients(szak_lin_reg1)
lm_val=fogyar_afacs_el_lm-fogyar_afacs_ut_lm
szlm_val=fogyar_afacs_el_szlm-fogyar_afacs_ut_szlm
paste("A csirkefogyasztás részesedése az áfacsökkentés előtt(lm):", as.numeric(fogyar_afacs_el_lm[1]),".")
paste("A csirkefogyasztás részesedése az áfacsökkentés után(lm):", as.numeric(fogyar_afacs_ut_lm[1]),".")
paste("A csirkefogyasztás részesedése az áfacsökkentés előtt(szlm):", as.numeric(fogyar_afacs_el_szlm[1]),".")
paste("A csirkefogyasztás részesedése az áfacsökkentés eután(szlm):", as.numeric(fogyar_afacs_ut_szlm[1]),".")
paste("A csirkefogyasztás részesedésének változása(lm):", as.numeric(lm_val[1]),"%.")
paste("A csirkefogyasztás részesedésének változása(szlm):", as.numeric(szlm_val[1]),"%.")
p_1=c(fogyar_afacs_el_szlm,100-fogyar_afacs_el_szlm)
p_1_labels=round(p_1,1)
p_1_labels<-paste(p_1_labels, "%", sep="")
pie(p_1, main="A sertés-, baromfihús-fogyasztás részesedése az áfacsökkentés előtt", col=(colors()[c(11,427)]), labels=p_1_labels)
p_2=c(fogyar_afacs_ut_szlm,100-fogyar_afacs_ut_szlm)
legend(0.9, 0.8, c("Baromfihús fogyasztás","Sertéshús fogyasztás"), cex=0.8,
       fill=(colors()[c(11,427)]))
p_2_labels=round(p_2,1)
p_2_labels<-paste(p_2_labels, "%", sep="")
pie(p_2, main="A sertés-, baromfihús-fogyasztás részesedése az áfacsökkentés után", col=(colors()[c(11,427)]), labels=p_2_labels)
legend(0.9, 0.8, c("Baromfihús fogyasztás","Sertéshús fogyasztás"), cex=0.8,
       fill=(colors()[c(11,427)]))