So I opened issue #55 and I tried using the array, but I ran into a data class error.
fit_MRE = metagenomeSeq::fitZig(obj = n_paired_merged_MRE, mod = mod, useCSSoffset = FALSE, control = zigControl())dupCor_MRE = duplicateCorrelation(fit_MRE$fit, mod, block = pData(n_paired_merged_MRE)$id)
I'm assuming I have to use one of these values as the object for duplicateCorrelation, but I'm not sure if I would use one of these components or use the fitted values matrix outputted by fitted(fitMRE$fit):
coefficients
stdev.unscaled
sigma
df.residual
cov.coefficients
pivot
rank
Amean
method
design
df
Hi,
So I opened issue #55 and I tried using the array, but I ran into a data class error.
fit_MRE = metagenomeSeq::fitZig(obj = n_paired_merged_MRE, mod = mod, useCSSoffset = FALSE, control = zigControl())
dupCor_MRE = duplicateCorrelation(fit_MRE$fit, mod, block = pData(n_paired_merged_MRE)$id)
I'm assuming I have to use one of these values as the object for duplicateCorrelation, but I'm not sure if I would use one of these components or use the fitted values matrix outputted by
fitted(fitMRE$fit)
: coefficients stdev.unscaled sigma df.residual cov.coefficients pivot rank Amean method design df