Learn R Programming

evolvability (version 2.0.1)

responsediffBeta: Calculate response differences along a set of selection gradients

Description

G1 and G2 need to be symmetric and positive definite.

Usage

responsediffBeta(G1, G2, Beta, means1 = 1, means2 = 1)

Value

An object of class

'responsdiffBeta', which is a list with the following components:

BetaThe matrix of selection gradients.d

Arguments

G1

A variance matrix.

G2

A variance matrix.

Beta

Either a vector or a matrix of unit length selection gradients stacked column wise.

means1

An optional vector of trait means of G1 (for internal mean standardization).

means2

An optional vector of trait means of G2 (for internal mean standardization).

Author

Geir H. Bolstad

Details

responsdiffBeta calculates the response difference along selection gradients between two additive-genetic variance matrices as described in Hansen and Houle (2008).

References

Hansen, T. F. & Houle, D. (2008) Measuring and comparing evolvability and constraint in multivariate characters. J. Evol. Biol. 21:1201-1219.

Examples

Run this code
G1 <- matrix(c(1, 1, 0, 1, 2, 2, 0, 2, 3), ncol = 3) / 10
G2 <- matrix(c(1, 2, 0,2, 1, 1,0, 1, 3), ncol = 3)
Beta <- randomBeta(5, 3)
X <- responsediffBeta(G1, G2, Beta)
summary(X)

Run the code above in your browser using DataLab