home *** CD-ROM | disk | FTP | other *** search
/ The World of Computer Software / World_Of_Computer_Software-02-385-Vol-1of3.iso / m / mcademo2.zip / CURVEFIT.MCD / CURVEFIT.MCD
Text File  |  1991-06-26  |  6KB  |  265 lines
  1. .MCD 30001 1 74
  2. .CMD PLOTFORMAT
  3. 0 0 1 0 0
  4. 0 0 1 0 0
  5. 0 1 0 0 NO-TRACE-STRING
  6. 0 2 1 0 NO-TRACE-STRING
  7. 0 3 2 0 NO-TRACE-STRING
  8. 0 4 3 0 NO-TRACE-STRING
  9. 0 1 4 0 NO-TRACE-STRING
  10. 0 2 5 0 NO-TRACE-STRING
  11. 0 3 6 0 NO-TRACE-STRING
  12. 0 4 0 0 NO-TRACE-STRING
  13. 0 1 1 0 NO-TRACE-STRING
  14. 0 2 2 0 NO-TRACE-STRING
  15. 0 3 3 0 NO-TRACE-STRING
  16. 0 4 4 0 NO-TRACE-STRING
  17. 0 1 5 0 NO-TRACE-STRING
  18. 0 2 6 0 NO-TRACE-STRING
  19. 0 3 0 0 NO-TRACE-STRING
  20. 0 4 1 0 NO-TRACE-STRING
  21. 0 1 21 15
  22.  
  23. .CMD FORMAT  rd=d ct=10 im=i et=3 zt=15 pr=3 mass length time charge
  24. .CMD SET ORIGIN 0
  25. .CMD SET TOL 0.001000000000000
  26. .CMD SET PRNCOLWIDTH 8
  27. .CMD SET PRNPRECISION 4
  28. .CMD PRINT_SETUP 1.200000 0
  29. .CMD DEFINE_FONTSTYLE_NAME fontID=0 name=Variables
  30. .CMD DEFINE_FONTSTYLE_NAME fontID=1 name=Constants
  31. .CMD DEFINE_FONTSTYLE_NAME fontID=2 name=Text
  32. .CMD DEFINE_FONTSTYLE_NAME fontID=3 name=Greek^Variables
  33. .CMD DEFINE_FONTSTYLE_NAME fontID=4 name=User^1
  34. .CMD DEFINE_FONTSTYLE_NAME fontID=5 name=User^2
  35. .CMD DEFINE_FONTSTYLE_NAME fontID=6 name=User^3
  36. .CMD DEFINE_FONTSTYLE_NAME fontID=7 name=User^4
  37. .CMD DEFINE_FONTSTYLE_NAME fontID=8 name=User^5
  38. .CMD DEFINE_FONTSTYLE_NAME fontID=9 name=User^6
  39. .CMD DEFINE_FONTSTYLE fontID=0 family=Tms^Rmn points=10 bold=0 italic=0 underline=0
  40. .CMD DEFINE_FONTSTYLE fontID=1 family=Tms^Rmn points=10 bold=0 italic=0 underline=0
  41. .CMD DEFINE_FONTSTYLE fontID=2 family=Helv points=10 bold=0 italic=0 underline=0
  42. .CMD DEFINE_FONTSTYLE fontID=3 family=Symbol points=10 bold=0 italic=0 underline=0
  43. .CMD DEFINE_FONTSTYLE fontID=4 family=Helv points=10 bold=0 italic=0 underline=0
  44. .CMD DEFINE_FONTSTYLE fontID=5 family=Courier points=10 bold=0 italic=0 underline=0
  45. .CMD DEFINE_FONTSTYLE fontID=6 family=System points=10 bold=0 italic=0 underline=0
  46. .CMD DEFINE_FONTSTYLE fontID=7 family=Script points=10 bold=0 italic=0 underline=0
  47. .CMD DEFINE_FONTSTYLE fontID=8 family=Terminal points=0 bold=0 italic=0 underline=0
  48. .CMD DEFINE_FONTSTYLE fontID=9 family=Modern points=10 bold=0 italic=0 underline=0
  49. .CMD UNITS U=1
  50. .TXT 2 1 0 0
  51. C a264,330,77
  52. {\rtf1\ansi \deff0
  53. {\fonttbl
  54. {\f0\fnil Helv;}
  55. }
  56. {\plain {}{\f0 \fs24 \b \i \ulnone }{}{\f0 
  57. \fs24 \b \i \ulnone POLYNOMIAL CURVE FITTING}{}}
  58. }
  59. .TXT 4 0 0 0
  60. C a547,574,218
  61. {\rtf1\ansi \deff0
  62. {\fonttbl
  63. {\f0\fnil Helv;}
  64. }
  65. {\plain {}{This document shows how to fit
  66.  a quadratic function to a set of data. 
  67. The same technique also works for other 
  68. types of curve fitting and for multiple 
  69. regression (regression with several independent
  70.  variables).}{}}
  71. }
  72. .TXT 8 0 0 0
  73. C a568,570,167
  74. {\rtf1\ansi \deff0
  75. {\fonttbl
  76. {\f0\fnil Helv;}
  77. }
  78. {\plain {}{(This document follows the presentation
  79.  in "Applied Linear Statistical Models",
  80.  by John Neter and William Wasserman, Richard
  81.  D. Irwin, 1974. Chapters 7 and 8.)}{}}
  82. }
  83. .TXT 5 0 0 0
  84. C a328,330,48
  85. {\rtf1\ansi \deff0
  86. {\fonttbl
  87. {\f0\fnil Helv;}
  88. }
  89. {\plain {}{First, read the data from external
  90.  files:}{}}
  91. }
  92. .EQN 4 0 0 0
  93. x:READPRN(data1)
  94. .EQN 0 18 0 0
  95. N:length(x)
  96. .EQN 0 14 0 0
  97. y:READPRN(data2)
  98. .EQN 0 18 0 0
  99. N=?_n_u_l_l_
  100. .TXT 4 -50 0 0
  101. C a208,210,33
  102. {\rtf1\ansi \deff0
  103. {\fonttbl
  104. {\f0\fnil Helv;}
  105. }
  106. {\plain {}{Compute sample statistics:}{}}
  107. }
  108. .EQN 4 0 0 0
  109. mean(x)=?_n_u_l_l_
  110. .EQN 0 14 0 0
  111. mean(y)=?_n_u_l_l_
  112. .EQN 0 18 0 0
  113. var(x)=?_n_u_l_l_
  114. .EQN 1 17 0 0
  115. var(y)=?_n_u_l_l_
  116. .EQN 3 -43 0 0
  117. stdev(x)=?_n_u_l_l_
  118. .EQN 0 21 0 0
  119. stdev(y)=?_n_u_l_l_
  120. .TXT 4 -27 0 0
  121. C a192,194,31
  122. {\rtf1\ansi \deff0
  123. {\fonttbl
  124. {\f0\fnil Helv;}
  125. }
  126. {\plain {}{Simple linear curve fit:}{}}
  127. }
  128. .EQN 4 0 0 0
  129. corr(x,y)=?_n_u_l_l_
  130. .EQN 0 17 0 0
  131. m:slope(x,y)
  132. .EQN 0 15 0 0
  133. b:intercept(x,y)
  134. .EQN 0 16 0 0
  135. linear(x):m*x+b
  136. .TXT 5 -48 0 0
  137. C a216,218,34
  138. {\rtf1\ansi \deff0
  139. {\fonttbl
  140. {\f0\fnil Helv;}
  141. }
  142. {\plain {}{Compute mean squared error:}{}}
  143. }
  144. .EQN 5 0 0 0
  145. SSE.L:{55}((((y-linear(x)))^(2))){49}
  146. .EQN 0 27 0 0
  147. MSE.L:(SSE.L)/(N-2)
  148. .EQN 0 18 0 0
  149. MSE.L=?_n_u_l_l_
  150. .TXT 5 -45 0 0
  151. C a296,298,44
  152. {\rtf1\ansi \deff0
  153. {\fonttbl
  154. {\f0\fnil Helv;}
  155. }
  156. {\plain {}{Quadratic fit using matrix operations}{}}
  157. }
  158. .TXT 4 0 0 0
  159. C a264,266,40
  160. {\rtf1\ansi \deff0
  161. {\fonttbl
  162. {\f0\fnil Helv;}
  163. }
  164. {\plain {}{Create second variable: x squared}{}}
  165. }
  166. .EQN 0 29 0 0
  167. x2:(((x)^(2))){49}
  168. .EQN 0 13 0 0
  169. i:0;N-1
  170. .EQN 4 -28 0 0
  171. (X)[(i,0):1
  172. .TXT 1 -14 0 0
  173. C a120,122,22
  174. {\rtf1\ansi \deff0
  175. {\fonttbl
  176. {\f0\fnil Helv;}
  177. }
  178. {\plain {}{Create X matrix}{}}
  179. }
  180. .EQN 0 25 0 0
  181. (X){52}(1):x
  182. .EQN 0 15 0 0
  183. (X){52}(2):x2
  184. .EQN 6 -33 0 0
  185. b:(((X){51}*X))^(-1)*((X){51}*y)
  186. .EQN 3 24 0 0
  187. b=?_n_u_l_l_
  188. .TXT 8 -31 0 0
  189. C a104,106,20
  190. {\rtf1\ansi \deff0
  191. {\fonttbl
  192. {\f0\fnil Helv;}
  193. }
  194. {\plain {}{Fitted curve:}{}}
  195. }
  196. .EQN 0 14 0 0
  197. quad(x):(b)[(0)+(b)[(1)*x+(b)[(2)*(x)^(2)
  198. .TXT 4 -14 0 0
  199. C a216,218,34
  200. {\rtf1\ansi \deff0
  201. {\fonttbl
  202. {\f0\fnil Helv;}
  203. }
  204. {\plain {}{Compute mean squared error:}{}}
  205. }
  206. .EQN 5 0 0 0
  207. SSE.Q:{55}((((y-quad(x)))^(2))){49}
  208. .EQN 1 25 0 0
  209. MSE.Q:(SSE.Q)/(N-3)
  210. .EQN 0 17 0 0
  211. MSE.Q=?_n_u_l_l_
  212. .TXT 4 -42 0 0
  213. C a336,338,49
  214. {\rtf1\ansi \deff0
  215. {\fonttbl
  216. {\f0\fnil Helv;}
  217. }
  218. {\plain {}{Now graph the two curves against
  219.  the data:}{}}
  220. }
  221. .EQN 1 0 0 0
  222. 7*(10)^(5)&0&(y)[(i),linear((x)[(i))@&&(x)[(i)
  223. 0 0 0 1 0
  224. 0 0 0 1 0
  225. 0 1 0 1 NO-TRACE-STRING
  226. 0 1 1 5 NO-TRACE-STRING
  227. 0 1 2 5 NO-TRACE-STRING
  228. 0 4 3 0 NO-TRACE-STRING
  229. 0 1 4 0 NO-TRACE-STRING
  230. 0 2 5 0 NO-TRACE-STRING
  231. 0 3 6 0 NO-TRACE-STRING
  232. 0 4 0 0 NO-TRACE-STRING
  233. 0 1 1 0 NO-TRACE-STRING
  234. 0 2 2 0 NO-TRACE-STRING
  235. 0 3 3 0 NO-TRACE-STRING
  236. 0 4 4 0 NO-TRACE-STRING
  237. 0 1 5 0 NO-TRACE-STRING
  238. 0 2 6 0 NO-TRACE-STRING
  239. 0 3 0 0 NO-TRACE-STRING
  240. 0 4 1 0 NO-TRACE-STRING
  241. 0 1 15 16
  242.  
  243. .EQN 0 38 0 0
  244. 7*(10)^(5)&0&(y)[(i),quad((x)[(i))@&&(x)[(i)
  245. 0 0 0 1 0
  246. 0 0 0 1 0
  247. 0 1 0 1 NO-TRACE-STRING
  248. 0 1 1 5 NO-TRACE-STRING
  249. 0 1 2 5 NO-TRACE-STRING
  250. 0 4 3 0 NO-TRACE-STRING
  251. 0 1 4 0 NO-TRACE-STRING
  252. 0 2 5 0 NO-TRACE-STRING
  253. 0 3 6 0 NO-TRACE-STRING
  254. 0 4 0 0 NO-TRACE-STRING
  255. 0 1 1 0 NO-TRACE-STRING
  256. 0 2 2 0 NO-TRACE-STRING
  257. 0 3 3 0 NO-TRACE-STRING
  258. 0 4 4 0 NO-TRACE-STRING
  259. 0 1 5 0 NO-TRACE-STRING
  260. 0 2 6 0 NO-TRACE-STRING
  261. 0 3 0 0 NO-TRACE-STRING
  262. 0 4 1 0 NO-TRACE-STRING
  263. 0 1 15 16
  264.  
  265.