data(Pheno.list)
data(Ancestry.list)
## Full NAM analysis
Parent.list <- Parents$Genotype[-1]
names(Parent.list) <- Parents$Family[-1]
ETest.nam <- purrr::imap(Parent.list, ~ EpiTest.fit(Ancestry = Ancestry.list[[.y]],
Pheno = Pheno.list[[.y]],
ParentName = c(P0 = 'B73',P1 = .x),
Parents = Parents,
Trait = 'GDD_DTA',
Weight = TRUE))
## Variance component plot
EpiTest.plot(ETest.nam,Title = 'Days to anthesis',Type = "var")
## Directional epistasis plot
EpiTest.plot(ETest.nam,Title = 'Days to anthesis',Type = "mean",Alpha = 5/100)
Run the code above in your browser using DataLab