Saludos a todos los miembros de la comunidad.

Estoy en estos momentos intentando resolver un problemita con R (una
metodología que tengo que repetir cerca de 100 veces con un conjunto de
datos) y me gustaría convertir el resultado de una regresión lineal en una
función que pueda ser integrada posteriormente, es decir:

Tengo los siguiente:
 >file <- commandArgs(trailingOnly=TRUE)
 >options(digits=16)
 >options(max.print=5.5E5)
 >x <- seq(50001,100000,1)
 >datos <- read.table(file)
 >attach(datos)
 >ajust <- lm(V1~x,data=datos)
 >ajust.

Hasta aquí no tengo problemas y se resolver integrales de manera numérica
en R de manera aislada, por ejemplo:

>## define the integrated function
>integrand <- function(x) {1/((x+1)*sqrt(x))}
>## integrate the function from 0 to infinity
>integrate(integrand, lower = 0, upper = Inf)

Cuando quiero intentar relacionar ambas cosas en una sola, no tengo idea de
como pasar el resultado que arroja "ajust" y convertirlo en una función que
pueda integrarse usando el ejemplo que os muestro arriba.

Alguna idea de como pudiera empezar para resolver este problemilla.!

Gracias de antemano.



-- 
*Daniel Arismendi
*

	[[alternative HTML version deleted]]

