- R-br - DINF UFPR

OFF TOPIC - Backup do histórico do firefox e decodificação da data de acesso
by Walmes Zeviani 17 Mai '11

17 Mai '11

Saudações, Peço desculpas pelo tema fora de tópico. Tenho o backup do histórico do navegador, mozilla firefox, usando o complemento FEBE ( https://addons.mozilla.org/pt-BR/firefox/addon/febe/) O arquivo com o backup do histórico de visitação é um arquivo de extensão *.json. Os registros de páginas visitadas apresentam essa estrutura (essa é apenas uma linha das 15000 que tenho) {"aURI":"https://addons.mozilla.org/","aTitle":"addons.mozilla.org ","aLastVisited":1305653472464928} Imagino que o último termo, o numérico com 16 dígitos seja uma forma de armazenar a data de acesso dessa página. Alguém saberia como decodificar esse número em data, no caso a data é 2011-05-17 14:31. 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

2 2

eixo log contour plot
by eder＠leg.ufpr.br 17 Mai '11

17 Mai '11

Pessoal, to me batendo aqui, Como pode ser feito para deixar apenas o eixo y em escala logaritima (base 10), e que os valores das variáveis continuem aparecendo no grafico. ## Exemplo x <- seq(pi/4, 5 * pi, length.out = 100) y <- seq(pi/4, 5 * pi, length.out = 100) r <- as.vector(sqrt(outer(x^2, y^2, "+"))) grid <- expand.grid(x=x, y=y) grid$z <- cos(r^2) * exp(-r/(pi^3)) contourplot(z~x*y, grid,col.regions=colorRampPalette(c("green","yellow","red")),region =TRUE) att

3 3

Dúvida básica
by Rodrigo Pereira 17 Mai '11

17 Mai '11

Caros colegas, Estou com um problema básico, mas que não estou conseguindo resolver (falta de prática com o R...) Estou com o seguinte banco de dados: resposta fator sujeito 1 0.00 R 1 2 2.96 R 1 3 10.28 R 1 4 12.42 R 1 5 15.60 R 1 6 14.23 R 1 7 10.55 R 1 8 8.00 R 1 9 6.12 R 1 10 3.24 R 1 11 2.82 R 1 12 1.18 R 1 13 0.42 R 1 14 0.00 R 1 15 0.00 R 1 16 0.00 T 2 17 4.01 T 2 18 19.82 T 2 19 17.83 T 2 20 12.74 T 2 21 7.37 T 2 22 7.73 T 2 23 4.70 T 2 24 4.11 T 2 25 3.11 T 2 26 3.19 T 2 27 2.25 T 2 28 0.44 T 2 29 0.16 T 2 30 0.00 T 2 Estou tentando utilizar a função unstack() para transformar este banco de dados longo em largo da seguinte forma: R 0 T 0 2,96 4,01 10,28 19,82 12,42 17,83 15,6 12,74 14,23 7,37 10,55 7,73 8 4,7 6,12 4,11 3,24 3,11 2,82 3,19 1,18 2,25 0,42 0,44 0 0,16 0 0 Alguém sabe como consigo fazer essa mudança com ou sem a função unstack()? Obrigado e abraço! Rodrigo Pereira (31)9326 3139

5 16

?Erros Correlacionados?
by Marcelo Cardoso mello 17 Mai '11

17 Mai '11

Estou com dificuldades em entender a idéia de erros correlacionados. Erros correlacionados com o que? Correlação entre quais variáveis? Entre os erros “puros” e as medias dos tratamentos? Entre os erros “puros” e os valores preditos dos tratamentos? Ou seriam os erros padronizados? Alguém poderia me mostrar um exemplo onde isto ocorre? A idéia é de correlação linear ou de qualquer tipo de dependência, ou seja, uma parabola me dá correlação linear nula mas ocorre relação e existe dependência (não linear). Abaixo segue um exemplo. Onde posso observer a dependência? t<-as.factor(c(rep(1,3),rep(2,3),rep(3,3), rep(4,3))) b<-as.factor(c(seq(1:3), seq(1:3), seq(1:3), seq(1:3))) r<-c(1,2,3,10,18,20,5,25,26,6,30,32) m<-aov(r~t+b) anova(m) e<-resid(m) p<-predict(m) a<-tapply(r,t, mean) b<-tapply(r,t, sd) d<-ave(r, t, FUN = mean) c<-ave(r, t, FUN = sd) ep<-e/(anova(m)$ "Mean Sq"[2]^0.5) ep1<-e/c plot(e~p, xlab= "Valores Preditos", ylab= "Resíduos") cor(e,p) plot(ep~p, xlab= "Valores Preditos", ylab= "Resíduos Padronizados") cor(ep,p) plot(ep1~p, xlab= "Valores Preditos", ylab= "Resíduos Padronizados") cor(ep1,p) plot(c~d, xlab= "Médias", ylab= "Desvios Padrões") cor(c,d) plot(t, e, xlab= "Tratamentos", ylab= "Resíduos") cor(a,b) cor(d,c)

2 1

Capturar dados via scanner
by Carlos Andrade 16 Mai '11

16 Mai '11

Prezados amigos da lista, Estamos realizando uma enquete com os alunos de Economia sobre a qualidade do ensino. Estamos testando um questionário. Observei que a impressora HP C4480 tudo-em-um, que tem um scanner, quando se troca o cartucho pede para imprimir uma página, com manchas quadradas e listras, e que a coloque no vidro do scanner e de lá a impressora lê e alinha os cartuchos. Então pesquisei na internet como programar o scanner para capturar dados de questionários. Assim pouparíamos o tempo e o aborrecimento da digitação. Nada encontrei de útil. Alguém da lista tem alguma dica ou solução para adaptarmos? De forma, também, a serem lidos os dados para o R. Agradeço 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

7 8

ajuste curva potencial
by Rodrigo Plei 16 Mai '11

16 Mai '11

Prezados, Como ajustar uma curva potencial (power curve) em uma dispersão x, y, mostrando a equação e o R2 no gráfico? Desde já, grato. Rodrigo -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777

2 1

Dados de contagem com análise de experimentos
by Edson Lira 16 Mai '11

16 Mai '11

Bom dia a todos. Estou tentando analisar o banco de dados anexo com a rotina abaixo: jat1<-read.csv2("...\\jg_09.csv") names(jat1) stopifnot(complete.cases(jat1$ngalhos) !=is.na(jat1$ngalhos))#eliminar os NA's dos dados jat<-jat1[complete.cases(jat1$ngalhos),] Gostaria de verificar o efeito da aração e tplan na contagem de galhos(ngalhos) de uma árvore. Tentei fazer a transformação para "normalidade" a transormação não foi possível. Tentei fazer um teste para verificar se os dados tem distribuição de "poisson" com o pacote vcd, não sei se fiz o correto. Alguém poderia dar algumas sugestões? []'s Edson Lira Estatístico Manaus-Amazonas

3 2

Colocar pontos no mapa
by Victor Eduardo 16 Mai '11

16 Mai '11

Alguém aqui já colocou pontos em um mapa através da latitude e longitude? Estou tentando fazer e preciso até então utilizar alguns pacotes, como por exemplo, o pacote rgdal, sendo que o mesmo vive dando erro ao instalar, pois não consigo usá-lo. Será que existe uma outra forma para fazer isso? Abraços, Victor Eduardo

4 6

Extração de data e hora de uma string de extrato bancario
by Walmes Zeviani 14 Mai '11

14 Mai '11

Saudações, Tenho um arquivo com uma coluna que é histórico de movimentação financeira. Os elementos dessa coluna são strings como o exemplo abaixo historico <- c("compra com cartão - 01/04 22:11 lanchonete", "compra com cartão - 02/04 12:51 restaurante") Como posso fazer para criar um vetor contendo só a data, padrão "dd/mm", e outro com a hora, padrão "hh:mm", que esteja no formato de data para que eu possa plotar os valores movimentados em função da data? 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 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

4 5

Formato horas
by eder＠leg.ufpr.br 13 Mai '11

13 Mai '11

Pessoal, gostaria de uma ajuda com datas, por padrão as horas '%H' variam de 0-23 (?strptime) porem tenho uma serie de dados que considera 24 horas, como no exemplo abaixo onde esta data pertence ao dia 20/02, desta forma alguem conhece como representar esta data adequadamente? z <- strptime("20/2/06 24:00:00.000", "%d/%m/%y %:%M:%OS") z Att

2 2

R 2.12 -> R 2.13
by Leonard Assis 13 Mai '11

13 Mai '11

Gostaria de saber se existe muita diferença entre R 2.12 e R 2.13 que justifiquem a migração. Leonard

6 6

Anova modelo para experimentos em faixas - Split Block
by Fernanda Balbino 13 Mai '11

13 Mai '11

Eder, Muito obrigada!! Era esse o modelo que estava procurando. Fernanda Balbino. ________________________________ De: Fernanda Balbino <ferbalbino(a)yahoo.com.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Quinta-feira, 12 de Maio de 2011 16:37:53 Assunto: Anova modelo para experimentos em faixas - Split Block Olá, Gostaria de saber se alguém conhece um exemplo de como inserir no R o modelo matemático da anova para experimentos em faixas. Tenho dois fatores, com 7 e 10 níveis. Muito obrigada pela atenção, Fernanda Balbino.

1 0

Anova modelo para experimentos em faixas - Split Block
by Fernanda Balbino 13 Mai '11

13 Mai '11

Olá, Gostaria de saber se alguém conhece um exemplo de como inserir no R o modelo matemático da anova para experimentos em faixas. Tenho dois fatores, com 7 e 10 níveis. Muito obrigada pela atenção, Fernanda Balbino.

3 2

Blop´s modelo fatorial
by Greice Laureano 12 Mai '11

12 Mai '11

Bom dia, Fiz uma análise fatorial de um Planejamento Fatorial Cruzado 23 em Blocos Casualizados. Utilizei a função aov() para executar a análise. Gostaria de saber se alguém sabe como faço para selecionar: 1° Resíduos Considerando os Blup's de blocos (Resíduos Condicionais). 2° Resíduos Desconsiderando os Blup's blocos (Resíduos Marginais). Obrigado! Att, Greice.

2 2

Error- entrada de dados no R- Urgente
by Carlos Magno 12 Mai '11

12 Mai '11

Auxilio R, Ao tentar inserir dados em formato txt, já configurada pr R ocorre o seguinte error: Erro em scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : linha 1 não tinha 10 elementos Gostaria de saber como sanar esse error pr eu proceguir com a minha analise , pois preciso com urgência Att Carlos Magno

5 4

Exportar dados do R para o formato Latex
by José Weligton Félix Gomes 12 Mai '11

12 Mai '11

Caros, Sou novato no uso do R, no entanto, tenho aprendido bastante com as notas de aulas ou materiais disponíveis na internet. Entretanto, ainda não consegui um tutorial que me ajudasse a exportar o dados do R para o formato Latex. Alguém poderia me ajudar? Agradeço de antemão. Weligton Gomes Mestrado em Economia - CAEN/UFC Universidade Federal do Ceará

7 9

Como calcular intervalos de confiança
by Claudio Carneiro 11 Mai '11

11 Mai '11

Olá a todos, Gostaria de saber como calcular intervalos de confiança para a média de uma amostra. Qualquer ajuda é bem vinda.Obrigado

5 5

Legenda de gráfico
by Danúbia Gomes 11 Mai '11

11 Mai '11

Olá, pessoal Estou tentando colocar uma legenda num mapa seguindo um modelo pré-determinado. No mapa, seguindo uma escala de zero a um, 0 é a cor mais clara e 1 a mais escura. Quando peço para inserir a legenda, ela aparece na vertical, com a cor mais clara em cima. Porém eu preciso inverter a ordem da legenda, de forma que a cor mais escura apareça primeiro. Vocês sabem se existe algum comando específico para isso? Muito obrigada, Danúbia Lacerda

2 2

Verificar independência de 2 variáveis normais
by Alexandro (Yahoo) 11 Mai '11

11 Mai '11

Obrigado, Muito Obrigado Paulo Justiniano, Walmes Zeviani e Leonard Assis pela resposta do tópico: calcular probabilidade com 2 variáveis normais ! Agora, caríssimos participantes da lista R-br: Suponha, que não se sabe, se as 2 variáveis normais (pelo shapiro.test) são independentes. Tem algum comando do software R que testa a independência de 2 variáveis normais ? Ou deve-se criar uma tabela de dupla entrada em intervalos dessas 2 variáveis para usar-se o comando chisq.test ou CrossTable (pacote gmodels)? Agradeço-lhe a todos pela atenção para comigo ! Alexandro Vieira Lopes

5 4

Fontes fa função splom do lattice
by Edson Lira 11 Mai '11

11 Mai '11

Boa tarde, pessoal tenho uma rotina do pacote lattice que gera um panel de scaterplot de um conjunto de variáveis. A função é a splom. splom(~data) Como faço para mudar o tamanho da fonte do nomes das variáveis e dos valores das variáveis? Edson Lira Estatístico Manaus-Amazonas

2 1

R-Serve - Acesso Simultâneo
by Marley Saraiva 11 Mai '11

11 Mai '11

Olá! Tenho um servidor que roda uma ferramenta em java que trabalha junto com o R-Serve. Quando um usuário acessa a ferramenta, os demais não podem usar porque o R-Serve não está aceitando vários acessos simultâneos. O R-Serve aceita acessos simultâneos? Marley Apolinario Saraiva

1 0

Ultima linha de arquivo
by Emerson Cotta Bodevan 11 Mai '11

11 Mai '11

Prezados, outra pergunta na linha da anterior. Tenho um arquivo com dados em que a ULTIMA linha tem lixo. Como faço para ler esse arquivo e ignorar a ultima linha. Pergunto isso porque meus dados sao compostos por 4 colunas e na ultima linha (o lixo) tem apenas um dado (um sinalizador de fim de arquivo. Por ex. -1) Agradeço qualquer ajuda, -- Emerson

2 2

Problemas com o pacote rgdal
by Victor Eduardo 11 Mai '11

11 Mai '11

Gente, estou com problema para instalar esse pacote, ele baixa normalmente, mas quando vou pedir usando library ou require ele diz que o pacote não existe. O que poderia ser? Abraços, Victor Eduardo

2 6

Primeira linha de arquivo texto
by Emerson Cotta Bodevan 11 Mai '11

11 Mai '11

Prezados, Tenho um arquivo texto com meus dados para análise. Esse arquivo é gerado de um programa em C++. A primeira linha é só lixo. Como faço para ler esse arquivo no R de tal maneira que ele ignore a primeira linha? (Por ex.: Se eu fosse usar a função read.table) Abs, -- Emerson

2 2

(sem assunto)
by Elias T. Krainski 11 Mai '11

11 Mai '11

Danúbia Suponha que seu mapa contenha a variável 'MESOREGIAO'. Então faça: bhm <- mg[mg$MESOREGIAO=="METROPOLITANA DE BELO HORIZONTE", ] O objeto 'bhm' é o mapa da Mesoregião Metropolitana de BH. Elias T. Krainski

2 1

Re: [R-br] Exportar dados do R para o formato Latex
by Jose Claudio Faria 10 Mai '11

10 Mai '11

Algumas dúvidas: 1. Você está em que ambiente: Windows ou Linux? 2. Usando quel Editor/GUI? Abs, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica - prof. Titular UESC/DCET/Brasil joseclaudio.faria(a)gmail.com ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

2 1

Mapas
by Danúbia Gomes 10 Mai '11

10 Mai '11

Olá pessoal, Estou fazendo um mapa de Minas Gerais e utilizei o script " http://leg.ufpr.br/~elias/ensino/ce223/aula17.R" Consegui traçar o mapa e colorir uma região específica, mas não sei como faço para traçar apenas as cidades dessa região. Por exemplo, se eu quiser um mapa apenas da Região Metropolitana de BH, não consigo. Com a malha municipal do IBGE, só desenha o mapa todo do estado. Alguém sabe como se faz? Um abraço, Danúbia

1 0

calcular probabilidade com 2 variáveis
by Alexandro (Yahoo) 10 Mai '11

10 Mai '11

#Saudações a todos os usuários do software R, #para se calcular probabilidades da distribuição normal com 1 variável pode-se fazer o seguinte: #Suponha X~N(100,100), logo: pnorm(95, 100, 10);# = P(X < 95) pnorm(110, 100, 10) - pnorm(90, 100, 10); # = P(90 < X < 110) pnorm(95, 100, 10, lower=F);# = P(X > 95) #Gostaria de saber como são calculadas probabilidades com 2 variáveis: #Exemplo: # X ~ N(100,20) e Y ~ N(80,10) em que X e Y são independentes. #P(X>80;Y>80) #P(x<100;Y>85) #P(85<X<105;75<Y<80) Muito Obrigado. Desejo-lhes uma excelente semana ! Alex

4 3

Re: [R-br] Microdados PNAD
by Elias T. Krainski 09 Mai '11

09 Mai '11

Régis, A função read.fwf() é idicada para ler microdados se vc tem MUITA memoria RAM ou se seu arquivo é pequeno (o que não é o seu caso). Com essa função vc lê todas as colunas dos dados, o que pode "estourar" a memória RAM. Eu acho apropriado ler apenas as colunas e linhas que interessam. Por exemplo, se quero estimar a proporção de mulheres que trabalham em segurança, basta eu ler a coluna sexo (e peso amostral) e apenas as linhas referentes aos indivíduos que trabalham com segurança. Eu uso uma função disponível em http://leg.ufpr.br/~elias/ensino/ce223/aula18.R No exemplo, eu considero dados da pnad 2006. Att. Elias T. Krainski

2 1

Redes neurais
by Marcelo Claro de Souza 09 Mai '11

09 Mai '11

Alguém teria um tutorial bem básico sobre redes neurais no R? Muito obrigado.

4 3

Rstudio Muitos dados
by eder＠leg.ufpr.br 09 Mai '11

09 Mai '11

Pessoal gostaria de comportilhar com vocês, Primeiramente gostaria de dizer que o Rstudio tem uma interface muito amigavel, porem estou usando o RStudio Desktop 0.93.89 - Windows XP/Vista/7 em um computador I7, 4gb Ram, 3Tb hd, Windows 7 64 bit e R 2.12.0, Quando gero por exemplo uma matrix de 500 x 500 e mando executa-la ele trava, travando toda a maquina, quando um procedimento destes é feito no RGui tudo roda sem problemas algum, desta forma gostaria de saber se alguem de vocês já teve problema similar com o Rstudio no Windows. Outra coisa que percebi é que quando se gera um objeto por exemplo uma matrix 500x500, avaliando o Rgui com o R executado no Rstudio o Rstudio aloca muito mais memória para o processo, acho isso muito estranho. Att

1 0

Microdados PNAD
by Régis B. de Oliveira 09 Mai '11

09 Mai '11

Boa Tarde a todos. Sou novato na lista e no uso do R. Trabalho com microdados do IBGE, mais especificamente os dados da Pesquisa Nacional por Amostra de Domicílios (PNAD). Alguem da lista tem experiência na manipulação destes dados com o R. Achei muito confusa a explicação que contém no site do IBGE (leitura em R). Preciso de um "empurrãozinho" inicial. Trabalho como o SAS e pretendo migrar. A dificuldade inicial é justamente ler os microdados. No SAS eu costumo ler todo o arquivo txt e salvar na extesão .sas7bdat. A partir daí, faço minhas análises (basicamente regressões e tabelas de frequencias). É isso. Desde já agradeço, -- Régis

5 5

Erros no output do Vim!
by Ivan Bezerra Allaman 08 Mai '11

08 Mai '11

Bom dia senhores!! Ao utilizar o Vim, estou encontrando alguns erros no output. Primeiro, é necessário eu executar várias vezes o comando ";l " para aparecer o output desejado. Segundo, o output está aparecendo pela metade. Por exemplo, quando rodo uma base de dados com 400 linhas por exemplo, na primeira execução só aparece as 50 primeiras, na segunda execução aparece os 150 e assim por diante até aparecer toda a base. Alguém já se deparou com este problema? Qual a solução? Só para complementar, já instalei a última atualização do "screen.vba". Abraços. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Reenvidando o email multcompBoxplot
by Edson Lira 07 Mai '11

07 Mai '11

Boa noite a todos, A rotina abaixo: require(multcompView) difat<-c(3.078444,3.567257,4.023913,3.7525714) names(difat) <- c("A-B", "A-C", "B-C", "A-D") (mcT4 <- multcompTs(difat)) mcT4.1 <- plot(mcT4, label.groups=0.05) (mcL4 <- multcompLetters(difat, Letters=LETTERS)) mcL4.1 <- plot(mcL4, label.groups=0.05) Mostra a seguinte mensagem: Erro em plot.window(...) : valores finitos são necessários para 'xlim' Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf 3: In min(x) : nenhum argumento não faltante para min; retornando Inf 4: In max(x) : nenhum argumento não faltante para max; retornando -Inf 5: In plot.window(...) : "label.groups" não é um parâmetro gráfico []'s Edson Lira Estatístico Manaus-Amazonas

1 0

Erro no multcompBoxplot
by Edson Lira 07 Mai '11

07 Mai '11

Boa noite, a rotina abaixo: require(multcompView) difat<-c(3.078444,3.567257,4.023913,3.7525714) names(difat) <- c("A-B", "A-C", "B-C", "A-D") (mcT4 <- multcompTs(difat)) mcT4.1 <- plot(mcT4, label.groups=0.05) (mcL4 <- multcompLetters(difat, Letters=LETTERS)) mcL4.1 <- plot(mcL4, label.groups=0.05) Está mostrando a mensagem abaixo: Erro em plot.window(...) : valores finitos são necessários para 'xlim' Além disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retornando Inf 2: In max(x) : nenhum argumento não faltante para max; retornando -Inf 3: In min(x) : nenhum argumento não faltante para min; retornando Inf 4: In max(x) : nenhum argumento não faltante para max; retornando -Inf 5: In plot.window(...) : "label.groups" não é um parâmetro gráfico O erro parece ser nos limites, não consegui identificar aonde inseri-los. Alguém teria alguma idéia? Edson Lira Estatístico Manaus-Amazonas

1 0

Re: [R-br] Microdados PNAD
by Régis B. de Oliveira 07 Mai '11

07 Mai '11

Edson, obrigado pela atenção. Os arquivos que utilizo estão disponíveis em (dados.zip): ftp://ftp.ibge.gov.br/Trabalho_e_Rendimento/Pesquisa_Nacional_por_Amostra_de_Domicilios_anual/microdados/2009/ Eu não preciso especificar para o R onde começa e onde termina cada variável do arquivo? Aproveito envio outra questão/dúvida. Você sugere o uso de softwares tipo o Tinn-R? Obrigado, Régis.

3 3

Erro no multcompBoxplot
by Edson Lira 07 Mai '11

07 Mai '11

Boa noite a todos, estou rodando um gráfico do pacote multcompView, tenho a variável altura (h), e a variável at que é a combinação dos níveis de dois fatores (ap,am, nap,nam). Quando rodo, aparece na tela gráfica, mas some rapidamente e me mostra a mensagem de erro. Erro em plot.window(...) : valres finitos são necesários para 'xlim' Alem disso: Warning messages: 1: In min(x) : nenhum argumento não faltante para min; retrnando Inf 2: In min(x) : nenhum argumento não faltante para max; retrnando -Inf Alguém poderia e dar uma dica? Edson Lira Estatístico Manaus-Amazonas

3 2

Modelo de efeito misto com correlacao AR1
by Bernardo Rangel Tura 06 Mai '11

06 Mai '11

Caros colegas, Estou precisando de um ajuda para um modelo de efeito misto com correlação AR1 A auto correlação se justifica pois os pacientes devem apreender e obter melhores resultado de um bloco para outro então usei os seguinte comando ------------------------------------------------------- base <- read.csv('jogov2.csv') require(nlme) modelo <- lme(valor~Gr+Idd+Escol,random=~1|Bloco,correlation=corAR1(, form=~1|Bloco) ,data=base) summary(modelo) ----------------------------------------------------- Que resultou em ----------------------------------------------------- Linear mixed-effects model fit by REML Data: base AIC BIC logLik 3928.017 3958.571 -1957.009 Random effects: Formula: ~1 | Bloco (Intercept) Residual StdDev: 2.480127 6.846336 Correlation Structure: AR(1) Formula: ~1 | Bloco Parameter estimate(s): Phi 0.1407996 Fixed effects: valor ~ Gr + Idd + Escol Value Std.Error DF t-value p-value (Intercept) -1.7057456 1.8249072 577 -0.934703 0.3503 GrB 1.7242603 0.6606029 577 2.610131 0.0093 Idd -0.0479666 0.0178547 577 -2.686502 0.0074 Escol 0.3441736 0.0926290 577 3.715615 0.0002 Correlation: (Intr) GrB Idd GrB -0.068 Idd -0.415 -0.101 Escol -0.602 -0.227 0.001 Standardized Within-Group Residuals: Min Q1 Med Q3 Max -3.16689925 -0.56838036 -0.09641692 0.55523377 3.66357685 Number of Observations: 585 Number of Groups: 5 ------------------------------------------------------- Agora vem o problema eu preciso saber se o efeito aleatório tem um média diferente diferente entre o grupo B e o grupo A e além disso se o valor de Phi é diferente para os dois grupos Alguém pode me ajudar ? Ps. Base em anexo -- []s Tura

4 4

erro em "lme"
by Samuel Carvalho 05 Mai '11

05 Mai '11

Novamente gostaria de saber se alguem ja viu este tipo de erro em na função augPred #Script plot(augPred(mod.lme, level=c(0,1))) "Erro em `[[<-.data.frame`(`*tmp*`, nm, value = c(28L, 28L, 28L, 28L, : replacement has 9558 rows, data has 29" ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

2 1

erro em "lmList"
by Samuel Carvalho 05 Mai '11

05 Mai '11

Caros Alguem ja conviveu com este tipo de erro na função lmList? 'Erro em na.fail.default(data) : missing values in object' Pensei que pudesse ser algo relacionado com celulas vazias (NA) mas não é. Valeu ==================================== Samuel P. C. Carvalho Mestre em Ciências Florestais [UFLA] Doutorando em Recursos Florestais [ESALQ/USP] =============================================

1 0

Extrair um vetor de diversos ajustes
by Gustavo Marcatti 05 Mai '11

05 Mai '11

Bom dia a todos! Quando desejo extrair os coeficientes de diversos ajustes, por exemplo, procedo da seguinte maneira: x<-1:60 y<-rep(sort(abs(rnorm(20))),3) caso <- rep(1:3,each = 20) dados <- data.frame(caso,x,y) exp1 <- lapply(split(dados,dados["caso"]),FUN=function(data){ coefs <- coef(lm(y ~ x, data=data)) unlist(coefs) }) do.call("rbind",exp1) Minha dúvida é como faço para extrair um vetor de varios casos, assim como no exemplo abaixo: exp2 <- lapply(split(dados,dados["caso"]),FUN=function(data){ preditos <- fitted(lm(y ~ x, data=data)) list(preditos) }) do.call("rbind",exp2) O resultado esperado seria um único vetor de tamanho 60 x 1 contendo todos os valores preditos para cada ajuste. Obs1: Esse exemplo é apenas para ilustrar, não é o que realmente estou precisando. Obs2: Os casos tem tamanhos variados, não exatamente 20 observações. Gustavo Marcatti Eng. Florestal UFV

3 3

Re: [R-br] Modelo de efeito misto com correlacao AR1
by Ivan Bezerra Allaman 05 Mai '11

05 Mai '11

Tura, Em síntese é o seguinte: Quando trata um efeito como aleatório, ou se ele tem esse perfil, quer dizer que o grupo que foi tratado como aleatório tem uma distribuição de probabilidade além do erro é claro. Isso é perfeito, pois agora podemos extrapolar nossa inferência para uma população, ao contrário de quando tratamos um grupo como efeito fixo, em que podemos somente fazer inferências para aquele grupo de observações por meio de comparações de médias. Suponhamos um exemplo bem esdrúxulo, que você avaliou o índice de infecção hospitalar por um novo vírus mortal chamado SAS em dez hospitais públicos. Se hospitais for tratado como efeito fixo, a sua inferência será restrita apenas ao grupo de hospitais em estudo. Se for tratado como aleatório, a inferência pode ser extrapolada para todos os hospitais daquela região, estado ou país. Os componentes da variância nada mais é do que a Esperança do Quadrado Médio daquela fonte de variação. Por exemplo, o componente da variância do erro é: E(QMerro) = sigma² erro, pois ele é aleatório. A componente da variância de hospital caso ele seja fixo, e supondo um delineamento simples é: E(QMhospital) = repetição Somatório (hospital)² + sigma² erro. Se ele for aleatório, o grupo passa a ter uma variância inerente a ela. E(QMhospital) = repetição sigma² hospital + sigma²erro. É claro que tudo isso irá depender do seu objetivo. Pode ser que o seu interesse não esteja no efeito aleatório. Para maiores detalhes veja: "Applied Mixed Models in Medicine". Helen Brown e Robin Prescott. Caso não tenha o livro me avise que eu lhe envio. Abraço. Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

1 0

Desdobramento da média
by rosanabarbosa＠inpa.gov.br 05 Mai '11

05 Mai '11

Boa noite colegas, Estou utilizando anova para analisar dois fatores. Quando anova mostra a interação destes dois fatores não consigo fazer o seu desdobramento. Eis minha dúvida que compartilho com vcs. No R, como fazer o desdobramento da média, se possível mostrando as letras qua a diferencia dentro de cada fator. Grata, Rosana

2 2

Re: [R-br] Coeficientes dos modelos polinomiais ajustados
by Wenceslau 05 Mai '11

05 Mai '11

Caro Walnes, Caro Gustavo, as dicas de remover o intercepto do modelo ( -1); e comando coefs resolveram a questão de conseguir os coeficientes. Muito obrigado. Saudações estatisticas Wenceslau Em 4/5/2011 12:00, r-br-request(a)listas.c3sl.ufpr.br escreveu: > 5. Re: Coeficientes dos modelos polinomiais ajustados

1 0

Modelo de efeito misto com correlacao AR1
by Ivan Bezerra Allaman 04 Mai '11

04 Mai '11

Tura, Não seria mais informativo você utilizar os componentes da variância, uma vez que se trata de efeitos aleatórios e não fixos? Allaman (S,f,P) M.Sc Ivan Bezerra Allaman Zootecnista Doutorando em Produção Animal/Aquicultura - UFLA email e msn - ivanalaman(a)yahoo.com.br Tel: (35)3826-6608/9900-2924

2 1

Atualização do R no Ubuntu 11.04
by Junior Beleti 04 Mai '11

04 Mai '11

Olá, estou utilizando o Ubunto 11.04 e gostaria de atualizar a versão do R, pois quero instalar o pacote "threads" e ele é dependente da versão 2.13.0. Através do arquivo INSTALL da versão 2.13.0 tentei o que está lá: ./configure make mas recebi a seguinte mensagem: make: *** Nenhum alvo indicado e nenhum arquivo make encontrado. Pare. Também tentei através dos comandos: sudo apt-get update sudo apt-get upgrade e não obtive sucesso. Já procurei nas discussões anteriores e não encontrei tópico com exatamente esse problema. Nos tópicos encontrei que pode-se ter problemas na atualização quanto a pacotes. Gostaria de saber como proceder para não ter problemas. Peço desculpas se for algo óbvio, pois estou iniciando minhas pesquisas com R. Obrigado. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

3 2

pacote em zip no Linux
by Carlos Andrade 04 Mai '11

04 Mai '11

Prezados Membros da Lista Já o uso o R em algumas tarefas. Estou iniciando a trabalhar com os microdados da Pnad. No site do IBGE, junto com os microdaddos vem um pacote do R chamado IBGEPesq, salvo no formato zip. Como eu faço para usar no R no sistema operacional Linux? Instalei o Debian 6.0 e baixei o R 2.13. -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Centro de Humanidades - Univesidade Federal de Campina Grande. UFCG/CH/UAECON

5 6

Coeficientes dos modelos polinomiais ajustados
by Wenceslau 04 Mai '11

04 Mai '11

Caros colegas, estou tentando fazer uma analise comparando modelos lineares da evolução do pH vs tempo, num experimento com diferentes doses de carbonato de cálcio (7 doses). Tenho este Tenho este script que faz a analise e também os gráficos. Mas não consegui obter os coeficientes de cada modelo individualmente (dose 0 vs tempo; dose 4 vs tempo, etc). Agradeço se puderem me ajudar ## Experimento pH vs tempo dados <- read.table("dosagem.txt", h=T) names(dados) attach(dados) B <- factor(Rep) D <- factor(Dose) T <- factor(Tempo) modelo <- lm(Leitura ~ D*T+B:D) anova(modelo) # Pode-se considerar o modelo fatorial diretamente fatorial <- lm(Leitura ~ D*(I(Tempo)+I(Tempo^2))) anova(fatorial) summary(fatorial) #plot(fatorial) #SELECIONANDO O MODELO QUADRATICO E GERANDO INTERVALOS DE CONFIANÇA DE 0.95 preditos <- predict(fatorial,interval="confidence") media <- preditos[,1] LI <- preditos[,2] LS <- preditos[,3] interaction.plot(Tempo,D,media,ylab="pH",xlab="Dias", ylim=c(4,8), pch=0:7, legend=T) #Plotando os pontos e as curvas ajustadas - Grafico 1 par(new=TRUE) interaction.plot(Tempo,D,Leitura,ylab="pH",xlab="Dias", type="p",pch=0:7,col=3:9, ylim=c(4,7.5),legend=T) par(new=TRUE) interaction.plot(Tempo,D,media,ylab="pH",xlab="Dias", ylim=c(4,7.5), pch=0:7, legend=F) #Plotando a curva ajustada e os IC par(new=TRUE) interaction.plot(Tempo,D,ylab="pH", media,ylim=c(4,7.5),legend=F) par(new=TRUE) interaction.plot(Tempo,D,ylab="pH",xlab="Dias",LI,col=2,ylim=c(4,7.5),legend=F) par(new=TRUE) interaction.plot(Tempo,D,ylab="pH",xlab="Dias", LS,col=3,ylim=c(4,7.5),legend=F) legend("topleft", c("IC - Superior", "IC - Inferior"), col=c("green", "red"), title="Legenda", lty=1, lwd=2)

2 1

Problemas com a execução do
by Marcelo Costa 03 Mai '11

03 Mai '11

Estou usando o pacote RgoogleMaps e tudo está funcionando bem, mas quando utilizo o comando PlotPolysOnStaticMap(TMap, poly, lwd=.25, col=col, add=T, border=NA) aparece a sequência Called from: FUN(data[x, , drop = FALSE], ...) Browse[1]> Pelo que pesquisei, me parece que a rotina entra no modo "debug". Só que tenho um loop que irá chamar este comando 2500 vezes. Até o momento tenho uma solução: manter o indicador pressionando a tecla <enter> até que todo o loop seja executado. Alguém teria outra alternativa para que o comando seja executado sem entrar em um modo de "espera"? Marcelo A Costa - UFMG

3 2

como trabalhar os dados desbalanceados no R
by Gilson Sanchez 03 Mai '11

03 Mai '11

Pessoal, mais uma vez consultando vc's tenho um amigo que esta avaliando 4 substratos, 5 cultivares e 4 repetições e muitas variáves de estudo. então dentro das avaliações ele perdeu alguns dados, no SAS esses dados perdidos são inseridos como ponto (.) no R são inseridos como NA, mas quando eu tento fazer uma anova eu não consigo, sai erro, então como eu posso trabalhar com esses tipos de dados dados<-read.table("dados.estatistica.txt", header=TRUE) dados attach(dados) anova <- lm(PTF~SUBST+CULT+REPET+SUBST:CULT, data=dados) > anova <- lm(PTF~SUBST+CULT+REPET+SUBST:CULT, data=dados.y) Erro em storage.mode(y) <- "double" : inválido mudar o modo de armazenamento de um fator Além disso: Mensagens de aviso perdidas: In model.response(mf, "numeric") : usando type="numeric" com um fator resposta será ignorada outro, eu utlizaria o erro tipo III na anova? -- MSc. Gilson Sánchez Chia Laboratório de Fisiologia Vegetal Embrapa Amazônia Ocidental Fone: (92) 3303-7841 --- 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.

3 8