Displaying 2 results from an estimated 2 matches for "vproducto".
Did you mean:
producto
2013 Dec 16
2
vectorizaciones
...forma de funcion
producto <- function(n1a,n2a,Ka,vec1a,P1a,P2a){
LLa = numeric(2*Ka)
for (i in 1:n1a){
for (j in 1:n2a){
LLa = LLa - vec1a%*%matrix(abs(P1a[i,]-P2a[j,]),
nrow=length(vec1a), ncol=2*Ka, byrow=TRUE)
}
}
return(sum(LLa))
}
producto(n1a,n2a,Ka,vec1a,P1a,P2a)
vproducto <- Vectorize(producto)
vproducto(n1a,n2a,Ka,vec1a,P1a,P2a)
# Error en P1a[i, ] : número incorreto de dimensiones
#...................
Y saludos
mm~
_______________________________
J. Miguel Marin
http://www.est.uc3m.es/jmmarin
Dep. de Estadistica
Universidad Carlos III de M...
2013 Dec 16
0
vectorizaciones
...rix(rnorm(n.val),nrow=3)
P1a <- matrix(rnorm(n.val),nrow=3)
LLa = numeric(2*Ka)
for (i in 1:n1a){
for (j in 1:n2a){
LLa = LLa - vec1a%*%matrix(abs(P1a[i,]-P2a[j,]), nrow=length(vec1a),
ncol=2*Ka, byrow=TRUE)
}
}
return(sum(LLa))
}
producto(n1a,n2a,Ka,vec1a,n.val)
vproducto <- Vectorize(producto)
sum(vproducto(3,3,10,as.vector(1:3),30))
vproducto <- Vectorize(producto, c("n1a","n2a","Ka","vec1a","n.val"))
sum(vproducto(3,3,10,as.vector(1:3),30))
#-----------------------------------------------
Saludos,
Carlos...