Dj *dZddlZddlmZdZddZy)zZAdditional functions prediction standard errors and confidence intervals A: josef pktd N)statsctj|}|jdk(r |dddf}|S|jdk(rtj|}|S|jdkDr t d|S)zM convert array_like to 2d from 1d or 0d not tested because not used Nrztoo many dimensions)npasarrayndim atleast_2d ValueError)xs jC:\Crop_Prediction\Backend\crop-ai-system\venv\Lib\site-packages\statsmodels/sandbox/regression/predstd.py atleast_2dcolr sp 1 A ! agJ H &&A+ MM!  H &&1*.// HcJ|j}|;|jj}|j}||jj}nt j |}|jd|jdk7r td|jj|j|}|d}nJt j|}|jdkDr&t||jdk7r td|j|z |t j||j j zj#dz}t j$|}t&j(j+|dz |j,}|||zz} |||zz } || | fS)acalculate standard deviation and confidence interval for prediction applies to WLS and OLS, not to general GLS, that is independently but not identically distributed observations Parameters ---------- res : regression result instance results of WLS or OLS regression required attributes see notes exog : array_like (optional) exogenous variables for points to predict weights : scalar or array_like (optional) weights as defined for WLS (inverse of variance of observation) alpha : float (default: alpha = 0.05) confidence level for two-sided hypothesis Returns ------- predstd : array_like, 1d standard error of prediction same length as rows of exog interval_l, interval_u : array_like lower und upper confidence bounds Notes ----- The result instance needs to have at least the following res.model.predict() : predicted values or res.fittedvalues : values used in estimation res.cov_params() : covariance matrix of parameter estimates If exog is 1d, then it is interpreted as one observation, i.e. a row vector. testing status: not compared with other packages References ---------- Greene p.111 for OLS, extended to WLS by analogy rzwrong shape of exogg?rz+weights and exog do not have matching shapeg@) cov_paramsmodelexog fittedvaluesweightsrr shaper predictparamsrsizelen mse_residdotTsumsqrtrtisfdf_resid) resrralphacovb predictedpredvarpredstdtppf interval_u interval_ls r wls_prediction_stdr+s^^ >> D |yy~~$$ ?ii''G}}T" ::a=DJJqM )23 3II%%cjj$7 ?Gjj)G||aCLDJJqM$A !NOO mmG#tbffT466.B.D.D'D&I&I!&LLGgggG 77;;uRx .DTG^+JTG^+J J **r)NNg?)__doc__numpyrscipyrr r+rr r0s  J+r