Prezados Edson e Felipe,
obrigado pela pronta ajuda.
Att.,
*Emerson*
Em 23 de agosto de 2016 13:57, Felipe
escreveu:
> Segue um exemplo simples com dados fictícios:
> # Exemplo com x1:Tempo e x2:temperatura e y:reposta
> x1 <-
Segue um exemplo simples com dados fictícios:
# Exemplo com x1:Tempo e x2:temperatura e y:reposta
x1 <- c(-1,-1,1,1,0,0,0,0,0,sqrt(2),-sqrt(2),0,0)
x2 <- c(-1,1,-1,1,0,0,0,0,0,0,0,sqrt(2),-sqrt(2))
y2 <- c(76.5,77,78,79.5,79.9,80.3,80,79.7,79.8,78.4,75.6,78.5,77)
mod <- lm(y2 ~
Veja se ajuda
https://www.r-project.org/conferences/useR-2011/TalkSlides/Invited/Gromping-Design_of_Experiments.pdfhttps://cran.r-project.org/doc/contrib/Vikneswaran-ED_companion.pdfhttps://cran.r-project.org/web/views/ExperimentalDesign.html
[ ]'s
Edson Lira
Estatístico
Manaus-Amazonas
Caros colegas, bom dia.
Estou estudando sobre superfície de resposta. Alguém poderia me indicar
pacotes sobre esse assunto no R?
Também, sendo possível, pacotes para planejamento e análises de
experimentos.
Agradeço qualquer sugestão e indicação.
Att.,
*Emerson*
as is
ou seja o argumento é interpretado como é
por exemplo
x1 * x2 é interpretado como uma interacao em uma formula (se forem
continuas, 3 parametros)
I(x1*x2) é um unico termo, com um parametro, dado pelo produto das
variaveis
mais em na sessao
23 Fórmulas e especificação de modelos
Olá Pessoal,
Estou precisando produzir um gráfico de superfície para o seguinte modelo :
lm(y=a+bx+cz+dx*z)
Alguém poderia me dizer que pacote poderia usar e qual o código?
Obrigada...
--
Evelyn da Fonseca Alecrim Bragion
Mestranda em Ecologia Aplicada - UFLA
Lavras - MG/Brasil
Evelyn,
Sempre poste um CMR, ajuda um monte...
uma opção pode ser:
x - rnorm(50)
z - rnorm(50)
y - x^2+z^2
m - lm(y~I(x^2)+I(z^2))
summary(m)
Ngrid-100
xx-seq(-2,2,l=Ngrid)
zz-seq(-2,2,l=Ngrid)
dat - expand.grid(xx=xx,zz=zz)
dat$ypre - predict(m,newdata=list(x=dat$xx,z=dat$zz))
Qual o significado deste I em I(x^2)+I(z^2)
Evelyn,
Sempre poste um CMR, ajuda um monte...
uma opção pode ser:
x - rnorm(50)
z - rnorm(50)
y - x^2+z^2
m - lm(y~I(x^2)+I(z^2))
summary(m)
Ngrid-100
xx-seq(-2,2,l=Ngrid)
zz-seq(-2,2,l=Ngrid)
dat - expand.grid(xx=xx,zz=zz)
dat$ypre -
Mauro,
?I()
Att
Em 10 de janeiro de 2012 21:54, Mauro Sznelwar sznel...@uol.com.br escreveu:
Qual o significado deste I em I(x^2)+I(z^2)
Evelyn,
Sempre poste um CMR, ajuda um monte...
uma opção pode ser:
x - rnorm(50)
z - rnorm(50)
y - x^2+z^2
m - lm(y~I(x^2)+I(z^2))
summary(m)