Ad Code

R Script for Diallel Analysis


install.packages("DiallelAnalysisR") 
library(DiallelAnalysisR) 
library(readxl) 

Griffing1Data2 <- read_excel("Griffing1Data2.xlsx") 
View(Griffing1Data2) 

***Diallel Analysis with Griffing's Aproach Method 2 & Model 1 

Griffing1Data2 <- Griffing(
  y = Yield,
  Rep = Rep,
  Cross1 = Cross1,
  Cross2 = Cross2,
  data = GriffingData2,
  Method = 2,
  Model = 1
) 

names(Griffing1Data2) 
Griffing1Data2 

Griffing1Data2Means <- Griffing1Data2$Means 
Griffing1Data2ANOVA <- Griffing1Data2$ANOVA 
Griffing1Data2Genetic.Components <- Griffing1Data2$Genetic.Components 
Griffing1Data2Effects <- Griffing1Data2$Effects 
Griffing1Data2StdErr <- as.matrix(Griffing1Data2$StdErr) 

sink("Griffing1Data2.txt") 
Griffing(
  y = Yield,
  Rep = Rep,
  Cross1 = Cross1,
  Cross2 = Cross2,
  data = GriffingData2,
  Method = 2,
  Model = 1
) 
names(Griffing1Data2) 
Griffing1Data2 

Griffing1Data2Means <- Griffing1Data2$Means 
Griffing1Data2ANOVA <- Griffing1Data2$ANOVA 
Griffing1Data2Genetic.Components <- Griffing1Data2$Genetic.Components 
Griffing1Data2Effects <- Griffing1Data2$Effects 
Griffing1Data2StdErr <- as.matrix(Griffing1Data2$StdErr) 
sink() 

****Diallel Analysis with Griffing's Aproach Method 2 & Model 2 

Griffing2Data2 <- Griffing(
  y = Yield,
  Rep = Rep,
  Cross1 = Cross1,
  Cross2 = Cross2,
  data = GriffingData2,
  Method = 2,
  Model = 2
) 

names(Griffing2Data2) 
Griffing2Data2 

Griffing2Data2Means <- Griffing2Data2$Means 
Griffing2Data2ANOVA <- Griffing2Data2$ANOVA 
Griffing2Data2Genetic.Components <- Griffing2Data2$Genetic.Components 

sink("Griffing2Data2.txt") 
Griffing2Data2 <- Griffing(
  y = Yield,
  Rep = Rep,
  Cross1 = Cross1,
  Cross2 = Cross2,
  data = GriffingData2,
  Method = 2,
  Model = 2
) 

names(Griffing2Data2) 
Griffing2Data2 

Griffing2Data2Means <- Griffing2Data2$Means 
Griffing2Data2ANOVA <- Griffing2Data2$ANOVA 
Griffing2Data2Genetic.Components <- Griffing2Data2$Genetic.Components 
sink()

Post a Comment

0 Comments

Close Menu