library(rJava)
library(xlsxjars)
setwd("/Users/charlotterunzel/Desktop/Results Data")
library(xlsx)
A <- read.xlsx("Cell Counts_Final copy.xlsx", sheetName="Sheet2")
library(ggplot2)
p3 = (3/69)
p2 = (2/69)
p5 = (5/69)
p4 = (4/69)
zz=0.0121
aa = ((mean(A[1:2,6])) - zz)^2 #2 - 32.51
bb = ((mean(A[3:6,6])) - zz)^2 #4 - 31.75
cc = ((mean(A[7:10,6])) - zz)^2 #4 - 30.63
dd = ((mean(A[11:13,6])) - zz)^2 #3 - 33.35
ee = ((mean(A[14:17,6])) - zz)^2 #4 - 31.1
ff = ((mean(A[18:22,6])) - zz)^2 #5 - 30.46
gg = ((mean(A[23:27,6])) - zz)^2 #5 - 29.1 
hh = ((mean(A[28:32,6])) - zz)^2 #5 - 32.4 
ii = ((mean(A[33:36,6])) - zz)^2 #4 - 33.95
jj = ((mean(A[37:41,6])) - zz)^2 #5 - 37.7
kk = ((mean(A[42:45,6])) - zz)^2 #4 - 29.7 
ll = ((mean(A[46:50,6])) - zz)^2 #5 - 26.365
mm = ((mean(A[51:55,6])) - zz)^2 #5 - 27.76
nn = ((mean(A[56:60,6])) - zz)^2 #5 - 4.96
oo = ((mean(A[61:64,6])) - zz)^2 #4 - 21.1
pp = ((mean(A[65:69,6])) - zz)^2 #5 - 15.5
qq = bb + cc + ee + ii + kk + oo 
rr = ff + gg + hh + jj + ll + mm + nn + pp

top = p2*(aa) +
  p3*(dd) + 
  p4*(qq) + 
  p5*(rr)

test = lm(Cover~Temp,A)
val = anova(test)["Residuals", "Mean Sq"]
f = sqrt(top/val)
library(pwr)
pwr.anova.test(k=5, f=f, sig.level=0.05, power=.8) #effect size is larger than .8

x=lm(data=A,Cover~Temp)
coefficients(x)
testa=lm(A$Cover~A$Temp)
summary(testa)

library(ggplot2)
plot1a = ggplot(data=A, aes(x=Temp, y=Cover)) + 
  geom_point() + 
  xlab("Temperature (ºC)") +
  ylab("% of Cells Aggregated") + geom_smooth(method = "loess") + 
  theme(axis.text=element_text(size=9)) 
  #annotate("text", x=15, y=0.035, label = "Residual Standard Error<0.01")
plot1a

plot1b = ggplot(data=A, aes(x=Temp, y=Cover)) + 
  geom_point() + 
  xlab("Temperature (ºC)") + 
  #ylab("% of Cells Aggregated") + geom_smooth(method = "loess")  
#annotate("text", x=15, y=0.035, label = "Residual Standard Error<0.01")
plot1b

test =loess(A$Cover~A$Temp)
summary(test)