- R-br - DINF UFPR

Ajuda para plotar um gráfico
by Daniele Portela de oliveira 05 Jul '11

05 Jul '11

Boa tarde colegas, estou recorrendo ao grupo pois estou com problemas para plotar um gráfico com dois tratamentos (Sol e Sombrite), e o tempo em função do ócio. Gostaria de colocar legenda também. Abaixo tem o modelo usando e também a minha rotina do grafico: mod.10<-glm(formula=ocio ~ tempo%in%trat + I(tempo^2)%in%trat+ I(tempo^3)%in%trat+ I(tempo^4)%in%trat, family=poisson(link=log), data = dados) summary(mod.10) pred.o4<-predict(mod.10) # gráfico pred<-cbind(pred.o4) trat1<-pred[1:280,] trat2<-pred[281:560,] grafico<-predict(mod.10) ordenax<- tempo[order(tempo)] ordenat1<-trat1[order(tempo)] ordenat2<-trat2[order(tempo)] plot(ordenax,ordenat1,col="red", xlim=c(8,18),type="l", xlab="Horas do Dia", ylab="Ócio") line(ordenax,ordenat2, col="blue",xlim=c(8,18),type="l") Daniele Portela de Oliveira Mestranda Produção Animal - UEM Cel: 44 9949 3525 67 9617 0965

2 2

Como trabalhar com uma base de dados muito grande
by Carlos Mendonça 05 Jul '11

05 Jul '11

Estou com problema para utilizar a minha base de dados(.mdb). Ela tem em torno de 2 milhões de registros, e > > quando tento executar aparece um erro de falta de memória ("Erro: não é > possível alocar vetor de tamanho 1.1 Gb"). > Como posso resolver essa questão? > Coloquei um exemplo reproduzivel, inclusive com a leitura de dados. > Criei um arquivo com parte da minha base de dados. O link é > http://www.datafilehost.com/download-529c5b03.html > > <http://www.datafilehost.com/download-529c5b03.htmlO> > O código que eu estou utilizando está no link > http://www.datafilehost.com/download-a5b8a3ca.html > Obrigado, > > Carlos Mendonça. >

3 8

Data entry e data management, quais as opções?
by Pedro Emmanuel Alvarenga Americano do Brasil 05 Jul '11

05 Jul '11

