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, $R$, de un tamaño dado de la distribución (21). Deje que $\hat{\mu}$ y $\hat{\sigma}^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 $$\hat{R}=\mu + (R-\hat{\mu})\hat{\sigma}^{-1}\sigma.$$ A continuación, la muestra modificada $\hat{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)^{\top} \barra invertida 1$ es de las componentes positivas.
Pregunta. Lo que significa $1$ 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 $1$?
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 $1$ 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.