-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy path.Rapp.history
51 lines (51 loc) · 2.15 KB
/
.Rapp.history
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
2*2
library(rptR)
browseVignettes(rptR)
browseVignettes("rptR")
library(devtools)
install_github("mastoffel/rptR", build_vignettes = TRUE)
library(rptR)
browseVignettes("rptR")
setwd("/Users/martin/Dropbox/projects/Rpackages/rptR")
md = read.table("ParalallusPolymorphData.txt", header=TRUE)#
#md = subset(md, Altitude >800) ##########################
require(lme4)#
require(lmerTest)#
require(rptR)#
require(colorspace)#
#
Alt = with(md, tapply(Altitude2, SiteID, mean, na.rm=TRUE))#
md$SiteSexID = paste(md$SiteID, md$Sex, sep="_")#
nboot = 1000#
#
# Variance components#
#
md$Trait = md$Pronotum#
summary(rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="F"), nboot=nboot))#
summary(rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="M"), nboot=nboot))#
md$Trait = md$Femur#
summary(rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="F"), nboot=nboot))#
summary(rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="M"), nboot=nboot))
md = read.table("ParalallusPolymorphData.txt", header=TRUE)
require(lme4)#
require(lmerTest)#
require(rptR)#
require(colorspace)
install.packages(c("lme4", "lmerTest", "rptR", "colorspace"))
Alt = with(md, tapply(Altitude2, SiteID, mean, na.rm=TRUE))#
md$SiteSexID = paste(md$SiteID, md$Sex, sep="_")#
nboot = 1000#
#
# Variance components#
#
md$Trait = md$Pronotum#
summary(rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="F"), nboot=nboot))#
summary(rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="M"), nboot=nboot))#
md$Trait = md$Femur#
summary(rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="F"), nboot=nboot))#
summary(rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="M"), nboot=nboot))
Alt = with(md, tapply(Altitude2, SiteID, mean, na.rm=TRUE))#
md$SiteSexID = paste(md$SiteID, md$Sex, sep="_")#
nboot = 1000
md$Trait = md$Pronotum
rptGaussian(Trait ~ I(Altitude/100) + (1|SiteID), grname="SiteID", data=subset(md, Sex=="F"), nboot=nboot)