Eric
2012-Apr-02 23:06 UTC
[R-es] graficar superficie de respuesta con los valores de las variables originales ...
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]]
Carlos Ortega
2012-Apr-03 09:34 UTC
[R-es] graficar superficie de respuesta con los valores de las variables originales ...
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 [[alternative HTML version deleted]]
Eric
2012-Apr-03 13:45 UTC
[R-es] graficar superficie de respuesta con los valores de las variables originales ...
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]]