Amigos de R, Desculpem o tópico um tanto fora do escopo da lista... :-( mas como imagino que muitos aqui trabalham com data management tambem, Ontem, a diretora da unidade que trabalho na FIOCRUZ veio me perguntar se acharia interessante a instituição se engajar nessa iniciativa aqui: http://www.project-redcap.org/ Eu lembro de ter lido umas mensagens antigas, provavelmente ainda no R_Stat, de alguns falando de softwares e linguagens de entrada de dados. Na minha experiencia (pra que o que eu faço), eu experimento e olho por ai, mas sempre acabo voltando para o EpiData - www.epidata.dk. Esse é bacana, é livre, é pequeno, é portatil, fácil fácil de usar, e faz o trabalho. O EpiData tem algumas desvantagens, nao permite um log de modificação de dados, ou seja, nao sei quem e quando os dados foram inseridos ou editados, nao permite entrada pela web, nao permite multiplos digitadores simultaneamente. Uma nova versão do Epidata está em teste e parece que muitas dessas coisas estão em solução. O redcap me parece ser iniciativa legal, e pelo que vi tem solução para essas questões, e melhor, exporta direto para R (entre outros). Alguem tem alguma experiencia com redcap, boa ou ruim? Alguem daria alguma dica de alguma coisa parecida? Abraço forte e que a força esteja com vocês, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.openoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 1

Re: [R-br] Levene.test
by Daniel C Bezerra 04 Jul '11

04 Jul '11

Eu gostaria de sugerir muito cuidado aos mais novos antes de usar o comando: dados<-na.omit(dados) Isto fara com que todas as linhas do banco de dados que tenham NAs sejam apagadas, o que obviamente pode levar a problemas em outras análises. O mais seguro é usar o argumento na.omit=T apenas na função de interesse. Ex.: mean(idade,na.omit=T) Para os modelos, não há necessidade de se preocupar pois por default o R usará na.omit apenas para as variáveis incluídas no modelo. O único problema é que se vc incluir duas variáveis que tenham números diferentes de observações, os modelos não estarão aninhados para comparações. Abs, D On Jul 4, 2011 8:56 AM, "Fernando Henrique Toledo" < fernandohtoledo(a)gmail.com> wrote:

2 1

Plotar média no box plot
by Edson Lira 04 Jul '11

04 Jul '11

Bom dia R-lista, como faço para mostrar no boxplot a média? Não vi nada no material que tenho. Edson Lira Estatístico Manaus-Amazonas

4 4

Coeficiente de Contingência
by Victor Eduardo 04 Jul '11

04 Jul '11

Gente, estou fazendo uma análise aqui em que tenho 2 variáveis qualitativas e preciso saber se elas estão ou não relacionadas, para isso, preciso do coeficiente de contingência, tem algum comando no R que faz esse cálculo para mim??? Abraços

4 3

Levene.test
by Victor Eduardo 04 Jul '11

04 Jul '11

Gente, alguém sabe como posso tirar os NAs para realizar o Levene.test? Tentei alguns recursos e está dando sempre erro.

2 4

Re: [R-br] Teste de Falta de Ajuste! Modelos lineares!
by Thiago De paula protásio 04 Jul '11

04 Jul '11

Obrigado Walmes e Emmanuel... Vou implementar as discas! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

OFF-TOPIC - Petição pública pela consolidação e ampliação do programa de Software Livre em todo o Governo Federal
by Daniel C Bezerra 03 Jul '11

03 Jul '11

Prezados amigos do R, Desculpem o assunto (um pouco) off-topic, mas acabo de assinar uma petição pública em prol da ampliação do programa de Software Livre em todo o Governo Federal. Como imagino que a esmagadora maioria dos membros desta lista deva ser a favor desta causa eu achei importante divulgar esta informação. Para saber mais sobre o assunto e assinar a petição basta dar uma olhada aqui: http://www.peticaopublica.com.br/?pi=SLGOV Um grande abraço a todos, D

1 0

Teste de Falta de Ajuste! Modelos lineares!
by Thiago De paula protásio 03 Jul '11

03 Jul '11

Membros da R-br; Novamente recorro a ajuda e conhecimento de vcs. Gostaria de saber se existe alguma função implementada no R para o teste de falta de ajuste de modelos lineares. Procurei nas "Rídiculas" e na lista e não encontrei mensagens a esse respeito. Muito obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

3 2

Carlos Andrade enviou-lhe o seguinte abaixo-assinado.
by Petição Pública 03 Jul '11

03 Jul '11

1 0

Quartis ponderados
by Carlos Mendonça 02 Jul '11

02 Jul '11

Como posso fazer para achar os quartis e a diferença interquartílica, por uma determinada variável e ponderados por uma outra. Por exemplo, nos dados a seguir, eu queria os quartis da variável valor, por código e ponderados pela variável quantidade. Data frame de entrada codigo valor quantidade 1.01 35867,14 125 1.01 5867,05 1045 1.01 62100,14 91 1.01 35000,99 21 1.02 401,87 2650 1.02 815,7 10345 Não sei se ajuda, mas no SAS eu utilizava o seguinte código PROC UNIVARIATE NOPRINT DATA = arq1 ; BY codigo ; VAR valor ; WEIGHT quantidade ; OUTPUT OUT = arq2 Q3=Q3 Q1=Q1 ; Obrigado, Carlos Mendonça.

7 23

Re: [R-br] Testes de Comparação Múltipla! Análise Multivariada!
by Cesar Rabak 01 Jul '11

01 Jul '11

Em 1/7/2011 18:55, Benilton Carvalho escreveu: > 2011/7/1 Cesar Rabak<cesar.rabak(a)gmail.com>: >> Em 1/7/2011 17:01, Benilton Carvalho escreveu: >>> >>> 2011/7/1 Cesar Rabak<cesar.rabak(a)gmail.com>: >>>> >>>> Em 30/6/2011 20:51, Benilton Carvalho escreveu: >>>>> >>>>> Para reflexão: se w é uma obs Bivariada com componentes (x1, y1), que >>>>> critérios usar p determinar se ela é maior / menor / igual q z = (x2, >>>>> y2)? [snipped] > Suponha que W represente o perfil de expressao genica (digamos, em > apenas 200 marcadores simultaneamente; os casos com que trabalho > envolvem alguns milhoes de marcadores simultaneamente com milhares de > pacientes) em pacientes com cancer... e que Z seja o analogo em > pessoas saudaveis. > > Se voce quiser avaliar a hipotese de que os pacientes apresentam um > perfil genico (no subconjunto acima) de maior expressao que pessoas > saudaveis, como voce definiria tal ordenacao? > > Atualmente, uma 'solucao' eh voce efetuar testes marginais e corrigir > por testes multiplos... Mas, se "soubessemos ordenar nesse dominio", > talvez pudessemos ter testes mais poderosos. > > Enfim... so' umas ideias pra quem gosta de refletir. > Benilton, Você passou a conversa de obs Bivariada para dois vetores de expressão gênica, que salvo melhor informação, são melhor pensados como seqüências de um dicionário formadas por um alfabeto de quatro símbolos. Nesse caso a melhor abordagem seria usar uma ordenação lexicográfica considerando esses símbolos com alguma "ordem" arbitrária, como por exemplo as próprias letras que sóem representar genes. Obviamente, essa ordenação não tenhuma relação com a expressão ou não, mas mantendo uma analogia com o dicionário, pode-se imaginar uma transformação entre essas seqüências e o que o nosso conhecimento já é capaz de reconhecer como "marcadores", por exemplo. []s -- Cesar Rabak GNU/Linux User 52247. Get counted: http://counter.li.org/

1 0

Como inserir frequências esperadas em um teste de Chi-quadrado
by Alberto Senra 01 Jul '11

01 Jul '11

Olá pessoal! Tenho que fazer um teste simples de chi-quadrado. Eu tenho as frequências esperadas mas não consigo inserir estas frequências para poder calcular o teste. Alguém pode me ajudar? Obrigado Alberto Senra Gonçalves, Ph.D Student Post Graduate Course in Animal Biodiversity Laboratory of Carcinology Biological Department, Universidade Federal de Santa Maria Avenida Roraima, number 1000, Cidade Universitária, Bairro Camobi, Santa Maria - RS, CEP: 97.105-900 Email: alberto.senra(a)gmail.com Fone: [0055] (55) 3212-5402 (Home fone) [0055] (55) 9966-7770 (celular fone) [cod. Brazil] (cod. city/Santa Maria)

4 3

Re: [R-br] Testes de Comparação Múltipla! Análise Multivariada!
by Thiago De paula protásio 01 Jul '11

01 Jul '11

Walmes; Obrigado pela resposta!Mas será que eu não posso utilizar o teste de Bonferroni para a comparação múltipla duas a duas? Vc não tem algum exemplo que possa ser reproduzido? Utilizando a comparação de modelos aninhados como vc sugeriu? Novamente agradeço-lhe! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

4 4

Testes de Comparação Múltipla! Análise Multivariada!
by Thiago De paula protásio 01 Jul '11

01 Jul '11

Pessoal; Recorro novamente a ajuda de vcs!Gostaria de saber se existe alguma função implementada no R para realizar a comparação múltipla multivariada....?? Muito obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

3 2

Re: [R-br] Testes de Comparação Múltipla! Análise Multivariada!
by Thiago De paula protásio 01 Jul '11

01 Jul '11

Obrigado Walmes pelo exemplo.... Quanto as comparações múltiplas....eu tb não vi em artigos, mas vi na tese do Prof. Paulo Trugilho (UFLA) defendida na UFV em 1995. Ele trabalhou com várias técnicas multivariadas e fez o teste de Roy para comparações múltiplas.... De fato, essa temática ainda é pouco explorada na estatística....Vou continuar pesquisando e qualquer novidade posto na lista......quem sabe poderemos colocar um exemplo nas "ridículas". Abraço! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Modelo misto
by Greice Laureano 30 Jun '11

30 Jun '11

Bom dia, Tenho que realizar uma analise conjunta de experimentos Fatorias e estou com problemas na hora de atribuir a parte aleatória.... Tenho 3 fatores e 5 experimentos. Preciso atribuir a parte aleatoria assim: Experimento + Sujeito(Experimento) + Experimento*Fatores Olhei o help do lme mas não consegui entender como devo colocar a fórmula fit <- lme(Desfecho ~ Fator1 + Fator2 + Fator3 + Fator1*Fator2*Fator3,random= "????", data=dados) Consegui fazer com lmer, mas eu preciso dos resíduos condicionais e marginais do modelo... Obrigado, Att Greice.

3 8

ARFIMA
by Ana Paula Sousa 30 Jun '11

30 Jun '11

Pessoal, Estou trabalhando com a previsão automática das ordens p,d,q do modelo ARFIMA através do código: arfima(x, drange = c(0, 0.5), estim = c("mle","ls"), auto.arima)## Ele calculou os parâmetros mas ocorreram os erros: *** Warning during (fracdf) fit: C fracdf() optimization failure *** Warning during (fdcov) fit: unable to compute correlation matrix; maybe change 'h' Alguém poderia ajudar com os códigos do R na estimativa desses parâmetros? Muito obrigada!

1 0

Integração R + Java
by Leandro Lima 30 Jun '11

30 Jun '11

Oi pessoal! Desculpa se alguém já postou algo na lista. Como não encontrei nada a respeito (na lista), vou perguntar: o que vcs recomendam para integrar R e Java? Obrigado! -- Leandro Lima http://www.ime.usp.br/~llima http://estudarcomputacao.com

2 2

Problema download.file
by Leandro Marino 30 Jun '11

30 Jun '11

Caros, Estou usando a função download.file para salvar uns arquivos que estou gerando no sidra. No entanto, para cada linha do arquivo está sendo criada uma linha em branco. Alguém faz ideia de como fazer o download vir correto? download.file('http://www.sidra.ibge.gov.br/download/2000BRMunicPopTot_b1.cs v', 'C:\\2000BRMunicPopTot_b1.csv',quiet = FALSE,mode="w",cacheOK=FALSE) Grato! Atenciosamente, Leandro Lins Marino Centro de Avaliação Fundação CESGRANRIO Rua Santa Alexandrina, 1011 - 2º andar Rio de Janeiro, RJ - CEP: 20261-903 ( (21) 2103-9600 R.:236 * <mailto:leandro@cesgranrio.org.br> leandro(a)cesgranrio.org.br P Antes de imprimir pense em sua responsabilidade e compromisso com o MEIO AMBIENTE Esta mensagem, incluindo seus anexos, pode conter informacoes privilegiadas e/ou de carater confidencial, nao podendo ser retransmitida sem autorizacao do remetente. Se voce nao e o destinatario ou pessoa autorizada a recebe-la, informamos que o seu uso, divulgacao, copia ou arquivamento sao proibidos. Portanto, se você recebeu esta mensagem por engano, por favor, nos informe respondendo imediatamente a este e-mail e em seguida apague-a.

2 2

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

29 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

gerar diversos subconjuntos
by Vinicius Brito Rocha 27 Jun '11

27 Jun '11

Senhores, estou com um problema muito peculiar. Preciso geara diversos subconjuntos, para ent'ao obter os histogrmas de cada um. Num primeiro momento estou mais preocupado com os subconjuntos do que com os histogramas. Note que tenho um data frame, com um output continuo e diversos fatores, cada um com diversos níveis. A idéia é obter todos os subconjuntos possíveis para a variável resposta, a partir dos fatores e então construir todos os histogramas possíveis. Pensei em alguma coisas com lapply, mas não obtive sucesso. Alguém tem alguma sugestão? abs -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

2 1

Re: [R-br] Dendrograma com UPGMA e Mahalanobis
by Ivan Bezerra Allaman 27 Jun '11

27 Jun '11

Bom dia Hugo! Após passar dias estudando para chegar em um gráfico semelhante ao que você estava querendo, creio consegui me aproximar uns 90%. Perceberás que o eixo x não está no topo como no da figura do artigo. Mandei um email para o criador da ggplot2 (pacote que usei para elaborar o gráfico) e ele me disse que não é possível ainda colocar o eixo x no topo do gráfico, mais que está trabalhando para que isso seja possível nas próximas versões. Creio que com o CMR que lhe enviarei será possível você manipulá-lo da maneira que lhe for mais conveniente. Para outros esclarecimentos consulte o site http://had.co.nz/ggplot2/. Segue o CMR abaixo: require(graphics) require(utils) require(ggplot2) require(NeatMap) dHN <- read.table('http://dl.dropbox.com/u/33619290/DadosHugoNeto.txt',h=T) dendHN <- hclust(dist(sqrt(dHN)),method='complete') clones=c("RB92579","RB768647","Q124","RB75126","RB925211","RB835486","RB946022", "RB912825","SP80-3280","RB855511","RB867515","SP83-5073","RB72454","SP80-1816", "RB825336","RB855156","RB855536","SP83-2847","RB855453","TUC77-42","SP70-2233", "RB931555","RB835054","RB855206","RB925345","RB8317","RB947501","RB93509", "RB813804","RB855113","RB83102","SP77-5181","RB845257","RB855046","RB845197", "RB835089","RB8495","RB825548","RB855463","RB925453","RB855595","RB945951", "RB946903","RB855035","RB92606","RB735200","RB945961","SP70-1143","RB855036", "IAC87-3396","CO62175","RB956911","TUC71-7","RB9350","RB966928","RB912850", "RB845210","L60-14","SP81-3250","RB928064","RB721012","SP79-1011","SP801842", "IAC52-326","RB935915","RB806043") ggplot.default() + opts(panel.background = theme_blank(), panel.grid.major = theme_blank(), panel.grid.minor = theme_blank(), legend.position ='none', axis.title.x = theme_blank(), axis.title.y = theme_blank(), axis.text.y = theme_text(angle=0,size=8))+ scale_x_continuous(breaks=seq(0:11),labels=seq(0,0.5,by=0.045))+ scale_y_continuous(breaks=seq(1:66),labels=clones) + geom_vline(xintercept=2.2,colour='red',linetype=2)+ geom_text(aes(x=2.6,y=-1,label='0.06',colour='red',size=5)) + geom_text(aes(x=3.1,y=67,label='Magic method',colour='red',size=5))+ geom_text(aes(x=5.3,y=30,label='88',colour='green',size=5))+ draw.dendrogram(dendHN, dendro.dir="right") 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

4 3

Re: [R-br] MANOVA! Comparações múltiplas multivariadas!
by Thiago De paula protásio 27 Jun '11

27 Jun '11

Entendi Ivan; Talvez seja uma boa opção mesmo.....usar a manova com caráter exploratório e partir para outras técnicas multivariadas.....ou um "fazer um mix" entre multivariada e univariada. Mas pelo que eu andei lendo a aproximação Bonferroni pode ser utilizada para esse tipo de comparação, mas este tópico de comparação múltipla multivariada tem poucas coisas implementadas nos softwares, mas como o R evolui a cada minuto...... Por isso recorri a ajuda de vcs! Muito obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Re: [R-br] MANOVA! Comparações múltiplas multivariadas!
by Ivan Bezerra Allaman 27 Jun '11

27 Jun '11

Thiago, A abordagem de comparações múltiplas sob o enfoque multivariada salve engano é diferente da abordagem univariada. Quando a hipótese de nulidade é rejeitada, você deve buscar métodos como cluster analisys, discriminante entre outras. 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

Re: [R-br] help - Ajuda com seleção de dados de uma coluna e criação de funções.
by Rafael Fernandes 27 Jun '11

27 Jun '11

Olá pessoal, Valeu pela ajuda, eu é que tava selecionando errado, agora já deu certo. Obrigado pela ajuda de todos. Em 27 de junho de 2011 11:59, Rafael Fernandes <rafadreux(a)gmail.com>escreveu: > Olá pessoal, agora eu consegui ler a tabela, com o read.csv, e trocando os > ; por , no notepad, e selecionando com o comando x2 = x1[1:49,4]. > Também consegui fazer o summary(x2), mas não estou conseguindo fazer a soma > total desse x1, é só sum(x1). > > > 2011/6/26 Rafael Fernandes <rafadreux(a)gmail.com> > >> Olá, segue em anexo parte da tabela que estou tentando trabalhar. >> Estou lendo com o read.table(). >> >> Em 23 de junho de 2011 12:01, Felipe Emanoel Barletta Mendes < >> felipe(a)leg.ufpr.br> escreveu: >> >> Rafael, pode enviar uma parte dos dados? >>> E como você está lendo eles no R? >>> >>> >>> Rafael Fernandes >>> > Cuntinua dando o mesmo erro de mudar os números. >>> > Não sei o que está acontecendo. >>> > >>> > Em 22 de junho de 2011 10:10, Felipe Emanoel Barletta Mendes < >>> > felipe(a)leg.ufpr.br> escreveu: >>> > >>> >> Rafael se dia é uma tabela de de 4 colunas e você quer selecionar os >>> >> valores de 2 a 49 desta coluna faça desta forma: >>> >> >>> >> >>> >> PAR<-dias[2:49,4] >>> >> PAR<-as.numeric(PAR) >>> >> >>> >> Seu objeto dias é um data frame? >>> >> is.data.frame(dias) >>> >> Se não for, transforme em data frame.Fica melhor para trabalhar com >>> ele. >>> >> as.data.frame(dias) >>> >> >>> >> >>> >> >>> >> >>> >> Rafael Fernandes >>> >> > Olá, eu acho que não adiantou, pq eu fiz o seguinte comando: >>> >> > >>> >> >> PAR <- as.numeric(dias$V4[2:49]) >>> >> > E deu esse resultado: >>> >> >> PAR >>> >> > [1] 1 1 2 1 1 1 1 1 1 5 7 13 41 14 24 32 40 53 62 71 73 >>> 46 >>> >> 80 >>> >> > 64 >>> >> > 66 >>> >> > [26] 57 60 43 26 42 50 34 10 11 28 35 17 8 2 2 3 1 1 1 1 1 >>> 1 >>> >> 1 >>> >> >> >>> >> > Que não é a quarta coluna dos valores 2 a 49. >>> >> > Sendo que o objeto "dias" é uma tabela de 4 colunas e eu quero >>> >> selecionar >>> >> > a >>> >> > quarta. >>> >> > Por favor me ajudem. >>> >> > >>> >> > >>> >> > >>> >> > Em 21 de junho de 2011 12:07, Felipe Emanoel Barletta Mendes < >>> >> > felipe(a)leg.ufpr.br> escreveu: >>> >> > >>> >> >> Olá rafael. >>> >> >> >>> >> >> O que acontece é que seus dados estão como fator. >>> >> >> Você precisa transformá-lo em numeric. >>> >> >> >>> >> >> dados$coluna<-as.numeric(dados$coluna) >>> >> >> >>> >> >> >>> >> >> >>> >> >> >>> >> >> Rafael Fernandes >>> >> >> > Olá Felipe, >>> >> >> > Eu não consegui fazer com que a função >>> >> >> > >>> >> >> >>soma = sum(c)*3.6 >>> >> >> > >>> >> >> > leia os dados, porque eu acho que eles não estão em uma coluna, >>> >> estão >>> >> >> só >>> >> >> > como uma sequência de números. >>> >> >> > Dá o seguinte erro: >>> >> >> > >>> >> >> >>Error in Summary.factor(c(1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 5L, >>> >> 7L, >>> >> >> : >>> >> >> > sum not meaningful for factors >>> >> >> > >>> >> >> > Como faço para mudar? >>> >> >> > >>> >> >> > Em 17 de junho de 2011 22:55, Felipe Emanoel Barletta Mendes < >>> >> >> > felipe(a)leg.ufpr.br> escreveu: >>> >> >> > >>> >> >> >> >>> >> >> >> >>> >> >> >> Rafael, >>> >> >> >> >>> >> >> >> Talvez isto lhe ajude: >>> >> >> >> >>> >> >> >> ===== criando matriz >>> >> >> >> x<-rnorm(200) >>> >> >> >> y <- rnorm(200) >>> >> >> >> dados <- data.frame(x,y) >>> >> >> >> head(dados) >>> >> >> >> >>> >> >> >> ===== selecionando a coluna y >>> >> >> >> dados2<-as.vector(dados[,2]) >>> >> >> >> soma <- function(c) >>> >> >> >> { >>> >> >> >> resp <- sum(c)*3.6 >>> >> >> >> resp >>> >> >> >> } >>> >> >> >> soma(dados2) >>> >> >> >> >>> >> >> >> >>> >> >> >> >>> >> >> >> >>> >> >> >> >>> >> >> >> Rafael Fernandes >>> >> >> >> > Olá, gostaria de saber como selecionar uma coluna de uma >>> tabela, >>> >> e >>> >> >> >> criar >>> >> >> >> > uma >>> >> >> >> > nova tabela com esses dados. E aí eu queria fazer uma função >>> que >>> >> >> some >>> >> >> >> essa >>> >> >> >> > coluna e multiplique por 3,6. >>> >> >> >> > >>> >> >> >> > Obrigado. >>> >> >> >> > >>> >> >> >> > -- >>> >> >> >> > Rafael Dreux Miranda Fernandes >>> >> >> >> > Estudante de Graduação em Engª Agronômica >>> >> >> >> > http://opinioesrafadreux.blogspot.com/ >>> >> >> >> > _______________________________________________ >>> >> >> >> > 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. >>> >> >> >> >>> >> >> >> >>> >> >> >> Felipe E. Barletta Mendes >>> >> >> >> (41)9189-5198 >>> >> >> >> (41)3025-2150 >>> >> >> >> (41)3328-7216 >>> >> >> >> http://www.leg.ufpr.br/doku.php/pessoais:felipe >>> >> >> >> >>> >> >> >> _______________________________________________ >>> >> >> >> 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. >>> >> >> >> >>> >> >> > >>> >> >> > >>> >> >> > >>> >> >> > -- >>> >> >> > Rafael Dreux Miranda Fernandes >>> >> >> > Estudante de Graduação em Engª Agronômica >>> >> >> > http://opinioesrafadreux.blogspot.com/ >>> >> >> > >>> >> >> >>> >> >> >>> >> >> Felipe E. Barletta Mendes >>> >> >> (41)9189-5198 >>> >> >> (41)3025-2150 >>> >> >> (41)3328-7216 >>> >> >> http://www.leg.ufpr.br/doku.php/pessoais:felipe >>> >> >> >>> >> >> >>> >> > >>> >> > >>> >> > -- >>> >> > Rafael Dreux Miranda Fernandes >>> >> > Estudante de Graduação em Engª Agronômica >>> >> > http://opinioesrafadreux.blogspot.com/ >>> >> > >>> >> >>> >> >>> >> Felipe E. Barletta Mendes >>> >> (41)9189-5198 >>> >> (41)3025-2150 >>> >> (41)3328-7216 >>> >> http://www.leg.ufpr.br/doku.php/pessoais:felipe >>> >> >>> >> >>> > >>> > >>> > -- >>> > Rafael Dreux Miranda Fernandes >>> > Estudante de Graduação em Engª Agronômica >>> > http://opinioesrafadreux.blogspot.com/ >>> > >>> >>> >>> Felipe E. Barletta Mendes >>> (41)9189-5198 >>> (41)3025-2150 >>> (41)3328-7216 >>> http://www.leg.ufpr.br/doku.php/pessoais:felipe >>> >>> >> >> >> -- >> Rafael Dreux Miranda Fernandes >> Agronomic Engineering Student of UFSCar - Brasil. >> http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4270671D0 >> >> > > > -- > Rafael Dreux Miranda Fernandes > Agronomic Engineering Student of UFSCar - Brasil. > http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4270671D0 > > -- Rafael Dreux Miranda Fernandes Agronomic Engineering Student of UFSCar - Brasil. http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4270671D0

1 0

Re: [R-br] Sites para colocar arquivos para ser lido no R
by Ivan Bezerra Allaman 27 Jun '11

27 Jun '11

Aeeeeeeeeeeeeeeeeeeeeeeeeeeee Beniltovisk!!! Valeuuuuuuuuuuuuuu agora ficou certinho. Muito obrigado! Gustavo, o link é esse: http://pastebin.com/kByUsY7d. Perceba que o arquivo não fica da maneira como foi colado antes. O R lê certinho, mais fica da maneira que está aí no link, e esta configuração não está correta. 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

Re: [R-br] Sites para colocar arquivos para ser lido no R
by Ivan Bezerra Allaman 27 Jun '11

27 Jun '11

Bom dia pessoal! Desde já agradeço a atenção de todos. Daniel, não deu certo com o dropbox. Fiz o meu cadastro, o download, criei a pasta, no entanto, quando clico com o mouse no botão direito em cima da pasta ou do arquivo txt, não aparece as opções que vc falou. Já procurei em tudo qualquer barra de ferramentas e não entrei nada que me gerasse um link. Gustavo, deu quase certo. O problema é que quando submeto o arquivo, ele modifica o arquivo. Veja as figuras em anexo antes de submeter o arquivo e depois de submetido. Gostaria que o arquivo ficasse do jeito que foi colado. Se puder me ajudar, fico grato. Quero aprender isso para poder ajudar o pessoal aqui do grupo mesmo, para que os códigos sejam mais otimizados, sem precisar fazer o downloado do arquivo. Obrigado a todos! 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

3 3

MANOVA! Comparações múltiplas multivariadas!
by Thiago De paula protásio 27 Jun '11

27 Jun '11

Prezados membros da R-br; Novamente recorro a experiência e conhecimento de vcs. Estou trabalhando com a análise de variância multivariada (manova). Obtive sucesso utilizando a função manova, mas como a hipótese de efeitos iguais para os tratamentos foi rejeitada gostaria de realizar a comparação múltipla das médias. Alguém conhece funções implementadas no R ou já trabalhou com comparações múltiplas multivariadas? Muito obrigado! Att., Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

1 0

Sites para colocar arquivos para ser lido no R
by Ivan Bezerra Allaman 26 Jun '11

26 Jun '11

Bom tarde senhores! Alguém conhece algum site onde eu possa disponibilizar um arquivo que possa ser lido pelo nosso grupo? Por exemplo, tem professores que disponibilizam arquivos via suas páginas pessoais, para que possa ser usado em um código mínimo reproduzível. Algo como: dados <- read.?('http://blablablabla/arquivo.?) Como eu não tenho site, gostaria de saber se alguém conhece algum site onde eu possa fazer algo semelhante ao exemplo. Já tentei usar o google docs por meio de um tutorial de um blog, mais o link gerado pelo google docs não está sendo lido pelo R. dweb <- read.csv('http://spreadsheets.google.com/spreadsheet/pub?key=0AiaQxeatudIYdHFSLWcxNTV…',header=TRUE) Erro em file(file, "rt") : não é possível abrir a conexão Obrigado mais uma vez. 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

4 5

odfWeave: tabelas com estilo ABNT
by Jakson Alves de Aquino 26 Jun '11

26 Jun '11

Olá a Todos! Preparei um script que modifica os estilos do odfWeave, tornando as tabelas mais próximas do sugerido pela ABNT. Os interessados deverão ler a seção 9.5 do "R para Cientistas Sociais" e examinar os arquivos: novo_estilo_odfWeave.R analisar_senado_2006.R exemploWeave2.odt O RparaCS (pdf e dados) está disponível em: http://www.lepem.ufc.br/aquino.php Encontrei as sugestões da ABNT para formatação de tabelas em: http://www2.unopar.br/tccs_adm/manual_adm/abnt_elementos_de_apoio_ao_texto.… Se alguém tiver alguma sugestão de melhoria, por favor, escreva. Até+, -- Jakson Aquino

1 0

Re: [R-br] Funcionamento (implementação) da função apply (lapply) em loops
by Junior Beleti 26 Jun '11

26 Jun '11

Muito obrigado pelas dicas sobre como visualizar o fonte das funções "apply". E obrigado Benilton pela explicação que foi esclarecedora. Vou continuar minhas pesquisas sobre programação paralela no R. Abraço a todos.. -- Carlos Roberto Beleti Junior Mestrado em Ciência da Computação Departamento de Informática Universidade Estadual de Maringá

2 1

Re: [R-br] Dendrograma com UPGMA e Mahalanobis
by Ivan Bezerra Allaman 24 Jun '11

24 Jun '11

Boa tarde Hugo! Ainda não estou conseguindo ler o seu arquivo. O que pode ser trivial para vc, para nós que estamos querendo ajudar as vezes não é. Por exemplo, salvei o seu arquivo com o nome de "clipboard", e mesmo assim a leitura do arquivo não sai. É preciso que vc passe detalhadamente o CMR, desde a leitura dos dados até a execução dos comandos. Não seria necessário primeiro um read.csv, read.table, ...., etc para ler o arquivo? Percebe-se que o seu CMR ainda está incompleto, pois "clone" em nenhum momento é requisitado pelas funções. Eu tenho um pouco de experiência com a "graphics", por isso é provável que eu possa ajudá-lo. Se vc quer utilizar a distância de Mahalanobis, procure por pacotes que tenham esta metodologia, pois ser digitar ?dist verás que não tem a metodologia de Mahalanobis. Instale a library "sos" para procurar por pacotes que tenha tal metodologia. install.packages("sos") library(sos) ???"Mahalanobis distance" 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

2 1

Ajuda com pacote Chron / Horas
by Rafael Garcia Cunha 22 Jun '11

22 Jun '11

Boa noite! Estou com um problema e gostaria de saber se algum de vocês poderia me ajudar. Estou gerando um Processo de Poisson não homogêneo com valores compreendidos entre 0 e 1. A ideia é que estou pensando em um banco e o 0 representa 10:00 h e o 1 representa 16:00 h. Gero os valores de 0 até 1 [que representam os instantes de ocorrências do processo de Poisson] e aplico a função "times" do pacote "chron" para que transforme esses valores em horas. Só que as horas geradas estão compreendidas ao longo de um dia, de 00:00 h até 23:59. Alguém sabe alguma maneira de colocar algo como "from = 10:00:00 to = 16:00:00" no pacote chron? Ou alguma outra maneira que eu possa resolver esse problema? Atenciosamente, Rafael Cunha # lambda(t) = 12 - t , 0 <= t <= 1 # max(lambda) = 12 library(chron) proc.5 <- function(T = 1){ t <- 0 I <- 0 S <- 0 U <- runif(1) while(t < T){ t <- t - 1/12*log(U) V <- runif(1) if(V <= 12 - t){ I <- I + 1 S[I] <- t } } return(S[-length(S)]) } x <- proc.5() times(x) -- YNWA

3 3

n misterioso
by Cristiano Strieder 21 Jun '11

21 Jun '11

Qndo chamo a função,mesmo sem pasar o valor de "n", ela executa normalmente. Porque? E debugando observo que "n" é um vetor de 15 posições. Como "n" adquire esse valor? dsad <- Vectorize(FUN= function(y,a,lambda){ poi <- function(y,n){ w <- y*log(a*n)-lfactorial(y)-a*n exp(w) } f1 <- function(n){ dexp(n,rate=lambda) * poi(y,n) } integrate(f1,0,Inf)$value }, "y") dsad(1,0.1,0.1)

2 1

Teste de friedman
by Fernanda Balbino 21 Jun '11

21 Jun '11

Olá, Gostaria de saber como deve ser a entrada dos dados para realizar o teste de friedman, estou olhando as informações do help mas não entendi como devo inserir os dados e o grupo e blocos. Em uma tentativa que fiz deu erro 'groups' faltante. Obrigada. Fernanda.

2 1

Re: [R-br] Dendrograma com UPGMA e Mahalanobis
by Ivan Bezerra Allaman 21 Jun '11

21 Jun '11

Neto, acho que seu exemplo não é reproduzível. Veja. require(graphics) require(utils) clones=c("RB92579","RB768647","Q124","RB75126","RB925211","RB835486","RB946022", "RB912825","SP80-3280","RB855511","RB867515","SP83-5073","RB72454","SP80-1816", "RB825336","RB855156","RB855536","SP83-2847","RB855453","TUC77-42","SP70-2233", "RB931555","RB835054","RB855206","RB925345","RB8317","RB947501","RB93509", "RB813804","RB855113","RB83102","SP77-5181","RB845257","RB855046","RB845197", "RB835089","RB8495","RB825548","RB855463","RB925453","RB855595","RB945951", "RB946903","RB855035","RB92606","RB735200","RB945961","SP70-1143","RB855036", "IAC87-3396","CO62175","RB956911","TUC71-7","RB9350","RB966928","RB912850", "RB845210","L60-14","SP81-3250","RB928064","RB721012","SP79-1011","SP801842", "IAC52-326","RB935915","RB806043")soca35 <- matrix(scan(file="clipboard",sep="\t"),ncol=66) # Esse daqui pessoal é da onde eu copio e colo minha matriz de Mahalanobis do Excel # Erro em scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : scan() esperava 'a real', obteve 'require(graphics)' Não estaria faltando o arquivo "clipboard"? De qualquer maneiro, nos envie um Código mínimo reproduzível (CMR) para que possamos ajudá-lo. 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

Aplicar teste CSR de Pearson a uma matriz
by Alexandre dos Santos 21 Jun '11

21 Jun '11

Boa tarde Pessoal, Tenho uma matriz que corresponde a dados de contagem em quadrats e gostaria de saber se tem como aplicar o teste de CSR de Pearson sobre os resultados da contagem, bom o spatstat faz por exemplo, mais queria fazer sem entrar com as coordenadas e criar um window, Alguem conhece algum pacote que faz diretamente sobre a matriz de resultados, Obrigado, > quadrat.ants [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 11 12 13 NA NA 16 17 18 19 20 [2,] 22 24 26 28 0 32 34 36 38 40 [3,] 33 36 39 42 45 48 51 54 57 60 [4,] 44 0 52 56 60 64 68 72 76 80 [5,] 55 60 65 70 75 80 85 0 95 100 [6,] 66 72 78 84 90 96 102 108 114 120 [7,] 77 84 91 98 105 112 119 126 133 140 [8,] 88 96 104 112 120 128 136 144 152 160 [9,] 99 108 117 126 135 144 153 162 171 NA [10,] 110 120 130 140 150 160 170 180 NA NA Alexandre dos Santos Ingenieur forestier, Msc. INRA- Biostatistique et Processus Spatiaux (BioSP) Domaine Saint-Paul Site Agroparc 84914 - Avignon - France Tél. : +33 (0)6 87 95 16 29

1 0

Granularidade dos eixos
by Vinicius Brito Rocha 21 Jun '11

21 Jun '11

Senhores, como fazer a granularidade dos eixos aumentarem? Quere intervalos menores, por exemplo : a escala está de 50 em 50, gostaria que fosse exibido de 10 em 10. Abs. -- Vinicius Brito Rocha. Estatístico e Atuário, Mestre em Pesquisa Operacional www.aplicademic.blogspot.com http://twitter.com/viniciusbritor "Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.

5 8

Lendo formato HTML
by Leonard Assis 21 Jun '11

21 Jun '11

Bom dia Tenho que acessar uma tabela em uma página asp. tenho duas perguntas simples: 1. consigo ler estes arquivos no R sem precisar converter para, por exemplo, CSV? 2. Em caso positivo, existe como saltar as 2 primeiras linhas, pois os dados começam na linha 3. Cabeçalho está na linha 2. uma base de dados que pode ser usada como exemplo seriam as disponíveis pra download no site das loterias da Caixa (Diferença é que estes estão em .zip e com o htm incluido). Já lidei com esses dados, na época eu optei por importar no excel e salvar como CSV. Portanto, aguardo sugestões de como ler uma tabela no formato html direto no R. Segue o link da caixa como exemplo de base de dados similar. http://www1.caixa.gov.br/loterias/_arquivos/loterias/D_megase.zip

3 3