# KFanalysis for spruce budworm, Choristoneura fumiferana

# Data are from Royama (1966) (Ecology, 77:87-93). D indicates foliage density of the previous year. It was read from Figure 2B of Morris et al. (1958) (Proceedings of the 10th International Congress of Entomology 1958 Vol. 4 Issue 137-149).

cat(file="Royama1996.txt",
"year D H1 H2 H3 H4 Hg Hr R X
1947 0.92 -0.023 -0.782 -0.62 -0.048 -1.473 2.256 0.783 0.407
1948 0.91 -0.023 -0.669 -0.791 -0.176 -1.659 2.063 0.404 1.19
1949 0.89 -0.143 -0.879 -0.433 -0.2 -1.655 2.617 0.962 1.594
1950 0.88 -0.2 -0.616 -0.702 -0.157 -1.675 1.78 0.105 2.556
1951 0.85 -0.126 -0.827 -0.601 -0.125 -1.679 1.306 -0.373 2.661
1952 0.82 -0.06 -0.582 -0.963 -0.102 -1.707 1.648 -0.059 2.288
1953 0.77 -0.038 -0.798 -1.285 -0.193 -2.314 2.846 0.532 2.228
1954 0.72 -0.086 -0.888 -0.905 -0.337 -2.216 1.643 -0.573 2.76
1955 0.63 -0.119 -1 -0.878 -0.398 -2.395 2.511 0.116 2.188
1956 0.49 -0.018 -1.491 -0.726 -0.283 -2.518 2.826 0.308 2.303
1957 0.28 -0.08 -1.22 -1.316 -0.288 -2.904 1.532 -1.372 2.612
1958 0.10 -0.041 -0.873 -1.724 -0.301 -2.939 1.978 -0.961 1.24
")

Royama1996 <- read.table("Royama1996.txt", header=TRUE)

source("B0_R_program_key-factor_key-stage_analysis.txt")

windows(width = 8*0.8, height =  6*0.8, record=TRUE)
KFanalysis(factor="D,X",stage="Hg,Hr",data=Royama1996)