Supongamos que tengo unas carteras que son el resultado de maximizar la utilidad esperada de diferentes aproximaciones de una función de utilidad, ¿cómo se prueban estas carteras fuera de muestra y cómo se comparan estos resultados?
Más concretamente, estoy tratando de comparar las carteras MV (media-varianza), MVS (media-varianza-estiramiento) y MVSK (media-varianza-estiramiento-curtosis) para un mismo conjunto de activos, basándome en la función de utilidad CARA. He maximizado estas tres ecuaciones y para cada una tengo una utilidad máxima esperada y la solución en forma de un conjunto de pesos. Así que ahora tengo tres carteras diferentes, y estas se basan en el mismo conjunto de muestra de 4 años (editar: 4 años de rendimientos diarios). ¿Puedo ahora probar estas carteras en un quinto año (un año fuera de la muestra, después del conjunto de muestra de cuatro años) y comparar algunas características? Si es así, ¿qué características debería comparar? ¿Puede comparar la utilidad alcanzada, o no es esta la intención correcta?
Edición: Algunas aclaraciones:
La aproximación de cuarto orden del valor esperado de la función de utilidad es la siguiente:
$E[U(W)] \approx U(\bar{W})+\frac{1}{2} U^{(2)}(\bar{W}) \sigma_{p}^{2}+\frac{1}{3 !} U^{(3)}(\bar{W}) s_{p}^{3}+\frac{1}{4 !} U^{(4)}(\bar{W}) \kappa_{p}^{4}$
donde $U(W)$ es la función de utilidad, $\bar{W}$ es la riqueza esperada, y $\sigma_{p}^{2}, s_{p}^{3}$ y $\kappa_{p}^{4}$ son la varianza, la asimetría y la curtosis de la cartera.
La función de utilidad que estoy utilizando es la función de utilidad CARA,
$U(W)=-\exp (-\gamma W)$ ,
donde $\gamma$ representa el coeficiente de aversión al riesgo constante.
El resultado es la siguiente ecuación,
$E[U(W)] \approx-\exp \left(-\gamma \mu_{p}\right)\left[1+\frac{\gamma^{2}}{2} \sigma_{p}+\frac{\gamma^{3}}{3 !} s_{p}^{3}+\frac{\gamma^{4}}{4 !} \kappa_{p}^{4}\right]$ ,
que es nuestra función objetivo a maximizar, en el caso de MVSK. Si consideramos MV, eliminamos los dos últimos términos del segundo operando de la última ecuación. En el caso de MVS, sólo eliminamos el último término.
Estos momentos están representados por $\begin{aligned} \mu_{p} &=\omega^{\prime} \mu \\ \sigma_{p}^{2} &=\omega^{\prime} \mathrm{M}_{2} \omega \\ s_{p}^{3} &=\omega^{\prime} \mathrm{M}_{3}(\omega \otimes \omega) \\ \kappa_{p}^{4} &=\omega^{\prime} \mathrm{M}_{4}(\omega \otimes \omega \otimes \omega) \end{aligned}$ donde $\mathrm{M}_{i}$ representa el tensor cumulante matricializado i, por lo que para $i = 2$ es la matriz de covarianza, para $i = 3$ el tensor de coscocidad matricializado y para $i = 4$ el tensor de cokurtosis matricializado. Con la matricización de un tensor, los cortes frontales del tensor se colocan lateralmente, por lo que para $i = n$ obtenemos un $n \times n^{(i-1)}$ matriz.
En cuanto a la optimización, puede resolverse utilizando múltiples intentos de SLSQP. Aunque no se garantiza un óptimo global, he comprobado que esto funciona bastante bien, ya que la mayoría de los intentos encuentran el mismo óptimo.
Ahora bien, encontremos o no un óptimo global para las optimizaciones MVS y MVSK, sigo sin entender cómo comparar mis carteras resultantes. Como mencionó Pontus Hultkrantz, la utilidad debería ser mi medida de rendimiento fuera de la muestra, pero ¿cómo puedo medir la utilidad realizada?
También debo señalar que los tensores de momento se están aproximando. Mi objetivo de investigación es demostrar que estas aproximaciones son lo suficientemente adecuadas para ser utilizadas en esta aplicación.
Gracias
0 votos
Si has utilizado la programación cuadrática para optimizar el modelo MV, ¿qué has utilizado para resolver los modelos MVS y MVSK?
0 votos
Utilicé la programación cuadrática por mínimos cuadrados secuenciales (SLSQP) de la función scipy.optimize.minimize de Python, para las tres optimizaciones.
2 votos
Pero la asimetría de la cartera no es convexa ni lineal porque es una función objetivo cúbica, no cuadrática. no se puede resolver con programación cuadrática, ya que de lo contrario se quedaría atascado en un óptimo local, no global, para los momentos superiores. ¿podría editar su pregunta mostrando las fórmulas que utilizó para los modelos MVS y MVSK CARA? sin ver lo que está haciendo, no se puede dar una respuesta adecuada en cuanto a la evaluación del rendimiento fuera de la muestra
0 votos
He añadido algunas aclaraciones. Por favor, indíqueme si debo proporcionar algunas más.
1 votos
Por qué la varianza, la asimetría y la curtosis reciben los mismos coeficientes de signo ( $+$ ¿Sabe que los inversores tienen una preferencia positiva por los momentos impares (media y asimetría) y negativa por los pares (varianza y curtosis)?
1 votos
Además, en su formulación, ¿está haciendo $\sigma_{p}^{2}, s_{p}^{3}$ y $\kappa_{p}^{4}$ todas las funciones de las ponderaciones de la cartera? ¿Cuáles son sus fórmulas? No veo las ponderaciones de la cartera en ninguna parte de su modelo
0 votos
He aclarado las fórmulas de los momentos de la cartera utilizados en la función objeto en mi post original. Soy consciente de que los inversores tienen preferencia positiva por los momentos impares y negativa por los pares, pero los coeficientes de signo han sido objeto de discusión. En Jondeau (2006) En la página de la web, se mencionan estas propiedades, así como la función del objeto (en la página 35), tal y como se describe en mi mensaje original.
0 votos
En un libro del mismo autor \textit{Financial Modeling Under Non-Gaussian Distributions} (2006), observo ahora que la función del objeto ha sido cambiada por la siguiente: $E[U(W)] \approx-\exp \left(-\lambda \mu_{p}\right)\left[1+\frac{\lambda^{2}}{2} \sigma_{p}^{2}-\frac{\lambda^{3}}{6} s_{p}^{3}+\frac{\lambda^{4}}{24} k_{p}^{4}\right]$ . ¿Cree que el primer documento ha cometido un error en la función del objeto?
1 votos
Gracias por las ediciones. Ahora veo dónde están las ponderaciones de la cartera y seguro que la asimetría de la cartera es una función cúbica, la curtosis: cuártica. No es cuadrática. Tu último comentario confirma las preferencias del momento. Ahora es una cuestión de si maximizan o minimizan el conjunto. ¿Qué piensas? ¿Y mencionaron qué algoritmos de optimización utilizaron para cada modelo?
1 votos
Maximizan la función objeto como se ha dicho, pero en la práctica minimizan (porque los optimizadores minimizan) la negación de la función objeto. En la bibliografía, Jondeau et al. sólo mencionan que el problema puede resolverse "con un paquete de optimización estándar".
0 votos
Hola a todos. ¿Tenéis código MATLAB y/o R para problemas de optimización de MVS y/o MVSK?