- R-br - DINF UFPR

Trocar a fonte de expression()
by Leonardo Monteiro 09 Set '14

09 Set '14

Boa noite, eu nao estou conseguindo trocar a fonte do texto que esta dentro do expression() gostaria de por para o estilo 'times new roman', e estava usando font = 6, entretanto sem sucesso. alguem teria uma solucao p isso? a<- 1:10 b<- 3/a plot(a,b) text(6, 2, expression(mg^-1), font = 6) obrigado, Leonardo

2 1

Classificação de pontos que pertencem a uma condição dada (x,y)
by ASANTOS 09 Set '14

09 Set '14

Boa tarde Pessoal, Fiz a representação de 48 parcelas com plantas no espaçamento de 2x4m usando a função expand.grid(), sendo as parcelas espaçadas em 7m (grid1), após criei pontos que presentam o centro das parcelas (grid2) e sobre esses pontos eu adicionei os códigos das parcelas (parc.names). Gostaria de saber se existe uma forma de criar uma classificação onde quero saber a qual código da parcela (parc.names) pertence cada um dos pontos que representam as plantas (grid1). Segue CRM: ##Plantas ---------------------------------------------------------------------- ## Coordenadas do grid no espaçamento 2x4 ## distâncias verticais e horizontais v <- 2; h <- 4 ## número de pontos vertical e horizontal npv <- 4; nph <- 5 ## distância entre vãos vertical e hotizontal V <- 7; H <- 7 ## número de linhas e colunas ncol <- 4; nrow <- 12 grid0 <- expand.grid(L=seq(2,length.out=nrow,by=npv*v+V), C=seq(2,length.out=ncol,by=nph*h+H), x=seq(2,length.out=npv,by=v), y=seq(2,length.out=nph,by=h)) grid1 <- with(grid0, data.frame(x=L+x, y=C+y)) plot(y~x, grid1, asp=1) ## Centro das parcelas ------------------------------------------------------ grid2 <- expand.grid(x1=seq(7,length.out=12,by=15), y1=seq(12,length.out=4,by=27)) points(grid2[,1],grid2[,2],col="red") #Código das parcelas parc.names<-c("2211","1221","3221","1122","2122","3222","1113","2113","3223","1124","2214","3114", "2111","1121","3211","1112","2212","3212","1213","2123","3213","1214","2124","3224", "2221","1211","3111","1212","2222","3112","1223","2213","3113","1224","2224","3214", "2121","1111","3121","1222","2112","3122","1123","2223","3123","1114","2114","3124") text(grid2$x1,grid2$y1,labels=parc.names,cex=1,col=2) # Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 2

fixar sigmasq na funcao likfit do pacote geoR
by Fernanda De Bastiani 09 Set '14

09 Set '14

Olá, estou trabalhando com o pacote geoR e utilizando a função likfit() para estimar os parametros. Percebi que a função disponibilizo a opção fix.nugget para podermos fixar o valor do nugget (efeito pepita). No entanto, gostaria saber se existe a opção de ficar o valor de sigmasq (sill), e usar a função likfit() para estimar somente os outros parâmetros. Além disso, alguém poderia por favor me fornecer algum material que possua algum exemplo pratico desta função (likfit()) utilizando a opção "realisations"? Grata pela atenção, Fernanda

2 2

obtencao do coef. de determinacao
by Leonardo Monteiro 09 Set '14

09 Set '14

Prezado Walmes, desculpa, mas minha eu nao consegui identificar a resposta da minha pergunta na sua resposta... se eu 'forcar' a reta de regressao a passar pelo zero, o R2 devera ser maior ou menor? obrigado, Leonardo. --

2 1

Sequencia de Imagens
by David Feitosa 09 Set '14

09 Set '14

Boa noite! Há algum tempo, eu lembro de ter visto um artigo que unia o recurso do R com Markdown/LaTex. No exemplo, era gerado um pdf com um animação, a partir de uma série de imagens. Recentemente, eu fiz uma implementação de um algoritmo que gera uma série de imagens e mostra uma "evolução" de um processo, mas não me recordo do pacote que gera o pdf com a animação. Alguém já usou este recurso ou sabe como gerar essa animação? Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

3 5

Re: [R-br] Sequencia de Imagens
by Robert Iquiapaza 09 Set '14

09 Set '14

library(animation) Att From: David Feitosa Sent: Monday, September 08, 2014 10:40 PM To: R-br(a)listas.c3sl.ufpr.br Subject: [R-br] Sequencia de Imagens Boa noite! Há algum tempo, eu lembro de ter visto um artigo que unia o recurso do R com Markdown/LaTex. No exemplo, era gerado um pdf com um animação, a partir de uma série de imagens. Recentemente, eu fiz uma implementação de um algoritmo que gera uma série de imagens e mostra uma "evolução" de um processo, mas não me recordo do pacote que gera o pdf com a animação. Alguém já usou este recurso ou sabe como gerar essa animação? Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("") -------------------------------------------------------------------------------- _______________________________________________ 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.

2 1

Leitura de número como caractere
by Vanessa Felix 09 Set '14

09 Set '14

