Displaying 1 result from an estimated 1 matches for "meanxj".
Did you mean:
meanx
2006 Jul 11
1
test regression against given slope for reduced major axis regression (RMA)
...r(i in 1:n){
sumXjYj <- sumXjYj + XjArray[i] * YjArray[i]
sumXj <- sumXj + XjArray[i]
sumYj <- sumYj + YjArray[i]
sumXjSquared <- sumXjSquared + XjArray[i]^2
sumYjSquared <- sumYjSquared + YjArray[i]^2
}
## Mean of X and Y vectors
meanyj <- sumYj / n
meanxj <- sumXj / n
## Create covariance
for(i in 1:n){
covariancexy <- covariancexy + ((XjArray[i] - meanxj) * (YjArray[i] - meanyj))
}
covariancexy <- covariancexy / n
## get variance of X and Y (SD)
varXj <- (n * sumXjSquared-sumXj^2)/(n*(n - 1))
varYj <- (n *...