Re: [R-br] RStudio para criar e testar pacotes

[Elias Carvalho via R-br]

Ola Valmes

Obrigado pelo retorno

Isso acontece quando tento usar alguma função do pacote

Em 17 de dezembro de 2016 10:37, Walmes Zeviani 
escreveu:

> Professor Elias,
>
> Eu uso os recursos do roxygen2 e devtools, porém sem usar o RStudio. Ainda
> não me deparei com os problemas que você mencionou. Isso está acontecendo
> quando roda o check()?
>
> À disposição.
> Walmes.
> ​
>



-- 


*In Jesu et Maria*
*Obrigado*
*Prof. Elias Carvalho*

*"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he
who has been able to understand the cause of things"*
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

[R-br] dose-resposta

[Maurício Lordêlo via R-br]

Caros,
Foi realizado um experimento no qual deseja-se saber qual composto extraído
de uma determinada espécie vegetal é mais eficiente para matar uma
determinado tipo de inseto.
Foram utilizadas 32 placas cada uma contendo 10 insetos vivos.
Em cada placa, foram realizadas 4 repetições das 8 diferentes doses do
composto de cada espécie.
O número de insetos mortos foi registrado após um determinado período de
tempo.
O objetivo é encontrar qual a dose, para cada um das três espécies, que
mata 50% dos insetos.
Alguma sugestão de como encontrar os modelos para as três espécies?



dose = rep(c(0,0.15625,0.31250,0.62500,1.2500, 2.5000, 5., 10.000),
each=4)
n_insetos = rep(10,32)
m_especie1 =
c(1,3,4,0,5,2,5,5,4,4,2,2,0,3,3,5,10,10,7,0,10,10,10,10,10,10,10,10,10,10,10,10)
m_especie2 = c(0,3,4,4,1,0,4,0,0,1,0,0,3,2,0,3,3,2,0,3,3,3,3,6,9,4,7,2,10,
10, 10, 10)
m_especie3 =
c(4,2,2,3,4,6,8,6,4,6,6,5,5,9,5,8,10,8,5,5,10,10,10,10,10,10,10,10,10,10,10,10)
pm_especie1 = m_especie1/n_insetos
pm_especie2 = m_especie2/n_insetos
pm_especie3 = m_especie3/n_insetos

dose.fator = as.factor(dose)

boxplot(pm_especie1~dose.fator, xlab= "Dose", ylab="Proporção de mortos
espécie 1",cex.axis=.85)
boxplot(pm_especie2~dose.fator, xlab= "Dose", ylab="Proporção de mortos
espécie 2",cex.axis=.85)
boxplot(pm_especie3~dose.fator, xlab= "Dose", ylab="Proporção de mortos
espécie 3",cex.axis=.85)

bartlett.test(pm_especie1~dose.fator)
bartlett.test(pm_especie2~dose.fator)
bartlett.test(pm_especie3~dose.fator)

modelo1 = glm(cbind(good=m_especie1,
bad=n_insetos-m_especie1)~dose,family="binomial", data=dados)
par(mfrow=c(2,2))
plot(modelo1)
shapiro.test(modelo1$residuals)

modelo2 = glm(cbind(good=m_especie2,
bad=n_insetos-m_especie2)~dose,family="binomial", data=dados)
par(mfrow=c(2,2))
plot(modelo2)
shapiro.test(modelo2$residuals)

modelo3 = glm(cbind(good=m_especie3,
bad=n_insetos-m_especie3)~dose,family="binomial", data=dados)
par(mfrow=c(2,2))
plot(modelo2)
shapiro.test(modelo3$residuals)
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

Re: [R-br] Modelar os Ys de um data frame

[Karina Rebuli via R-br]

Oi, André.

Eu gosto de usar o lapply nessas situações. Facilita até para comparar os
modelos, aplicar predcit etc., pois retorna uma list (com a estutura
abaixo, em que cada elemento terá o nome da respectiva coluna do data.frame
passado como objeto).

models <- lapply( dados[, 1:4], function( y ){
  lm( y ~ X1*X2 )
})

Abs,
Karina


On Tue, Dec 20, 2016 at 2:37 PM Andre Oliveira via R-br <
r-br@listas.c3sl.ufpr.br> wrote:

>
> Pessoal como eu poderia modelar os Ys com  um for ou de alguma outra
> forma? (Tenho 20 Ys)
>
> Exemplo.
> library(mvtnorm)
> Sigma <- matrix(c(1,1,-.4,-.9,1,4,1.6,-1.2,-.4,1.6,4,
> 1.8,-.9,-1.2,1.8,9),ncol=4)
> Y <- rmvnorm(10, c(2,4,6,8), Sigma)
> X1=as.factor(c(rep("A",5),rep("B",5)))
> X2=as.factor(c(rep("1",2),rep("3",3),rep("2",2),rep("0",3)))
> dados = data.frame(Y,X1,X2)
>
> for (i in 1:4) lm(dados[,i] ~X1*X2))
>
> No verdade desejo usar ajustar com a função!
>
> require(ExpDes.pt)
> fat2.dic(X1, X2, dados[i], quali = c(TRUE, TRUE), mcomp = "sk", fac.names
> = c("Doses", "Substrato"), sigT = 0.05, sigF = 0.05)
>
> obg
>
>
>
>
> André Oliveira Souza. Graduação em Matemática, mestrado em estatística
> aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito
> Santo.  IFES
> ___
> R-br mailing list
> R-br@listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.

[R-br] Modelar os Ys de um data frame

[Andre Oliveira via R-br]

Pessoal como eu poderia modelar os Ys com  um for ou de alguma outra forma? 
(Tenho 20 Ys) 

Exemplo. 
 library(mvtnorm) 
 Sigma <- matrix(c(1,1,-.4,-.9,1,4,1.6,-1.2,-.4,1.6,4, 
1.8,-.9,-1.2,1.8,9),ncol=4)
 Y <- rmvnorm(10, c(2,4,6,8), Sigma) 
 X1=as.factor(c(rep("A",5),rep("B",5)))
 X2=as.factor(c(rep("1",2),rep("3",3),rep("2",2),rep("0",3)))
 dados = data.frame(Y,X1,X2)
 for (i in 1:4) lm(dados[,i] ~X1*X2))
 No verdade desejo usar ajustar com a função! 

 require(ExpDes.pt)
 fat2.dic(X1, X2, dados[i], quali = c(TRUE, TRUE), mcomp = "sk", fac.names = 
c("Doses", "Substrato"), sigT = 0.05, sigF = 0.05)
 obg


 André Oliveira Souza. Graduação em Matemática, mestrado em estatística 
aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. 
 IFES___
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo 
m�nimo reproduz�vel.