Março 2011 - R-br

Testando a criação de um novo tópico.
by Pedro Rafael 24 Mar '11

24 Mar '11

Esse tópico é um simples teste. -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 1

Dúvida
by Pedro Rafael 24 Mar '11

24 Mar '11

Testando envio -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

1 0

Ler linhas
by Victor Eduardo 24 Mar '11

24 Mar '11

Como faço para ler linhas de arquivos que começam somente pelos numeros 82 ou 83? Por exemplo, vou colocar alguns dados aqui: V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 8345 2010 5 20 15 0 -38.20 -12.70 11 -9999 1016 1017 -9999 -9999.0 309.2 301.4 5122 2010 5 20 15 0 -43.80 -21.20 1000 -9999 887 998 284 2.5 292.1 289.1 8222 2010 5 20 15 0 -43.30 -22.60 33 -9999 1012 1016 131 0.4 294.8 292.6 8343 2010 5 20 15 0 -37.00 -10.90 5 -9999 1012 1013 -9999 -9999.0 301.6 296.9 4107 2010 5 20 15 0 -38.20 -9.40 253 -9999 984 1012 109 3.8 300.6 293.9 No caso, o meu interesse está na V1, quero pegar os numeros que comecam por 82 ou 83 e ao pegar, queria que os dados viessem dessa maneira, por exemplo: 8345 2010 5 20 15 0 -38.20 -12.70 11 -9999 1016 1017 -9999 -9999.0 309.2 301.4 8222 2010 5 20 15 0 -43.30 -22.60 33 -9999 1012 1016 131 0.4 294.8 292.6 8343 2010 5 20 15 0 -37.00 -10.90 5 -9999 1012 1013 -9999 -9999.0 301.6 296.9

5 7

Email teste (Como beber um litro de jurubeba numa golada só? )
by Ivan Bezerra Allaman 24 Mar '11

24 Mar '11

Senhores, este é um email teste visto que alguns membros não estão conseguindo postar suas dúvidas via email. Por favor, o mantenedor do grupo apague-a após recebe-la. Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Dúvida no rotina rcorr no pacote Hmisc
by Edson Lira 24 Mar '11

24 Mar '11

Como faço para colocar os nomes das variáveis na rotina abaixo. Faço a rotina para qa correlação:require(Hmisc) rcorr(cbind(corre$prevalência,corre$precipitação,corre$temp_max, + corre$temp_min) + ,type=c("pearson")) Esta saída é gerada. Como faço para colocar os nomes da variável. [,1] [,2] [,3] [,4] [1,] 1.00 -0.04 -0.15 -0.28 [2,] -0.04 1.00 -0.80 -0.41 [3,] -0.15 -0.80 1.00 0.67 [4,] -0.28 -0.41 0.67 1.00 n= 63 P [,1] [,2] [,3] [,4] [1,] 0.7841 0.2516 0.0249 [2,] 0.7841 0.0000 0.0010 [3,] 0.2516 0.0000 0.0000 [4,] 0.0249 0.0010 0.0000 Edson Lira Estatístico Manaus-Amazonas

2 2

Ajuda com modelo Log-linear (agora com dados)
by Diogo Borges Provete 24 Mar '11

24 Mar '11

Olá a todos novamente e desculpem-me pela submissão anterior, O problema que tenho é o seguinte: gostaria de uma ajuda para ajustar um modelo log-linear para dados com 3 categorias de pigmentação na superfície de órgãos da cavidade abdominal de anfíbios. Os dados em anexo referem-se ao coração. Os dados brutos consistem das categorias de pigmentação (0, 1 e 2) de 10 indivíduos divididos nos grupos experimentais (nas colunas). Eu gostaria de saber se há diferença nas categorias de pigmentação de acordo com os tratamentos experimentais. A tabela de contingência tem a informação da frequência das categorias de pigm distribuídas nos grupos experimentais. Eu usei a função loglm do pacote MASSque permite que se insira a fórmula parecida com o lm e glm. Eu inclui somente os dados pq a função diz que se a fórmula não for especificada pode-se colocar somente a origem dos dados. Mas não soube interpretar os resultados. Agradeço desde já Atenciosamente, Diogo Borges Provete ============================== Biólogo Mestre em Biologia Animal (UNESP) Doutorando PPG Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás - UFG Goiânia-GO CP: 131 74001-970 Brazil Tel. +55 62 3521-1732 Cel. +55 62 8231-5775 Skype: diogoprovete MSN: diogoprov(a)yahoo.com.br Personal web page Traduza conosco: <<<American Journal Experts>>> <<<D-Lang Soluções linguisticas>>> <<<Perfil no ProZ.com como tradutor freelancer>>> ==============================

1 0

Estimativa do intervalo de confianca em um modelo ZINB
by Humberto Hazin (hotmail) 24 Mar '11

24 Mar '11

