Kleber,
Outros já deram algumas respostas na direção certa, mas há ainda mais uma
observação a fazer:
"a capacidade de representação dos números estão na faixa de 1.xxe-308 a
1.xxe+308" diz respeito à faixa de *magnitudes* dos números resta saber
qual é *mantissa* deles (denominado "signifcando"
Matematicamente a resposta seria sim, *mas* seria bom tomar cuidado pois se
o número de pontos for muito grande esses testes podem ser muito
'sensíveis' e acabar rejeitando a hipótese nula mesmo que não faça sentido
e se o número de casos for pequeno eles podem deixar passar óbvias não
Ola Cesar.
Então qdo li sobre o teste de homocedasticidade uma das pressuposições é
que ele serve para grandes amostras, mas o que seria uma grande amostra?
Qual o limiar? Minha análise tem 83 pontos isso é grande?
Att.
Wagner
Em 22/09/2015 18:29, "Cesar Rabak" escreveu:
Olá pessoal da lista, a dúvida é a seguinte!
É possível utilizar testes como Shapiro-Wilk e Breusch Pagan para
verificar, respectivamente a normalidade e homocedasticidade de resíduos
oriundos de uma validação cruzada da krigagem?
Att
Wagner Wolff
___
Alguém poderia passar o código anterior, pois eu apaguei o que chegou antes.
vc pode fazer assim:
SMR2 <- function(a,w=1)
sum(w*a$Hosp_Death)/sum(w*a$SAPS3Pro2)
SMR2(a) ### SMR observada nos dados
require(boot)
plot(boo <- boot(a,SMR2,999,stype='i'))
boot.ci(boo)
.
---
Este email
Mauro, segue o script,
dput(a[sample(1:nrow(a),100),c("Hosp_Death","SAPS3Pro2")])
structure(list(Hosp_Death = c(0L, 1L, 0L, 1L, 0L, 0L, 0L, 1L,
0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
1L, 1L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L,
0L, 0L, 0L, 0L, 0L, 1L,
Tentei rodar e não consegui.
> dput(a[sample(1:nrow(a),100),c("Hosp_Death","SAPS3Pro2")])
Error in dput(a[sample(1:nrow(a), 100), c("Hosp_Death", "SAPS3Pro2")]) :
object 'a' not found
Amigos de R,
Eu estou envolvido num projeto que o coordenador deseja estimar SMR em diversos
vc pode fazer assim:
SMR2 <- function(a,w=1)
sum(w*a$Hosp_Death)/sum(w*a$SAPS3Pro2)
SMR2(a) ### SMR observada nos dados
require(boot)
plot(boo <- boot(a,SMR2,999,stype='i'))
boot.ci(boo)
___
R-br mailing list
R-br@listas.c3sl.ufpr.br
> 12345678912345678+1
[1] 12345678912345680
> 1.2345678912345678+0
[1] 1.2345678912345679
Prezados,
pq exatamente o fato acima apresentado acontece?
se a capacidade de representação dos números estão
na faixa de 1.xxe-308 a 1.xxe+308
agradeço antecipadamente pela atenção
cleber
---
Este
On 09/22/2015 07:54 AM, Cleber N.Borges wrote:
> 12345678912345678+1
[1] 12345678912345680
> 1.2345678912345678+0
[1] 1.2345678912345679
Prezados,
pq exatamente o fato acima apresentado acontece?
se a capacidade de representação dos números estão
na faixa de 1.xxe-308 a 1.xxe+308
agradeço
O objeto a é o banco de dados que precisa ser criado a partir da saida do
dput().
Pedro Brasil
Em 22 de setembro de 2015 11:19, ASANTOS
escreveu:
> Prezados,
>
>Também tentei rodar a rotina e o objeto a não foi encontrado, sedo
> o erro
>
> Erro em
Esse é o grande Deus da informática: "floating point"
Faça
options(digits = 22)
k = 123456789123456789
sprintf('%1.2f', k)
Sugiro a leitura:
http://www.burns-stat.com/pages/Tutor/R_inferno.pdf
saudações
Em Ter, 2015-09-22 às 07:54 -0300, Cleber N.Borges escreveu:
> > 12345678912345678+1
>
Prezados,
O problema foi resolvido. Abaixo a solução:
fs=as.numeric()
for(i in 2:ncol(y)){
m1 <- capscale((y[,1:i]) ~ x, dist="bray")
fs[i]=anova(m1)$F[1]
}
Nicolay
2015-09-21 16:28 GMT-04:00 Nicolay Cunha :
> Prezados,
>
>
>
> Gostaria de ajuda para automatizar a
Prezados,
Também tentei rodar a rotina e o objeto a não foi encontrado,
sedo o erro
Erro em NROW(data) : objeto 'a' não encontrado
--
==
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação,
14 matches
Mail list logo