Diantres, y ahora
debería hacer una superficie de selección para ver que pasa…
Bueno, hay varias
opciones, pero como sabemos que las superficies paramétricas (esas que resultan
de usar los parámetros del modelo de selección de Lande y Arnold), pueden dar una idea engañosa de la realidad
y sólo son capaces de representar un rango limitado de formas (no lo digo
yo…) vamos a usar una superficie no paramétrica, cuyo origen se rastrea a Dolph
Schluter y sus splines cúbicos (1,2).
Hace un tiempo, leyendo
el paper de Mcglothlin et al. (3) encontré cómo hacerlo en R con el paquete fields. Desde entonces lo he aplicado,
publicado, y la gente parece no quejarse al respecto.
Necesitan dos
medidas fenotípicas y una medida de fitness. la rutina es en extremo sencilla.
#cargamos el paquete
library(fields)
#no deben existir datos faltantes.
datos<- na.omit(datos)
#unimos las dos variables fenotípicas
sup<-cbind(datos$X1, datos$X2)
#ajustamos el modelo
#el primer término son las variables
predictoras,
#el segundo es el fitness
fit<-Tps(sup,
datos$fitness)
Y ya. La explicación
de lo que el programa hizo, es un poquito oscura (tomen el curso de Métodos en
Ecología Evolutiva que damos con Mariano Ordano y Andrea Cocucci!! pero no
prometo que la entiendan de todas formas). Se ha calculado una aproximación a una superficie no paramétrica, que es
lo que necesitamos para visualizar la selección correlacionada. Esto último se puede
hacer con image y persp.
image(predict.surface(fit))
persp(predict.surface(fit))
A lo cual, por
supuesto, les podemos manipular los parámetros gráficos para que se vean mejor.
image(predict.surface(fit),
col = gray(seq(0, 1, by=0.02)),
xlab="nectario",
ylab="número de flores")
contour(predict.surface(fit),
add=TRUE, col="red")
points(x=datos$nectario, y=datos$flores,
pch=20)
persp(predict.surface(fit),
border="slateblue4", col= "aquamarine2", xlab=
"nectario", ylab= "número de flores", zlab=
"polinarios", theta= -45, phi= 20, r=5, d=1, ticktype=
"detailed", shade= 0.1 )
1-Schluter
D. 1988. Estimating the form of natural selection on a quantitative trait. Evolution 42:
849-861.
2-Schluter
D, Nychka D. 1994. Exploring fitness surfaces. The American Naturalist 143:
597-616.
3-Mcglothlin JW,
Parker PG, Nolan VJ, Ketterson ED. 2005. Correlational selection leads to
genetic integration of body size and an attractive plumage trait in dark-eyed
juncos. Evolution 59: 658-671.