Claro, incluyo el codigo que uso para generar las superficies de respuesta
que aparecen con los ejes codificados ... en el punto #1. se codifican los
datos de acuerdo a lo sugerido por la literatura de sup. de respuesta y por
la misma literatura de R ... luego el grafico aparece con las unidades
codificadas en los ejes ... ahora, acabo de leer anoche que no es forzoso
codificar los datos, pero que si ahorra mucho tiempo de calculo, voy a
probar y les cuento como me va, pero sigue en pie la pregunta de si es
posible, habiendo codificado los datos, obtener un grafico con las unidades
originales de las variables.

Saludos y muchas gracias por su tiempo.

Eric.

pd. ayer conteste la encuesta carlos, respecto al uso de R.



# AJUSTE DE UNA SUPERFICIE DE RESPUESTA PARA HEXANO
library("rsm")
setwd("/home/neo/Dropbox/Tesis/4tesis/1.objesp/medicionhc/gravimetria/")

# 0. PREPARANDO LOS DATOS PARA SU CODIFICACION
dathchexgrav <-
read.table("/home/neo/Dropbox/Tesis/4tesis/1.objesp/medicionhc/gravimetria/hchex.csv",
header=TRUE,sep="\t",dec=".",as.is=FALSE,na.strings = "NA")
summary(dathchexgrav)

# 1. CODIFICACION DE DATOS
dathchexgrav.coded <- coded.data(dathchexgrav, x1 ~ (tiempo - 2)/1, x2 ~
(ratio - 7.5)/2.5)
dathchexgrav.coded

# 2. AJUSTANDO EL MODELO LINEAL con interaccion
dathchexgrav.coded.rsm1 <- rsm(masadw ~ FO(x1, x2)+TWI(x1,x2), data =
dathchexgrav.coded)
summary(dathchexgrav.coded.rsm1)
            Call:
            rsm(formula = masadw ~ FO(x1, x2) + TWI(x1, x2), data =
dathchexgrav.coded)

            Residuals:
               Min     1Q Median     3Q    Max
            -6.450 -2.388 -2.033  3.133  8.550

            Coefficients:
                    Estimate Std. Error t value Pr(>|t|)
            (Intercept) 12.45000    1.07727  11.557 3.54e-09 ***
            x1           2.25000    1.39076   1.618    0.125
            x2           1.25000    1.39076   0.899    0.382
            x1:x2        0.08333    1.39076   0.060    0.953
            ---
            Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

            Residual standard error: 4.818 on 16 degrees of freedom
            Multiple R-squared: 0.1765,    Adjusted R-squared: 0.02207
            F-statistic: 1.143 on 3 and 16 DF,  p-value: 0.3618

            Analysis of Variance Table

            Response: masadw
                    Df Sum Sq Mean Sq F value Pr(>F)
            FO(x1, x2)   2  79.50  39.750  1.7126 0.2119
            TWI(x1, x2)  1   0.08   0.083  0.0036 0.9530
            Residuals   16 371.37  23.210
            Lack of fit  1   4.03   4.033  0.1647 0.6906
            Pure error  15 367.33  24.489

            Stationary point of response surface:
             x1  x2
            -15 -27

            Stationary point in original units:
            tiempo  ratio
               -13    -60

            Eigenanalysis:
            $values
            [1]  0.04166667 -0.04166667

            $vectors
                  [,1]       [,2]
            [1,] 0.7071068 -0.7071068
            [2,] 0.7071068  0.7071068



#3. PLOTEANDO LA SR
par(cex=1.3)
contour(dathchexgrav.coded.rsm1,~x1+x2+x1*x2,col="black",xlab=c("Time","Solvent:biomass
ratio"),labcex=1.3)
savePlot(filename = "graficos/3dathchexgrav_coded_rsm1_cp.png",type =
c("png"),device = dev.cur())










On Tue, Apr 3, 2012 at 5:34 AM, Carlos Ortega <cof@qualityexcellence.es>wrote:

> Hola,
>
> De forma directa con funciones del propio paquete rsm no veo solución...
>
> ¿Has probado a crear el contour directamente sobre tus datos
> no-codificados?.
> No sé si será posible si estás definiendo un modelo entre tus variables
> complicado. Otra alternativa es hacerlo con contourplot() de lattice.
> Si puedes proporcionar el código que utilizas y un conjunto de datos,
> quizás podríamos ir un poco más lejos...
>
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>
> El 3 de abril de 2012 01:06, Eric <ericconchamunoz@gmail.com> escribió:
>
>> estimados, ajusté una superficie de respuesta con el paquete "rsm" de R y
>> hago un grafico 2D de la superficie de respuesta con contour(), pero como
>> el argumento de esta función es la superficie de respuesta ajustada y ésta
>> se ha ajustado codificando las variables, de acuerdo a la recomendación de
>> la documentación del paquete rsm, obtengo un gráfico con los ejes en
>> unidades codificadas ... la pregunta es: hay alguna forma de obtener el
>> gráfico de superficie de respuesta con las unidades originales habiendo
>> ajustado la superficie de respuesta con las unidades codificadas ?
>>
>> saludos y muchas gracias por vuestro tiempo y ayuda,
>>
>> eric.
>>
>>
>>
>>
>> --
>> Nota: las tildes se han omitido para evitar conflictos con algunos
>> lectores
>> de correo.
>>
>> * Un economista es un experto que sabrá mañana por qué las cosas que
>> predijo
>> ayer no han sucedido hoy (Laurence Peter).
>>
>> * SATYÂT NÂSTI PARO DHARMAH
>>
>>        [[alternative HTML version deleted]]
>>
>>
>> _______________________________________________
>> R-help-es mailing list
>> R-help-es@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-help-es
>>
>>
>
>
> --
> Saludos,
> Carlos Ortega
> www.qualityexcellence.es
>



-- 
Nota: las tildes se han omitido para evitar conflictos con algunos lectores
de correo.

* Un economista es un experto que sabrá mañana por qué las cosas que predijo
ayer no han sucedido hoy (Laurence Peter).

* SATYÂT NÂSTI PARO DHARMAH

	[[alternative HTML version deleted]]