Pessoal boa noite, Gostaria de saber porque estranhamente o R leu 4 colunas dos meus dados (anos 1991,1992,1994,1995) como caracter, sendo que nas outras ele leu normalmente como numérico. Poderiam me ajudar tbm em como converter de caractere para numérico. Pq já aconteceu comigo outras vezes. Obgda, # 2 dado_conc<-read.csv("PIB_estadual_pcorrentes_1939_2010.csv",sep=";",dec=",",header=TRUE, as.is=TRUE) str(dado_conc) 'data.frame': 151 obs. of 54 variables: $ X : chr "" "" "PIB POR UF (R$MIL)" "1939-1994: EST.SEC.XX (PART.NO PIB TOTAL)" ... $ X.1 : int NA NA NA NA NA NA 12 27 13 16 ... $ X.2 : chr "PIB BRASIL (R$MIL)" "" "" "" ... $ X1939: num 1.69e-08 NA NA NA NA ... $ X1947: num 6.49e-08 NA NA NA NA ... $ X1948: num 7.54e-08 NA NA NA NA ... $ X1949: num 8.79e-08 NA NA NA NA ... $ X1950: num 1.02e-07 NA NA NA NA ... $ X1951: num 1.27e-07 NA NA NA NA ... $ X1952: num 1.49e-07 NA NA NA NA ... $ X1953: num 1.78e-07 NA NA NA NA ... $ X1954: num 2.44e-07 NA NA NA NA ... $ X1955: num 2.96e-07 NA NA NA NA ... $ X1956: num 3.74e-07 NA NA NA NA ... $ X1957: num 4.54e-07 NA NA NA NA ... $ X1958: num 5.65e-07 NA NA NA NA ... $ X1959: num 8.43e-07 NA NA NA NA ... $ X1960: num 1.16e-06 NA NA NA NA ... $ X1961: num 1.69e-06 NA NA NA NA ... $ X1962: num 2.71e-06 NA NA NA NA ... $ X1963: num 4.86e-06 NA NA NA NA ... $ X1964: num 9.53e-06 NA NA NA NA ... $ X1965: num 1.55e-05 NA NA NA NA ... $ X1966: num 2.28e-05 NA NA NA NA ... $ X1970: num 7.07e-05 NA NA NA NA ... $ X1975: num 0.000382 NA NA NA NA ... $ X1980: num 0.00455 NA NA NA NA ... $ X1985: num 0.476 NA NA NA NA ... $ X1986: num 1.27 NA NA NA NA ... $ X1987: num 4.04 NA NA NA NA ... $ X1988: num 29.4 NA NA NA NA ... $ X1989: num 426 NA NA NA NA ... $ X1990: num 11549 NA NA NA NA ... $ X1991: chr "60285,99927" "" "" "" ... $ X1992: chr "640958,7676" "" "" "" ... $ X1993: num 14097114 NA NA NA NA ... $ X1994: chr "349204679" "" "" "" ... $ X1995: chr "705640892,1" "" "" "" ... $ X1996: num 8.44e+08 NA NA NA NA ... $ X1997: num 9.39e+08 NA NA NA NA ... $ X1998: num 9.79e+08 NA NA NA NA ... $ X1999: num 1.06e+09 NA NA NA NA ... $ X2000: num 1.18e+09 NA NA NA NA ... $ X2001: num 1.3e+09 NA NA NA NA ... $ X2002: num 1.48e+09 NA NA NA NA ... $ X2003: num 1.7e+09 NA NA NA NA ... $ X2004: num 1.94e+09 NA NA NA NA ... $ X2005: num 2.15e+09 NA NA NA NA ... $ X2006: num 2.37e+09 NA NA NA NA ... $ X2007: num 2.66e+09 NA NA NA NA ... $ X2008: num 3.03e+09 NA NA NA NA ... $ X2009: num 3.24e+09 NA NA NA NA ... $ X2010: num 3.77e+09 NA NA NA NA ... $ X2011: num 4.14e+09 NA NA NA NA ... # 3 anos<-dado_conc[1,4:54] str(anos) pib_bra<-dado_conc[2,4:54] str(pib_bra) pib_anos_bra<-rbind(anos,pib_bra) fix(pib_anos_bra) -- Vanessa Felix - Estatística (32) 8854-4881 / 9142-8428

6 10

obtencao de coef. de determinacao
by Leonardo Monteiro 08 Set '14

08 Set '14

Pessoal, mais outra questao sendo: a<- 1:10 b<- c(3,6,5,6,8,10,8,11,10,12) quando eu faco summary(lm(b~a)) , me resulta em um Rsquared = 0,8743 Entretanto, quando eu forço a passagem pela origem do grafico: summary(lm(b ~a-1)), o valor de Rsquared é maior (~0.95) O valor do Rsquared quando eu forco a passagem pela origem nao deveria ser menor?? obrigado, Leonardo.

3 2

Salvando figura em 300 dpi
by Leonardo Monteiro 08 Set '14

08 Set '14

Prezados, mais uma vez gostaria da ajuda de voces.. Bem, eu preciso salvar graficos em 300 dpi, entretanto nao estou conseguindo (infelizmente estou conseguindo salva-los em apenas 96 dpi) Alguem poderia me ajudar? mais uma vez, obrigado, Leonardo

3 2

Excluindo valores NaN
by Leonardo Monteiro 08 Set '14

08 Set '14

Pessoal, bom dia Preciso excluir valores NaN e na de um cjto de dados, p fazer o grafico de dispersao. alguem poderia me ajudar? obrigado, Leonardo

3 2

Dúvida com tabela do Censo 2010
by Carlos Andrade 06 Set '14

06 Set '14

Prezados da lista, Quero fazer uma análise da distribuição de renda por bairro no município de Campina Grande, Paraíba. Baixei do site do IBGE a planilha Resultados do Universo do Censo Demográfico 2010 e peguei a Tabela 4.12.7.1 - Pessoas de 10 anos ou mais de idade, por classes de rendimento nominal mensal, segundo as mesorregiões, microrregiões, os municípios e os bairros - Paraíba-2010. As classes de rendimento nominal mensal, em salário mínimo, são: Até 1/2, Mais de 1/2 a 1, Mais de 1 a 2, Mais de 2 a 5, Mais de 5 a 10, Mais de 10 a 20, Mais de 20 e Sem rendimento. A dúvida: para calcular as estatísticas descritivas média, mediana, quartis e percentis a classe Sem rendimento é incluída? Agradeço antecipadamente qualquer ajuda. -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

4 5

ler aquivo .spt e transformar para .csv no r
by Alisson Lucrécio 06 Set '14

06 Set '14

Caros colegas, Boa noite. Tenho um arquivo .spt (.SPT) e que quero transformar para .csv (.CSV), algum colega saberia como fazer isso no R? Obrigado. -- Alisson Lucrecio da Costa

5 13

Sugestão de mateiral GAM
by Sérgio Henrique almeida da silva ju 05 Set '14

05 Set '14

Olá Alguém sabe de um bom material de modelos aditivos generalizados no R com variáveis resposta binária? Abraços -- Sérgio Henrique Almeida da Silva Junior Instituto Benjamin Constant - IBC 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) 968463637

3 4

Instalar RStudio por meio de repositórios *.deb
by walmes . 05 Set '14

05 Set '14

Saudações, Aqui na UFPR todo sistema coletivo de computadores é mantido pelo c3sl (centro de computação científica e softaware livre, é quem disponibiliza o nosso serviço de lista de discussão). O sistema, para funcionar com milhares de usuários em rede, com segurança, etc, obviamente é Linux. O c3sl, por medidas de segurança e manutenção, só instala aplicativos disponíveis nos repositórios debian. Sendo assim, R, texlive, Emacs, libreOffice, firefox, estão todos disponíveis. Quando existe atualização, o sistema automaticamente incorpora. Mas o RStudio não está na lista. Ele é instalado, mesmo o server, por meio de um *.deb e para ser atualizado deve-se periodicamente baixar e instalar o novo deb, o que vai contra a política de manutenação do c3sl. Dito isso, alguém saberia uma forma de incorporar as atualizações do RStudio automaticamente, olhando para um repositório, tal como acontece com o R e Emacs? Grato. 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 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

Problema com representação de parcelas usando expand.grid()
by ASANTOS 04 Set '14

04 Set '14

Boa noite Pessoal, Estou tentando representar as parcelas de um experimento em campo e estou tendo dois problemas sendo: 1) _Representação do centro das parcelas_: ## Coordenadas do grid no espaçamento das plantas de 2x4m ## distâncias verticais e horizontais entre plantas v <- 2; h <- 4 ## número de plantas vertical e horizontal npv <- 4; nph <- 5 ## distância entre vãos vertical e horizontal (distância entre parcelas) V <- 7; H <- 7 ## número de parcelas em linhas e colunas ncol <- 4; nrow <- 12 grid0 <- expand.grid(L=seq(2,length.out=nrow,by=npv*v+V), C=seq(2,length.out=ncol,by=nph*h+H), x=seq(2,length.out=npv,by=v), y=seq(2,length.out=nph,by=h)) grid1 <- with(grid0, data.frame(x=L+x, y=C+y)) plot(y~x, grid1, asp=1) ### Representação das plantas Até aí tudo bem, porém quando tento representar por pontos o centro das parcelas, começa o problema, sendo: ### O centro de cada parcela vai se localizar a cada 13 m na horizontal e a cada 23 na vertical grid2 <- expand.grid(seq(7,length.out=13*4,by=13), seq(12,length.out=23*12,by=23)) points(grid2[,1],grid2[,2],col="red") ## Vejo que os pontos caem errado e não sei por que? 2) _Criar uma coluna para identificar a parcela a qual pertence cada planta_: pos<-1:48 text(grid2[,1],grid2[,2]+0.5,labels=pos, cex=1,col=2) e depois que o centro cair no local certo, gostaria de fazer grid1$pos para identificar a qual centro da parcela cada ponto pertence. Espero que alguém possa dar um help, Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

1 0

Popularização da Estatística
by jose gomes de souza filho Gomes 04 Set '14

04 Set '14

Achei muito legal o blog do professor Dani Gamerman, professor Titular do Departamento de Métodos Estatísticos - UFRJ. O blog se destina à divulgação e popularização da Estatística para a sociedade em geral. http://statpop.blogspot.com.br/p/sobre-o-blog.html Gomes.

5 11

Contraste
by Daniel Batista Lemes 04 Set '14

04 Set '14

Boa tarde, Pessoal, sou novo n alista, mas já tenho minhas dúvidas. Sou usuário SAS e estou migrando para o R e consegui fazer contraste das médias, porém o parâmetro retornado é o estimate e gostaria que tivesse os mesmos dados da saída do contrats do SAS que traz a soma de quadrados, alguém sabe como saber ou algum package que faça isso? Att Daniel Lemes -- @lemes_daniel

6 9

identify(...)
by ana paula coelho madeira 04 Set '14

04 Set '14

Prezados colegas, boa tarde. Trabalhando com os gráficos de pontos alacanca, de pontos influentes e dois de resíduos (glm), tento identificar qual (ou quais) pontos são influentes usando a função identify(...). No entanto o R fica rodando, rodando... até que eu paro a computação. Alguém poderia me dizer o que pode estar ocorrendo? Desde já agradeço. Atenciosamente, Ana Paula

4 4

Limites API Google maps
by Leonard de Assis 04 Set '14

04 Set '14

Alguém sabe qual o limite diário desta API? Construí um código hoje à tarde e rodava. Fui gerar o pdf hoje a noite e deu um erro louco gerando um gráfico. Quando rodei o código novamente, onde era pra ter as coordenadas estava tudo cheio de NA

5 5

Índice Moran bivariado
by Jacqueline Cantú 04 Set '14

04 Set '14

Olá pessoal,Alguém sabe como calcular o índice de Moran bivariado no R? Att, Jacqueline

1 0

Instalar versão mais recente do R no Debian Squeeze
by Nei 03 Set '14

03 Set '14

Oi Pessoal, Estou planejando instalar o R no Debian Squeeze. Pelo que vi a versão que está disponível no repositório é a 2.11.1-6 Pelo site http://cran.r-project.org/ a última versão é 3.1.1. Como faço para instalar a versão mais recente e não ter problemas de dependência com as bibliotecas instaladas? Obrigado, Nei

3 2

Codificação
by Giselle Davi 03 Set '14

03 Set '14

Prezados Colaboradores, Tenho um arquivo com uma matriz formada por letras ( SNP: bases nitrogenadas: A, C, T, G) e preciso codificá-la para números, onde as letras A e T recebam 1 e as letras C e G recebam 0, formando assim uma matriz de genótipos Gn x m = {0;1;2} onde n é o número de indivíduos e m o número de SNPs. Anexado a este e-mail encontra-se parte dessa matriz com dimensão 6 x 30 ( Facilitando, portanto, o envio) que necessita dessa codificação. Matematicamente, a sequência de SNP é atribuída ( A qual eu não precisei construir, pois a sequência foi originada por genotipagem ) como: Dado um valor para PA (probabilidade de ocorrer a base A) e para PC (probabilidade de C), de modo que PA= PT, PC = PG e PA+PC = 0; 5, então as bases são geradas do seguinte modo: Se y existe (E virado para direita) [0;PA), a base é A; Se y existe [PA;0:5), a base é C; Se y existe [0:5; (0:5+PC)), a base é G; Se y existe [(0:5+PC);1], a base é T. Cada indivíduo possui um par de cromossomos homólogos, que se mantêm juntos por ligações entre as bases A e T e ligações entre as bases C e G. A codificação das sequências ( SNPs) deve ser feita, portanto, do seguinte modo: Contando o alelo de interesse (no caso, A) e seu correspondente no pareamento (T) em cada fita numa posição do genótipo, e atribuir 1 quando um dos dois (A ou T) ocorre e 0 se não ocorrem. Esta codificação resulta em: 0 - CC, CG, GC, GG)P(0) = 4*PC elevado a 2 1 - AC, AG, TC, TG, CA, GA, CT, GT)P(1) = 8*PA*PC 2 - AA, AT, TA, TT )P(2) = 4*PA elevado a 2 Sendo assim, preciso utilizar estas informações matemáticas no R de modo a proceder a codificação. Alguém poderia me ajudar ? Desde já agradeço. Att, Giselle

4 3

percentage prediction error
by Humberto Hazin 02 Set '14

02 Set '14

Prezados, Estou com uma dúvida que pode parecer besteira, porém gostaria da ajuda de vocês! A dúvida é a seguinte como calcular o percentual de "Bias" (percentage prediction error) por linha de um modelo ZINB (package- pscl)? O meu objetivo é fazer um boxplot com esse percentual por ano Eu iniciei da seguinte forma: dois bancos de dados um de try.csv (10% retirados do banco original) e outro original.csv try<-read.csv(try.csv, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) <https://www.dropbox.com/s/v7vruztus11nnbi/try.csv> https://www.dropbox.com/s/v7vruztus11nnbi/try.csv original<-read.csv(original.csv, header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) <https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv> https://www.dropbox.com/s/7htqc2bqi45rkku/original.csv # Criar Fatores para o banco original original$area<-rep(1,nrow(original)) original$area[original$lat<=-15&original$lat>=-45]<-2 original$A=as.factor(original$area) summary(original) original$S<-factor(original$id.group.fleet) original$Q<-factor(original$Quarter) original$Y<-factor(original$Year) original$logEf<-log(original$effort1) summary(original) # Criar Fatores para o banco try try$area<-rep(1,nrow(try)) try$area[try$lat<=-15&try$lat>=-45]<-2 try$A=as.factor(try$area) summary(try) try$S<-factor(try$id.group.fleet) try$Q<-factor(try$Quarter) try$Y<-factor(try$Year) try$logEf<-try$effort1*0 summary(try) zeroinfl.fit_original<-zeroinfl(BSH~Y+S+Q+A+offset(logEf)|Y+S+Q+A,dist="negb in",data=original) pred<-predict(zeroinfl.fit_original,newdata=try,type="response") Daqui em diante não tenho mais ideia de como gerar esse percentage prediction error Alguem poderia me ajudar! Desde já Agradeço Humberto

3 2

instalar pacote AED (Zuur et al.)
by Ana Maria Nievas 02 Set '14

02 Set '14

Bom dia! Alguém sabe como instalar o pacote AED do Zuur? A dificuldade está no fato de ele não aparecer no CRAN. Obrigada!! -- Ana Maria Nievas Bióloga e Mestre em Ciências, com ênfase em Ecologia Aplicada (16) 98843-7744

3 3

Re: [R-br] Mensagem de erro: "testes-F ANOVA sobre um ajuste essencialmente perfeito, são incertos"
by andrebvs＠bol.com.br 02 Set '14

02 Set '14

1 0

Leitura de dados da PNAD 2008 com o pacote dicionariosIBGE
by Leonardo Ferreira Fontenelle 02 Set '14

02 Set '14

Bom dia, pessoal. Estou começando a namorar os microdados da Pesquisa Nacional por Amostra de Domicílios de 2008, mas não estou conseguindo utilizar o pacote dicionariosIBGE para ler os microdados. Na saída abaixo, substituí o prompt ">" por "." para evitar que clientes de e-mail interpretem como citação de e-mails prévios. . library("dicionariosIBGE") . data(dicPNAD2008) . dicdom2008 inicio cod tamanho desc 1 1 V0101 4 ANO DA PESQUISA 2 5 UF 2 UNIDADE DA FEDERACAO 3 5 V0102 8 NUMERO DE CONTROLE 4 13 V0103 3 NUMERO DE SERIE 5 16 V0104 2 TIPO DE ENTREVISTA [...] . amostra <- le.pesquisa(dicdom2008, "~/Dropbox/Cobertura/2008/dados/DOM2008.txt", codigos = c("V0101", "UF", "V0102", "V0103", "V4105", "V4106", "V4606", "V4604", "V4605", "V4606", "V4607", "V4608", "V4609", "V4610", "V4611", "SUBAMO", "V0223"), nlines = 150591) |=====================================================================| 100% . summary(amostra$V0101) Mode NA's logical 1015101 . summary(amostra$UF) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 0.00 3.00 12.00 19.69 31.00 99.00 264343 Como os microdados da PNAD 2008 estão reponderados utilizando informações posteriores à publicação original do arquivo, considerei a possibilidade de que a organização dos dados tivesse mudado. Mas, conferindo o dicionário de dados atualmente disponibilizado, vi que os números de início, código e tamanho continuam os mesmos, pelo menos para as primeiras variáveis (são setenta variáveis, só neste banco de dados). O que será que estou fazendo de errado? Grato, Leonardo Ferreira Fontenelle http://lattes.cnpq.br/9234772336296638

4 9

Mensagem de erro: "testes-F ANOVA sobre um ajuste essencialmente perfeito, são incertos"
by andrebvs＠bol.com.br 02 Set '14

02 Set '14

3 2

Fit não linear com função que aprensenta singularidades.
by Michelle Baú Graczyk 02 Set '14

02 Set '14

Desculpe Fernando, seguem os dados: > dput(file) structure(list(x = 0:372, y = c(1.552918816, 1.177870257, 0.728582521, 0.754977859, 1.014857186, 0.952781237, 0.92487764, 0.930448803, 0.825743421, 1.061320976, 0.778891815, 0.987829557, 0.839076467, 1.035174095, 1.026961486, 1.180269616, 0.897432732, 0.750135342, 0.713731659, 0.986891474, 0.814066753, 0.838472048, 0.819105964, 0.913483671, 1.134245404, 0.790984967, 0.785299927, 0.795018406, 0.839324923, 0.735915178, 0.593269981, 0.923184631, 1.018998131, 0.803540789, 0.80479661, 0.72999522, 0.699068733, 0.618802888, 0.815272405, 0.803152178, 0.701841573, 0.734495312, 0.73549551, 0.903936404, 0.99491774, 0.990528339, 0.752979853, 0.708288706, 0.822200524, 0.787724769, 0.852584706, 0.741692593, 0.833094394, 0.835140976, 0.843551076, 0.789762592, 0.71986185, 0.641050114, 0.800506751, 0.884314696, 0.987518987, 0.768463801, 0.866236287, 1.065961223, 0.806744447, 0.766739735, 0.751533707, 0.696793602, 0.936172067, 0.998768819, 0.790549371, 1.112424217, 0.727549674, 0.91716513, 0.785049082, 0.835587721, 0.673062803, 0.849652599, 0.849677285, 0.792995714, 0.817868458, 0.790538222, 0.718977121, 0.969363333, 0.821173251, 0.682762966, 0.547801611, 0.864380029, 0.637068436, 0.716060144, 0.925191396, 0.82309242, 1.217327091, 0.950335691, 0.825029108, 0.858161444, 0.620884509, 0.743739175, 0.782014247, 0.815065357, 0.746548647, 0.778370215, 0.747050338, 0.85151045, 0.853198681, 0.741675073, 0.980756505, 0.78396288, 0.87737065, 0.76027349, 0.814334321, 0.955472852, 0.797915475, 0.64547774, 0.703556084, 1.160420157, 0.697551714, 0.680614453, 0.87038201, 0.885639002, 0.931712588, 0.634435751, 0.771353703, 0.723837157, 0.864745547, 0.683628583, 0.677707032, 0.625280281, 0.780408834, 0.740849273, 0.816243934, 0.743649986, 0.767603759, 0.788182662, 0.638866562, 0.749819197, 0.61128548, 0.699447788, 0.742390183, 0.654970059, 0.621842501, 0.841482993, 0.71000242, 0.719206466, 0.651004308, 0.657988171, 0.659939193, 0.590479621, 0.916040704, 0.735814044, 0.780472541, 0.690410176, 0.852647617, 0.663802217, 0.774593992, 0.679669203, 0.841208257, 0.644216344, 0.790331175, 0.655288593, 0.598600651, 0.850198089, 0.639451072, 0.790941168, 0.56980038, 0.751766237, 0.800986145, 0.661843231, 0.699207295, 0.729910808, 0.713261821, 0.658346522, 0.698184004, 0.725423458, 0.754622693, 0.681435475, 0.471715731, 0.705996852, 0.570707406, 0.660579447, 0.738280295, 0.68915117, 0.568583579, 0.61153553, 1.02975742, 0.661527086, 0.567946511, 0.543164549, 0.702431658, 0.509399922, 0.511034002, 0.683268639, 0.743713692, 0.583393828, 0.599969552, 0.722769271, 0.562765552, 0.596352596, 0.640454455, 0.610732823, 0.734572557, 0.644839078, 0.714663372, 0.679286961, 0.669174296, 0.639958338, 0.732628702, 0.521852221, 0.591436817, 0.622606983, 0.592701398, 0.584728486, 0.808503552, 0.557689709, 0.683618231, 0.671007461, 0.612336643, 0.867552629, 0.71857975, 0.581004025, 0.57724771, 0.65777316, 0.53544089, 0.81379122, 0.607679673, 0.561732705, 0.583045033, 0.642250988, 0.612925931, 0.619729026, 0.655734541, 0.674082112, 1.199682684, 0.565012811, 0.59800977, 0.839234937, 0.806546956, 0.58329986, 0.673788265, 0.627473389, 0.69975597, 0.795051056, 0.606857855, 0.767389545, 0.691856322, 0.637890254, 0.615875558, 0.527046717, 0.69097876, 0.543992738, 0.721710941, 0.541890412, 0.623967124, 0.765809615, 0.697307239, 0.648002123, 0.695423905, 0.58107251, 0.663865924, 0.676050654, 0.576674348, 0.75346084, 0.684342896, 0.692429683, 0.749339007, 0.681089865, 0.694518471, 0.61196794, 0.704916225, 0.752186703, 0.973376864, 0.884314696, 0.892780539, 0.644331016, 0.810050036, 0.659620659, 0.875652955, 0.739896856, 1.050720157, 0.701284935, 0.939293702, 1.023670231, 0.841587313, 0.925278197, 0.780156396, 0.887787516, 0.823984315, 0.985433384, 0.894676614, 0.83659827, 0.96598687, 0.787544001, 0.764039361, 0.898974437, 0.861125405, 0.933638924, 0.86870652, 0.861905018, 0.992322483, 1.381073587, 0.899413218, 1.055144597, 1.048128881, 0.947446586, 0.816658028, 0.85940532, 0.972548675, 0.875777979, 0.974942459, 0.790028568, 0.71600679, 0.792859541, 0.889052096, 0.966600048, 0.89341283, 0.958572986, 1.020383755, 0.76941144, 1.110825175, 0.843357567, 1.120730793, 0.889411244, 0.999906782, 1.02455496, 1.051479065, 1.095087992, 0.923559705, 1.128647961, 1.313131831, 0.915153587, 1.058557692, 0.853607998, 0.949029701, 0.969001, 1.091106314, 1.033087696, 0.957625347, 0.90042775, 1.315533579, 0.86301272, 0.970012346, 0.996177543, 0.872745533, 1.169665612, 1.130164184, 1.12352832, 1.059315803, 1.008628249, 1.296573626, 1.045684131, 1.193302444, 1.181863084, 1.096267365, 1.240134937, 1.16928576, 1.233877333, 1.092117661, 1.15841578, 1.500965885, 1.229769834, 1.438586533, 1.346249834, 1.459758706, 1.542678738, 1.583001188, 1.512468952, 1.330386034, 1.215865019, 1.642284388, 1.588310357, 1.602909577, 1.479793712, 1.527826283 )), .Names = c("x", "y"), class = "data.frame", row.names = c(NA, -373L)) Boa Tarde, Eu estou fazendo um fit nâo linear com um polinômio. Porém, a fiunção que eu estou usando : y= - {[a*(x^(2*n)) + b]/[(x-186)(x+186)]} apresenta duas assíntotas em +- 186 e por isso dá um erro no fit: Erro em numericDeriv(form[[3l]], names(ind), env): obtido valor faltante ou infinito quando avaliando o modelo. Pelo que eu entendi, o motivo é que a derivada também possui assíntotas nestes mesmo postos e daí na função vai para o infinito. Eu fiz D( (-(a*(x^(2*n)) + b)/((x-186)(x+186))), "x") e a mensagem de erro foi: Erro em x^(2*n)argumento não numérico para operador binário. Aí eu pensei que era por causa do n e deixei então só o x^2 e o erro foi o mesmo. Eu procurei na internet e não achei nada que explicasse o que está querendo dizer. O que seria um operador binário? Voltando às assintotas, eu achei uma função chamada ‘NLSstAsymptotic(xy)’ que é para assintotas verticais em gráficos mas eu tentei emplementá-la e nâo deu certo. Tem outra função que eu possa usar para que o fit seja feito mesmo que a função usada para o fit apresente singularidades? Abaixo está o código: file <- read.table(file="/Users/bau/Desktop/teste/DISparam.txt", col.names=c('x','y')) plot(x=file$x, y=file$y, xlab='Valores de x', ylab='Valores de y', lwd=0.5,) dados <- data.frame(x=file$x,y=file$y) #chute inicial guess <- list(a=1, b=3, n=3) # funcao usada para o fit # func <- function(x,n, a, b){ y <- (-(a*(x^(2*n)) +b)/((x-186)*(x+186))) return(y) } #executa o fit e armazena resultado na variavel 'fit' fit <- nls(y ~ func(x,n,a,b), data=dados, start=guess, lower=-200, upper=200) fit summary(fit) #gera valores de y a partir do resultado do fit, para plotar #o resultado no mesmo grafico que os dados yfit <- predict(fit) #plota a linha do fit lines(dados$x, yfit, col='blue',lwd=2) #carrega os parametros param <- coef(fit) Obrigada, Michelle

5 7

Uma pequena dúvida sobre métodos
by Jose Claudio Faria 02 Set '14

02 Set '14

Olá pessoal, Gostaria de fazer uma consulta ao pessoal mais ligado ao desenvolvimento de soluções para o R. Estamos trabalhando para fazer o pacote fdth (http://cran.r-project.org/web/packages/fdth/index.html) mais genérico, flexível e por conseguinte mais útil aos usuários. Entre as novidades: - os métodos tabulares e gráficos foram estendidos para as variáveis categóricas com as novas opções (barchart, poligon e dotchart para qualquer coluna da fdt assim como o gráfico de pareto) para as categóricas; - permite a criação de tabelas de distribuição de frequências (fdt) para variáveis contínuas e categóricas com um mínimo de argumentos baseado em uma fdt prévia existente na literatura (ou seja, você não conhece - de fato - o vetor da variável, apenas as frequências); - cálculo de medidas de tendência central (média e mediana) e posição (quantil genérico) a partir de uma fdt. Para gerar uma tabela de uma variável contínua (por exemplo) à partir do vetor de frequências é bem simples: > tb <- make.fdt(x=vetor_freq, start, end) Uma vez criado 'tb' irá pertencer a classe 'fdt'. Foram implementando métodos para as funções genéricas (mean, median, quantile) para a classe 'fdt'. por exemplo: > mean(tb) > median(tb) > quantile(tb, i=10, probs=seq(0, 1, .1) para calcular o 10 percentil. Minha questão é a seguinte, as funções 'var' e 'sd' não são genéricas. Assim não posso implementar o métodos: var.fdt e sd.fdt para chamar assim: > var(tb) > sd(tb) Poderia se criasse nova funções 'var' e 'sd' genéricas no pacote. Contudo, fico um pouco apreensivo com esta abordagem. Qual seria uma opção razoável para este caso? A mais simples seria escrever as funções 'var.fdt' e 'sd.fdt' e o usuário teria que usar assim: > var.fdt(x=tb) > sd.fdt(x=tb) Mas foge do padrão do restante do pacote que realmente é orientado a objetos S3. Nunca tinha me deparado com uma caso desses. Sugestões são muito bem vindas... 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 ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

4 4

Pacote formulário para criação de banco de dados
by Fátima Lima Paula 01 Set '14

01 Set '14

Prezados, existe algum pacote que nos permita criar um formulário para ser preenchido para montar um banco de dados? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

3 4

Re: [R-br] IMPORTAR ARQUIVO PNG TRANSPARENTE
by andrebvs＠bol.com.br 01 Set '14

01 Set '14

1 0

IMPORTAR ARQUIVO PNG TRANSPARENTE
by andrebvs＠bol.com.br 01 Set '14

01 Set '14

3 2

Arquivos com extensão json.
by Edson Lira 31 Ago '14

31 Ago '14

Caros amigos, uso o aplicativo ODKAggregate para gerenciar dados coletados via smartfones, ele me permite exportar os dados coletados no formato csv, porém, sua conversão tem problemas com a acentuação, já tentei usar a opcão fileEncod e não consegui corrigir o problema. Quando exporto no formato json, o problema da acentuação se resolve, consigo abrir no notepad, e ver uma estrutura tipo de programação. No link abaixo, tem o conjunto de dados exportado pelo ODKAggregate. http://www.datafilehost.com/d/d60c05e4 Alguém poderia dar uma olhada e verificar se há uma saída analisar esses dados com o uso do R? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

5 8

Acerca do uso de Excel em estatística
by Cesar Rabak 31 Ago '14

31 Ago '14

Na "trédi" CAPTURAR CÓDIGO DE COR a seguinte questão foi colocada por André Rutz: > Walmes, fiquei curioso com seu comentário sobre a qualidade do ms > Excel. Pode falar mais sobre o porque de vc considerar ele de > qualidade nula? Sua opinião sobre isto é muito valiosa. Embora a questão tenha sido endereçada para o Walmes e eu não tenho procuração para falar em nome dele ;-), coloco aqui à guisa de subsídio os seguintes links: http://trust.guidestar.org/2014/07/09/friends-dont-let-friends-use-excel-to… http://www.dogoodconsulting.org/resources/documents/friendsdonotletfriendsu… Este link, embora de um documento de 2001, ainda está acurado até a mais recente versão: http://www.stat.uiowa.edu/~jcryer/JSMTalk2001.pdf Neste link mais uma plétora de links sobre os problemas do Excel em estatística: http://biostat.mc.vanderbilt.edu/wiki/Main/ExcelProblems http://www.pucrs.br/famat/viali/tic_literatura/artigos/planilhas/pottel.pdf Note-se que o assunto é rebarbativo, com defensores do seu uso (uma busca rápida na Internet sobre livros e cursos baseados nessa ferramenta) havendo inclusive um verdadeiro ecossistema em volta do Excel com os "add ons", "plugins" e demais acréscimos. Meus 0,0199999..... -- Cesar Rabak

2 1

Utilização de softwares desconhecidos
by Fernando Souza 31 Ago '14

31 Ago '14

Caros colegas, Estou corrigindo um artigo onde a autora utilizou um software que desconheço: BARBOSA, J.C.; MALDONADO JR., W. *AgroEstat:* sistema para análises estatísticas de ensaios agronômicos. Versão 1.0. Jaboticabal: Departamento de Ciências Exatas, 2010. Fiz uma busca na internet para obter informações sobre o software e obtive alguns trabalhos publicados que utilizaram esse software mas nenhum link ou referência que me direcione para um local onde posso baixar ou comprá-lo ou algum manual do software. Minha dúvida é: Como proceder com softwares desconhecidos. Há quem diga que o resultado é o mesmo que softwares consagrados e que não há diferenças, pois Anova é tudo igual. Ou que quem realiza a análise acaba validando o software, pois cabe ao analista assegurar os resultados. Há revistas que questionam o uso de softwares não conhecidos. Eu creio que não e que softwares mais conhecidos são validados pelo seu uso pela comunidade e caso não conhecemos seu funcionamento, podemos recorrer a um colega ou buscar alguma comunidade ou informações na net. Qual é a opinião de vocês? Atenciosamente --- Este email está limpo de vírus e malwares porque a proteção do avast! Antivírus está ativa. http://www.avast.com

5 4

ler shapefile
by Aline Lipsky 29 Ago '14

29 Ago '14

Caros colegas, Tentei ler um shapefile com os seguintes comandos: readShapeSpatial(file.choose())readShapePoly(file.choose())readShapeLines(file.choose()) Aparece o seguinte erro: Error in getinfo.shape(filen) : Error opening SHP file. Alguém saberia como resolver? Agradeço desde a ajuda, Att., Aline

5 5

Read .spt e export como .csv
by Alisson Lucrécio 29 Ago '14

29 Ago '14

Caros Colegas, Bom dia. Tenho um arquivo .spt (.SPT) e que quero transformar para .csv (.CSV), algum colega saberia como fazer isso no R? Obrigado. -- Alisson Lucrecio da Costa

2 1

DENSIDADES EMPÍRICAS
by andrebvs＠bol.com.br 29 Ago '14

29 Ago '14

2 1

Re: [R-br] CAPTURAR CÓDIGO DE COR
by andrebvs＠bol.com.br 29 Ago '14

29 Ago '14

1 0

CAPTURAR CÓDIGO DE COR
by andrebvs＠bol.com.br 28 Ago '14

28 Ago '14

5 8

modelo BCRM
by Jacqueline Cantú 27 Ago '14

27 Ago '14

Olá pessoal,Alguém sabe dizer se já existe algum pacote ou função implementada no R sobre o modelo espacial gaussiano bivariado de co-regionalização (BCRM)? Att, Jacqueline

2 1

Fit não linear com função que aprensenta singularidades.
by Michelle Baú Graczyk 27 Ago '14

27 Ago '14

Boa Tarde, Eu estou fazendo um fit nâo linear com um polinômio. Porém, a fiunção que eu estou usando : y= - {[a*(x^(2*n)) + b]/[(x-186)(x+186)]} apresenta duas assíntotas em +- 186 e por isso dá um erro no fit: Erro em numericDeriv(form[[3l]], names(ind), env): obtido valor faltante ou infinito quando avaliando o modelo. Pelo que eu entendi, o motivo é que a derivada também possui assíntotas nestes mesmo postos e daí na função vai para o infinito. Eu fiz D( (-(a*(x^(2*n)) + b)/((x-186)(x+186))), "x") e a mensagem de erro foi: Erro em x^(2*n)argumento não numérico para operador binário. Aí eu pensei que era por causa do n e deixei então só o x^2 e o erro foi o mesmo. Eu procurei na internet e não achei nada que explicasse o que está querendo dizer. O que seria um operador binário? Voltando às assintotas, eu achei uma função chamada ‘NLSstAsymptotic(xy)’ que é para assintotas verticais em gráficos mas eu tentei emplementá-la e nâo deu certo. Tem outra função que eu possa usar para que o fit seja feito mesmo que a função usada para o fit apresente singularidades? Abaixo está o código: file <- read.table(file="/Users/bau/Desktop/teste/DISparam.txt", col.names=c('x','y')) plot(x=file$x, y=file$y, xlab='Valores de x', ylab='Valores de y', lwd=0.5,) dados <- data.frame(x=file$x,y=file$y) #chute inicial guess <- list(a=1, b=3, n=3) # funcao usada para o fit # func <- function(x,n, a, b){ y <- (-(a*(x^(2*n)) +b)/((x-186)*(x+186))) return(y) } #executa o fit e armazena resultado na variavel 'fit' fit <- nls(y ~ func(x,n,a,b), data=dados, start=guess, lower=-200, upper=200) fit summary(fit) #gera valores de y a partir do resultado do fit, para plotar #o resultado no mesmo grafico que os dados yfit <- predict(fit) #plota a linha do fit lines(dados$x, yfit, col='blue',lwd=2) #carrega os parametros param <- coef(fit) Obrigada, Michelle

2 2

Re: [R-br] CAPTURAR CÓDIGO DE COR
by andrebvs＠bol.com.br 27 Ago '14

27 Ago '14

1 0

Re: [R-br] CAPTURAR CÓDIGO DE COR
by andrebvs＠bol.com.br 27 Ago '14

27 Ago '14

1 0

Dados em uma tabela
by jose gomes de souza filho Gomes 27 Ago '14

27 Ago '14

Boa Noite, Tenho o seguinte banco de dados em txt( em anexo): JOSÉ A................................. 2 B.....................................0 C...................................0 D.................................... 6 E.......................0 F............................. 1 G................................ 13 MARIA A................................. 3 B.....................................0 C...................................0 D.................................... 1 E.......................0 F............................. 1 G................................ 6 PEDRO A.................................0 B..................................... C...................................0 D...................................0 E.......................0 F.............................0 G................................ 3 Gostaria que ficasse assim: Nome A B C D E G F JOSÉ 2 0 0 6 0 1 13 MARIA 3 0 0 1 0 1 6 PEDRO 0 0 0 0 0 0 3 Como fazer isso? Att, Gomes

4 4

Sugestão de site para download de livros grátis.
by andrebvs＠bol.com.br 27 Ago '14

27 Ago '14

9 13

Números grandes estourando!
by Jodavid Ferreira 25 Ago '14

25 Ago '14

Alguém sabe algum tipo númerico de armazenamento maior que o "as.double" no R. Estou fazendo uma simulação, porém os valores são tão altos que estão estourando -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Bacharelando em Estatística - Universidade Federal da Paraíba - www.de.ufpb.br Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://blog.jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

1 0

R lendo arquivos .spt
by Alisson Lucrécio 24 Ago '14

24 Ago '14

Caros colegas da r-br, Boa noite. Algum de vocês saberia se R consegue ler arquivos .spt? Obrigado. -- Alisson Lucrecio da Costa

1 0

kml - importação para análise de dados
by Leonard de Assis 24 Ago '14

24 Ago '14

1. Qual o melhor jeito de ler um arquivo kml contendo polígonos no R? 2. Como identificar em quais destes polígonos está contida uma determinada coordenada? É que eu tenho dados geoespaciais para vários polígonos parte armazenados em arquivos kml e parte eu tenho a referência (x,y) com os dados da ocorrência, preciso de uma maneira saudável para unir essas informações

3 5

Re: [R-br] Dúvida Simples - Histograma
by andrebvs＠bol.com.br 23 Ago '14

23 Ago '14

4 4