PSTricks---plotstyle LSM(最小二乘法)

PSTricks---plotstyle LSM(最小二乘法)

我正在探索 PSTricksplotstyle=LSM来进行回归分析y=m x + c。这是我的代码:

\documentclass[12pt]{article}
\usepackage{filecontents}
\begin{filecontents*}{data.dat}
0.0     0.0
144.321  0.84
159.407     0.8925

\end{filecontents*}

\usepackage{pst-plot,pstricks-add}

\begin{document}
\readdata{\data}{data.dat}

\psset{xunit=0.6mm,yunit=30mm,xlabelFactor=$$,labelFontSize=\footnotesize,mathLabel=false}
\begin{pspicture}(0,0)(150,3)
\psaxes[axesstyle=frame,Dx=50,Dy=0.5, ticksize=0 4pt](150,3)
\listplot [xStart=0, xEnd=150,xunit=1,plotstyle=values,rot=90] {\data}
\listplot [xStart=0,xEnd=150,xunit=1,linewidth=1.5pt,linecolor=red,showpoints, dotstyle=square] {\data}
\listplot[PstDebug=1,xStart=0, xEnd=150,xunit=1,plotstyle=LSM,linecolor=blue, linestyle=dashed]{\data}
\end{pspicture}

\end{document}

该参数xEnd=150使得线性回归分析仅考虑数据文件的前两行。

经过编译,我得到线方程为

y = 0.291017 x + 9.64966e-08

我发现的梯度0.291017似乎是错误的,因为梯度应该是0.84/144.321 = 0.0058...

我在这里做错了什么?

答案1

这是不同轴单元的问题。尝试使用此版本:

http://texnik.dante.de/tex/generic/pst-plot/pst-plot.tex

有了这个文件,我就可以正常工作。

顺便说一下,现在有一个 pst-fit,见http://texnik.dante.de/doc/generic/pst-fit/pst-fit-doc.pdf

相关内容