R Script for RBD, Genetic Variability and Correlation - Path Analysis

install.packages("variability") 
library(variability) 

#################Importing data into R############# 
traits <- read.table("clipboard", h = TRUE) 
attach(traits) 
str(traits) 

library(readxl) 
traits <- read_excel("traits.xlsx") 
View(traits) 

#################Estimation of genetic parameters############# 
genvar <- gen.var(traits[3:13], traits$Genotypes, traits$Replication) 
genvar 

sink("genvar.txt") 
print(genvar) 

#################Genotypic and phenotypic correlations############# 
gencor <- geno.corr(traits[3:13], traits$Genotypes, traits$Replication) 
gencor 

sink("gencor.txt") 
print(gencor) 

phecor <- pheno.corr(traits[3:13], traits$Genotypes, traits$Replication) 
phecor 

sink("phecor.txt") 
print(phecor) 

#################Genotypic and phenotypic path coefficients############# 
genpath <- geno.path(traits[13], traits[3:12], traits$Genotypes, traits$Replication) 
genpath 

sink("genpath.txt") 
print(genpath) 

phepath <- pheno.path(traits[13], traits[3:12], traits$Genotypes, traits$Replication) 
phepath 

sink("phepath.txt") 
print(phepath)