library(bivpois)
dataku <- read.csv("E://datapaperundip.csv",sep=",",header=TRUE)
y1=dataku[,1]
y2=dataku[,2]
x1=dataku[,3]
x2=dataku[,4]
x3=dataku[,5]
x4=dataku[,6]
x5=dataku[,7]
x6=dataku[,8]
x7=dataku[,9]
x8=dataku[,10]
#varians konstan
model.a<-lm.bp(y1~x1+x2+x3+x4+x5+x6+x7+x8, y2~x1+x2+x3+x4+x5+x6+x7+x8,data=dataku)
model.a$coef
model.a$parameters
model.a$AIC
model.a$loglikelihood


#standard error varians konstan
n=length(y1)
bootrep <- 200
results.a<-matrix(NA,bootrep,19)
for (i in 1:bootrep)
{
bootx1<-rpois(n,model$lambda1)
bootx2<-rpois(n,model$lambda2)
bootx3<-rpois(n,model$lambda3)
bootx<-bootx1+bootx3
booty<-bootx2+bootx3
data = cbind(bootx,booty,x1,x2,x3,x4,x5,x6,x7,x8)
data=data.frame(data)
testtempt<-lm.bp (bootx~x1+x2+x3+x4+x5+x6+x7+x8,booty~x1+x2+x3+x4+x5+x6+x7+x8,data=data)
betafound<-c(testtempt$beta1,testtempt$beta2,testtempt$beta3)
results.a[i,]<-betafound
}
write.table(results.a, file = "E://BAHAN BELAjAR/PENELITIAN_KLATEN/Jurnal POISSON/covkonstan.csv", sep = ",", col.names = NA, qmethod = "double")