Outubro 2012 - R-br

[off-topic] Concurso público 03 vagas para Estatística UTFPR
by Walmes Zeviani 22 Out '12

22 Out '12

Saudações, Encaminho link do edital da UTFPR sobre divulgação de abertura de 03 vagas para Professor de Estatística dedicação exclusiva. O concurso tem vagas para outras áreas do conhecimento também. http://www.utfpr.edu.br/concursos/campi/ct/arquivos-gerais/Edital_068_2012_… Aos interessados, boa sorte. 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

script verificar dias completos em um ano
by luc.souza＠usp.br 22 Out '12

22 Out '12

Boa tarde! Existe algum script que verifique em uma série histórica climática, por exemplo, se em um determinado ano tenha todos os dias completos e insira esses dias faltosos automaticamente? Obrigado Lucas Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutorando Engenharia de Sistemas Agrícolas - ESALQ/USP tel: (19) 8160-7068

4 8

Re: [R-br] Blocos com algumas repetições do material
by Eder David Borges da Silva 22 Out '12

22 Out '12

> Boa tarde Pessoal, > > Tenho experimentos em blocos(Locais) e dentro deles gostaria de por > algumas de repetições de alguns tratamentos para estimar a variância dentro > do bloco, que desconfio ser grande. Por exemplo, tenho 5 Cultivares > gostaria de repetir 3 delas dentro do bloco, tendo ao todo 8 parcelas, o > ideal seria ter uma repetição de tudo, porem na pratica isso não > é possível. > Minhas duvidas são: > 1) Estou definindo de maneira correta essas repetições dentro do modelo? > 2) Como o método de calculo dos dois modelos abaixo é igual para os > efeitos fixos, por que o intercepto esta com valores diferentes? > 3) Neste dois casos(m0 e m1, abaixo) o BLUP seria a média mais o efeito da > cultivar? > Obrigado pela ajuda > > require(lme4) > ### Simulando dados > dados <- > expand.grid(Variedade=factor(c(letters[1:3],letters[1:5])),Local=factor(1:5)) > dados$Rep <- factor(rep(c(2,2,2,1,1,1,1,1),t=5)) > dados$Resp <- sort(rnorm(nrow(dados))) > dados > > ### Modelos > # Sem efeito dentro do local > m0 <- lmer(Resp~1+(1|Local)+(1|Variedade),dados) > summary(m0) > ranef(m0) > # Com efeito dentro do local > m1 <- lmer(Resp~1+(0+Rep|Local)+(1|Variedade),dados) > summary(m1) > ranef(m1) >

4 6

Carlos Antonio Soares de Andrade enviou-lhe o seguinte abaixo-assinado.
by Petição Pública 21 Out '12

21 Out '12

1 0

read.csv
by Edson Lira 21 Out '12

21 Out '12

