miércoles, 14 de agosto de 2013

Software para ajuste de rectas

La página principal del R está acá. Para probar, un archivo con los datos en formato
x y ex ey
y el programa para correr en R.

Desde una sesión de R hacer:

md <- read.table( "pru.dat" )

fit <- deming( md$V1, md$V2, md$V3, md$V4 )

print( fit )

y salen los parámetros de la recta con su error:

Call:

deming(x = md$V1, y = md$V2, xstd = md$V3, ystd = md$V4)

Coef se(coef) z p

Intercept -0.02009639 0.08497616 -61.6346400 0.0000000

Slope 2.01167055 0.04715615 0.5702868 0.1270237

Scale= 0.0009688619

redondeando, si la recta es y = a x + b, daría:

a=2.01 +/- 0.05

b=-0.02 +/- 0.08

No hay comentarios:

Publicar un comentario