# Revisions dat<-read.table("genital_morpho_master_with_pca.txt",header=T) contacts<-rbind(dat[dat$SPOT=="CZW",],dat[dat$SPOT=="CZE",]) model<-aov(contacts$PC1.y~ contacts$cycle + contacts$yeartype + contacts$SPOT + contacts$cycle*contacts$SPOT + contacts$yeartype*contacts$SPOT) summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$cycle 1 19.7 19.66 12.471 0.000479 *** contacts$yeartype 1 23.7 23.74 15.058 0.000129 *** contacts$SPOT 1 55.6 55.56 35.245 8.15e-09 *** contacts$cycle:contacts$SPOT 1 3.6 3.62 2.295 0.130863 contacts$yeartype:contacts$SPOT 1 16.8 16.85 10.690 0.001205 ** Residuals 295 465.0 1.58 dat<-read.table("wing_contemporary_morpho_master_with_pca.txt",header=T) contacts<-rbind(dat[dat$SPOT=="CZW",],dat[dat$SPOT=="CZE",]) model<-aov(contacts$PC1.y~ contacts$cycle + contacts$yeartype + contacts$SPOT + contacts$cycle*contacts$SPOT + contacts$yeartype*contacts$SPOT) summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$cycle 1 0.0013 0.00133 0.763 0.3830 contacts$yeartype 1 0.1172 0.11719 67.383 5.80e-15 *** contacts$SPOT 1 0.0660 0.06604 37.972 2.19e-09 *** contacts$cycle:contacts$SPOT 1 0.0035 0.00353 2.029 0.1553 contacts$yeartype:contacts$SPOT 1 0.0053 0.00532 3.058 0.0813 . Residuals 315 0.5478 0.00174 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 model<-aov(contacts$PC2.y~ contacts$cycle + contacts$yeartype + contacts$SPOT + contacts$cycle*contacts$SPOT + contacts$yeartype*contacts$SPOT) summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$cycle 1 0.00144 0.001436 2.111 0.1473 contacts$yeartype 1 0.00345 0.003452 5.075 0.0250 * contacts$SPOT 1 0.00578 0.005784 8.504 0.0038 ** contacts$cycle:contacts$SPOT 1 0.00040 0.000398 0.585 0.4449 contacts$yeartype:contacts$SPOT 1 0.00066 0.000662 0.973 0.3248 Residuals 315 0.21424 0.000680 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 model<-aov(contacts$PC3.y~ contacts$cycle + contacts$yeartype + contacts$SPOT + contacts$cycle*contacts$SPOT + contacts$yeartype*contacts$SPOT) summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$cycle 1 0.00375 0.003748 11.784 0.000677 *** contacts$yeartype 1 0.00249 0.002491 7.831 0.005453 ** contacts$SPOT 1 0.00003 0.000028 0.089 0.765217 contacts$cycle:contacts$SPOT 1 0.00044 0.000437 1.373 0.242114 contacts$yeartype:contacts$SPOT 1 0.00076 0.000755 2.374 0.124345 Residuals 315 0.10020 0.000318 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 model<-aov(contacts$PC4.y~ contacts$cycle + contacts$yeartype + contacts$SPOT + contacts$cycle*contacts$SPOT + contacts$yeartype*contacts$SPOT) summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$cycle 1 0.00183 0.0018305 6.252 0.01291 * contacts$yeartype 1 0.00008 0.0000848 0.290 0.59074 contacts$SPOT 1 0.00150 0.0015032 5.134 0.02414 * contacts$cycle:contacts$SPOT 1 0.00003 0.0000270 0.092 0.76166 contacts$yeartype:contacts$SPOT 1 0.00202 0.0020222 6.907 0.00901 ** Residuals 315 0.09223 0.0002928 model<-aov(contacts$PC5.y~ contacts$cycle + contacts$yeartype + contacts$SPOT + contacts$cycle*contacts$SPOT + contacts$yeartype*contacts$SPOT) > > summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$cycle 1 0.00060 0.000605 1.305 0.254 contacts$yeartype 1 0.01715 0.017152 37.009 3.41e-09 *** contacts$SPOT 1 0.00928 0.009281 20.024 1.07e-05 *** contacts$cycle:contacts$SPOT 1 0.00042 0.000423 0.913 0.340 contacts$yeartype:contacts$SPOT 1 0.01484 0.014835 32.010 3.45e-08 *** Residuals 315 0.14599 0.000463 # Historic dat<-read.table("wing_historic_morpho_master_with_pca.txt",header=T) contacts<-rbind(dat[dat$SPOT2=="CZW",],dat[dat$SPOT2=="CZE",]) model<-aov(contacts$PC1~ contacts$SPOT2 + contacts$yeartype + t*contacts$yeartype) Df Sum Sq Mean Sq F value Pr(>F) contacts$SPOT2 1 0.00814 0.008145 2.859 0.0938 . contacts$yeartype 1 0.00537 0.005370 1.885 0.1726 contacts$SPOT2:contacts$yeartype 1 0.00516 0.005157 1.811 0.1813 Residuals 106 0.30193 0.002848 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 model<-aov(contacts$PC2~ contacts$SPOT2 + contacts$yeartype + contacts$SPOT2*contacts$yeartype) Df Sum Sq Mean Sq F value Pr(>F) contacts$SPOT2 1 0.00005 0.000048 0.056 0.813 contacts$yeartype 1 0.01929 0.019285 22.594 6.33e-06 *** contacts$SPOT2:contacts$yeartype 1 0.00167 0.001675 1.962 0.164 Residuals 106 0.09048 0.000854 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > model<-aov(contacts$PC3~ contacts$SPOT2 + contacts$yeartype + contacts$SPOT2*contacts$yeartype) > > summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$SPOT2 1 0.00018 0.0001796 0.199 0.657 contacts$yeartype 1 0.00062 0.0006181 0.684 0.410 contacts$SPOT2:contacts$yeartype 1 0.00123 0.0012257 1.357 0.247 Residuals 106 0.09576 0.0009034 model<-aov(contacts$PC4~ contacts$SPOT2 + contacts$yeartype + contacts$SPOT2*contacts$yeartype) summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$SPOT2 1 0.00175 0.001753 3.229 0.0752 . contacts$yeartype 1 0.00101 0.001006 1.853 0.1763 contacts$SPOT2:contacts$yeartype 1 0.00173 0.001726 3.179 0.0774 . Residuals 106 0.05755 0.000543 > model<-aov(contacts$PC5~ contacts$SPOT2 + contacts$yeartype + contacts$SPOT2*contacts$yeartype) > > summary(model) Df Sum Sq Mean Sq F value Pr(>F) contacts$SPOT2 1 0.00034 0.0003441 0.675 0.4130 contacts$yeartype 1 0.00256 0.0025604 5.026 0.0271 * contacts$SPOT2:contacts$yeartype 1 0.00009 0.0000904 0.177 0.6745 Residuals 106 0.05400 0.0005094 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1