Skip to content

Instantly share code, notes, and snippets.

@erikseulean

erikseulean/post_impact.r

Last active January 25, 2022 16:02
Show Gist options
  • Save erikseulean/639fd875c14582066b11b152d829f31b to your computer and use it in GitHub Desktop.
Save erikseulean/639fd875c14582066b11b152d829f31b to your computer and use it in GitHub Desktop.
Download ZIP
interval
Raw
post_impact.r
npost = nrow(filter(df, Impact==1))
mm = mean(filter(df, Impact==1)$AbundanceRate)
ll = mm - qnorm(1 - alpha/2) * sd(filter(df, Impact==1)$AbundanceRate)/sqrt(npost)
uu = mm + qnorm(1 - alpha/2) * sd(filter(df, Impact==1)$AbundanceRate)/sqrt(npost)
print(paste("lower", round(ll, 2)))
print(paste("upper", round(uu, 2)))
Raw
quilt.r
library(fields)
par(mrow=c(2,2))
quilt.plot(df$XPos[df$Impact==0], df$YPos[df$Impact==0], df$Nhat[df$Impact==0]/df$Area[df$Impact==0], pch=20,
main="Impact = 0", xlab="X co-ordinate", ylab="Y co-ordinate",
nrow=20, ncol=20, asp=1)
quilt.plot(df$XPos[df$Impact==1], df$YPos[df$Impact==1], df$Nhat[df$Impact==1]/df$Area[df$Impact==1], pch=20,
main="Impact = 1", xlab="X co-ordinate", ylab="Y co-ordinate",
nrow=20, ncol=20, asp=1)
Raw
standard_deviation.r
alpha = 0.02
quant = qt(p=alpha/2, df=n-1, lower.tail=F)
n = nrow(filter(df, Impact==0))
m = mean(filter(df, Impact==0)$AbundanceRate)
s = sd(filter(df, Impact==0)$AbundanceRate)
lower = m - quant * s/sqrt(n)
upper = m + quant * s/sqrt(n)
print(paste("lower", round(lower, 2)))
print(paste("upper", round(upper, 2)))
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment