- R-br - DINF UFPR

Ler dados do epi info
by Emerson Cotta Bodevan 11 Mai '12

11 Mai '12

Prezados, como faço para importar dados do epi info para o R.? -- *Emerson*

4 4

Simulação de um DQL no R
by Edmundo Caetano 11 Mai '12

11 Mai '12

Prezados, Pretendo simular um experimento em quadrados latinos no R. Como devo proceder??? Edmundo Caetano

4 8

Problemas com comando witch
by Tito Conte 11 Mai '12

11 Mai '12

tenho dois vetores: p=c(2,3,4,5,6,7,8,9) r=c(3,4,5) quando uso which(p==r) ou which(r==p) ele me retorna as posições onde r ocorre em p como fazer para saber onde p ocorre em r? obrigado (já dei um google nisso e não encontrei ) Tito Conte

4 7

Escalas
by Eric Keven 10 Mai '12

10 Mai '12

Boa noite, Como eu faço para alterar a escala dos gráficos e fazer com que o programacomece num valor que eu fixar para a escala?Obrigado! Att, Eric Keven Silva Coordenador do Núcleo de Estudos de Processos da Indústria de Alimentos (PROALI) Graduando em Engenharia de Alimentos Universidade Federal de Lavras

5 4

Comparações múltiplas usando MMC (mean–mean multiple comparisons) plots
by ASANTOS 10 Mai '12

10 Mai '12

*Bom dia pessoal, Estou tendo problemas em visualizar os resultados de uma comparação múltipla usando GLM, com o uso MMC plots do pacote HH. Bom, tenho quatro grupos funcionais de cupins que são humívoros (humi), ceifadores(ceifa), intermediários (inter) e xilófagos (xilo) que são meus tratamentos e como variável resposta o número de ocorrência de espécies de cada grupo (acorr). Para coleta da variável resposta, fizemos 6 transectos que são as repetições. Para ver se existe diferença entre os grupos funcionais fiz:* # mod2<-glm(ocorr~grupo, family="quasipoisson", data=trat1_90antes)## Um plot(mod2) prévio mostrou que quasipoisson foi a distribuição mais adequada aos dados que poisson summary(mod2) glm(formula = ocorr ~ grupo, family = "quasipoisson", data = trat1_90antes) Deviance Residuals: Min 1Q Median 3Q Max -1.9148 -0.4946 0.1213 0.4178 1.3813 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.5041 0.1624 9.260 1.13e-08 *** grupohumi -1.0986 0.3249 -3.382 0.00296 ** grupointer -0.8979 0.3019 -2.974 0.00750 ** grupoxilo -0.6568 0.2780 -2.363 0.02838 * --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Dispersion parameter for quasipoisson family taken to be 0.7123377) Null deviance: 28.663 on 23 degrees of freedom Residual deviance: 16.889 on 20 degrees of freedom AIC: NA Number of Fisher Scoring iterations: 5 anova(mod2, test="Chi") ## Anova do modelo Analysis of Deviance Table Model: quasipoisson, link: log Response: ocorr Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev Pr(>Chi) NULL 23 28.663 grupo 3 11.774 20 16.889 0.0008834 *** --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 plot(mod2) ##Teste de comparação múltipla summary(glht(mod2, linfct=mcp(grupo="Tukey"))) Simultaneous Tests for General Linear Hypotheses Multiple Comparisons of Means: Tukey Contrasts Fit: glm(formula = ocorr ~ grupo, family = "quasipoisson", data = trat1_90antes) Linear Hypotheses: Estimate Std. Error z value Pr(>|z|) humi - ceifa == 0 -1.0986 0.3249 -3.382 0.00405 ** inter - ceifa == 0 -0.8979 0.3019 -2.974 0.01491 * xilo - ceifa == 0 -0.6568 0.2780 -2.363 0.08255 . inter - humi == 0 0.2007 0.3794 0.529 0.95122 xilo - humi == 0 0.4418 0.3606 1.225 0.60632 xilo - inter == 0 0.2412 0.3401 0.709 0.89176 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 (Adjusted p values reported -- single-step method) > ##Resultados > tapply(trat1_90antes$ocorr, trat1_90antes$grupo,mean, na.rm = TRUE)##Media dos tratamentos ceifa humi inter xilo 4.500000 1.500000 1.833333 2.333333 > tapply(trat1_90antes$ocorr,trat1_90antes$grupo,sd,na.rm=TRUE)/sqrt(tapply(trat1_90antes$ocorr,trat1_90antes$grupo,length))###Erros padrão da media ceifa humi inter xilo 0.6708204 0.5000000 0.5426274 0.3333333 *Aqui tornou-se difícil saber a diferença entre os grupos funcionais em termos de atribui a, b, c, as médias, então resolvi fazer uma análise gráfica usando o pacote HH:* > ##Visualizando as comparações > contrast.1<- + if.R(r=glht(mod2, linfct=mcp(grupo="Tukey")), + s=multicomp(mod2, plot=FALSE)) > plot(contrast.1) *Que me deu um gráfico pouco informativo, porém seguindo o help do mmc {HH}, encontrei que o gráfico de comparações deveria ser feito usando as funções glht.mmc e multicomp.mmc, que requerem o uso de aov(), então fiz:* > > ## Contraste por mmc > contrast.aov <- aov(ocorr ~ grupo, data=trat1_90antes) > contrast.mmc <- + if.R(r=glht.mmc(contrast.aov, linfct=mcp(grupo="Tukey")), + s=multicomp.mmc(contrast.aov, plot=FALSE)) > plot(contrast.mmc) > # *Pois bem, agora tenho a comparação entre os grupos, mas não utilizando a distribuição de erros quasipoisson e não sei o que fazer, pois o pacote HH não me permite utilizar os objetos criados usando glm(), para se fazer as comparações múltiplas. Alguém do grupo teria alguma sugestão, Obrigado, Alexandre* -- Alexandre dos Santos Engenheiro Florestal, Dr. Universidade Federal de Lavras Departamento de Entomologia Laboratório de Entomologia Florestal Caixa Postal 3037 37200-000 - Lavras/MG Fone: +55 (35) 9223-0304

2 2

Re: [R-br] Teste de médias em GLMM com covariável!
by Ivan Bezerra Allaman 10 Mai '12

10 Mai '12

Valeu Walmes, era exatamente isso que eu procurava. Agora entendi. Muito obrigado! \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature} ________________________________ De: Walmes Zeviani <walmeszeviani(a)gmail.com> Para: Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> Enviadas: Quinta-feira, 10 de Maio de 2012 0:14 Assunto: Re: [R-br] Teste de médias em GLMM com covariável! Ivan, Bem, de fato não estou seguro sobre ter entendido o que você quer, mas vamos lá. Partindo do fato que o nosso modelo aponta interação entre F e x nos teremos retas não paralelas (na escala do preditor linear), portando diferenças de médias não são constantes a medida que variamos o valor de x. Podemos comparar as médias em um valor fixo de x. Por exemplo, podemos comparar o peso final de suínos do sexo macho vs fêmea (F) supondo que os animais tem mesma idade (x=100 dias). Outra possibilidade, que deve ser a que você quer, é comparar nas médias amostrais de x nos níveis de F. Assim, comparamos o peso final dos machos com idade média de 100 dias com as fêmeas com idade média final de 95 dias. Se isso é uma comparação justa ou não é outra história e depende do contexto. Eu fiz o CMR desse último caso. É isso que você procura? Eu prefiro sempre que possível evitar as comparações (principalmente seguidas de letras) e dar um gráfico com curvas ajustadas acompanhadas das bandas de confiança. Assim o leitor pode ver como são as diferenças em cada nível de x, entender a relação de x com F e apontar para quais intervalos em x existe ou se pronuncia a diferença entre os níveis de F. É sem dúvida mais informativo. Espero ter ajudado. pred <- expand.grid(F=levels(da$F), x=seq(0,1,l=30)) pred$eta <- predict(ac0, newdata=pred) pred$y <- predict(ac0, newdata=pred, type="response") require(lattice); require(latticeExtra) # retas não paralelas. em que ponto vou compara-las? # em cada lugar dá uma conclusão diferente! xyplot(eta~x, groups=F, data=pred, type="l") # você pode comparar todas em um valor fixo, como x=0.5 # ou pode comparar ela no valor médio de x para cada nível de F xm <- with(da, tapply(x, F, mean)); xm mpm==0.5 mpm2 <- mpm for(i in 1:length(xm)) mpm2[i,mpm2[i,]==0.5] <- xm[i] mpm2 mpm2%*%betas exp(mpm2%*%betas) # diferença nas médias populacionais marginais duas à duas nlevels(da$F) comp <- t(combn(1:nlevels(da$F), 2)) rownames(comp) <- apply(comp, 1, paste, collapse="-") dmpm <- t(apply(comp, 1, function(vs) mpm2[vs[1],]-mpm2[vs[2],])) rownames(dmpm) <- rownames(comp) apply(dmpm, 1, function(x) x%*%betas) require(multcomp) summary(glht(ac0, linfct=mpm2)) # médias populacionais marginais summary(glht(ac0, linfct=dmpm)) # contraste entre as médias populacionais marginais xyplot(y~x, groups=F, data=da, type="p")+ as.layer(xyplot(y~x, groups=F, data=pred, type="l")) À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

Teste de médias em GLMM com covariável!
by Ivan Bezerra Allaman 09 Mai '12

09 Mai '12

Boa tarde senhores! Estou a alguns dias pesquisando em como fazer um teste de comparações em modelo com covariáveis e cujos slopes são diferentes significativamente, já que, as comparações devem ser feitas com as médias ajustadas para as covariáveis. Já procurei no nablle e pela internet, mais nada esclarecedor. Se alguém já trabalhou com este tipo de problema e que queira me dar uma dica, principalmente no que diz respeito a algum material que explique filosoficamente e matricialmente como funciona as comparações múltiplas neste casos de covariável e ficarei extremamente grato. Segue um CRM do que eu vos falo. set.seed(1) da <- data.frame(Trat=gl(5,10),covi=sample(rpois(100,30),50,replace=TRUE),vari=sample(rpois(100,60),50,replace=TRUE),animal=seq(1:50)) da library(lme4) ancova <- glmer(vari ~ Trat*covi + (1|animal%in%Trat),data=da,family=poisson) Anova(ancova)#A interação mostra que as diferenças de slopes entre Trat são significativas ancova2 <- update(ancova,~.-Trat:covi) anova(ancova2,ancova) library(multcomp) summary(glht(ancova,linfct=mcp(Trat='Tukey'),interaction_average=TRUE,covariate_average=TRUE))#Acredito que isto é feito sem correção para as covariáveis. Em algumas buscas, encontrei que a função effect da library ``effects'' faz tal ajutes, mais não estou seguro quanto ao seu uso. Por isso, com o entendimento de como funciona isto matricialmente e acho que consigo resolver via 'glht'. \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

4 5

Bandas de confiança em plot de dispersão
by Marcelo Cardoso mello 09 Mai '12

09 Mai '12

Pessoal, gostaria de adicionar bandas de confiança na CMR que segue. Agradeço qualquer dica. x<-structure(list(week = c(9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L, 9L, 10L, 11L, 12L), kid = c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 13L, 13L, 13L, 13L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, 16L, 17L, 17L, 17L, 17L), treatment = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L ), gain = c(1.2, 1, 1.1, 1.3, 1.2, 1.1, 1.4, 1.5, 1.3, 1.4, 1.4, 1.6, 1.1, 1.1, 1.2, 1.3, 1.2, 1.3, 1.2, 1.3, 1.1, 1.1, 1.1, 1.2, 1.1, 1.2, 1.3, 1.5, 1.3, 1.3, 1.3, 1.4, 1.2, 1.5, 1.9, 2.1, 1.3, 1.2, 1.4, 1.7, 1.5, 1.7, 1.6, 1.7, 1.4, 1.5, 1.7, 1.8, 1.2, 1.2, 1.4, 1.6, 1, 1.1, 1.4, 1.5, 1.4, 1.8, 2.1, 2.1, 1.1, 1.3, 1.4, 1.8, 1.2, 1.5, 1.7, 1.9)), .Names = c("week", "kid", "treatment", "gain"), class = "data.frame", row.names = c(NA, -68L)) y<-aggregate(.~week+treatment, x, FUN=mean) s1<-subset(y, treatment==1) s2<-subset(y, treatment==2) plot(gain~week, data=s1, col=1, type="o", ylim=c(1,2), lty=1) points(gain~week, data=s2, col=2, type="o", lty=2) legend(9.5,1.7, c("t1","t2"), col=c(1,2), lty=c(1,2))

3 2

Dúvida urgente sobre array
by Victor Eduardo 09 Mai '12

09 Mai '12

Galera, queria tirar uma dúvida sobre os arrays, isso está me dando dor de cabeça para realizar uma função para gerar um modelo. Um exemplo seria: > ar_K<-array(1,dim=c(1,2,2)) > ar_K , , 1 [,1] [,2] [1,] 1 1 , , 2 [,1] [,2] [1,] 1 1 > ar_K[,,1] [1] 1 1 > is.matrix(ar_K[,,1]) [1] FALSE > is.vector(ar_K[,,1]) [1] TRUE Queria selecionar uma matriz desse array, mas sempre vem como vetor. Como posso selecionar só a matriz dentro do array? Abraços!!!

2 3

Label=T em barplot
by Diego Leal 09 Mai '12

09 Mai '12

Boa noite pessoal, estou usando a funcao barplot() para gerar um grafico. Existe algum parametro equivalente ao label=T da funcao hist() que eu possa usar e que mostre os valores plotados? esp=c(0,12,30,5,2,1,0) barplot(rev(esp),horiz=TRUE,names.arg=(length(esp)-1):0) gostaria que os valores 0, 12, 30 ... aparececem rotulando as barras sem ter que faze isso manualmente. Att, Diego Leal Togni

3 2

Duvida sobre a função query() do pacote seqinr
by Augusto Ribas 09 Mai '12

09 Mai '12

Ola pessoal, eu estou tentando procurar e extrair algumas sequências de DNA pelo R Em um livro que estou lendo fui guiado até o pacote seqinr, ele tem uma função chamada query() pra fazer isso. A principio tudo parece bem simples e intuitivo. Por exemplo, eu tenho interesse na sequência 16s rRNA de varias especies de sapinhos. Ai se procuro pra uma especie tenho: > library(seqinr) > s <- choosebank("genbank") > query("SP","SP=Rhinella rubescens AND K=16S RIBOSOMAL RNA",s$socket) > SP$nelem [1] 1 > getName(SP) [1] "GU907196.RR1" Blz eu achei 1 sequência especifica pro que eu quero, e o nomezinho pra acessar ela. Procuro outra > query("SP","SP=Leptodactylus syphax AND K=16S RIBOSOMAL RNA",s$socket) > print(SP$nelem) [1] 2 > print(getName(SP)) [1] "JF789923" "JF789924" Blz, 2 sequencias,ai eu fui ao site do genbank, que é a base da dados da uma olhada e estao la blz So que muitas especies que procuro o R para. Por exemplo: > query("SP","SP=Leptodactylus marmoratus AND K=16S RIBOSOMAL RNA",s$socket) O R fica parado e tenho que fechar ele. Eu não entendo o que esta dando errado, principalmente pq algumas especies funcionam e outras não. Essa especie em especifico num era pra achar nada, não achei registro dela procurando pelo site do genbank. Mas a hora que vou usar o R num vai. Outro exemplo é "Leptodactylus syphax" que devia acha 2 sequencias mas ele para também. E a maior parte das especies num vai. Eu uso windows 7 como SO. Estou com o R 2.15.0, internet residencial. Alguem tem experiencia com isso? Sabe dizer o que esto fazendo de errado? E alguem conhece alguma outra forma de fazer o mesmo tipo de busca usando o R? O task view indica o Biocondutor que tem o pacote SRAdb que deve fazer a mesma coisa, mas não consegui usar. Eu vi que o livro que é de 2011 a função mudou a ordem dos argumentos, e teve um update no pacote, mas não sei se estou fazendo algo errado, se tenho que fazer mais alguma coisa. Os manuais que vi não tem um faq com problemas comum e soluções. Minha lista de especies caso alguem queria ver é: lista.spp<-c("Rhinella rubescens", "Rhinella fernandezae", "Elosia rustica", "Crossodactylus gaudichaudii", "Leptodactylus pustulatus", "Leptodactylus syphax", "Dendropsophus cachimbo", "Leptodactylus marmoratus", "Physalaemus soaresi", "Hylodes nasus", "Scinax fuscomarginatus", "Leptodactylus petersii", "Chiasmocleis capixaba", "Ceratophrys aurita", "Pleurodema diplolister", "Cystignatus gigas", "Scinax trilineatus", "Elosis nasus", "Dryadophis bifossatus" ) Se alguem puder dar uma luz, não sei nem por onde começar a procurar pra entender o que ta acontecendo. Vou quebrando a cabeça por enquanto, Boa noite a todos. -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

2 4

Criar Buffer
by Ricardo Dobrovolski 09 Mai '12

09 Mai '12

Olá, Como posso criar um buffer a partir de uma linha? Essa linha foi criada a partir de um par de pontos. Eu preciso desse buffer no formato raster. Abraço, RD.

1 0

Como ler formato dbc no R
by Bernardo Rangel Tura 09 Mai '12

09 Mai '12

Primeiro de tudo bom dia a todos Na minha atividade costumo a usar bancos de dados do Ministério da Saúde. Recentemente precisei trabalhar com todas as internações do país disponíveis para análise. O que dava uma quantidade boa de arquivos pois são 17 anos de arquivos mensais por estado = 5.486 arquivos. Mas isto não chega a ser um problema. O problema é que o arquivo vem num formato de dados .dbc, que até onde consegui descobrir é de DBF compactado utilizado pelo foxpro. Como eu uso linux tenho que emular um windows (wine) utilizar o programa tabwin para transformar dbc em dbf. Como isto o conjunto de dados sai de 10 Gb para 35 Gb. No intuito de simplificar minha vida gostaria de saber de vocês; é possível um mortal que utiliza R 2.15 - atualizado e com praticamente todo e qualquer pacote instalado - rodando Ubuntu 12.04 64 bits num phenom 2 X6 ler diretamente estes arquivos .dbc ? Se for possível o que devo fazer? Antes de mais nada envio o menor arquivo que trabalho como Minimo Exemplo Reprodutível. Ele esta em 2 formatos .dbc e .dbf corresponde ao extrato das internações de Roraima de dezembro de 1997 Certo de sua atenção e genialidade já me despeço agradecendo []s Tura

1 0

Tukey para parcelas subdividad
by Gilson Sanchez 08 Mai '12

08 Mai '12

Pessoal, alguém tem um script para análise de dados em DBC com parcelas subdividas, alias preciso do teste de Tukey (5% de prob). Obrigado -- Bach. Agr.Gilson Sánchez Chia, MSc. Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 / 3303-7919 http://lattes.cnpq.br/3684553005529634 --- Antes de imprimir, pense na sua responsabilidade com o Meio Ambiente. Please consider the environment before printing this email. Renew, Reduce and Recycle. The Planet Knows.*** Salva un Árbol.* * *** * PIENSA**** EN EL MUNDO PARA TUS HIJOS** *

5 6

Contraste
by Adriano Carvalho Costa 08 Mai '12

08 Mai '12

Boa noite, Gostaria de saber se alguém tem um scrip para contraste. Obrigado!

3 2

Criando Arquivo Shape com Linha
by Ricardo Dobrovolski 07 Mai '12

07 Mai '12

Caros R Users, Eu preciso criar um conjunto de arquivos shape, sendo que cada arquivo deve conter como objeto apenas uma linha que conecta dois pontos conhecidos. Desde já, agradeço. Atenciosamente, Ricardo Dobrovolski

1 0

windows.options() - Linux
by Luís Gustavo 07 Mai '12

07 Mai '12

