Tengo una duda, quiero coger numero por numero en una lista de 100 datos, os paso el script a ver si me podeis ayudar.Quiero que Rho, vaya cogiendo datos uno a uno para que me calcule W, ya que me da mal.
rm(list=ls())
r=matrix(c(0.095, 0.070, 0.090, 0.075), nrow=4)
x=matrix(c(0.038, 0.0085, 0.0089, 0.0066, 0.0085, 0.0331, 0.0156, 0.0039, 0.0089, 0.0156, 0.0334, 0.007, 0.0066, 0.0039, 0.007, 0.0240),nrow=4)
varcov=solve(x)
t<-matrix(c(1,1,1,1),nrow=4)
rho <- seq(0.05, 0.15, 0.001)
lrho=length(rho)
a=t(t)%*%solve(x)%*%t
ba=t(t)%*%solve(x)%*%r
c=t(r)%*%solve(x)%*%r
d=a%*%c-ba^2
g=varcov%*%((t%*%c-r%*%ba))%*%solve(d)
h=varcov%*%((r%*%a-t%*%ba))%*%solve(d)
i=seq(1,lrho)
w=g+h*rho[i]