我正在探索 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