Foro de debate

Duda

 
Imagen de David Tarin Bernad
Duda
de David Tarin Bernad - sábado, 11 de febrero de 2017, 21:31
 

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]