Coverage for cvx/risk/factor/factor.py: 100%

1# Copyright 2023 Stanford University Convex Optimization Group 

2# 

3# Licensed under the Apache License, Version 2.0 (the "License"); 

4# you may not use this file except in compliance with the License. 

5# You may obtain a copy of the License at 

6# 

7# http://www.apache.org/licenses/LICENSE-2.0 

8# 

9# Unless required by applicable law or agreed to in writing, software 

10# distributed under the License is distributed on an "AS IS" BASIS, 

11# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. 

12# See the License for the specific language governing permissions and 

13# limitations under the License. 

14"""Factor risk model 

15""" 

16from __future__ import annotations 

17 

18from dataclasses import dataclass 

19 

20import cvxpy as cvx 

21import numpy as np 

22 

23from ..bounds import Bounds 

24from ..linalg import cholesky 

25from ..model import Model 

26 

27 

28@dataclass 

29class FactorModel(Model): 

30 """Factor risk model""" 

31 

32 assets: int = 0 

33 k: int = 0 

34 

35 def __post_init__(self): 

36 self.parameter["exposure"] = cvx.Parameter( 

37 shape=(self.k, self.assets), 

38 name="exposure", 

39 value=np.zeros((self.k, self.assets)), 

40 ) 

41 

42 self.parameter["idiosyncratic_risk"] = cvx.Parameter( 

43 shape=self.assets, name="idiosyncratic risk", value=np.zeros(self.assets) 

44 ) 

45 

46 self.parameter["chol"] = cvx.Parameter( 

47 shape=(self.k, self.k), 

48 name="cholesky of covariance", 

49 value=np.zeros((self.k, self.k)), 

50 ) 

51 

52 self.bounds_assets = Bounds(m=self.assets, name="assets") 

53 self.bounds_factors = Bounds(m=self.k, name="factors") 

54 

55 def estimate(self, weights, **kwargs): 

56 """ 

57 Compute the total variance 

58 """ 

59 var_residual = cvx.norm2( 

60 cvx.multiply(self.parameter["idiosyncratic_risk"], weights) 

61 ) 

62 

63 y = kwargs.get("y", self.parameter["exposure"] @ weights) 

64 

65 return cvx.norm2( 

66 cvx.vstack([cvx.norm2(self.parameter["chol"] @ y), var_residual]) 

67 ) 

68 

69 def update(self, **kwargs): 

70 exposure = kwargs["exposure"] 

71 k, assets = exposure.shape 

72 

73 self.parameter["exposure"].value[:k, :assets] = kwargs["exposure"] 

74 self.parameter["idiosyncratic_risk"].value[:assets] = kwargs[ 

75 "idiosyncratic_risk" 

76 ] 

77 self.parameter["chol"].value[:k, :k] = cholesky(kwargs["cov"]) 

78 self.bounds_assets.update(**kwargs) 

79 self.bounds_factors.update(**kwargs) 

80 

81 def constraints(self, weights, **kwargs): 

82 y = kwargs.get("y", self.parameter["exposure"] @ weights) 

83 

84 return ( 

85 self.bounds_assets.constraints(weights) 

86 + self.bounds_factors.constraints(y) 

87 + [y == self.parameter["exposure"] @ weights] 

88 )