Re: [R-br] precisão da operação

2015-09-22 Por tôpico Cesar Rabak
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"

Re: [R-br] Validação cruzada da krigagem

2015-09-22 Por tôpico Cesar Rabak
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

Re: [R-br] Validação cruzada da krigagem

2015-09-22 Por tôpico Wagner Wolff
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:

[R-br] Validação cruzada da krigagem

2015-09-22 Por tôpico Wagner Wolff
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 ___

[R-br] RES: Intervalos de confiança com boot.ci

2015-09-22 Por tôpico Mauro Sznelwar
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

Re: [R-br] RES: Intervalos de confiança com boot.ci

2015-09-22 Por tôpico ASANTOS
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,

[R-br] RES: Intervalos de confiança com boot.ci

2015-09-22 Por tôpico Mauro Sznelwar
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

Re: [R-br] Intervalos de confiança com boot.ci

2015-09-22 Por tôpico Elias Teixeira Krainski
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

[R-br] precisão da operação

2015-09-22 Por tôpico Cleber N.Borges
> 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

Re: [R-br] precisão da operação

2015-09-22 Por tôpico Bernardo Rangel Tura
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

Re: [R-br] Intervalos de confiança com boot.ci

2015-09-22 Por tôpico Pedro Emmanuel Alvarenga Americano do Brasil
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

Re: [R-br] precisão da operação

2015-09-22 Por tôpico salah
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 >

Re: [R-br] Automatizar a inclusão de váriaveis

2015-09-22 Por tôpico Nicolay Cunha
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

Re: [R-br] Intervalos de confiança com boot.ci

2015-09-22 Por tôpico ASANTOS
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,