IUBio GIL .. BIOSCI/Bionet News .. Biosequences .. Software .. FTP

Simple graph reading software

Jari Oksanen jari.oksanen at helsinki.fi
Fri Apr 7 00:30:40 EST 2000

Kurt Sims wrote:
> I'm looking for a nice simple program that will let me plot the points
> for a scattergram, give it a best fit line and then read an unknown.  I
> would love to hear some suggestions.
> All I need it to do is accept the standards for the X axis and Y axis
> and then give me the reading when I give it an unknown reading off of
> one of the axis.
> The program can be either PC or Mac based, or even on the web.

Any decent statistical program should be able to do that. In R (or
S-plus) you can do like this:

> conc <- c(0.5,1,1.5,2)
> reading <- c(15.2,30.2,44.6,59)
> standard <- lm(conc~reading)
> summary(standard)

lm(formula = conc ~ reading)

        1         2         3         4 
 0.006096 -0.008255 -0.002032  0.004191 

              Estimate Std. Error t value Pr(>|t|)    
(Intercept) -0.0273051  0.0099380  -2.748    0.111    
reading      0.0342901  0.0002444 140.296 5.08e-05 ***
Signif. codes:  0  `***'  0.001  `**'  0.01  `*'  0.05  `.'  0.1  ` '  1 

Residual standard error: 0.007969 on 2 degrees of freedom
Multiple R-Squared: 0.9999,     Adjusted R-squared: 0.9998 
F-statistic: 1.968e+04 on 1 and 2 degrees of freedom,   p-value:

> predict(standard,data.frame(reading=38.3))
[1] 1.286005

That is, for reading=38.3 the `standard' predicts a concentration of

With the following command, you get more information than you want
(i.e., the confidence intervals of your prediction: usually they are
wider than you want to believe):

> predict(standard,data.frame(reading=38.3),interval="prediction")
          fit      lwr      upr
[1,] 1.286005 1.247655 1.324354

Yes, of course you can omit the intercept and and get the standard line
go through the origin. And of course, you can get all your readings in
one run, just put them in a vector in a data.frame.

cheers, j.oksanen
J.Oksanen -- UHLA (University of Helsinki at Lahti)
C-h i m emacs RET m text RET m formatted RET
"GNU Emacs is evil." Linus Torvalds

More information about the Bio-soft mailing list

Send comments to us at archive@iubioarchive.bio.net