Estoy leyendo el papel por Zhao et al (2008) y tiene un problema con las definiciones utilizadas en el texto en la página 1535.
En primer lugar, generamos una muestra, RR, de un tamaño dado de la distribución (21). Deje que ˆμ^μ y ˆσ2^σ2 ser la media de la muestra del vector y la varianza de la muestra matriz de covarianza. Entonces, podemos modificar el generado de la muestra a ˆR=μ+(R−ˆμ)ˆσ−1σ.^R=μ+(R−^μ)^σ−1σ. A continuación, la muestra modificada ˆR^R tiene el mismo primer y segundo momentos de la distribución original. A fin de probar, si la muestra modificada es el arbitraje libre, usamos el Matlab barra diagonal inversa de la función de examinar si la solución a (1+R)⊤\barrainvertida1(1+R)⊤\barrainvertida1 es de las componentes positivas.
Pregunta. Lo que significa 11 en la última línea? Es una identidad de la matriz o un vector columna de unos? Y lo que es de dimensiones de este 11?
He tratado de examinar la solución de la Tabla 1.
library(pracma)
n <- 5
R<-matrix(
c(
0.0025, 0.0377, 0.0110, 0.0769, 0.0047,
0.0025, 0.0431, 0.0001, 0.0045, 0.0562,
0.0025, 0.0469, 0.0643, 0.0400, 0.0370,
0.0025, 0.0504, 0.0422, 0.0169, 0.0333,
0.0025, 0.0596, 0.0038, 0.1896, 0.0663), ncol=n)
#I<-ones(n)
I <-diag(n)
mldivide(t(I+R),I)
Agregar después de JejeBelfort del comentario. Si el 11 es un vector columna de unos, a continuación, lo que es de ++? De la unión de operación o producto de Kronecker de operador?
I <- rep(1, n)
cbind(I, R)
De referencia.
Yonggan Zhao, William T. Ziemba (2008) el Cálculo de riesgo neutral probabilidades y óptima de la cartera de pólizas en un dinámico modelo de inversión con el inconveniente del riesgo de control. European Journal of Operational Research 185 (2008) 1525-1540.