Novembro 2014 - R-br

Seria um caso de R inferno?
by Jose Claudio Faria 03 Nov '14

03 Nov '14

Olá, Estava preparando um material didático para aulas de introdução à estatística, na parte das variáveis aleatórias contínuas e me deparei com um problema que ainda não encontrei resposta. Seria um caso de R inferno? # Função de densidade: chisq fchisq <- function(x, n) (1/(2^(n/2)*gamma(n/2)))*x^(n/2-1)*exp(-x/2) par(mfrow=c(1, 2)) plot(fchisq(x=1:100, n=10)) curve(dchisq(x, 10), col='red', lw=2, add=TRUE, n=1e3) ## OK: a função fchisq está correta! plot(fchisq(x=seq(0, 50, len=100), n=10)) curve(dchisq(x, 10), col='red', lw=2, add=TRUE, n=1e3) ## ? Gastei um bom tempo para perceber que o uso de "seq" não é correto nesses casos. Ab, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9966.9100 - VIVO 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

3 5

Re: [R-br] Função likfitBGCCM e predict
by Jacqueline.cantu 03 Nov '14

03 Nov '14

Ao trocar y1por y2 todos os parâmetros da função mudam.. Por se tratar do modelo bivariado BGCCM os parâmetros não deveriam ser os mesmos? Enviado por Samsung Mobile "walmes ." <walmeszeviani(a)gmail.com> wrote: Não é só uma questão de trocar y1 por y2 em todas as ocorrências, e vice versa? À disposição. Walmes. Em 03/11/2014 16:23, "Jacqueline Cantú" <jacqueline.cantu(a)hotmail.com> escreveu: Caros, Usando os comando que seguem (os quais estão disponíveis em: http://www.leg.ufpr.br/geoR/tutorials/CCM.R) Gostaria de saber como proceder para obter a predição espacial para a segunda variável (y2), uma vez que essa rotina me fornece apenas os valores da predição da primeira variável (y1). require(geoR) set.seed(30) all.coords <- round(cbind(runif(200), runif(200)), dig=3) set.seed(111); ind1 <- sample(1:200, 120) coords1 <- all.coords[ind1,] set.seed(222); ind2 <- sample(1:200, 80) coords2 <- all.coords[ind2,] ## Model parameters: mu1 <- 10; mu2 <- 50 sigma01 <- 2; sigma1 <- 1.5; sigma02 <- 8; sigma2 <- 6 phi0 <- 0.2; phi1 <- 0.15; phi2 <- 0.25 ## Model parameters (reparametrised) sigma <- sigma01 nu1 <- sigma1/sigma eta <- sigma02/sigma; nu2 <- sigma2/sigma ## Simulating model components S0 <- grf(grid=all.coords, cov.pars=c(sigma, phi0))$data S1 <- grf(grid=coords1, cov.pars=c(sigma1, phi1))$data S2 <- grf(grid=coords2, cov.pars=c(sigma2, phi2))$data ## Y1 and Y2 data y1 <- as.geodata(cbind(coords1,10+S0[ind1]+S1)) y2 <- as.geodata(cbind(coords2,50+S0[ind2]+S2)) ## maximum likelihood estimation fit12 <- likfitBGCCM(y1, y2, ini.s=c(2,1.5,8,6), ini.phi=c(.2,.15,.25), control=list(trace=T)) fit12 ## Prediction locs <- cbind(c(0.2, 0.5, 0.7, 0.2), c(0.3, 0.5, 0.2, 0.8)) pred12 <- predict(fit12, loc=locs) Att, Jacqueline. _______________________________________________ R-br mailing list R-br(a)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.

1 0

separador decimal e milhar
by Cleber N.Borges 03 Nov '14

03 Nov '14

Prezados, Boa Tarde... Tenho uma tabela de valores monetários em que há: 1) separação decimal por vírgula e 2) separação de milhar por ponto. como ler isso direto no R? Pensei na função read.table, porém só vi a opção *dec=","* e nada para a questão da separação de casa de milhar. Agradeço antecipo pela ajuda Cleber Exemplo de linhas: 1 2.195,46 666,66 0,00 1.436,83 46,03 20,94 25,00 200.000,00 2 2.238,43 666,66 0,00 1.479,95 45,88 20,94 25,00 199.333,34 3 2.233,33 666,66 0,00 1.475,00 45,73 20,94 25,00 198.666,68

3 2

dados Facebook
by Vinicius Brito Rocha 03 Nov '14

03 Nov '14

Senhores, alguém já extraiu informação do facebook com o R? O pacote Rfacebook está com problemas. Abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário * *M.Sc. Engenharia de Produção/PO* *Ph.D. Estudante Sistemas Computacionais*

1 0

Dúvida teórica sobre hipótese nula e alternativa
by Fernando Antonio de souza 02 Nov '14

02 Nov '14

Caros colegas Não sei se é uma questão ingênua, me ocorreu agora e queria uma opinião. A hipótese nula representa a hipótese comumente aceita e que o pesquisador deseja provar ser falsa. Esta hipótese não necessariamente expressa igualdade de dois tratamentos ou da media da população a um valor pre-determinado. Por exemplo se é sabido q dois tratamentos são diferentes e você deseja demonstrar q são iguais, a hipótese nula não é de igualdade mas de diferença. Estou certo? Como procedo na Anova uma vez q a hipótese nula no R é de igualdade? Penso q conduzo a analise da mesma forma só q meu interesse estará em aceitar a hipótese nula ao invés de rejeita-la . Não deveria o software permitir a definição da hipótese nula pelo usuário?

3 2