Boa tarde pessoal, Estou com a seguinte dúvida: qual o comando que é similar ao windows.options() no Linux, pois este não existe para o Linux. O objetivo final é criar uma função que me gera a figura da legenda, porém esta função tem rodar tanto no Windows quanto no Linux. Pois gostaria que a figura da legenda se ajuste para os mais variados números de faixas. Abaixo segue um exemplo. (Acredito que função não a mais eficiente, qualquer sugestão é bem-vinda) Obrigado. brks <- c(0,5,10,15,20) r <- 3 cols <- brewer.pal(4,"Blues") gera_figura_legenda=function(brks, cols,r) gera_figura_legenda=function(brks, cols,r) { brks = round(brks,r) dest.fig = tempfile(pattern="legenda",fileext=".png",tmpdir=getwd()) windows.options(width=2.5, height=1.7, pointsize=12, reset = TRUE) CairoPNG(filename = dest.fig, width = 290, height = 290, bg="#FFFFFF00", pointsize = 10) par(bg="transparent", mar=c(0, 0, 0, 0)) brks = round(brks,r) plot(1,1,col="white",axes=F,col.axis="white",xlab="",ylab="") aux1 = 2:(length(brks)-1) aux2 = expand.grid(aux1,aux1) d = aux2[,1]-aux2[,2] aux3 = aux2[d==-1,] legend( "center", title = "Legenda", c(paste("<",brks[2]), paste(brks[aux3[,1]],brks[aux3[,2]],sep="-"), paste(">",brks[(length(brks)-1)])), fill=c(cols), cex=1.3, bg="#FFFFFF") dev.off() return (dest.fig) } -- Luís Gustavo Silva e Silva

5 5

Favor desconsidera o email anterior com inscrição RBRAS.
by Carlos Pereira Consultoria Estatística 07 Mai '12

07 Mai '12

Favor desconsidera o email anterior com inscrição RBRAS. Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MG Cel: 035 - 9119 - 2245 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

1 0

Dúvida no uso de 2 versões do R
by Marco Pollo 06 Mai '12

06 Mai '12

Caros, Utilizo atualmente o R versão 2.14.1 na plataforma windows 64bits. No entanto, preciso utilizar um pacote que só roda em versões anteriores à versão 2.10("Erro: package ‘esamp2’ was built before R 2.10.0: please re-install it"). Posso utilizar duas versões diferentes sem prejuízo do uso que faço da versão mais atual? atencipo meus agradecimentos. Marco Pollo Almeida

2 6

Dúvida básica regressão
by Peterson Gardin 04 Mai '12

04 Mai '12

Pessoal Estou orientando um TCC de uma aluna que não é da minha área. Eu quando trabalho com regressão só uso a polinomial, sou agrônomo. Ela está fazendo um TCC que usa doses (0,1, 1, 10 ,100, 1000) de fungicidas (5 marcas de fungicidas) e 5 repetições. Devo extrair o log10 de doses e analisar por regressão linear?? as doses ficariam -1, 0,1, 2, 3 ou devo usar um modelo com log?? Não conheço esses modelos. Obrigado desde já. João Peterson Pereira Gardin Doutor em Fisiologia Vegetal, Eng. Agr., Pesquisador Fones: (49) 3566-0054 e 9933-4290 MSN: petersongardin(a)yahoo.com.br - SKYPE: petersongardin Endereço: Rua João Zardo, 1660, Bairro Campo Experimental, Videira - SC - CEP 89560000

3 3

Loop varrendo quantidade enorme de arquivos
by Thiago Veloso 03 Mai '12

03 Mai '12

Boa noite pessoal, Preciso de ajuda para fazer um loop que percorra uma lista enorme de arquivos. Organizei uma série (+50GB) de imagens de satélite que está organizada em anos (2001 a 2011) e, por ano, a cada oito dias (julianos). lai200101.tif lai200109.tif lai200117.tif lai200125.tif . . . lai2011361.tif onde a nomenclatura segue o padrão: lai{ano}{diajuliano}.tif. Percebam que o incremento entre os dias é sempre oito. Preciso fazer duas operações nesse banco de dados: 1) um loop para abrir cada um dos arquivos, nomeando a variável com o ano e mês, e multiplicar a mesma por 0.001. Algo assim: library(raster) lai{ano}{diajuliano}<-raster('lai{ano}{diajuliano}.tif') #carregando imagem lai{ano}{diajuliano}<-lai{ano}{diajuliano}*0.001 #calibrando imagem 2) Ao longo de todos os arquivos, somar cada quatro para produzir um mês de dados, por exemplo: lai_jan_2001 <- lai200101.tif + lai200109.tif + lai200117.tif + lai200125.tif lai_feb_2001 <- lai200133.tif + lai200141.tif + lai200149.tif + lai200157.tif . . . lai_dec_2011 <- lai2011337.tif + lai2011345.tif + lai2011353.tif + lai2011361.tif Um verdadeiro desafio de programação, concordam? Alguém pode me dar uma ideia de como começar? Agradece desde já, Thiago.

4 5

Predição de dados - Modelagem
by Sérgio Henrique almeida da silva ju 03 Mai '12

03 Mai '12