Olá a todos, Tenho 2 problemas se alguém puder me ajudar (Um exemplo dos dados encontra-se em anexo): 1- como faço para estimar os intervalos de confiança dos coeficientes de um modelo ZINB. Estou usando o package (pscl). Já tentei por bootstrap mais ele não estima os intervalos. Segue abaixo o script que usei data(bum2) F1a<-formula(BUM ~ Year + Quarter + area + Strategy+offset(log(effort1))) options(contrasts=c("contr.treatment","contr.poly")) ZINB<-zeroinfl(F1a,dist="negbin",link = "logit",data=bum2) boot.zinb <- function(data, indices){ require(pscl) data <- data[indices,] try(mod <- zeroinfl(F1a, data = data, link = "logit", dist = "negbin", trace = FALSE EM = FALSE)) if (exists("mod")) { coef(mod) } else { NA } } system.time(zinb.boot.out <- boot(data = bum2, statistic = boot.zinb, R = 99,sim="permutation")) zinb.boot.out boot.ci(zinb.boot.out,type="bca") 2- como faço para estimar os predicts e os intervalos de confiança apenas da variável YEAR de um modelo ZINB. Abaixo o código usado: x<-seq(1980,2008) newdata=cbind.data.frame(Year=as.factor(x),Quarter=as.factor(rep(1,length(x) )),area=as.factor(rep(1,length(x))),Strategy=as.factor(rep(4,length(x))),eff ort1=rep(1000,length(x))) newdata fm_zip <- zeroinfl(F1a,dist="negbin",link = "logit",data=bum2) fit <- predict(fm_zip) Pearson <- resid(fm_zip, type = "pearson") VarComp <- resid(fm_zip, type = "response") / Pearson year <- bum2$Year bootstrap <- replicate(999, { yStar <- pmax(round(fit + sample(Pearson) * VarComp, 0), 0) predict(zeroinfl(yStar ~ year | 1), newdata = newdata) }) newdata0 <- newdata newdata0$fit <- predict(fm_zip, newdata = newdata, type = "response") newdata0<- t(apply(bootstrap, 1, quantile, c(0.025, 0.975))) Desde Ja agradeço a vossa atenção Humberto

1 0

Re: [R-br] [Dúvida] Ler dados direto no HD
by beniltoncarvalho 24 Mar '11

24 Mar '11

Eu nao mencionei converter p Excel :) . Ha' vezes que conversao para txt e' um mal necessario, mas se vc estiver em Windows pode ser que haja problemas tambem (depende do filesystem que vc estiver usando). De qualquer forma, voce pode usar o pacote RODBC com os comandos odbcConnectDbase e sqlFetch. Verifique o arquivo de ajuda para os detalhes. -- View this message in context: http://r-br.2285057.n4.nabble.com/Duvida-Ler-dados-direto-no-HD-tp3400886p3… Sent from the R-br mailing list archive at Nabble.com.

1 0

Ajuda com modelo Log-linear
by Diogo Borges Provete 23 Mar '11

23 Mar '11

Olá a todos, gostaria de uma ajuda para ajustar um modelo log-linear para dados com 3 categorias de pigmentação na superfície de órgãos da cavidade abdominal de anfíbios. Os dados seguem em anexo, a planilha é referente ao coração, a parte da esquerda são os dados originais, que consistem das categorias de pigmentação (0, 1 e 2) de 10 indivíduos divididos nos grupos experimentais (nas colunas). Eu gostaria de saber se há diferença nas categorias de pigmentação de acordo com os tratamentos experimentais. O pacote MASS tem uma função loglm que permite que se insira a fórmula parecida com o lm e glm, mas não sei como usar isso neste caso. Agradeço desde já Atenciosamente, Diogo Borges Provete ============================== Biólogo Mestre em Biologia Animal (UNESP) Doutorando PPG Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás - UFG Goiânia-GO CP: 131 74001-970 Brazil Tel. +55 62 3521-1732 Cel. +55 62 8231-5775 Skype: diogoprovete MSN: diogoprov(a)yahoo.com.br Personal web page Traduza conosco: <<<American Journal Experts>>> <<<D-Lang Soluções linguisticas>>> <<<Perfil no ProZ.com como tradutor freelancer>>> ==============================

2 1

visita Edzer pebesma ao LEG e mini-curso
by Paulo Justiniano 23 Mar '11

23 Mar '11

Caros Na proxima semana, entre 29-31/03/2011 estaremos recebendo no LEG o Edzer Pebesma Edzer Pebesma Institute for Geoinformatics (ifgi) University of Münster Germany. Phone: +49 251 O Edzer é o autor do pacote gstat, co-autor do pacote pacote sp e co-autor do livro: Analysis of Spatial Data with R (ASDAR) (com Roger Bivand e Virgilio G. Rubio) Durante a visita serah ministrado um curso de "Spatial temporal analysis" 1 hr na terca a tarde (29) 4 hrs na quata de manha (30) e se necessario 1 hr extra na quinta (31) (na segunda 28/03) rodaramos uma intro sobre recursos antes do inicio do curso Paginas relevantes: Instituto: http://ifgi.uni-muenster.de Pagina do Edzer no Instituto: http://ifgi.uni-muenster.de/en/staff/7-mitarbeiterinnen/3-pebesma-edzer-pro… Lista de Publicações: http://ifgi.uni-muenster.de/~epebe_01/publ.html Projeto/Comunidade http://www.52north.org/geostatistics Paulo Justiniano Ribeiro Jr LEG (Laboratorio de Estatistica e Geoinformacao) Universidade Federal do Parana Caixa Postal 19.081 CEP 81.531-990 Curitiba, PR - Brasil Tel: (+55) 41 3361 3573 VOIP: (+55) (41) (3361 3600) 1053 1066 Fax: (+55) 41 3361 3141 e-mail: paulojus AT ufpr br http://www.leg.ufpr.br/~paulojus

1 0