14 Questions 3 Answers 0 Followers
Questions related from Moshood Bakare
A linear mixed model was fitted using lmer function of lme4 package in R. Does anyone have idea of how to extract studentized conditional residual for individual data point?
01 January 2020 2,139 3 View
Dear All, I am using ASREML-R to fit unstructured (UN) and factor analytic (FA) model to explore complex structure of genotype by environment interaction in multienvironment yield data. There are...
11 November 2019 781 4 View
I need to filter a vcf file for genotype quality and depth. How do I choose my filter parameter value? vcftools --gzvcf snps.txt.gz --quality ?? --depth ?? --recode --recode-INFO-all --out...
11 November 2018 3,242 0 View
A genotype data for Genome Wide Association Analysis (GWAS) is provided as 0, 1, and 2 format where individuals along the row and each column is a SNP. The minor allele frequency (MAF) is...
04 April 2018 9,638 3 View
Assuming 12 varieties are evaluated in 4 replicates in Randomized Complete Block Design (RCBD) where 20 variables are measured, is it appropriate to run correlation and principal component...
05 May 2016 2,641 3 View
These statistics measures association between ordinal variables: gamma, Kendall’s tau-, Stuart’s tau-, and Somers’. If I am to report correlation between ordinal and nominal variables, what is...
05 May 2016 3,033 2 View
A GSTAT library of R was used to generate a grid size of (i) 10 x 10 m and (ii) 10 x 5 m running from an Easting of 299677 m to 301297 m and a Northing of 5737278 m to 5738128 m. This resulted in...
11 November 2014 413 6 View
A number of variogram models were fitted where estimates like nugget, range, and sill were obtained from gstat of R. Could anyone guide me on how to obtain standard error for each of these...
04 April 2014 7,152 2 View
I used gstat to fit a Matern variogram model to my crop yield by setting the kappa = 0.5. To my surprise, its parameter estimates are exactly the same with that of exponential. Does it mean that...
04 April 2014 7,422 3 View
The four variogram models below are fitted to the same dataset. I got a warning message while cross validating a fitted linear variogram model as Warning message: In sqrt(ret[[var.name]]) : NaNs...
04 April 2014 5,696 2 View
I am surprised to get different results for skewness and kurtosis from different application or library. Take for example, with moments library package in R we have the following...
03 March 2014 3,188 14 View
An exponential model is fitted to empirical semivariogram using gstat of R package. The variogram parameters obtained is model psill range 1 Nug 0.04281109 0.00000 2 Exp 0.76810071...
03 March 2014 219 7 View
I am using R version of 2.15.3 in R studio environment. I used a function "skew" but got an error message that the name function cannot be found.
03 March 2014 2,927 7 View
I tried to compare the result obtained by fitting the spherical variogram model using fit.variogram and nls functions. The large difference in the results is a great concern for me knowing that...
09 September 2013 7,875 1 View
