search

cran / sommer

:exclamation: This is a read-only mirror of the CRAN R package repository. sommer — Solving Mixed Model Equations in R. Homepage: https://github.com/covaruber/sommer Report bugs for this package: https://github.com/covaruber/sommer/issues
15 stars 8 forks source link

Error in mmer(Y = y, Z = ETA, method = "EMMA") : unused argument (Y = y, Z=ETA) #3

Open jeremysutherland opened 4 years ago

jeremysutherland commented 4 years ago

library(vroom) library(textshape) library(plyr) library(dplyr) library(sommer) #Version 4.1

Sys.setenv("VROOM_CONNECTION_SIZE" = "13200000") options(expressions = 5e5)

READ IN TABLES

pheno <- read.table("diversity.csv", sep = ",", header = TRUE) geno <- vroom("Y_filt.csv", col_names = TRUE)

CLEAN UP TABLES

geno <- textshape::column_to_rownames(geno, loc = 1) pheno <- textshape::column_to_rownames(pheno, loc = 1) geno <- subset(geno, rownames(geno) %in% rownames(pheno)) pheno <- subset(pheno, rownames(pheno) %in% rownames(geno))

SORT TABLES

geno <- geno[ order(row.names(geno)), ] pheno <- pheno[ order(row.names(pheno)), ]

Narrow Sense Heritability

y <- pheno$Shannon Za <- diag(length(y)) Zd <- diag(length(y)) Ze <- diag(length(y)) A <- A.mat(geno) D <- D.mat(geno) E <- E.mat(geno)

ETA <- list(add=list(Z=Za,K=A),dom=list(Z=Zd,K=D),epi=list(Z=Ze,K=E))

ans.ADE <- mmer(Y=y, Z=ETA, method = "EMMA")

Here is where the error comes in. It simply won't accept y and ETA as inputs.

gaborcsardi commented 4 years ago

Hi, this is a read only mirror of CRAN, please see the package authors in the DESCRIPTION file. Thanks!