Junho 2011 - R-br - DINF UFPR

Rótulo de Dados
by Enéias Façanha 30 Jun '11

30 Jun '11

Olá, Como faço para inserir os rótulos de dados junto aos respectivos pontos no exemplo abaixo? x<-c(5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22) y<-c(2,3,5,7,2,5,4,9,1,7,2,1,0,0,1,0,0,1) plot(x,y,type="p",xaxs="i",yaxs="i",xlim=c(0,23),ylim=c(0,10)) grid() axis(1,1:23) axis(2,1:10) Obrigado -- Enéias Façanha

2 2

[OFF-TOPIC] logotipo do grupo!!
by Ivan Bezerra Allaman 29 Jun '11

29 Jun '11

Boa noite senhores!! Gostaria de saber quando será aberto a enquete para logotipo do grupo. O que está faltando? 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

quadrado latino
by Daniela Recchia 29 Jun '11

29 Jun '11

Boa tarde pessoal, tenho o seguinte conjunto de dados: > dds V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 1 E, 29.0 F, 14.5 D, 20.5 A, 22.5 B, 16.0 C, 6.5 2 B, 17.5 A, 29.5 E, 12.0 C, 9.0 D, 33.0 F, 12.5 3 F, 17.0 B, 30.0 C, 13.0 D, 29.0 A, 27.0 E, 12.0 4 A, 31.5 D, 31.5 F, 24.0 E, 19.5 C, 10.5 B, 21.0 5 D, 25.0 C, 13.0 B, 31.0 F, 26.0 E, 19.5 A, NA 6 C, 12.2 E, 13.0 A, 34.0 B, 20.0 F, NA D, NA Existem 6 variedades de nabo que foram "arranjados" em um quadrado lation. A resposta do experimento eh o peso do nabo. Existem 3 valores que nao estao disponiveis, pois a planta foi danificada por vandalos. Desejo sabe ser se a sequencia dos fatores influencia no resultado. Meus fatores aqui sao A;B;C;D;E;F; Pensei em reescrever meus dados de tal forma que eu tenha: Fatores Peso E 29 B 17.5 F 17 ... e depois ajustar um modelo aov peso~fatores e por fim fazer um HSD.test. Porem nao sei se essa ideia funciona aqui para esse experimento de quadrado latino, pois nao sei se estaria alterando a sequencia quando analiso de tal forma. Depois refaco a analise para uns certos valores que devem ser substituidos onde se tem NA, mas a analise sera a mesma, ne? Alguem teria uma sugestao? Boa semana a todos, abracos -- 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

3 3

multcomp lmer contraste interação..
by Eder David Borges da Silva 29 Jun '11

29 Jun '11

Pessoal da lista: Minha duvida: require(lme4) require(multcomp) ## Exemplo no lme4 moodificado fm2 <- lmer(Reaction ~ Days/Subject + (0+Days|Subject), sleepstudy) summary(fm2) # Tenho interesse neste contraste por exemplo: summary(glht(fm2,c('Days:Subject3-Days:Subject1 = 0')) ) ## Funcionou... ## Agora com meus dados querro coisa semelhante ### Agora com meu conjunto de dados: DADOSS <- expand.grid(sitio=1,ponto=as.character(c('SS')),serie=2,camara=as.character(LETTERS[1:5]),tempo=c(0,15,30,45),armazen=c('VIDRO','HEAD','PP')) DADOSS <- rbind(DADOSS,expand.grid(sitio=1,ponto=as.character(c('SL')),serie=2,camara=as.character(LETTERS[6:10]),tempo=c(0,15,30,45),armazen=c('VIDRO','HEAD','PP'))) DADOSS <- rbind(DADOSS,expand.grid(sitio=1,ponto=as.character(c('SL','SS')),serie=2,camara='Z',tempo=1:2,armazen=c('VIDRO','HEAD','PP'))) DADOSS$ID <- paste(DADOSS[,1],DADOSS[,2],DADOSS[,3],DADOSS[,4],DADOSS[,5],DADOSS[,6],sep='',collapse = NULL) DADOSS$camaraM <- paste(DADOSS$sitio,DADOSS$ponto,DADOSS$camara,sep='') DADOSS$resp <- rnorm(nrow(DADOSS)) model <- lmer(resp~armazen/tempo+(1|camaraM)+(1|ponto),DADOSS) summary(model) anova(model) ### Eu tenho interesse apenas em este contraste (e nos outros que envolvem a interação ): summary(glht(model,c('armazenVIDRO:tempo-armazenPP:tempo = 0'))) ## PQ não funciona????????????(Apesar de ter pouco cabelo, já arranquei alguns para tentar descobrir....) # OBS: tempo é númerico mesmo.... ## Alguma dica? ## Tem formas mais inteligente de comparar interação com multcomp? Att

2 3

Somar uma coluna de um data frame a cada mudança de valor de uma outra coluna
by Carlos Mendonça 29 Jun '11

29 Jun '11

Como posso fazer para somar uma coluna de um data frame a cada mudança de valor de uma outra coluna e gerando um outro data frame com esse resultado. Segue exemplo. Data frame de entrada codigo valor 1.01 35867,14 1.01 5867,05 1.01 62100,14 1.01 35000,99 1.02 401,87 1.02 815,7 soma o valor a cada mudança de codigo data frame de saída codigo valor 1.01 138835,32 1.02 1217,57 Obrigado, Carlos Mendonça. -- Um abraço, Mendonça

4 4

Seemi
by Alejandro C. Frery 29 Jun '11

29 Jun '11

Bom dia pessoal, Recomendo o Seemi - Simpósio de Estatística Espacial e Modelagem de Imagens, que este ano será realizado em Valparaíso, Chile. Os trabalhos podem ser enviados até 16 de setembro. http://www.uc.cl/agronomia/simposio/ Alejandro -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008

1 0

Re: [R-br] multcomp lmer contraste interação..
by Ivan Bezerra Allaman 29 Jun '11

29 Jun '11

Eder, É provável que o seu modelo esteja superparemetrizado ou algo do tipo. Procure ver mais sobre reparametrizações para ver como se adequa ao seu caso. Veja que quando eu tiro o intercepto do seu modelo, dá certo. model <- lmer(resp ~ 0 + armazen/tempo+(1|camaraM)+(1|ponto),DADOSS) summary(glht(model,c('armazenVIDRO:tempo - armazenPP:tempo = 0'))) 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

Numeração das observações em arquivo de saída
by Carlos Mendonça 28 Jun '11

28 Jun '11

Como posso fazer para retirar a numeração das observações de um data frame no arquivo de saída, tenho utilizado o comando select=c(codigo, valor), porém a numeração continua sendo gravada no arquivo conforme exemplo abaixo. codigo valor 1 1.01 35867,14 2 1.02 401,87 3 1.03 815,7 Obrigado, Carlos

2 2

Gravação de uma planilha Excel com os resultados
by Carlos Mendonça 28 Jun '11

28 Jun '11

Tenho utilizado os comandos abaixo para importar dados de um arquivo excel, porém, não consigo exportar os resultados utilizando write.xls, o programa dá uma mensagem de erro, dizendo que o número de argumentos informados não corresponde ao número de argumentos da função. Porém eu informo todos os argumentos que o help diz que precisa. Tentei, também, só informando o nome dos arquivos, sem os outros argumentos, e não funcionou. O que posso fazer? # lê arquivo de dados em excel library(RODBC) require(RODBC) xlscon <- odbcConnectExcel("D:\\Dados_13.xls") entra <- sqlFetch(xlscon,"Dados") # grava arquivo de resultados em excel write.xls(final, "D:\\Arquivo Final do Programa 2.xls", colNames = TRUE, sheet = 1, from = 1, rowNames = NA, naStrings = "") Obrigado, Carlos

2 2

Ajuda com um erro...
by Lincoln Thiago Neves 27 Jun '11

27 Jun '11

pessoal estou tentando realizar uma Análise de Correspondência, mas quando tento plotar o output dá o seguinte erro: Erro em 1:K : argumento de comprimento zero alguém pode me ajudar???? -- Att; Lincoln T. Neves MSN: lincolnthiago(a)hotmail.com E-mail: lincolnthiago.neves(a)gmail.com Web: https://sites.google.com/site/consultoriaeaulasdeestatistica/home Cel: 8883-3632

2 1