Olá Rodei um modelo: reg1<-lm(est ~ e1 + f1 + g1 + i1 + k1 + m1 + o1 + r1 + s1 + t1, data = dados) summary (reg1)$coefficients Estimate Std. Error t value Pr(>|t|) (Intercept) 33.811287 0.7983445 42.351753 1.557462e-202 e1 -3.300681 1.0234877 -3.224934 1.313683e-03 f1 -1.602262 0.6537985 -2.450697 1.448143e-02 g1 -1.337623 0.6721036 -1.990203 4.692499e-02 i1 -6.224414 0.9810556 -6.344608 3.820429e-10 k1 -1.279617 0.5544328 -2.307975 2.126718e-02 m1 -1.569079 0.5133349 -3.056638 2.316421e-03 o1 2.984603 1.6371984 1.822994 6.869551e-02 r1 1.340099 0.6338281 2.114295 3.481391e-02 s1 1.598583 0.4717877 3.388352 7.392205e-04 t1 6.286149 0.9227351 6.812518 1.946057e-11 E gostaria de criar uma variável em um outro banco (dados2) com a estimação da minha variável gerada pelo meu modelo, estou usando o seguinte comando, porém dá um erro. dados2$estim=predict(idfu.reg1, dados2, se.fit = F) Error in `$<-.data.frame`(`*tmp*`, "esqui", value = list(fit = c(37.4586673781906, : replacement has 4 rows, data has 5565 Como posso fazer isso? -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

4 4

Sobre distância euclidiana
by Marcio Leles 03 Mai '12

03 Mai '12

Calculo de distância euclidiana Bom dia. Tenho um arquivo que ajusto um modelo de regressão gerando três betas para cada árvore. Tenho casos de ter a medição de mil árvores. Após fazer os ajustes passo a ter beta 0, beta 1 e beta 2 para cada árvore. Gostaria de calcular a distância euclidiana do beta 0 com o beta 0, beta 1 com o beta 1 e beta 2 com o beta 2, entre todas as árvores. O meu objetivo é ter uma matriz de distância euclidiana entre todos os betas de todas as árvores para saber por exemplo qual árvore é mais próxima da árvore de número 21. Abaixo um pedaço do arquivo. Obrigado pela ajuda. ARV y x1 x2 coefb_0 coefb_1 coefb_2 1 2,631605 0,011628 0,000135 2,443489 -7,54249 6,537685 1 1 0,151163 0,02285 2,443489 -7,54249 6,537685 1 0,830123 0,267442 0,071525 2,443489 -7,54249 6,537685 2 1,531406 0,007576 5,74E-05 1,341076 -2,17855 0,89864 2 1 0,098485 0,009699 1,341076 -2,17855 0,89864 2 0,832656 0,174242 0,03036 1,341076 -2,17855 0,89864 3 1,645648 0,007194 5,18E-05 1,397167 -2,65584 1,382927 3 1 0,093525 0,008747 1,397167 -2,65584 1,382927 3 0,808183 0,165468 0,02738 1,397167 -2,65584 1,382927 4 2,08642 0,007042 4,96E-05 1,652506 -3,76968 2,395018 4 1 0,091549 0,008381 1,652506 -3,76968 2,395018 4 0,840278 0,161972 0,026235 1,652506 -3,76968 2,395018 -- Atenciosamente; Marcio Leles

3 2

Write.table Posicional
by Vinicius Brito Rocha 03 Mai '12

03 Mai '12

Senhores, Estou precisando exportar uma tabela onde não há delimitador e as variáveis são determinadas através de posições específicas.. Ex: Var1 (ocupa a posição 1-10 ) Var2 (ocupa a posição 11-31) e assim por diante. A função write table permite fazer algo assim ou apenas exporta com delimitador? Abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

3 3

Problemas com R 2.15
by Tito Conte 03 Mai '12

03 Mai '12

Acabei de instalar o R para Ubuntu e não estou conseguindo instalar nenhum pacote novo, fiz isso para o pacote ncdf, e tentei pra mais alguns e nada todos dão o mesmo erro (achava que era problema com o NetCDF) a mensagem que aparece segue abaixo installing *source* package ‘ncdf’ ... ** package ‘ncdf’ successfully unpacked and MD5 sums checked checking for nc-config... /usr/bin/nc-config configure: creating ./config.status config.status: creating src/Makevars ** libs gcc -std=gnu99 -I/usr/share/R/include -DNDEBUG -I/usr/include -fpic -O3 -pipe -g -c ncdf.c -o ncdf.o ncdf.c:3:20: fatal error: netcdf.h: No such file or directory compilation terminated. make: *** [ncdf.o] Error 1 ERROR: compilation failed for package ‘ncdf’ * removing ‘/home/tito/R/i686-pc-linux-gnu-library/2.15/ncdf’ The downloaded source packages are in ‘/tmp/Rtmpn8hiQ4/downloaded_packages’ Warning message: In install.packages("ncdf") : installation of package ‘ncdf’ had non-zero exit status Obrigado Tito

4 4

Loop varrendo quantidade enorme de arquivos
by Thiago Veloso 02 Mai '12

02 Mai '12

Boa noite pessoal, Preciso de ajuda para fazer um loop que percorra uma lista enorme de arquivos. Organizei uma série (+50GB) de imagens de satélite que está organizada em anos (2001 a 2011) e, por ano, a cada oito dias (julianos). lai200101.tif lai200109.tif lai200117.tif lai200125.tif . . . lai2011361.tif onde a nomenclatura segue o padrão: lai{ano}{diajuliano}.tif. Percebam que o incremento entre os dias é sempre oito. Preciso fazer duas operações nesse banco de dados: 1) um loop para abrir cada um dos arquivos, nomeando a variável com o ano e mês, e multiplicar a mesma por 0.001. Algo assim: library(raster) lai{ano}{diajuliano}<-raster('lai{ano}{diajuliano}.tif') #carregando imagem lai{ano}{diajuliano}<-lai{ano}{diajuliano}*0.001 #calibrando imagem 2) Ao longo de todos os arquivos, somar cada quatro para produzir um mês de dados, por exemplo: lai_jan_2001 <- lai200101.tif + lai200109.tif + lai200117.tif + lai200125.tif lai_feb_2001 <- lai200133.tif + lai200141.tif + lai200149.tif + lai200157.tif . . . lai_dec_2011 <- lai2011337.tif + lai2011345.tif + lai2011353.tif + lai2011361.tif Um verdadeiro desafio de programação, concordam? Alguém pode me dar uma ideia de como começar? Agradece desde já, Thiago.

1 0

precisão numérica e which.max com empates
by Pedro Emmanuel Alvarenga Americano do Brasil 02 Mai '12

02 Mai '12

Amigos de R, Estou tentando comparar alguns métodos de estimação e um dos passos necessários é encontrar valores maximos ou minimos de vetores numéricos. O problema numero 1 é que quando comparo visualmente os vetores ocorre que ocasionalmente acontece 2 ou 3 empates nos valores máximos ou mínimos como no vetor abaixo. Reparem que o valor .7 aparece na posição 21 e 23. Mas quando utilizo a função which.max ou max eu so encontro a posição 23. POssivlemente ha alguma coisa relacionado com precisão nume´rica apesar de não fazer muito sentido no momento. > test.values$Youden [1] 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 0.55 0.60 0.65 0.60 0.65 0.60 0.65 [21] 0.70 0.65 0.70 0.65 0.60 0.65 0.60 0.65 0.60 0.55 0.50 0.45 0.40 0.35 0.30 0.25 0.20 0.15 0.10 0.05> max(test.values$Youden)[1] 0.7 > which.max(test.values$Youden)[1] 23 > which(test.values$Youden==max(test.values$Youden)) [1] 23 > test.values$Youden[21]==test.values$Youden[23] [1] FALSE O problema 2 é que quando de fato ocorrem empates a função which.max retorna somente a primeira posição. Então, a questão seria como retornar todas as posições com valores máximos ou encontrar uma mediana dos valores empatados? Qualquer luz é bem vinda. Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

3 4

Retirar coordenadas de cidade do google maps com R
by Augusto Ribas 02 Mai '12

02 Mai '12

Opa a algum tempo postaram uma função aqui pra calcula a distancia entre cidades pelo googe maps. Eu gostaria de pegar a coordenada da cidade. Mas eu não acho onde esta essa informação na pagina do google depois que abro ela no R E também nao sei procurar pra falar a verdade, pq fica gigante depois que le a pagina. Alguém pode me dar uma luz? #Função anteriormente postada distancia <- function(cidade1, uf1, cidade2, uf2){ link <- paste( "http://maps.google.com/maps?saddr=", cidade1, "+-+", uf1, ",+Brasil&daddr=", cidade2, "+-+", uf2, ",+Brasil&hl=pt-BR&ie=UTF8", sep="") link <- iconv(link, 'latin1', 'UTF-8') Encoding(link) <- 'bytes' pagina <- readLines(url(link)) linha <- grep('<div class="altroute-rcol altroute-info"> <span>', pagina) distancia <- pagina[linha] distancia <- strsplit(distancia, '<div class="altroute-rcol altroute-info"> <span>')[[1]][2] distancia <- strsplit(distancia, '</span>')[[1]][1] return(list(distancia = distancia, link=link)) } distancia("Dourados", "MS", "Campo.Grande", "MS") distancia("Brasília", "DF", "Goiânia", "GO") distancia("Brasília", "DF", "Natal", "RN") distancia("Brasília", "DF", "Florianópolis", "SC") distancia("Brasília", "DF", "Porto.Alegre", "RS") #O que eu gostaria por exemplo para a cidade de Dourados MS uf1<-c("MS") cidade1<-c("Dourados") link<-paste( "http://maps.google.com/maps?saddr=%20 ",cidade1,"+-+",uf1,",+Brasil&hl=pt-BR&ie=UTF8",sep="") Encoding(link) <- 'bytes' pagina <- readLines(url(link)) #é saber onde em pagina esta as coordenadas da cidade. #que linha olhar, se é que ela esta nessa pagina, e colocar como o resultado acima algo tipo #coordenadas são tal -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056

1 1

Ajuda com Correlação de Matrizes com tamanhos diferentes
by Débora Martins 02 Mai '12

02 Mai '12

Olá pessoal, Estou analisando alguns dados para um artigo, entretanto não consigo rodar a correlação do jeito que desejo. Meu problema é um pouco parecido com o descrito nesse link http://stats.stackexchange.com/questions/24980/correlation-between-matrices…, pois eu desejo obter um único valor que corresponda a correlação, mas como eu divido os valores do arquivo em algumas dimensões que me interesso, não consigo obter esse valor único. Eu li o arquivo,separo cada coluna na dimensão que desejo analisar e tento fazer a correlação entre as dimensões mas só obtenho a correlação dos itens. teste<-read.table("1_analise.txt", header=T) a<-teste[1] b<-teste[2] c<-teste[3] d<-teste[4] e<-teste[5] estrutura<-cbind(a,b,c,d,e) a<-teste[6] b<-teste[7] c<-teste[8] d<-teste[9] e<-teste[10] f<-teste[11] g<-teste[12] cultura<-cbind(a,b,c,d,e,f,g) correlacao1<-cor(estrutura,cultura) Com a última linha do código, gostaria de obter algo como a tabela 3 deste link http://www.tlainc.com/articl292.htm. -- *Débora de Vasconcelos Martins** *Graduanda em Ciência da Computação Universidade Estadual do Ceará Bolsista do Programa de Educação Tutorial - PET http://petcomp.uece.br

2 3

Diferença entre krigagem ordinária e krigagem externa a deriva
by Mauro Sznelwar 02 Mai '12

02 Mai '12

Quem sabe explicar em poucas palavras a diferença entre elas?

2 1

Duvida Regressão
by Sérgio Henrique almeida da silva ju 30 Abr '12

30 Abr '12

Olá Estou trabalhando com a taxa de uma doença e gostaria de rodar um modelo de regressão. Meus dados não apresentaram distribuição normal e por teste de aderência deu uma distribuição lognomal -> fitdistr(tx_esqui, "lognormal") meanlog sdlog 6.25725753 1.90754231 (0.11630491) (0.08223999) ->ks.test(tx_esqui, "plnorm", sdlog=1.90754231, meanlog=6.25725753,alternative=c("two.sided")) One-sample Kolmogorov-Smirnov test data: tx_esqui D = 0.0701, p-value = 0.1419 alternative hypothesis: two-sided Tem como eu rodar um modelo de regressão com essa variável? Como fazê-lo? Qual a melhor forma de trabalhar com essa variável? Abraços Segue a minha variável dependente: structure(list(tx_esqui = c(3L, 3L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 9L, 10L, 10L, 19L, 22L, 22L, 23L, 23L, 25L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L, 34L, 36L, 39L, 41L, 42L, 45L, 46L, 49L, 49L, 56L, 57L, 58L, 58L, 59L, 61L, 62L, 62L, 68L, 70L, 74L, 76L, 80L, 82L, 82L, 83L, 83L, 86L, 101L, 102L, 107L, 107L, 112L, 113L, 118L, 120L, 121L, 123L, 126L, 135L, 138L, 140L, 153L, 156L, 158L, 179L, 191L, 200L, 201L, 204L, 207L, 209L, 212L, 215L, 223L, 227L, 228L, 233L, 247L, 248L, 249L, 255L, 257L, 258L, 258L, 259L, 264L, 274L, 277L, 279L, 287L, 297L, 309L, 310L, 313L, 314L, 326L, 328L, 329L, 331L, 332L, 332L, 340L, 354L, 358L, 365L, 372L, 379L, 383L, 395L, 396L, 418L, 423L, 447L, 449L, 451L, 453L, 493L, 532L, 535L, 553L, 554L, 565L, 573L, 581L, 584L, 599L, 603L, 604L, 611L, 672L, 703L, 708L, 727L, 756L, 799L, 809L, 812L, 815L, 816L, 862L, 864L, 866L, 932L, 946L, 951L, 964L, 969L, 983L, 985L, 989L, 1032L, 1102L, 1102L, 1112L, 1112L, 1199L, 1218L, 1245L, 1252L, 1283L, 1302L, 1327L, 1327L, 1328L, 1340L, 1367L, 1382L, 1413L, 1441L, 1468L, 1487L, 1493L, 1512L, 1523L, 1567L, 1576L, 1585L, 1598L, 1620L, 1623L, 1634L, 1664L, 1704L, 1813L, 1816L, 1824L, 1842L, 1892L, 1939L, 1977L, 2010L, 2035L, 2060L, 2100L, 2152L, 2184L, 2210L, 2248L, 2282L, 2288L, 2358L, 2462L, 2473L, 2527L, 2541L, 2555L, 2576L, 2603L, 2705L, 2711L, 2895L, 2916L, 2958L, 2963L, 3088L, 3110L, 3304L, 3427L, 3489L, 3535L, 3629L, 3638L, 3682L, 3704L, 3796L, 3932L, 4160L, 4293L, 4294L, 4366L, 4516L, 4715L, 4746L, 5152L, 5263L, 5307L, 5823L, 6025L, 6306L, 6313L, 6330L, 6649L, 7114L, 7158L, 7340L, 7488L, 7665L, 8055L, 8394L, 8773L, 8850L, 8918L, 9273L, 9356L, 10110L, 10281L, 10310L, 10519L, 11111L, 13525L, 14465L, 14955L, 16629L)), .Names = "tx_esqui", class = "data.frame", row.names = c(NA, -269L)) -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

7 16

Contagem das linhas de uma matriz
by Diego Leal 30 Abr '12

30 Abr '12

Boa tarde pessoal, a questão é a seguinte: queria saber se existe algum comando equivalente ao comando table() para tabelar os valores das linhas de uma matriz. Por exemplo: matriz=matrix(rep(1:20),10,10,byrow=TRUE) funcaoDesconhecida(mat) # essa funcao acima deveria me retornar falando que a linha (1,2,...,10) aparece 5 vezes e a linha (11,12,13,...,20) aparece 5 vezes. alguma ideia? Att, Diego Leal Togni

2 3

Tutorial do WinBUGS no R
by Gilbert Queiroz 30 Abr '12

30 Abr '12

Alguém tem, ou sabe onde tem, tutorial de uso do WinBUGS no R? Abs. Gilbert.

2 1

Pacote GPCM
by Felipe Buchbinder 30 Abr '12

30 Abr '12

Prezados, Alguém sabe onde eu posso encontrar o pacote gpcm para instalá-lo? Eu li sobre ele em um artigo, mas não estou conseguindo encontrá-lo nem no site do CRAN nem através do comando "Install packages" do R. Desde já, obrigado, Felipe

2 4

Ajuda em função para criação do Filtro De Kalman
by Victor Eduardo 30 Abr '12

30 Abr '12

Galera, estou criando uma função para realizar o Filtro de Kalman sendo que em certo momento, preciso transpor uma matriz, só que essa matriz está dentro de um array, logo não consigo realizar esta operação. Está dentro de um for, vou colar aqui a parte que está dando erro: ar_F[,,t]<-ar_Z[,,t]%*%ar_P[,,t]%*%t(ar_Z[,,t])+ar_H[,,t] Acredito que o erro esteja justamente em t(ar_Z[,,t]), pois é um array, e não uma matriz. Como poderia fazer isso? Abraços

2 3

sobre distribuções tempered stable
by Marcelo Zeuli 30 Abr '12

30 Abr '12

Bom dia, Gostaria de saber se a pacote do R que implemente distribuições do tipo tempered stable. Grato, Marcelo Zeuli

2 1

Perda de Parcela.
by Carlos Pereira 30 Abr '12

30 Abr '12

Boa tarde queria saber se poderiam me ajudar, estou com uma duvida de como fazer ANOVA, quando existe a perda de parcela, como posso fazer para que o R estime essa parcela perdida para realizar ANOVA. Em um delineamento em Blocos ao Acaso. Desde já obrigado pela atenção! Carlos Pereira da Silva Bacharel em Estatística - UEPB Campina Grande Cel: 083 - 9672- 1966 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

4 3

Manipulação de matrizes
by Diego Leal 29 Abr '12

29 Abr '12

Boa tarde, estou com um problema ao manipular matrizes. Tenho uma matriz e quero encontrar a linha de um vetor dentro desta matriz. Por exemplo: Suponha esta matriz: mat=matrix(1:25, 5, 5, byrow=TRUE) [,1] [,2] [,3] [,4] [,5] [1,] 1 2 3 4 5 [2,] 6 7 8 9 10 [3,] 11 12 13 14 15 [4,] 16 17 18 19 20 [5,] 21 22 23 24 25 preciso saber se existe alguma funcao q eu passe um vetor c(6,7,8,9,10) e ele me retorne a linha 2. Ex: quando eu rodar o seguinte comando, a funcao me retorne 2 funcao_desconhecida(mat,c(6,7,8,9,10)) Att, Diego Leal Togni

2 2

R_HOME
by Leandro Marino 29 Abr '12

29 Abr '12

*Caros,* * * *estou tentando configurar o StatET no Ubuntu no entanto ele solicita o caminho para o R_HOME, qual é este caminho no linux? Lembrando que fiz a instalação via apt-get e não compilando das fontes.* * * *Um abraço,* *Leandro* * * * *Atenciosamente, Leandro Marino Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

4 6

uso do knitr com Lyx
by Cleber N.Borges 28 Abr '12

28 Abr '12

Olá a todos, Alguém tem experiência de usar o knitr com Lyx? Estou tentando usá-lo, porém é só abrir um documento em branco no Lyx e pedir para carregar o módulo knitr, digitar uma letra que seja e compilar que aparece a seguinte mensagem de erro: --------------------------------------------- File ended while scanning use of \@xdblarg. <*> pls_detail.tex *** (job aborted, no legal \end found) --------------------------------------------- Outro ponto é que consigo compilar os exemplos baixados da internet, então eu sei que tenho os pacotes necessários e que o R é encontrado no path sem problemas. Cao alguém puder ajudar, agradeço antecipadamente. Obrigado. Cleber

2 1

off toppic: Extensão DWG
by Leonard de Assis 27 Abr '12

27 Abr '12

Alguém sabe qual software (Gratuíto e multiplataforma de preferência) que consiga abrir a extensão DWG 1. Autocad nem CorelCAD contam 2. Preciso abrir isto em uma máquina linux (Ubuntu 11.10 e em outra windows -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

4 3

Shape file + Seta
by Vinicius Brito Rocha 27 Abr '12

27 Abr '12

Senhores estou plotando um shape file de um Estado e colorindo dois municípios de interesse, queria plotar uma seta apontando de um para o outro. Tentei usar a função arrows em conjunto com a latitude e longitude do shape das cidades, mas não consegui. Alguém tem alguma sugestão? Abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 2

alguem sabe como descarregar uma coleção de SST da http://podaac.jpl.nasa.gov/dataset/JPL-L2P-MODIS_A (NOAA)
by carlos pombo sonderblohm 27 Abr '12

27 Abr '12

olá pessoal, sei que devo fazer os post mas específicos e directos, mas estou mesmo perdido. Gostaria de criar e visulaizar uma colecção de SST para minha área de estudo (TOP = 38, bottom= 35, left= -10, right =-6; sur de Portugal), do parámetro SST do dataset http://podaac.jpl.nasa.gov/dataset/JPL-L2P-MODIS_A e possivel construir o query da informação no R, p.e.: SST Mensais desde 2006 ate 2012? em que formato e armazenada a collecção? logo posso usar no meu GIS (QGIS)? alguma dica, desculpem ser tão impreciso, obrigado Carlos -- Carlos A. Pombo Sonderblohm PhD Student on Marine Science (Fisheries) Faculdade de Ciências e Tecnología Universidade do Algarve, Campus de Gambelas 8005-139 Faro Portugal Tef. 289 800 905 ext. 7605

2 1

ANOVA/MANOVA
by T Branquinho 26 Abr '12

26 Abr '12

Caros colegas, Estou aprendendo agora a usar o R. Então gostaria de saber se é possível, e se for possível como fazer? Utilizando outro software para fazer regressão linear múltipla e redes neurais para predição de atividade. O programa me informa o coeficiente de correlação, mean absolute error, root mean squared error, relative absolute error, root relative squared error. É possível, utilizando o valor real e o valor predito, para efetuar a MANOVA pelo R (teste de F?)? Agradeço antecipadamente,TBO

2 1

unzip
by Vinicius Brito Rocha 26 Abr '12

26 Abr '12

Senhores, estou tendo uma mensagem de erro no função unzip. Erro 1. Alguém sabe o que pode ser? -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 2

Função LENGTH
by Carlos Pereira 26 Abr '12

26 Abr '12

Boa tarde tenho uma duvida com respeito ao a função LENGTH essa função dar a dimensão de um vetor em um vetor com 3000 valores onde esses valores são caractere por EXEMPLO ( A, B, C) teria como eu saber quantos A eu teria em no meu vector? Carlos Pereira da Silva Bacharel em Estatística - UEPB Campina Grande Cel: 083 - 9672- 1966 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

4 5

Calcular Raiz elevada a um determinado número t.
by Jodavid Ferreira 26 Abr '12

26 Abr '12

Gostaria de saber com calcular uma raiz , de um determinado número t no R. exemplo: raiz elevada a 4 de 16. -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Responsável pelo Ministério Jovem e de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com

4 3

Séries Temporais no R
by Sérgio Henrique almeida da silva ju 25 Abr '12

25 Abr '12

Olá Estou com um conjunto de dados de Hanseníase para os anos de 2001 a 2009 e gostaria de fazer algumas análises de séries temporais. Vocês indicariam referências e tutorial sobre o tema no R? Abraços Em 18 de abril de 2012 12:16, Benilton Carvalho <beniltoncarvalho(a)gmail.com>escreveu: > talvez > > f(dados1$linhaa, "A39|G00") > > mas e' apenas achometro. > > boa sorte. > > 2012/4/18 Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com>: > > Se por exemplo eu busca além do código A39 querer também o G00 > > > > dados1$causaa = f(dados1$linhaa, "A39") > > dados1$causab = f(dados1$linhab, "A39") > > dados1$causac = f(dados1$linhac, "A39") > > dados1$causad = f(dados1$linhad, "A39") > > dados1$causaii = f(dados1$linhaii, "A39") > > dados1$causaa1 = f(dados1$linhaa, "G00") > > dados1$causab1 = f(dados1$linhab, "G00") > > dados1$causac1 = f(dados1$linhac, "G00") > > dados1$causad1 = f(dados1$linhad, "G00") > > dados1$causaii1 = f(dados1$linhaii, "G00") > > E depois agregar essas colunas em uma só! > > > > Tem como eu buscar o código A39 e G00 em um comando só? como por exemplo: > > > > dados1$causaa = f(dados1$linhaa, "A39" | "G00") > > > > Abraços > > > > > > > > Em 18 de abril de 2012 11:56, Benilton Carvalho < > beniltoncarvalho(a)gmail.com> > > escreveu: > > > >> Como sugerido no guia de postagem, sinta-se a vontade para colocar seu > >> arquivo em algum servico como http://www.datafilehost.com/ , assim > >> quem esta' tentando te ajudar tem a chance de entender mais > >> rapidamente qual eh o problema. > >> > >> O que vc quer dizer com "mais de um codigo"? > >> > >> Em teoria, vc precisa apenas de um wrapper para: > >> > >> with(dados1, gsub('NA', '', > >> > >> > paste(dados1$causaa,dados1$causab,dados1$causac,dados1$causad,dados1$causaii, > >> sep=''))) > >> > >> mas sem entender o que exatamente vc quer dizer, minha ajuda acaba por > >> aqui. > >> > >> b > >> > >> 2012/4/18 Sérgio Henrique almeida da silva ju < > sergio.edfisica(a)gmail.com>: > >> > > >> > > >> > ---------- Mensagem encaminhada ---------- > >> > De: Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com> > >> > Data: 18 de abril de 2012 11:22 > >> > Assunto: Re: [R-br] Duvida String > >> > Para: r-br(a)listas.c3sl.ufpr.br > >> > > >> > > >> > Uma alternativa que encontrei foi > >> > > >> > > >> > f = function(v, dat, v1=sapply(gregexpr(dat,v),"[",i=1)){ > >> > n = ifelse(v1 > 1, nchar(v), NA) > >> > substr(v, v1,n) > >> > } > >> > dados1$causaa = f(dados1$linhaa, "A39") > >> > dados1$causab = f(dados1$linhab, "A39") > >> > dados1$causac = f(dados1$linhac, "A39") > >> > dados1$causad = f(dados1$linhad, "A39") > >> > dados1$causaii = f(dados1$linhaii, "A39") > >> > > >> > dados1$causa<- with(dados1, gsub('NA', '', > >> > > >> > > paste(dados1$causaa,dados1$causab,dados1$causac,dados1$causad,dados1$causaii, > >> > sep=''))) > >> > > >> > > >> > Como faço para buscar mais de um código? Tem que fazer um monte de > >> > comando > >> > como esses? > >> > Existe uma forma mais simples de fazer isso? > >> > > >> > Segue o banco anexo > >> > > >> > > >> > Desculpe enviar o arquivo anexo, mas é uma tentativa para vocês me > >> > ajudarem, > >> > já que o dput não é uma boa saida! > >> > > >> > Abraços > >> > > >> > > >> > > >> > > >> > -- > >> > Sérgio Henrique Almeida da Silva Junior > >> > Doutorando em Epidemiologia em Saúde Pública > >> > Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > >> > http://lattes.cnpq.br/1611345552843383 > >> > Tel: (21) 94429486/78101651 id: 123*20942 > >> > > >> > > >> > > >> > > >> > > >> > -- > >> > Sérgio Henrique Almeida da Silva Junior > >> > Doutorando em Epidemiologia em Saúde Pública > >> > Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > >> > http://lattes.cnpq.br/1611345552843383 > >> > Tel: (21) 94429486/78101651 id: 123*20942 > >> > > >> > > >> > > >> > _______________________________________________ > >> > 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. > >> _______________________________________________ > >> 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. > > > > > > > > > > -- > > Sérgio Henrique Almeida da Silva Junior > > Doutorando em Epidemiologia em Saúde Pública > > Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > > http://lattes.cnpq.br/1611345552843383 > > Tel: (21) 94429486/78101651 id: 123*20942 > > > > > > > > _______________________________________________ > > 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. > _______________________________________________ > 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. > -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

2 1

Divergencia entre o valor do 1 quartil
by Daniel Marcelino 25 Abr '12

25 Abr '12

Ola caros, um colega me perguntou no exemplo abaixo porque o R calcula o primeiro quartil como sendo 1.75 e no boxplot ele apresenta como sendo 1.5. q8<-c(1,1,3,2,3,2,3,2) summary(q8) boxplot(q8) summary(q8) Min. 1st Qu. Median Mean 3rd Qu. Max. 1.000 1.750 2.000 2.125 3.000 3.000

4 3

Re: [R-br] Divergencia entre o valor do 1 quartil
by Ivan Bezerra Allaman 25 Abr '12

25 Abr '12

Complementando... quartil1 <- NULL for(i in 1:9){ quartil1[i] <- quantile(q8,0.25,type=i) } quartil1 quartil3 <- NULL for(i in 1:9){ quartil3[i] <- quantile(q8,0.75,type=i) } quartil3 (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0