- R-br - DINF UFPR

teste de hipoteses
by Daniela Recchia 13 Abr '11

13 Abr '11

Boa tarde pessoal, desejo fazer um teste de hipoteses bem simples, porem nao consigo visualizar como deveria fazer isso no R, uma vez que os testes que ja fiz sempre foram utilizando um modelo (lm) e a funcao linearHypothesis. Nesse caso tenho os dados da seguinte forma: > cov V1 V2 V3 V4 V5 V6 45 65 81 55 31 0 41 80 73 51 36 0 42 79 74 40 45 0 43 77 76 62 * 0 Onde as colunas V1...V6 sao niveis de agua adicionados a uma planta, e cada linha eh uma caixa onde o experimento foi realizado, nesse caso 4 caixas para cada nivel de agua. Os valores sao o numero de semestes que germinaram. Quero testar a hipotese de que o nivel de agua nao influencia na germinacao das sementes. A principio pensei em usar anova, mas nao tenho claro como devo escrever meu modelo. Alguem poderia me dar uma dica? Abracos e bom dia! Daniela -- Daniela Rodrigues Recchia Master Student of Statistics - Technische Universität Dortmund. "Not everything that can be counted counts; and not everything that counts can be counted.” George Gallup

2 1

Re: [R-br] Res: Blocos incompletos no R
by Julia Maria 13 Abr '11

13 Abr '11

Walmes obrigado pela ajuda. A proposito, nao sei se deu tempo, mas naquele topico antigo vc pediu pra eu lhe mandar um .xls com os dados, ta lembrado? Voce chegou a dar uma olhada? Estou perguntando porque nao sei se recebi alguma mensagem sua, o modo com que essas mensagens chegam aqui estao meio bagunçadas. POr exemplo essa duvida que vc me respondeu hoje, eu usei o topico que eu criei la atras pra responder, mas ela apareceu sozinha na lista.

1 0

Blocos incompletos no R
by julia_est 13 Abr '11

13 Abr '11

Pessoal, tomei conhecimento dessa lista através do site do leg. Por acaso alguém saberia me informar se no R existem pacotes ou funções já prontas para analisar os diferentes tipos de Blocos incompletos? Tipo I, Tipo II, Tipo III... Andei conversando com alguns colegas e eles apontaram como saída fazer tudo a mão, montar aquelas tabelas auxiliares e por aí vai. Eu queria saber se não tem como utilizar os comandos normais para os outros delineamentos, criando modelos via LM ou aov. Alguém poderia me ajudar? Grata Julia -- View this message in context: http://r-br.2285057.n4.nabble.com/Blocos-incompletos-no-R-tp3402407p3402407… Sent from the R-br mailing list archive at Nabble.com.

7 12

loop em nchar
by Daniel Marcelino 13 Abr '11

13 Abr '11

Caros, etou tentando obter o número de caracteres para cada célula para descobrir se trata-se de CNPJ ou CPF, mas acho que meu looping está com problemas pois não para de computar e não me retorna nada. O que está errado? for (i in CPF_CNPJ) { receitas2002$len = nchar(CPF_CNPJ) } Daniel

3 3

Importação correta de palavras acentuadas com a gdata::read.xls()
by Walmes Zeviani 13 Abr '11

13 Abr '11

Saudações, Tem pouco tempo que importo arquivos de extensão *xls. Estou com um problema que ao importar um arquivo possuo nomes de cidades brasileiras com acentos no nomes. Ao importar, as letras acentuadas são representadas de outra forma. Quero saber como importar corretamente ou se é possível devolver os nomes certos após importação. Veja o CMR que preparei require(gdata) # contém a função read.xls() para importar planilhas teste <- read.xls("http://www.leg.ufpr.br/~walmes/foo/teste.xls", sheet=1, method="tab") head(teste) # abaixo o print desse comando Sigla Codigo mun a2000 1 AC 1200013 Acrel\xe2ndia 1042.38 # â trocado por \xe2 2 AC 1200054 Assis Brasil 429.22 3 AC 1200104 Brasil\xe9ia 2243.25 4 AC 1200138 Bujari 786.68 5 AC 1200179 Capixaba 603.84 6 AC 1200203 Cruzeiro do Sul 9363.32 Os demais acentos possuem o seu "\xe?" correspondente. Alguém tem uma solução? Outro pacote/função? Meu sistema é Ubuntu 10.04. Os dados são um subconjunto dos dados do ipeadata, baixados no formato.xls. Outra coisa que percebi é que se o nome das colunas tiver acentos, a importação falha, ex "municipio" (sem acento) importa, mas com "município" (com acento) no nome da coluna não importa. Como evitar isso sem ter que sair renomeando as colunas manualmente? Planilhas de sites coleta de dados do Brasil sempre tem acentos por causa do nomo idioma. Minhas buscas no nabble não retornaram resultado. # sheet 2 está escrito "Município" como nome da coluna teste <- read.xls("http://www.leg.ufpr.br/~walmes/foo/teste.xls", sheet=2, method="tab") Erro em make.names(col.names, unique = TRUE) : string multibyte inválida 3 Certo sua atenção, agradeç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 ==========================================================================

8 12

Anova multivariada.
by Edson Lira 13 Abr '11

13 Abr '11

Boa noite a todos. Estou replicando um exemplo de um livro que verifica a dependência entre o tempo de viagem(y) ~ número de entregas(x1) +número de milhas percorridas. Quando peço a anova no R. anova(mod.cam), o R me dá a saída abaixo. Analysis of Variance Table Response: cam$tem_via Df Sum Sq Mean Sq F value Pr(>F) cam$n_entrega 1 9.0449 9.0449 27.535 0.001190 ** cam$milhas 1 12.5556 12.5556 38.222 0.000453 *** Residuals 7 2.2994 0.3285 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Eu gostaria da ANOVA com o Df da Regressão e do Resíduo, ou seja a anova reduzida. Alguém poderia me indicar qual a função ou rotina para isso? Em anexo o banco de dados em csv e a rotina. []' s Edson Lira Estatístico Manaus-Amazonas

2 2

estimadores
by Killa p'unchay Morales 12 Abr '11

12 Abr '11

--- Em ter, 12/4/11, Killa p'unchay Morales <killa.punchay(a)yahoo.com.br> escreveu: De: Killa p'unchay Morales <killa.punchay(a)yahoo.com.br> Assunto: help Para: "RRRR" <r-br-request(a)listas.c3sl.ufpr.br> Data: Terça-feira, 12 de Abril de 2011, 3:04 olá amigos estou com dificultade de encontrar os 4 parâmetros desta função de verossimilhança....primeiro to tentando derivar uma somatoria da acumulada da Gama que chamei de XX com respeito de alpha e derivar a somatoria do log da acumulada da Gama com respeito de alpha = YY, fez analiticamente as derivadas com respeito dos 4 parâmetros......não sei utilizar o optim para maximizar essa função. estou anexando o que fez. Obrigada alguem me ajude.....

2 1

Problemas write.table
by Alan Rodrigo Panosso 12 Abr '11

12 Abr '11

Caros amigos, Como salvar a lista gerada pelo código abaixo, utilizando a função write.table? Estou encontrando o erro abaixo! set.seed(200) dados<-rnorm(8836,36,6) NL=94 NC=94 Xmin=0 Xmax=64 Ymin=0 Ymax=64 Zmin=min(dados) Zmax=max(dados) ma<-matrix(dados,ncol=NC) mc<-matrix(c(NL, Xmin,Ymin,Zmin,NC,Xmax,Ymax,Zmax),ncol=2) lac<-list(mc,ma) write.table(lac,"C://r//Saida.txt",sep="\t",col.names=F,row.names=F) Erro em data.frame(c(94, 0, 0, 12.7834627990016, 94, 64, 64, 61.968478487737 : arguments imply differing number of rows: 4, 94 Abraços a todos. Alan Rodrigo Panosso

3 6

Correlação de Spearman - "Não é possível calcular valores p exatos quando há empates"
by Edson Ferraz Evaristo de Paula 12 Abr '11

12 Abr '11

Boa noite! O que seria essa "warning message" que aparece quando rodo cor.test (pelo método de spearman)? "Não é possível calcular valores p exatos quando há empates" Obrigado -- Edson F. Evaristo de Paula (41) 9669-3646 (42) 8813-7268

3 2

Como responder a um tópico da lista!
by Ivan Bezerra Allaman 11 Abr '11

11 Abr '11

Boa noite senhores! Me parece que muitos (inclusive eu) estavam tendo problemas ao responder um determinado tópico da lista, gerando inúmeros tópicos repetidos. Eu particulamente, não estou a contento da maneira que está, em só responder os tópicos por email, quando poderíamos responder tais tópicos diretamente pela lista, como acontece no R internacional. Mais já foi justificado o porque disto. Logo, para evitarmos estes inúmeros tópicos repetidos, o que polui a lista, ao responder tal tópico, basta fazer o seguinte: Suponhamos que eu queira responder o tópico: [R-br] Anova duas vias com medidas repetitivasBasta enviarmos para a lista desta maneira: Para: r-br(a)listas.c3sl.ufpr.br Cc: Assunto: Re: [R-br] Anova duas vias com medidas repetitivas Pronto, basta colocar "Re:" na frente do título do tópico que sua resposta irá para o mesmo tópico da pessoa que a enviou, evitando deste modo outro tópico com o mesmo nome. 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