Advertisement
Guest User

Untitled

a guest
Jun 26th, 2019
59
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.17 KB | None | 0 0
  1. data1={{2014.,0.015},{2015.,0.005},{2016.,0.0},{2017.,0.01},{2018.,0.02},{2019.,0.014}};
  2.  
  3. ListPlot[data1
  4. ,Frame->True
  5. ,PlotRange->{{2013.,2022.},{-0.01,0.03}}
  6. ,PlotStyle->Directive[Orange,PointSize[Large]]
  7. ]
  8.  
  9. Mean[data1]
  10. StandardDeviation[data1]
  11.  
  12. data2=Standardize[data1];
  13.  
  14. ListPlot[data2
  15. ,Frame->True
  16. ,PlotRange->{{-2.,2.},{-2.,2.}}
  17. ,PlotStyle->Directive[Orange,PointSize[Large]]
  18. ]
  19.  
  20. lmFit[data_List,degree_Integer]:=LinearModelFit[data,Table[x^i,{i,degree}],x]
  21.  
  22. lmFitPlot[data_List,degree_Integer,{xmin_,xmax_,ymin_,ymax_}]:=Module[{lmf,ss},
  23.  
  24. lmf=lmFit[data,degree];
  25. ss=Total[lmf["FitResiduals"]^2]; (* Sum of squared residuals *)
  26.  
  27. Show[
  28. {Plot[lmf[x],{x,xmin,xmax}
  29. ,PlotRange->{{xmin,xmax},{ymin,ymax}}]
  30. ,ListPlot[data,PlotStyle->Directive[Orange,PointSize[Large]]]
  31. }
  32. ,Frame->True
  33. ,FrameLabel->{{"",""},{"Year",Row[{"Sum squared residuals= ",ss}]}}
  34. ,ImageSize->Medium
  35. ]
  36. ]
  37.  
  38. lmFitPlot[data1,5,{2013.,2022.,-0.01,0.03}]
  39.  
  40. lmFitPlot[data2,5,{-2.,2.,-2.,2.}]
  41.  
  42. MatrixRank[lmFit[data1,5]["DesignMatrix"]]
  43. MatrixRank[lmFit[data2,5]["DesignMatrix"]]
  44.  
  45. FindGeometricTransform[data1,data2]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement