Hola!!
Estoy trabajando con una función del tipo
h <- function(x) matrix(c(1,x[1],x[2]),ncol=1)
b0=c(1,2,1)
hb0 = function(x) sum(h(x) *b0)
hb0 <-Vectorize(hb0)
hb0(c(1,3)) #Si lo hace bien
Pero me gustaría que como argumento introducir una matiz de nx2, de tal
manera que mi función sea evaluada por renglón, es decir, que el resultado
sea un vector de n entradas. Tal como sucede al utilizar la función
"dmvnorm" de la librería "mvtnorm"
mu1=c(2,-2)
S1=4*diag(2)
dmvnorm(a,mu1,S1) #Esto si lo hace bien
#Sin embargo la función que yo declaré no lo hace
a=matrix(c(1,2,3,4,5,6,7,8),nrow=4,ncol=2)
hb0(a) #pero no evalúa por renglón
#Otra función que me gustaría optimizar es
b1=0.8
V=diag(2)
v<-function(x,y){ exp(-b1*t(x-y)%*%solve(V)%*%(x-y))}
v(c(1,3),c(2,1)) #De esta forma si funciona
De nuevo me gustaría que cada matriz la evalúe por renglón
v(a,-a)
Agradeceré mucho su ayuda. De antemano muchas gracias por tomarse el tiempo
de leer este correo.
[[alternative HTML version deleted]]