Hola Jorge. Perfecto ya resolvi la primera cuestion con tu indicacion Gracias!
en cuanto a lo segundo estoy calculando una "suma de cuadrados de
residuales penalizada (PRSS)",para elegir uno de los splines de entre
los muchos posibles, segun varian los grados de libertad o (df).
El asunto es que esto requiere la "integral de la derivada 2da del
spline al cuadrado" como parte del calculo...
Usando predict(smooth.spline("mis datos"),deriv=2) consigo los valores
de la derivada 2 pero no consigo ver la formula resultante y menos
calcular su integral....
cualquier consejo o pista que me ponga en camino es bienvenido.
El 5 de julio de 2009 19:47, Jorge Ivan Velez
<jorgeivanvelez en gmail.com> escribió:>
> Buenas noches Nahuel,
> Bienvenido a R-help-es!
> La extracción de lambda del objeto que contiene el smooth.spline() ajustado
puede obtenerse fácilmente al analizar la estructura de dicho objeto utilizando
la función str() del R. Consideremos el ejemplo que aparece en ?smooth.spline():
> attach(cars)
> plot(speed, dist, main = "data(cars) & smoothing splines")
> cars.spl <- smooth.spline(speed, dist)
> cars.spl
> str(cars.spl)
> El resultado de str(cars.spl) es una lista de 15 elementos y el valor de
lambda puede obtenerse como:
> cars.spl$lambda
> # [1] 0.1112206
> En cuanto a tu segunda pregunta, podrías darnos un poco más de información
al respecto?
> Espero sea de utilidad,
> Jorge Ivan Velez
>
> 2009/7/5 Nahuel farias <nahuelfarias en gmail.com>
>>
>> hola soy nuevo usuario de R y necesito crear un objeto p tal que,
>>
>> p=lambda*integral[f´´(x)^2 dx ]
>> donde "lambda" es uno de los parametros que resultan de la
funcion
>> "smooth.spline()" y la integral es sobre la derivada 2 de esa
misma
>> funcion...
>> dos cuestiones:
>> 1) como extraigo lambda de los resultados de smooth.spline() para
usarlo
>> como objeto cuando lo requiera y
>> 2) como puedo obtener con R la integral del segundo termino?
>>
>> gracias
>>
>> --
>> Lic. Nahuel E. Farias
>> Laboratorio de Invertebrados
>> Departamento de Biología
>> Facultad de Ciencias Exactas y Naturales,
>> Universidad Nacional de Mar Del Plata, CC 1245, CP 7600 Mar del Plata,
>> Argentina
>> Tel. + 54 - 0223 - 4752426 (Int. 466)
>>
>> [[alternative HTML version deleted]]
>>
>>
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es en r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>
--
Lic. Nahuel E. Farias
Laboratorio de Invertebrados
Departamento de Biología
Facultad de Ciencias Exactas y Naturales,
Universidad Nacional de Mar Del Plata, CC 1245, CP 7600 Mar del Plata, Argentina
Tel. + 54 - 0223 - 4752426 (Int. 466)