Hi, I have an equation Xcp and would like to look at errors affecting it. I
am applying errors of -.58 to .1 to R and -.45 to .47 to Xc. I would like
to look at all combinations. I set up a matrix function as follows and it
does not work. On the last line, when it run, it does not seem like the two
X’s are not the same.
Example: it is NOT 1+(162/1), but using another “R” value.
R<-c(1:5)
Xc<-c(16:20)
Xcp<-R+((Xc^2)/R)
test=R+(Xc^2/R)
x=c(-.58, -.44, -.31, -.17, -.04, .1)
y=c(-.45,-.27,-.08,.10,.29,.47)
actual=Xcp
X=((R%o%x+R))
Y=(Xc%o%y+Xc)
num.x.col <- length(X[1,])
num.y.col <- length(Y[1,])
num.rows <- length(X[,1])
Z <- matrix(nrow=num.rows, ncol=num.x.col*num.y.col)
for( i in 1:num.rows) {
Z[i,] <- as.vector(X[i,] + ((Y[i,])^2 %*% t(X[i,])^-1 ))
}
--
M. Keith Cox, Ph.D.
Alaska NOAA Fisheries, National Marine Fisheries Service
Auke Bay Laboratories
17109 Pt. Lena Loop Rd.
Juneau, AK 99801
Keith.Cox@noaa.gov
marlinkcox@gmail.com
U.S. (907) 789-6603
[[alternative HTML version deleted]]