Caros amigos, estou tentando ler um arquivo csv com a rotina abaixo. h<-read.csv("/Users/edsonlira/Edinho/Edinho/Consulta/Ciesa/BSC_2307.csv") está me dando esse erro. Erro em read.table(file = file, header = header, sep = sep, quote = quote, : more columns than column names O arquivo com algumas linhas está postado no endereço abaixo. Estou tentando ler no MacBookPro http://www.datafilehost.com/download-996309b6.html Boa noite a todos. Edson Lira Estatístico Manaus-Amazonas

2 1

dúvidas sobre utilização de disitribuição beta
by Nicolay Cunha 21 Out '12

21 Out '12

Olá a todos, Tenho uma variável dependente (Y) que varia de 0 a 1 e quero testar em relação a duas variáveis categóricas (X1 com três níveis e X2 com 7 níveis). A categórica X2 transformo em números para ganhar graus de liberdade. Os dados não atendem nenhum tipo de premissa de linearidade, portanto tive que buscar outras alternativas para analisa-los. Ao estudar sobre a melhor forma de ajustar esses dados, encontrei o artigo ( http://cran.r-project.org/web/packages/betareg/vignettes/betareg.pdf) sobre o tema, e os exemplos dados aparentam encaixar com o meu caso. No entanto meu Y tem valores inteiros (0 e 1) e a distribuiçãoa aparenta não aceitar isso, para tal transformei todos os 0s em 0.0001 e 1s em 0.9999. Ainda tenho dúvidas sobre a interpretação dos dados (inclusive o papel do phi), utilização e escolha da melhor forma de gerar o modelo da forma que estou propondo, logo gostaria de sugestões e orientações sobre como proceder com essas análises. Abaixo um CMR com um exemplo do que estou fazendo onde aparece uma mensagem de aviso perdida ao estimar o modelo. *"In betareg.fit(X, Y, Z, weights, offset, link, link.phi, type, control) : no valid starting value for precision parameter found, using 1 instead"* ############################################################################# X1<-factor(rep(c("A", "B", "C"),50)) ## categoria 1 X2<- rep(sample(1:7), len=150 ) ## categoria 2 Y<-rep(sample(0:100)/100, len=150) ## minha variável dependente em proporção Y[Y=='1']<-0.99999 # transformo 1 em 0.99999 pois a distribuição beta não aceita 0 e 1, somente o intervalo entre eles Y[Y=='0']<-0.00001 # transformo 0 em 0.00001 library(betareg) beta1<-betareg(Y~X1+X2, link="loglog") coef(beta1) summary(beta1) beta2<-betareg(Y~X1+X2|X2, link="loglog") ## adicionando X2 como um regressor adicional coef(beta2) summary(beta2) AIC(beta1, beta2, k = log(nrow(data.frame(X1,X2,Y)))) ############################################################################## Agradeço a ajuda! Nicolay -- Nicolay Leme da Cunha Biólogo, Mestre, Doutorando em Ecologia e Conservação Universidade Federal de Mato Grosso do Sul, 79070-900 Campo Grande, MS, Brasil E-mail: nicolaycunha(a)gmail.com lattes.cnpq.br/5916316648872099 -- Nicolay Leme da Cunha Biólogo, Mestre, Doutorando em Ecologia e Conservação Universidade Federal de Mato Grosso do Sul, 79070-900 Campo Grande, MS, Brasil E-mail: nicolaycunha(a)gmail.com lattes.cnpq.br/5916316648872099

4 5

ggplot2 ou Lattice
by Leonard de Assis 20 Out '12

20 Out '12

Qual a mais tranquila de usar segundo a opinião de vocês Eu atualmente uso a ggplot2 e acho meio sacal a forma de lidar com ela. Mas já acostumei com algumas coisas. Estou lendo ultimamente sobre elas, mas não cheguei a nenhuma conclusão -- []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com

5 6

problemas com planilhas Excel xlsx
by Rodrigo Plei 20 Out '12

20 Out '12

Prezados, Ha tempos venho tendo problemas com a digitação/edição/uso de dados digitados em planilhas do Excel 2010 (formato . xlsx). Quando salvo e uso os dados em .txt ou .csv, o R usa os dados em problemas. No entanto, quando abro esses mesmos arquivos novamente no Excel 2010, o programa bagunça todos os valores que originalmente eram números fracionados Por exemplo, se tenho em uma coluna qualquer uma valore de, digamos, 5.677899 na planilha original, quando abro a mesma planilha em formato .txt ou . csv esse mesmo valor se transforma em 5.677.899. E esse arquivo, se usado posteriormente no R, gera resultados errôneos. Alguém já passou por isso e teria alguma dica para compartilhar? Obrigado desde logo. R. -- =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 Google Scholar: http://scholar.google.com.br/citations?user=zX_EezEAAAAJ&hl=pt-BR

5 5

mclapply
by clara luz 18 Out '12

18 Out '12

Olá galera do R! Tudo bom? Gostaria de saber se alguem conhece ou já trabalhou com a função mclapply, "a Parallel version of lapply". Segundo sua descrição, deve-se inserir como argumento um vetor X, sendo que para cada elemento de X será aplicada a função que eu desejo (FUN): mclapply(X, FUN, ..., mc.preschedule = TRUE, mc.set.seed = TRUE, mc.silent = FALSE, mc.cores = getOption("cores"), mc.cleanup = TRUE) Porém eu gostaria de aplicar FUN a uma matriz inteira. Ou seja, eu precisaria que cada elemento de X fosse uma matriz... X<-(matrizA, matrizB, matrizC, ...) Não sei se me fiz entender, mas será que isso é possível? Abraços eRRantes! Clara Luz B. Sant'Anna Bióloga, Bacharel e Licenciada Mestranda em Ecologia - Unicamp Laboratório de Biogeografia da Conservação - UFG (19) 9233.1062

3 3

ordenar colunas - seleção expandida
by marcelo claro de souza 18 Out '12

18 Out '12

Olá pessoal, Preciso reordenar uma coluna de um conjunto de dados, mas gostaria que todas as outras fossem reordenadas juntamente a primeira em estilo seleção expandida, assim como no excel. Tentei o comando sort , mas só consegui ordenar uma coluna de vez. Alguma dica? Muito obrigado. Abraço Marcelo

6 11