Skip to contents

Two measures of femur length femur.length for each of 25 walking sticks (Timema cristinae). Note that specimen is not coded as a factor.

Usage

WalkingStickFemurs

Format

A data frame with 50 observations on the following 2 variables.

specimen

a integer denoting specimen number.

femur.length

a numeric vector of femur length

Source

Nosil, P. and B.J. Crespi. 2006. Experimental evidence that predation promotes divergence in adaptive radiation. Proceedings of the National Academy of Sciences (USA) 103: 9090-9095.

Examples


demo(sec15.6)
#> 
#> 
#> 	demo(sec15.6)
#> 	---- ~~~~~~~
#> 
#> > data(WalkingStickFemurs)
#> 
#> > str(WalkingStickFemurs)
#> 'data.frame':	50 obs. of  2 variables:
#>  $ specimen    : int  1 1 2 2 3 3 4 4 5 5 ...
#>  $ femur.length: num  0.26 0.26 0.23 0.19 0.25 0.23 0.26 0.26 0.23 0.22 ...
#> 
#> > WalkingStickFemurs$specimenF <- factor(WalkingStickFemurs$specimen)
#> 
#> > fm <- aov(femur.length ~ specimenF, data = WalkingStickFemurs)
#> 
#> > summary(fm)
#>             Df  Sum Sq  Mean Sq F value   Pr(>F)    
#> specimenF   24 0.05913 0.002464   6.921 4.08e-06 ***
#> Residuals   25 0.00890 0.000356                     
#> ---
#> Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
#> 
#> > if (require(ICC)){
#> +   ICCest(specimenF,
#> +          femur.length,
#> +          data = WalkingStickFemurs,
#> +          alpha = 0.05,
#> +          CI.type = "Smith")
#> + }
#> Loading required package: ICC
#> $ICC
#> [1] 0.7475028
#> 
#> $LowerCI
#> [1] 0.5727477
#> 
#> $UpperCI
#> [1] 0.9222579
#> 
#> $N
#> [1] 25
#> 
#> $k
#> [1] 2
#> 
#> $varw
#> [1] 0.000356
#> 
#> $vara
#> [1] 0.001053917
#>