Ochlero<-read.table("culicid_DT.txt",header=TRUE)
head(Ochlero)
with(Ochlero,plot(Temperature_C,ln.Development_time.,col=localities,pch=as.numeric(localities)+7))
model<-lm(Ochlero$ln.Development_time. ~ Ochlero$Temperature_C)
par(mfrow=c(2,4))
llmodel<-lm(log(log(Ochlero$ln.Development_time.)) ~ Ochlero$Temperature_C)
plot(model)
plot(llmodel)

Sar <- subset(Ochlero, localities=="Sarmiento")
BAires <- Ochlero[Ochlero$localities=="BuenosAires",]
regrSar<-lm(log(log(Sar$ln.Development_time.)) ~ Sar$Temperature_C)
regrBA<-lm(log(log(BAires$ln.Development_time.)) ~ BAires$Temperature_C)
summary(regrSar)[[4]][2]; summary(regrBA)[[4]][2]
summary(regrSar)[[4]][4]; summary(regrBA)[[4]][4]

Z<-abs(summary(regrSar)[[4]][2]-summary(regrBA)[[4]][2])/ sqrt(summary(regrSar)[[4]][4]^2+ summary(regrBA)[[4]][4]^2)
#Z
qt(0.975,83)