@@ -33,14 +33,14 @@ ylim1.1 <- range(samples2[,"phi1"])*c(1, 1.15)
3333par(mar = pm.1 )
3434plot(samples2 [," phi1" ], type = " l" , col = " gray" , xlab = " " , ylab = " " , bty = " l" , xpd = TRUE , ylim = ylim1.1 ) # plot full samples of variable phi1
3535points(samples2.max [," phi1.i" ], samples2.max [," phi1" ], col = " blue" , xpd = TRUE ) # add circles for annual maxima of phi1
36- mtext(bquote(phi1 ), side = 2 , line = 0.75 )
36+ mtext(bquote(Driver ~ ( phi1 ) ), side = 2 , line = 0.75 )
3737text(x = 25 , y = 0.85 , " A" , font = 2 )
3838
3939par(mar = pm.1 )
4040ylim1.2 <- range(samples2 [," phi2" ])* c(1 , 1.15 )
4141plot(samples2 [," phi2" ], type = " l" , col = " gray" , xlab = " " , ylab = " " , bty = " l" , xpd = TRUE , ylim = ylim1.2 ) # plot full samples of variable phi2
4242points(samples2.max [," phi2.i" ], samples2.max [," phi2" ], col = " red" , xpd = TRUE ) # add circles for annual maxima of phi2
43- mtext(bquote(gamma ), side = 2 , line = 0.75 )
43+ mtext(bquote(Driver ~ ( gamma ) ), side = 2 , line = 0.75 )
4444text(x = 25 , y = 0.85 , " B" , font = 2 )
4545
4646 # ==========
@@ -50,7 +50,7 @@ text(x=25, y=0.85, "B", font=2)
5050par(mar = pm.2 )
5151plot(samples2 [," X" ], type = " l" , col = " gray" , xlab = " " , ylab = " " , bty = " l" , xpd = TRUE ) # plot full samples of variable phi1
5252points(samples2.max [," X.i" ], samples2.max [," X" ], col = " black" , xpd = TRUE ) # add circles for annual maxima of phi1
53- mtext(bquote(X ), side = 2 , line = 0.75 )
53+ mtext(bquote(Population ~ ( X ) ), side = 2 , line = 0.75 )
5454mtext(bquote(Time ), side = 1 , line = 0.75 )
5555text(x = 25 , y = 115 , " C" , font = 2 )
5656
@@ -86,7 +86,7 @@ mtext(bquote(Density), side=1, line=0.75)
8686# part 4: difference in xi's between population and the environment
8787par(mar = pm.3 )
8888colorDens(vals = list (xis [,1 ], xis [,2 ], xis [,3 ], (xis [,3 ]- (xis [,1 ]+ xis [,2 ]))), cols = c(" blue" ," red" ," black" ," green" ), revxy = FALSE , yaxt = " s" , bty = " l" )
89- mtext(bquote(xi ), side = 1 , line = 0.75 )
89+ mtext(bquote(Tailedness ~ ( xi ) ), side = 1 , line = 0.75 )
9090mtext(bquote(Density ), side = 2 , line = 0.75 )
9191text(y = 2.5 , x = - 0.75 , " G" , font = 2 )
9292
0 commit comments