- R-br - DINF UFPR

R lê o Censo Demográfico, com o Revolution e tals?
by Richard Santos 25 Nov '12

25 Nov '12

Caros, O R lê o Censo Demográfico, ou uma base no mínimo com o mesmo tamanho, com os pacotes para grandes bases, com o Revolution, e tals? Ou ninguém teve essa experiência? Obrigado. Richard

6 10

Seminários dos Usuários de R no Brasil
by Fabio Mathias Corrêa 25 Nov '12

25 Nov '12

Nobres colegas, Em diversos países os usuários do R promovem seminários de 2 ou 3 dias para apresentarem seus trabalhos relacionados ao R. No RBras que ocorreu em São Carlos realizamos um encontro dos usuários do R bem descontraído de forma a conhecermos as pessoas da lista e demais usuários que não participam da lista. Não lembro de terem feito este encontro nos eventos posteriores. Ao conversar com o Walmes hoje, acreditamos que temos potencial para expandirmos e promovermos algo mais elaborado. Claro que este tipo de evento depende da colaboração de todos! Por isso, venho através da lista saber do interesse dos usuários sobre este tipo de evento e assim podemos começarmos a planejar. Valeu! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076

33 39

Formulário para descrição dos usuários do R no Brasil
by Fabio Mathias Corrêa 24 Nov '12

24 Nov '12

Caros colegas, Atendendo a sugestão de um dos contribuintes da lista, cujo o nome eu esqueci, criei um formulário simples com perguntas objetivas para que os usuários do R no Brasil possam preencher e assim identificarmos o perfil dos usuários do R no Brasil. Caso alguém da lista conheça algum usuário do R que não participa da lista, pode enviar o link de acesso ao formulário para preenchimento. Não é obrigatório a participação no grupo R_br para preencher o formulário. https://docs.google.com/spreadsheet/viewform?formkey=dFpZSEh1QWhuNWpsNVlHam… Agradeço a contribuição de todos! Muito obrigado pela sugestão! Valeu! Fábio Mathias Corrêa Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e da Terra - DCET Campus Soane Nazaré de Andrade, km 16 Rodovia Ilhéus-Itabuna CEP 45662-900. Ilhéus-Bahia Tel.: 73-3680-5076

3 3

distancia euclidiana
by Natalia Martins 24 Nov '12

24 Nov '12

Prezados, precisei construir uma matriz de ponderacao utilizando a distancia euclidiana, entrei com a lat e a long em UTM. Por favor, qual seria a unidade da distancia que eu obtive? Att -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

2 2

Seminários dos Usuários de R no Brasil
by Ary Casado 23 Nov '12

23 Nov '12

Sugestões para organizar os interessados em seminários e congressos. 1) Fazer um cadastro dos usuários da lista interessados em fazer parte de grupos de discussão; 2) Classificar os interessados por cidade ou região próxima; 3) Identificar os interessados que possam colaborar na coordenação de reuniões por cidade ou região; A partir daí acredito que a "coisa" crescerá. Arykerne Nascimento Eng. Químico - UFPE Especialista em Informática - UFPE Mestrando em Ciências Nucleares - UFPE

1 0

Criar data.frame utilizando a PNAD
by Alexandre Loures 23 Nov '12

23 Nov '12

Boa noite pessoal! Sou iniciante no R e estou tentando criar um data.frame a partir dos dados da PNAD. No Windows obtive sucesso, porém, no Ubuntu 12.10 está ocorrendo o erro abaixo. library(IBGEPesq) data(dicPNAD2009) microdados <- file.path("/home/alexandre","Dados","PES2011.txt") dados<-le.pesquisa(dicionario=dicpes2009,pathname.in=microdados,pathname.out="dados.Rdata",codigos=c("V4803","V4704","V4805","V4707","V4729"), nomes=c("ANOEST","CONDATV","CONDOCP","HTRAB","PESO")) Erro em waitbar() : não foi possível encontrar a função "windows" Obrigado! -- *Alexandre Rodrigues Loures* Mestrando em Economia do Desenvolvimento (Bolsista Capes) Faculdade de Administração, Contabilidade e Economia - FACE Pontifícia Universidade Católica do Rio Grande do Sul - PUCRS Tel.: +55 32 9102-0690 (MG) ou +55 51 8241-6579 (RS) e-Mail: alexandre(a)loures.pro.br Skype: alexandre_loures Site: http://loures.pro.br

1 0

RES: Simulação de variáveis correlacionadas usando Cópulas
by Gabriel Bruno de Lemos 23 Nov '12

23 Nov '12

Mauro, estranho pq aqui em meu computar esta função está rodando normalmente. Segue o código desta função e também o output fornecido pela rotina que passei... > rMvdc function (n, mvdc) { dim <- mvdc@copula@dimension u <- rCopula(n, mvdc@copula) x <- u for (i in 1:dim) { qdf.expr <- asCall(paste0("q", mvdc@margins[i]), mvdc@paramMargins[[i]]) x[, i] <- eval(qdf.expr, list(x = u[, i])) } x } > require(copula) > > myCop.clayton <- archmCopula(family = "clayton", dim = 3, param = 2) > > myMvd <- mvdc(copula = myCop.clayton, margins = c("exp", "norm", + "chisq"), paramMargins = list(list(rate= 2), list(mean = 0, + sd = 1), list(df=5))) > > > myCop.t <- ellipCopula(family = "t", dim = 3, dispstr = "un", + param = c(0.8, 0.5,0.2), df = 8) > > myMvd <- mvdc(copula = myCop.t, margins = c("exp", "norm","chisq"), + paramMargins = list(list(rate= 2), list(mean = 0,sd = 1), list(df=5))) > > x.samp <- rMvdc(10000,myMvd) > cor(x.samp) [,1] [,2] [,3] [1,] 1.0000000 0.7194154 0.4536219 [2,] 0.7194154 1.0000000 0.1652259 [3,] 0.4536219 0.1652259 1.0000000 De: Mauro Sznelwar [mailto:sznelwar@uol.com.br] Enviada em: quarta-feira, 21 de novembro de 2012 23:08 Para: gblemos(a)usp.br Assunto: Re: [R-br] Simulação de variáveis correlacionadas usando Cópulas Tentei rodar o seu script e diz que este comando rMvdc não existe! Boa noite a todos! Já tem algum tempo que estou tentando aprender um pouco mais sobre cópulas... e confesso que não está sendo a coisa mais simples do mundo. Sinto falta de alguma aplicação (rotina no R) mais simples mostrando uma utilização aplicada deste tipo de função (principalmente para casos acima de 2 dimensões). Utilizo-a no meu trabalho para precificar portfólios... mas em um software caixa preta... Hoje consegui rodar um pequeno exemplo de uma das possíveis utilidades destas cópulas. Gerar variáveis aleatórias correlacionadas... com as marginais possuindo QUALQUER função de distribuição de probabilidade (não ficando restrito ao caso da normal multivariada). Parece que está funcionando legal... se alguém tiver mais alguns exemplos ou possíveis correções na rotina abaixo será bem vinda! Abs require(copula) ## Primeiro define-se a cópula a ser utilizada t neste caso. O comando dispstr=un especifica que desejo informar 3 correlações diferentes na matriz de correlações. myCop.t <- ellipCopula(family = "t", dim = 3, dispstr = "un", param = c(0.8, 0.5,0.2), df = 8) ## Definida a cópula, especifico quais marginais desejo que minha função possua. Neste caso defini exponencial, normal e qui-quadrado. myMvd <- mvdc(copula = myCop.t, margins = c("exp", "norm","chisq"), paramMargins = list(list(rate= 2), list(mean = 0,sd = 1), list(df=5))) x.samp <- rMvdc(10000,myMvd) cor(x.samp) par(mfrow=c(1,3)) hist(x.samp[,1]) hist(x.samp[,2]) hist(x.samp[,3]) Gabriel Bruno de Lemos Mestrando em Estatística e Experimentação Agronômica Esalq / Usp skype: gb_lemos msn: gb_lemos(a)hotmail.com (19) 8212-6999 .

2 1

duvida no pacote ff (read.csv.ffdf)
by Vinicius Brito Rocha 23 Nov '12

23 Nov '12

Pessoal, estou tentando ler um txt com teste=read.table.ffdf(filename, first.rows=2, next.rows=10, VERBOSE=TRUE,colClasses=c("character", "character", "character","character", "character","character","character","character", "character","POSIXct", "POSIXct","character") ) o filename possui o caminho com o nome do arquivo Estou a seguinte mensagem de erro: Erro em read.table.ffdf(filename, first.rows = 2, next.rows = 10, VERBOSE = TRUE, : only ffdf objects can be used for appending (and skipping the first.row chunk) Alguém sabe o que pode ser? -- *Vinicius Brito Rocha.* *Estatístico e Atuário (IM / UFRJ)** Mestre em Pesquisa Operacional (COPPE / UFRJ)* 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.

1 0

Re: [R-br] Movimento Valorização do profissional de Estatística, já!
by Richard Santos 23 Nov '12

23 Nov '12

E esclarecer que: um software de geoprocessamento não tem estatística, tem estatística básica e conglomerados. Um bioestatístico não é estatístico, é bioestatístico. Um econometrista não é estatístico, é econometrista. Um economista que sabe modelos matemáticos e econômicos, e que de modelos estatísticos só teve uma disciplina de regressão e Séries Temporais, só sabe regressão e Séries Temporais, pelo menos em que se formou. Pela grade de formação dele, não sabe mais, não sabe inferência, amostra, etc. Quem sabe amostra por cotas é porque não sabe amostra probabilística, e eu não tenho coragem de discutir a validade da aleatoriedade, sem fazer um doutorado antes sobre isso. Um estatístico necessariamente mexe com bases, e muitos com BIG DATA, outro profissional que trabalha com BIG DATA necessariamente não é um estatístico. Se as pessoas, a ciência, valoriza tanto a estatística, não vamos deixar acontecer um sequestro mal feito, e valorizá-la também. Abçs. Richard Em 23/11/12, Richard Santos<jamesrichardsantos(a)gmail.com> escreveu: > Caros, > > Devido aos comentários recentes e conclusões, que tal, quem concordar, > inclusive tendo em vista o ANO INTERNACIONAL DA ESTATÍSTICA, em 2013, > lançarmos o movimento "valorização do profissional de Estatística já"? > > Esse movimento contemplaria: > > - Só repetindo, de novo, acolher e valorizar os outros profissionais > que se especializam em alguma técnica estatística, mas não os > confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos > preciso. > - Parabenizar as instituições que fazem concurso para Estatístico e > vagas privadas que contemplam o estatístico, quando o trabalho for em > sua maioria fazer estatística; > - Motivar o trabalho em equipe, economista, com estatístico, com > engenheiro, com TI, com cientista social, etc; > - Motivar os estatísticos a fazerem o concurso de Analista Geral > existentes atualmente no IPEA e no IBGE, e para nos fortalecer lá > dentro, fazendo frente aos economistas; > - Definir estatístico, na regulamentação: graduação, outro > profissional que tenha pós em estatística? Matemática com pós em > estatística? Economista com pós em estatística? Manter a definição de > estatístico e técnico em estatística; > - conscientizar nos cursos de Estatística básica, nos mini-cursos, etc; > - Fazer marketing à sociedade, à sociedade civil, aos empregadores e > intelectuais brasileiros com relação à Estatística e ao profissional > de estatística; > - Só repetindo, de novo, acolher e valorizar os outros profissionais > que se especializam em alguma técnica estatística, mas não os > confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos > preciso. > > Porque se um dia no Brasil a estatística virar uma coisa que qualquer > um pode fazer ou se virar totalmente "promíscua", eu só recomendo > livros, ensino quais botões apertar, e não mudo de país não, eu vou > fazer outra coisa, aqui no Brasil mesmo, rs. > > Abraços, > Att. > Richard > Estatístico - CONRE 1 >

2 2

Concursos na Ufal
by Alejandro C. Frery 23 Nov '12

23 Nov '12

Boa tarde, Há vários editais para concursos na Ufal. Talvez nesta lista haja candidatos interessados. Ver http://www.copeve.ufal.br/index.php?opcao=concurso&idConcurso=153 Alejandro -- Alejandro C. Frery Maceió, AL - Brazil http://sites.google.com/site/acfrery/ http://www.researcherid.com/rid/A-8855-2008

1 0

R ler bases SAS
by Richard Santos 23 Nov '12

23 Nov '12

Um colega disse que "Especificamente com o formato sas7bdat tem o pacote sas7bdat. Eu já li bases do sas com esse pacote no R open source. Funciona! http://cran.r-project.org/web/packages/sas7bdat/index.html" Se ler o Censo Demográfico, que nem é um BIGDATA, eu já fico feliz. Rsrs. Lê? Rsrsrs. Abraços, Richard

1 0

Seminários dos Usuários de R no Brasil
by Ary Casado 23 Nov '12

23 Nov '12

Pelo que vejo na lista de usuários temos representantes de vários estados brasileiros e de várias instituições. Acredito que seja a hora de organizar os usuários e desenvolvedores em torno deste foco de interesse: a utilização e desenvolvimento do R tanto no meio científico como no âmbito da aplicabilidade prática. Não sei quantos interessados existem no Recife, mas acredito que já exista o interesse de iniciar encontros e seminários com o objetivo maior gerar um grupo de usuários e desenvolvedores. Arykerne Nascimento Engenheiro Químico - UFPE Especialista em Informática - UFPE Mestrando em Ciências Nucleares - UFPE

2 1

Movimento Valorização do profissional de Estatística, já!
by Richard Santos 23 Nov '12

23 Nov '12

Caros, Devido aos comentários recentes e conclusões, que tal, quem concordar, inclusive tendo em vista o ANO INTERNACIONAL DA ESTATÍSTICA, em 2013, lançarmos o movimento "valorização do profissional de Estatística já"? Esse movimento contemplaria: - Só repetindo, de novo, acolher e valorizar os outros profissionais que se especializam em alguma técnica estatística, mas não os confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos preciso. - Parabenizar as instituições que fazem concurso para Estatístico e vagas privadas que contemplam o estatístico, quando o trabalho for em sua maioria fazer estatística; - Motivar o trabalho em equipe, economista, com estatístico, com engenheiro, com TI, com cientista social, etc; - Motivar os estatísticos a fazerem o concurso de Analista Geral existentes atualmente no IPEA e no IBGE, e para nos fortalecer lá dentro, fazendo frente aos economistas; - Definir estatístico, na regulamentação: graduação, outro profissional que tenha pós em estatística? Matemática com pós em estatística? Economista com pós em estatística? Manter a definição de estatístico e técnico em estatística; - conscientizar nos cursos de Estatística básica, nos mini-cursos, etc; - Fazer marketing à sociedade, à sociedade civil, aos empregadores e intelectuais brasileiros com relação à Estatística e ao profissional de estatística; - Só repetindo, de novo, acolher e valorizar os outros profissionais que se especializam em alguma técnica estatística, mas não os confundir com estatístico, e, mesmo assim, corrigi-los quando acharmos preciso. Porque se um dia no Brasil a estatística virar uma coisa que qualquer um pode fazer ou se virar totalmente "promíscua", eu só recomendo livros, ensino quais botões apertar, e não mudo de país não, eu vou fazer outra coisa, aqui no Brasil mesmo, rs. Abraços, Att. Richard Estatístico - CONRE 1

1 0

Comando install.packages(pacote)
by Richard Santos 23 Nov '12

23 Nov '12

Pessoal, Eu dei o comando install.packages('pacote'), e o R me deu as mensagens: --- Please select a CRAN mirror for use in this session --- Aviso: unable to access index for repository http://cran.fiocruz.br/src/contrib Aviso: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/src/contrib Aviso: unable to access index for repository http://cran.fiocruz.br/bin/windows/contrib/2.15 Aviso: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.15 Mensagens de aviso perdidas: 1: In open.connection(con, "r") : não foi possível conectar a 'cran.r-project.org' na porta 80. 2: package ‘pacote’ is not available (for R version 2.15.2) What the hell? Em 23/11/12, Richard Santos<jamesrichardsantos(a)gmail.com> escreveu: > Excelentes dicas, Flávio, obrigado pela atenção. Gostei do que foi > falado sobre a linguagem, e sobre a instalação de pacotes local. Vou > testar, realmente, para aprender o R aos poucos, também. Só uma > confirmação: essas opções do Revolution e dataframe, etc, são para ler > bases grandes também, né? > Grato. > Richard > > Em 22/11/12, Flavio Barros<flaviomargarito(a)gmail.com> escreveu: >> - sempre prefiro a programação, mas a programação do R é mais extensa >> do que o SAS e o SPSS. Por exemplo, para um comando de igualdade no >> SAS ou no SPSS se digita X4=1 e no R se digita dados$x4 != 1; >> >> Os programas em R são bem compactos na verdade (em termos de sintaxe), >> mesmo se comparados a outras linguagens como python, ruby ou mesmo o >> matlab. Vide ( >> http://blog.revolutionanalytics.com/2012/11/which-programming-language-is-t… >> ) >> >> Quando a usar X4 = 1 ou dados$x4 = 1 no R é possível das duas formas, >> basta >> 'fixar' o conjunto de dados antes com o attach(dados). Particularmente >> eu >> prefiro escrever x <- 4 (padrão do google também! >> http://google-styleguide.googlecode.com/svn/trunk/google-r-style.html#assig… >> que diga-se de passagem é um usuário ativo do R.) >> >> - para grandes bases, tem-se algumas opções como Revolution, mas mesmo >> o Revolution parece que apresenta um problema que encontrei similar ao >> Minitab, que as variáveis envolvidas tem que ter o mesmo tamanho ou >> formato. Eu não queria encontrar nenhuma limitação quanto a ler >> grandes bases; >> >> Na verdade isso não é necessário. No R existem algumas estruturas de >> dados >> como o matrix, o data.frame e o list. De fato, se você quiser ler os seus >> dados em uma matriz, todos os dados tem que ser do mesmo tipo, mas em um >> data.frame não! O formato XDF da revolution também não impõe esse tipo de >> limitação. Caso as variáveis tenham tamanhos diferentes ainda é possível >> ler como um list de forma que nem o tamanho das variáveis precisa ser >> igual. Mande um e-mail para o suporte da revolution que acredito que o >> problema foi sua implementação. Eles respondem rápido. Já recebei uma >> resposta no domingo depois de 20 minutos! >> >> - para cada procedimento que se quer fazer, pode ser necessário >> instalar um pacote diferente, eu tentei pedir para o administrador do >> meu trabalho que uma pasta minha fiique livre para eu instalar >> pacotes, mas não sei se isso vai ser feito. >> >> Acho que já foi falado aqui, mas não custa lembrar: não é necessário >> privilégio de administrador para instalar pacotes. Caso você não seja >> administrador basta dar um install.packages('pacote') que ele vai criar >> um >> biblioteca local de pacotes. Esse não é um problema. >> >> No mais, é normal para aprender qualquer software ter que gastar >> tempo, mas as questões que citei acima se quiserem podem considerar >> como gosto pessoal mesmo. Mas acho importante eu citar essas questões >> para os empregadores saberem, pelo menos. Pois creio que essas >> questões não permitem ainda dizer para todos que o R é uma opção fácil >> de programa livre, porque não é tão fácil assim. Mas ele é ótimo, não >> que não seja ótimo. Apenas acho que ainda há argumentos para que uma >> empresa compre tanto o SPSS, quando achar que precisa, quanto o SAS, >> quando for necessário. >> >> Na minha opinião existem casos e casos para migrar do SAS ou SPSS para o >> R. >> Uma questão importante é o treinamento: se a empresa já tem profissionais >> bem treinados usando SAS, muito código pronto que seria difícil migrar a >> melhor solução é continuar a usar essas ferramentas. Agora no caso de um >> consultor, uma startup ou uma empresa pública que vai optar por um pacote >> estatístico acredito que o R é sim uma melhor opção. >> >> O R é muito utilizado no Google por exemplo. Se é bom para o Google dá >> para >> usar em qualquer empresa. >> >> Não advogo todo mundo usar o R, mas a linguagem tem muitas >> características >> interessantes, é bem intuitiva para quem programa em outras linguagens >> (python, C, Java por ex) e tem custo muito menor que o SAS ou SPSS. E é >> open source. >> >> No seu caso Richard que acredito que deve ter uma boa experiência em SAS, >> trocar pode não ser a melhor opção para você. Se você tem acesso a uma >> licença do SAS ou do SPSS e tem experiência com essas ferramentas não >> vejo >> porque trocar. No entanto é interessante ir aprendendo o R em paralelo. >> >> Flavio Barros >> [image: Facebook] >> <http://www.facebook.com/flavio.barros.1650?ref=tn_tnmn> [image: >> LinkedIn] <http://www.linkedin.com/profile/view?id=61839390&trk=tab_pro> >> Contact me: [image: Google Talk] flaviomargarito(a)gmail.com >> “"Sing like no one's listening, love like you've never been hurt, dance >> like nobody's watching, and live like its heaven on earth." - Mark >> Twain<http://www.quotesdaddy.com/quote/1040448/mark-twain/sing-like-no-ones-liste…> >> ” Get this email app! >> <http://www.wisestamp.com/apps/quotes?utm_source=extension&utm_medium=email&…> >> >> Create your free signature: >> <http://r1.wisestamp.com/r/landing?promo=33&dest=http%3A%2F%2Fwww.wisestamp.…> >> CLICK >> HERE!<http://r1.wisestamp.com/r/landing?promo=33&dest=http%3A%2F%2Fwww.wisestamp.…> >> >> >> >> 2012/11/22 Richard Santos <jamesrichardsantos(a)gmail.com> >> >>> ** >>> >>> >>> Caros, >>> >>> Apenas gostaria de informar que eu, como estatístico, particularmente, >>> não vou divulgar para as empresas e instituições que o R é uma opção >>> com os mesmos níveis de facilidade do que o SAS e o SPSS. Eu não tenho >>> dúvidas de que alguém que fique expert no R consiga fazer tudo e ainda >>> mais rápido do que o SAS. >>> >>> Apenas que eu pessoalmente não vou optar pelo R agora, pelo menos por >>> enquanto, ainda que eu tenha que trazer notebook para o trabalho, por >>> alguns motivos que não tive adesão por questão de gosto pessoal: >>> >>> - sempre prefiro a programação, mas a programação do R é mais extensa >>> do que o SAS e o SPSS. Por exemplo, para um comando de igualdade no >>> SAS ou no SPSS se digita X4=1 e no R se digita dados$x4 != 1; >>> >>> - para grandes bases, tem-se algumas opções como Revolution, mas mesmo >>> o Revolution parece que apresenta um problema que encontrei similar ao >>> Minitab, que as variáveis envolvidas tem que ter o mesmo tamanho ou >>> formato. Eu não queria encontrar nenhuma limitação quanto a ler >>> grandes bases; >>> >>> - para cada procedimento que se quer fazer, pode ser necessário >>> instalar um pacote diferente, eu tentei pedir para o administrador do >>> meu trabalho que uma pasta minha fiique livre para eu instalar >>> pacotes, mas não sei se isso vai ser feito. >>> >>> No mais, é normal para aprender qualquer software ter que gastar >>> tempo, mas as questões que citei acima se quiserem podem considerar >>> como gosto pessoal mesmo. Mas acho importante eu citar essas questões >>> para os empregadores saberem, pelo menos. Pois creio que essas >>> questões não permitem ainda dizer para todos que o R é uma opção fácil >>> de programa livre, porque não é tão fácil assim. Mas ele é ótimo, não >>> que não seja ótimo. Apenas acho que ainda há argumentos para que uma >>> empresa compre tanto o SPSS, quando achar que precisa, quanto o SAS, >>> quando for necessário. >>> >>> Obrigado. >>> >>> Abraços, >>> Att. >>> Richard >>> >>> >> >

5 13

pacote fora do ar
by Rodrigo Plei 23 Nov '12

23 Nov '12

Prezados, O pacote "PBSmapping" foi retirado recentemente do CRAN. Alguém teria o arquivo zip do mesmo para compartilhar? Agradeço desde já. 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 Google Scholar: http://scholar.google.com.br/citations?user=zX_EezEAAAAJ&hl=pt-BR

2 2

Função Bartlett
by tiago souza marçal 23 Nov '12

23 Nov '12

Alguém sabe como selecionar o valor p após a realização do teste de Bartlett? Att. Tiago.

2 2

Re: [R-br] Criando Rotina
by Giovane Pretto Azambuja 23 Nov '12

23 Nov '12

Boa tarde, pessoal! Consegui realizar. Obrigado pela ajuda de vocês. Grande abraço Giovane P. Azambuja Enviado do Email do Windows *De:* Luis Iván Ortiz Valencia *Enviado:* 22 de novembro de 2012 12:33 *Para:* r-br(a)listas.c3sl.ufpr.br *Assunto:* Re: [R-br] Criando Rotina Como sao varios objetos, pode criar uma lista contendo esse objetos, depois pode chamar ( com o for ) um a um os elementos da lista(objetos) e fazer o summary. Ivan Em 22 de novembro de 2012 12:27, Alisson Lucrecio <alissonluc(a)yahoo.com.br>escreveu: Use o pacote plyr e ver do for. # exemplo para fazer uma aov list.dados <- lapply(data[,-c(1)], function(y){ dm <- data[,c(1)] y <- y dfm <- cbind(dm, y) return(list( dm = dm, y = y, dfm = dfm)) }) list.av <- llply(list.dados, function(list){ with(list, aov(y ~ trar, data=dfm)) }) list.anova <- lapply(list.av, anova) Alisson Lucrécio da Costa ------------------------------ *From:* Giovane Pretto Azambuja <gibuja(a)gmail.com> *To:* "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> *Sent:* Thursday, November 22, 2012 11:47 AM *Subject:* [R-br] Criando Rotina Olá, pessoal Sou um pouco novato em R, mas vamos tentar. Tenho de fazer a análise estatística, pela função summary(), de uma pancada de objetos. Gostaria de saber se existe uma maneira em que o R faça essas análises, e as armazene em um objeto diferente para cada análise, de maneira automática. Tipo uma rotina. Creio que consigo fazer isso com o comando for, porém como são vários objetos não me veio em mente como organiza-lo. OBS: Não precisa formular a rotina inteira, apenas dê uma luz. Agradeço quem puder ajudar Giovane P. Azambuja Enviado do Email do Windows _______________________________________________ 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. _______________________________________________ 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. -- Luis Iván Ortiz Valencia Doutorando Saúde Pública - Epidemiologia, IESC, UFRJ Estatístico Msc. Spatial Analyst Msc. _______________________________________________ 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.

1 0

Path Analysis
by Malcon do Prado Costa 23 Nov '12

23 Nov '12

Caros colegas usuários do programa R, gostaria de saber se existe um pacote ou rotina para elaboração de uma análise de caminhos (Path Analysis) no R. Estou ajudando um amigo que trabalha com palmeiras do RS e ele quer fazer relações entre variáveis ambientais e biométricas, a partir da análise de caminhos. Os senhores conhecem algum pacote que trabalha com esta análise? Grato -- Malcon do Prado Costa "Aprenda como se você fosse viver para sempre. Viva como se você fosse morrer amanhã.(...) A vida é um todo indivisível.(...) É melhor que fale por nós a nossa vida, que as nossas palavras.(...) A minha fé, nas densas trevas, resplandece mais viva." (Mahatma Gandhi) -- Malcon do Prado Costa - APRENDIZ de jardineiro, mateiro e curandeiro. "Aprenda como se você fosse viver para sempre. Viva como se você fosse morrer amanhã.(...) A vida é um todo indivisível.(...) É melhor que fale por nós a nossa vida, que as nossas palavras.(...) A minha fé, nas densas trevas, resplandece mais viva." (Mahatma Gandhi)

5 4

uma ajuda
by alanarocha＠sapo.pt 23 Nov '12

23 Nov '12

bom dia a todos, oEu estou a fazer um trabalho onde tenho uma variavel condicional cm dois niveis "H" e "C" para relacionar com duas variaveis ou seja verificar se no H ha relaçao entre pthi e dose e em C ver se se ha relacao entre pthi e dose. eu fiz analise de regressão e uma anova que nao sei se será necessário. mostro meus dados e codigo. ao grafico que fiz nao cnsigo cola-lo, existe algum mil para onde possa envialo em anexo se for necessario. Doente Grupo Pthi Dose_Cinacalcet 1 54020637 H 221.60 1080 2 54020698 H 317.00 2790 3 54020681 H 412.60 930 4 54131227 H 462.00 900 5 54030338 H 474.90 930 6 54220370 H 547.00 930 7 54110496 H 547.50 930 8 54110160 H 589.00 930 9 54140258 H 709.00 930 10 54110468 H 804.90 930 11 54130380 H 842.00 930 12 54300217 H 855.40 930 13 54110542 H 911.50 930 14 54300021 H 979.30 930 15 54020732 H 1000.00 930 16 54300160 H 1054.00 840 17 54300040 H 1178.00 900 18 54131239 H 1409.00 930 19 54140376 H 1579.00 420 20 54020303 H 2794.00 1860 21 54030135 H 204.00 420 22 54140295 H 236.00 2790 23 54140011 H 253.00 930 24 54210252 H 334.00 930 25 54130307 H 358.00 390 26 54300160 H 360.80 930 27 54210085 H 426.00 930 28 54030451 H 494.50 1860 29 54020681 H 535.90 900 30 54030344 H 594.60 1860 31 54110542 H 615.70 900 32 54110255 H 659.00 1860 33 54110468 H 666.10 900 34 54020637 H 705.70 480 35 54300040 H 765.20 930 36 54300217 H 875.80 930 37 54110496 H 953.20 900 38 54131239 H 982.00 930 39 54300021 H 1275.00 930 40 54110160 H 1297.00 900 41 54110070 H 1312.00 1440 42 54130380 H 1330.00 930 43 54020732 H 1388.00 900 44 54020303 H 2060.00 2790 45 54110386 H 2432.00 900 46 54210217 H 63.50 1860 47 54110255 H 86.19 1710 48 54210104 H 236.00 2040 49 54210333 H 249.00 1860 50 54300040 H 256.00 390 51 54030135 H 265.90 60 52 54140295 H 322.70 2790 53 54030451 H 354.60 1860 54 54110542 H 436.00 930 55 54060515 H 458.80 510 56 54110468 H 460.40 930 57 54210252 H 480.10 930 58 54300160 H 480.90 900 59 54210085 H 627.80 930 60 54110496 H 700.50 930 61 54110160 H 700.90 930 62 54020732 H 739.90 870 63 54030344 H 760.70 1860 64 54020681 H 808.60 870 65 54300021 H 830.10 900 66 54210233 H 883.00 930 67 54140326 H 1033.00 390 68 54300217 H 1074.00 900 69 54210181 H 1138.00 2820 70 54130380 H 1145.00 900 71 54110070 H 1174.00 1350 72 54020637 H 1308.00 870 73 54220660 H 1401.00 1860 74 54300248 H 1462.00 1860 75 54020303 H 1536.00 2700 76 54131031 C 76.50 780 77 54210217 C 159.30 390 78 54210252 C 161.10 780 79 54210181 C 268.00 2070 80 54300021 C 280.90 390 81 54030451 C 364.10 720 82 54300160 C 574.30 390 83 54110468 C 603.90 390 84 54030135 C 651.80 390 85 54300255 C 727.10 360 86 54110542 C 738.60 390 87 54110255 C 745.00 390 88 54210233 C 764.50 1170 89 54020681 C 771.20 420 90 54300217 C 772.30 840 91 54140011 C 788.70 540 92 54030344 C 807.50 780 93 54210085 C 813.60 900 94 54110160 C 856.30 780 95 54300268 C 865.10 840 96 54110070 C 882.70 810 97 54140326 C 883.90 540 98 54140295 C 916.30 600 99 54110496 C 923.70 780 100 54210333 C 984.70 780 101 54300068 C 1044.00 420 102 54210104 C 1084.00 810 103 54020637 C 1287.00 2970 104 54020732 C 1332.00 390 105 54130380 C 1334.00 600 106 54210180 C 1401.00 420 107 54300248 C 1452.00 390 108 54110386 C 2664.00 780 109 54020303 C 2980.00 1620 110 54300068 C 176.70 390 111 54210217 C 217.00 390 112 54210104 C 221.60 1080 113 54110160 C 369.70 720 114 54300160 C 373.00 420 115 54030451 C 484.40 840 116 54300156 C 512.10 390 117 54030135 C 574.50 390 118 54210252 C 613.00 780 119 54210181 C 613.30 1080 120 54130380 C 626.30 840 121 54060515 C 706.10 840 122 54140326 C 777.30 720 123 54030344 C 793.20 780 124 54110468 C 819.00 360 125 54300255 C 822.60 390 126 54110255 C 855.20 390 127 54140011 C 861.50 960 128 54300021 C 948.70 960 129 54110496 C 984.50 660 130 54300268 C 1069.00 780 131 54140295 C 1130.00 720 132 54110070 C 1177.00 780 133 54210180 C 1250.00 1170 134 54300217 C 1318.00 780 135 54210333 C 1382.00 870 136 54300248 C 1404.00 570 137 54110542 C 1446.00 360 138 54210085 C 2230.00 1170 139 54110386 C 2577.00 720 140 54020303 C 2924.00 1440 141 54210180 C 116.80 1080 142 54210181 C 229.70 960 143 54210217 C 267.90 270 144 54210252 C 285.80 690 145 54300068 C 311.80 390 146 54300156 C 470.20 420 147 54110255 C 485.90 390 148 54030451 C 526.40 780 149 54110374 C 628.10 720 150 54110070 C 687.10 780 151 54210333 C 702.70 1170 152 54030344 C 722.20 840 153 54140326 C 765.30 840 154 54130380 C 811.80 780 155 54300217 C 845.20 720 156 54030135 C 853.80 420 157 54300255 C 949.30 570 158 54210104 C 978.10 780 159 54140295 C 1076.00 780 160 54300268 C 1240.00 780 161 54210085 C 1391.00 1260 162 54300248 C 1738.00 780 grupo<-dados$Grupo pthi<-dados$Pthi dose<-dados$Dose_Cinacalcet grupo pthi dose #'pthi ~ dose+grupo' indica: modele pthi como função estatística das variáveis dose e grupo (efeito aditivo dos modelos lineares); #a variável grupo entrou alterando o intercepto da regressão entre (pthi) e (dose). class(grupo) [1] "factor" class(pthi) [1] "numeric" class(dose) [1] "integer" modelo <- lm(pthi~dose+grupo) modelo Call: lm(formula = pthi ~ dose + grupo) Coefficients: (Intercept) dose grupoH 760.1821 0.1729 -159.5969 summary(modelo3) Call: lm(formula = pthi ~ dose + grupo) Residuals: Min 1Q Median 3Q Max -858.70 -322.92 -64.76 219.93 1939.70 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 760.18206 84.59910 8.986 7.13e-16 *** dose 0.17291 0.08286 2.087 0.0385 * grupoH -159.59687 92.60311 -1.723 0.0868 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 541.9 on 159 degrees of freedom Multiple R-squared: 0.0325, Adjusted R-squared: 0.02033 F-statistic: 2.671 on 2 and 159 DF, p-value: 0.07231 anova(modelo3) Analysis of Variance Table Response: pthi Df Sum Sq Mean Sq F value Pr(>F) dose 1 696390 696390 2.3711 0.12559 grupo 1 872380 872380 2.9703 0.08675 . Residuals 159 46698728 293703 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 x11() par(col="green") plot(pthi~dose+grupo,x1lab="dose",x2lab="grupo",ylab="pthi", main="grafico de dspersão") #Adicionar a reta da regressão no gráfico: abline(modelo3) #recta da regressão text("y=760.1821+0.1729dose-159.5969grupo, p=0.07231") ue preciso obter um grafico dedispersão onde se veja uma curva do H e outra do C mesmo não havend relação a minha onclusao O modelo entre pthi~dose entre grupos Ao relacionar pthi~dose sendo o grupo H ou C a estimativa da recta dos mínimos quadrados/regressão/relação foi Pthi estimado=760.1821+0.1729*dose-159.5969*grupo a tabela ANOVA revela cada pvalue superior a 5 e 10% o que traduz a ineficiência do modelo de regressão ou seja não há relação. Estes modelos de regressão não têm qualquer significância estatística. o que posso escrever sobre a anova para quem nao sabe o que é perceber melhor? Agradeço desde ja a vossa ajuda. Ana Rocha

1 0

uma ajuda
by alanarocha＠sapo.pt 23 Nov '12

23 Nov '12

bom dia a todos, oEu estou a fazer um trabalho onde tenho uma variavel condicional cm dois niveis "H" e "C" para relacionar com duas variaveis ou seja verificar se no H ha relaçao entre pthi e dose e em C ver se se ha relacao entre pthi e dose. eu fiz analise de regressão e uma anova que nao sei se será necessário. mostro meus dados e codigo. ao grafico que fiz nao cnsigo cola-lo, existe algum mil para onde possa envialo em anexo se for necessario. Doente Grupo Pthi Dose_Cinacalcet 1 54020637 H 221.60 1080 2 54020698 H 317.00 2790 3 54020681 H 412.60 930 4 54131227 H 462.00 900 5 54030338 H 474.90 930 6 54220370 H 547.00 930 7 54110496 H 547.50 930 8 54110160 H 589.00 930 9 54140258 H 709.00 930 10 54110468 H 804.90 930 11 54130380 H 842.00 930 12 54300217 H 855.40 930 13 54110542 H 911.50 930 14 54300021 H 979.30 930 15 54020732 H 1000.00 930 16 54300160 H 1054.00 840 17 54300040 H 1178.00 900 18 54131239 H 1409.00 930 19 54140376 H 1579.00 420 20 54020303 H 2794.00 1860 21 54030135 H 204.00 420 22 54140295 H 236.00 2790 23 54140011 H 253.00 930 24 54210252 H 334.00 930 25 54130307 H 358.00 390 26 54300160 H 360.80 930 27 54210085 H 426.00 930 28 54030451 H 494.50 1860 29 54020681 H 535.90 900 30 54030344 H 594.60 1860 31 54110542 H 615.70 900 32 54110255 H 659.00 1860 33 54110468 H 666.10 900 34 54020637 H 705.70 480 35 54300040 H 765.20 930 36 54300217 H 875.80 930 37 54110496 H 953.20 900 38 54131239 H 982.00 930 39 54300021 H 1275.00 930 40 54110160 H 1297.00 900 41 54110070 H 1312.00 1440 42 54130380 H 1330.00 930 43 54020732 H 1388.00 900 44 54020303 H 2060.00 2790 45 54110386 H 2432.00 900 46 54210217 H 63.50 1860 47 54110255 H 86.19 1710 48 54210104 H 236.00 2040 49 54210333 H 249.00 1860 50 54300040 H 256.00 390 51 54030135 H 265.90 60 52 54140295 H 322.70 2790 53 54030451 H 354.60 1860 54 54110542 H 436.00 930 55 54060515 H 458.80 510 56 54110468 H 460.40 930 57 54210252 H 480.10 930 58 54300160 H 480.90 900 59 54210085 H 627.80 930 60 54110496 H 700.50 930 61 54110160 H 700.90 930 62 54020732 H 739.90 870 63 54030344 H 760.70 1860 64 54020681 H 808.60 870 65 54300021 H 830.10 900 66 54210233 H 883.00 930 67 54140326 H 1033.00 390 68 54300217 H 1074.00 900 69 54210181 H 1138.00 2820 70 54130380 H 1145.00 900 71 54110070 H 1174.00 1350 72 54020637 H 1308.00 870 73 54220660 H 1401.00 1860 74 54300248 H 1462.00 1860 75 54020303 H 1536.00 2700 76 54131031 C 76.50 780 77 54210217 C 159.30 390 78 54210252 C 161.10 780 79 54210181 C 268.00 2070 80 54300021 C 280.90 390 81 54030451 C 364.10 720 82 54300160 C 574.30 390 83 54110468 C 603.90 390 84 54030135 C 651.80 390 85 54300255 C 727.10 360 86 54110542 C 738.60 390 87 54110255 C 745.00 390 88 54210233 C 764.50 1170 89 54020681 C 771.20 420 90 54300217 C 772.30 840 91 54140011 C 788.70 540 92 54030344 C 807.50 780 93 54210085 C 813.60 900 94 54110160 C 856.30 780 95 54300268 C 865.10 840 96 54110070 C 882.70 810 97 54140326 C 883.90 540 98 54140295 C 916.30 600 99 54110496 C 923.70 780 100 54210333 C 984.70 780 101 54300068 C 1044.00 420 102 54210104 C 1084.00 810 103 54020637 C 1287.00 2970 104 54020732 C 1332.00 390 105 54130380 C 1334.00 600 106 54210180 C 1401.00 420 107 54300248 C 1452.00 390 108 54110386 C 2664.00 780 109 54020303 C 2980.00 1620 110 54300068 C 176.70 390 111 54210217 C 217.00 390 112 54210104 C 221.60 1080 113 54110160 C 369.70 720 114 54300160 C 373.00 420 115 54030451 C 484.40 840 116 54300156 C 512.10 390 117 54030135 C 574.50 390 118 54210252 C 613.00 780 119 54210181 C 613.30 1080 120 54130380 C 626.30 840 121 54060515 C 706.10 840 122 54140326 C 777.30 720 123 54030344 C 793.20 780 124 54110468 C 819.00 360 125 54300255 C 822.60 390 126 54110255 C 855.20 390 127 54140011 C 861.50 960 128 54300021 C 948.70 960 129 54110496 C 984.50 660 130 54300268 C 1069.00 780 131 54140295 C 1130.00 720 132 54110070 C 1177.00 780 133 54210180 C 1250.00 1170 134 54300217 C 1318.00 780 135 54210333 C 1382.00 870 136 54300248 C 1404.00 570 137 54110542 C 1446.00 360 138 54210085 C 2230.00 1170 139 54110386 C 2577.00 720 140 54020303 C 2924.00 1440 141 54210180 C 116.80 1080 142 54210181 C 229.70 960 143 54210217 C 267.90 270 144 54210252 C 285.80 690 145 54300068 C 311.80 390 146 54300156 C 470.20 420 147 54110255 C 485.90 390 148 54030451 C 526.40 780 149 54110374 C 628.10 720 150 54110070 C 687.10 780 151 54210333 C 702.70 1170 152 54030344 C 722.20 840 153 54140326 C 765.30 840 154 54130380 C 811.80 780 155 54300217 C 845.20 720 156 54030135 C 853.80 420 157 54300255 C 949.30 570 158 54210104 C 978.10 780 159 54140295 C 1076.00 780 160 54300268 C 1240.00 780 161 54210085 C 1391.00 1260 162 54300248 C 1738.00 780 grupo<-dados$Grupo pthi<-dados$Pthi dose<-dados$Dose_Cinacalcet grupo pthi dose #'pthi ~ dose+grupo' indica: modele pthi como função estatística das variáveis dose e grupo (efeito aditivo dos modelos lineares); #a variável grupo entrou alterando o intercepto da regressão entre (pthi) e (dose). class(grupo) [1] "factor" class(pthi) [1] "numeric" class(dose) [1] "integer" modelo <- lm(pthi~dose+grupo) modelo Call: lm(formula = pthi ~ dose + grupo) Coefficients: (Intercept) dose grupoH 760.1821 0.1729 -159.5969 summary(modelo3) Call: lm(formula = pthi ~ dose + grupo) Residuals: Min 1Q Median 3Q Max -858.70 -322.92 -64.76 219.93 1939.70 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 760.18206 84.59910 8.986 7.13e-16 *** dose 0.17291 0.08286 2.087 0.0385 * grupoH -159.59687 92.60311 -1.723 0.0868 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 541.9 on 159 degrees of freedom Multiple R-squared: 0.0325, Adjusted R-squared: 0.02033 F-statistic: 2.671 on 2 and 159 DF, p-value: 0.07231 anova(modelo3) Analysis of Variance Table Response: pthi Df Sum Sq Mean Sq F value Pr(>F) dose 1 696390 696390 2.3711 0.12559 grupo 1 872380 872380 2.9703 0.08675 . Residuals 159 46698728 293703 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 x11() par(col="green") plot(pthi~dose+grupo,x1lab="dose",x2lab="grupo",ylab="pthi", main="grafico de dspersão") #Adicionar a reta da regressão no gráfico: abline(modelo3) #recta da regressão text("y=760.1821+0.1729dose-159.5969grupo, p=0.07231") ue preciso obter um grafico dedispersão onde se veja uma curva do H e outra do C mesmo não havend relação a minha onclusao O modelo entre pthi~dose entre grupos Ao relacionar pthi~dose sendo o grupo H ou C a estimativa da recta dos mínimos quadrados/regressão/relação foi Pthi estimado=760.1821+0.1729*dose-159.5969*grupo a tabela ANOVA revela cada pvalue superior a 5 e 10% o que traduz a ineficiência do modelo de regressão ou seja não há relação. Estes modelos de regressão não têm qualquer significância estatística. o que posso escrever sobre a anova para quem nao sabe o que é perceber melhor? Agradeço desde ja a vossa ajuda. Ana Rocha

1 0

uma ajuda
by alanarocha＠sapo.pt 23 Nov '12

23 Nov '12

bom dia a todos, oEu estou a fazer um trabalho onde tenho uma variavel condicional cm dois niveis "H" e "C" para relacionar com duas variaveis ou seja verificar se no H ha relaçao entre pthi e dose e em C ver se se ha relacao entre pthi e dose. eu fiz analise de regressão e uma anova que nao sei se será necessário. mostro meus dados e codigo. ao grafico que fiz nao cnsigo cola-lo, existe algum mil para onde possa envialo em anexo se for necessario. Doente Grupo Pthi Dose_Cinacalcet 1 54020637 H 221.60 1080 2 54020698 H 317.00 2790 3 54020681 H 412.60 930 4 54131227 H 462.00 900 5 54030338 H 474.90 930 6 54220370 H 547.00 930 7 54110496 H 547.50 930 8 54110160 H 589.00 930 9 54140258 H 709.00 930 10 54110468 H 804.90 930 11 54130380 H 842.00 930 12 54300217 H 855.40 930 13 54110542 H 911.50 930 14 54300021 H 979.30 930 15 54020732 H 1000.00 930 16 54300160 H 1054.00 840 17 54300040 H 1178.00 900 18 54131239 H 1409.00 930 19 54140376 H 1579.00 420 20 54020303 H 2794.00 1860 21 54030135 H 204.00 420 22 54140295 H 236.00 2790 23 54140011 H 253.00 930 24 54210252 H 334.00 930 25 54130307 H 358.00 390 26 54300160 H 360.80 930 27 54210085 H 426.00 930 28 54030451 H 494.50 1860 29 54020681 H 535.90 900 30 54030344 H 594.60 1860 31 54110542 H 615.70 900 32 54110255 H 659.00 1860 33 54110468 H 666.10 900 34 54020637 H 705.70 480 35 54300040 H 765.20 930 36 54300217 H 875.80 930 37 54110496 H 953.20 900 38 54131239 H 982.00 930 39 54300021 H 1275.00 930 40 54110160 H 1297.00 900 41 54110070 H 1312.00 1440 42 54130380 H 1330.00 930 43 54020732 H 1388.00 900 44 54020303 H 2060.00 2790 45 54110386 H 2432.00 900 46 54210217 H 63.50 1860 47 54110255 H 86.19 1710 48 54210104 H 236.00 2040 49 54210333 H 249.00 1860 50 54300040 H 256.00 390 51 54030135 H 265.90 60 52 54140295 H 322.70 2790 53 54030451 H 354.60 1860 54 54110542 H 436.00 930 55 54060515 H 458.80 510 56 54110468 H 460.40 930 57 54210252 H 480.10 930 58 54300160 H 480.90 900 59 54210085 H 627.80 930 60 54110496 H 700.50 930 61 54110160 H 700.90 930 62 54020732 H 739.90 870 63 54030344 H 760.70 1860 64 54020681 H 808.60 870 65 54300021 H 830.10 900 66 54210233 H 883.00 930 67 54140326 H 1033.00 390 68 54300217 H 1074.00 900 69 54210181 H 1138.00 2820 70 54130380 H 1145.00 900 71 54110070 H 1174.00 1350 72 54020637 H 1308.00 870 73 54220660 H 1401.00 1860 74 54300248 H 1462.00 1860 75 54020303 H 1536.00 2700 76 54131031 C 76.50 780 77 54210217 C 159.30 390 78 54210252 C 161.10 780 79 54210181 C 268.00 2070 80 54300021 C 280.90 390 81 54030451 C 364.10 720 82 54300160 C 574.30 390 83 54110468 C 603.90 390 84 54030135 C 651.80 390 85 54300255 C 727.10 360 86 54110542 C 738.60 390 87 54110255 C 745.00 390 88 54210233 C 764.50 1170 89 54020681 C 771.20 420 90 54300217 C 772.30 840 91 54140011 C 788.70 540 92 54030344 C 807.50 780 93 54210085 C 813.60 900 94 54110160 C 856.30 780 95 54300268 C 865.10 840 96 54110070 C 882.70 810 97 54140326 C 883.90 540 98 54140295 C 916.30 600 99 54110496 C 923.70 780 100 54210333 C 984.70 780 101 54300068 C 1044.00 420 102 54210104 C 1084.00 810 103 54020637 C 1287.00 2970 104 54020732 C 1332.00 390 105 54130380 C 1334.00 600 106 54210180 C 1401.00 420 107 54300248 C 1452.00 390 108 54110386 C 2664.00 780 109 54020303 C 2980.00 1620 110 54300068 C 176.70 390 111 54210217 C 217.00 390 112 54210104 C 221.60 1080 113 54110160 C 369.70 720 114 54300160 C 373.00 420 115 54030451 C 484.40 840 116 54300156 C 512.10 390 117 54030135 C 574.50 390 118 54210252 C 613.00 780 119 54210181 C 613.30 1080 120 54130380 C 626.30 840 121 54060515 C 706.10 840 122 54140326 C 777.30 720 123 54030344 C 793.20 780 124 54110468 C 819.00 360 125 54300255 C 822.60 390 126 54110255 C 855.20 390 127 54140011 C 861.50 960 128 54300021 C 948.70 960 129 54110496 C 984.50 660 130 54300268 C 1069.00 780 131 54140295 C 1130.00 720 132 54110070 C 1177.00 780 133 54210180 C 1250.00 1170 134 54300217 C 1318.00 780 135 54210333 C 1382.00 870 136 54300248 C 1404.00 570 137 54110542 C 1446.00 360 138 54210085 C 2230.00 1170 139 54110386 C 2577.00 720 140 54020303 C 2924.00 1440 141 54210180 C 116.80 1080 142 54210181 C 229.70 960 143 54210217 C 267.90 270 144 54210252 C 285.80 690 145 54300068 C 311.80 390 146 54300156 C 470.20 420 147 54110255 C 485.90 390 148 54030451 C 526.40 780 149 54110374 C 628.10 720 150 54110070 C 687.10 780 151 54210333 C 702.70 1170 152 54030344 C 722.20 840 153 54140326 C 765.30 840 154 54130380 C 811.80 780 155 54300217 C 845.20 720 156 54030135 C 853.80 420 157 54300255 C 949.30 570 158 54210104 C 978.10 780 159 54140295 C 1076.00 780 160 54300268 C 1240.00 780 161 54210085 C 1391.00 1260 162 54300248 C 1738.00 780 grupo<-dados$Grupo pthi<-dados$Pthi dose<-dados$Dose_Cinacalcet grupo pthi dose #'pthi ~ dose+grupo' indica: modele pthi como função estatística das variáveis dose e grupo (efeito aditivo dos modelos lineares); #a variável grupo entrou alterando o intercepto da regressão entre (pthi) e (dose). class(grupo) [1] "factor" class(pthi) [1] "numeric" class(dose) [1] "integer" modelo <- lm(pthi~dose+grupo) modelo Call: lm(formula = pthi ~ dose + grupo) Coefficients: (Intercept) dose grupoH 760.1821 0.1729 -159.5969 summary(modelo3) Call: lm(formula = pthi ~ dose + grupo) Residuals: Min 1Q Median 3Q Max -858.70 -322.92 -64.76 219.93 1939.70 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 760.18206 84.59910 8.986 7.13e-16 *** dose 0.17291 0.08286 2.087 0.0385 * grupoH -159.59687 92.60311 -1.723 0.0868 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 541.9 on 159 degrees of freedom Multiple R-squared: 0.0325, Adjusted R-squared: 0.02033 F-statistic: 2.671 on 2 and 159 DF, p-value: 0.07231 anova(modelo3) Analysis of Variance Table Response: pthi Df Sum Sq Mean Sq F value Pr(>F) dose 1 696390 696390 2.3711 0.12559 grupo 1 872380 872380 2.9703 0.08675 . Residuals 159 46698728 293703 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 x11() par(col="green") plot(pthi~dose+grupo,x1lab="dose",x2lab="grupo",ylab="pthi", main="grafico de dspersão") #Adicionar a reta da regressão no gráfico: abline(modelo3) #recta da regressão text("y=760.1821+0.1729dose-159.5969grupo, p=0.07231") ue preciso obter um grafico dedispersão onde se veja uma curva do H e outra do C mesmo não havend relação a minha onclusao O modelo entre pthi~dose entre grupos Ao relacionar pthi~dose sendo o grupo H ou C a estimativa da recta dos mínimos quadrados/regressão/relação foi Pthi estimado=760.1821+0.1729*dose-159.5969*grupo a tabela ANOVA revela cada pvalue superior a 5 e 10% o que traduz a ineficiência do modelo de regressão ou seja não há relação. Estes modelos de regressão não têm qualquer significância estatística. o que posso escrever sobre a anova para quem nao sabe o que é perceber melhor? Agradeço desde ja a vossa ajuda. Ana Roca

1 0

Re: [R-br] [STAT-MATH] R como opção de software livre
by Richard Santos 23 Nov '12

23 Nov '12

Ótima dúvida, abaixo. Encaminhando para o grupo do R também. Ou seja, como parametrizar o R para propor às empresas? (Revolution, etc). Grato. Richard Em 23/11/12, Jose Fernandes<josefer75(a)yahoo.com.br> escreveu: > Meus caros, > > Aproveito para fazer uma breve apresentação. Sou Prof de estatística na > Universidade de Cabo Verde (aqui vai o mapa para quem nunca ouviu falar > neste país > http://maps.google.pt/maps?hl=pt-PT&q=cabo+verde&bav=on.2,or.r_gc.r_pw.r_qf…) > Sempre trabalhei com SAS, mas aqui em cabo verde as empresas não utilizam > SAS (por ser caro, entre outras razões...). Assim, pensei em utilizar o R > por ser open source. Mas como ainda conheço pouco de R gostaria de saber > como que faço para parametrizar o R de modo que posso propor às > empresas? Acabei de ver o Revolutions. Penso que isso poderá ser importante > para mim e para as empresas aqui em Cabo verde. > > Cptos, > > José Moniz Fernandes > PhD Statistics > > > ________________________________ > De: Leonard de Assis <assis.leonard(a)gmail.com> > Para: STAT-MATH(a)yahoogrupos.com.br > Enviadas: Quinta-feira, 22 de Novembro de 2012 23:15 > Assunto: Re: [STAT-MATH] R como opção de software livre > > > > Outra dica boa: > https://docs.google.com/document/edit?id=1esDVxyWvH8AsX-VJa-8oqWaHLs4stGlIb… > > R Coding conventions > > []s > Leonard de Assis http://about.me/ldeassis > Em 22/11/2012 20:31, Luciano Kalatalo escreveu: > > >>Muito boa esta dica dos padrões do R no Google >> >> >>Valeu >> >> >> >> >>Luciano >> >> >> >> >> >>2012/11/22 Flavio Barros <flaviomargarito(a)gmail.com> >> >> >>>- sempre prefiro a programação, mas a programação do R é mais extensa >>>do que o SAS e o SPSS. Por exemplo, para um comando de igualdade no >>>SAS ou no SPSS se digita X4=1 e no R se digita dados$x4 != 1; >>> >>> > Os programas em R são bem compactos na verdade (em termos de sintaxe), mesmo > se comparados a outras linguagens como python, ruby ou mesmo o matlab. Vide > (http://blog.revolutionanalytics.com/2012/11/which-programming-language-is-t…) > >>> >>> >>>Quando a usar X4 = 1 ou dados$x4 = 1 no R é possível das duas formas, >>> basta 'fixar' o conjunto de dados antes com o attach(dados). >>> Particularmente eu prefiro escrever x <- 4 (padrão do google >>> também! http://google-styleguide.googlecode.com/svn/trunk/google-r-style.html#assig… >>> que diga-se de passagem é um usuário ativo do R.) >>> >>> >>>- para grandes bases, tem-se algumas opções como Revolution, mas mesmo >>>o Revolution parece que apresenta um problema que encontrei similar ao >>>Minitab, que as variáveis envolvidas tem que ter o mesmo tamanho ou >>>formato. Eu não queria encontrar nenhuma limitação quanto a ler >>>grandes bases; >>> >>> >>>Na verdade isso não é necessário. No R existem algumas estruturas de dados >>> como o matrix, o data.frame e o list. De fato, se você quiser ler os seus >>> dados em uma matriz, todos os dados tem que ser do mesmo tipo, mas em um >>> data.frame não! O formato XDF da revolution também não impõe esse tipo de >>> limitação. Caso as variáveis tenham tamanhos diferentes ainda é possível >>> ler como um list de forma que nem o tamanho das variáveis precisa ser >>> igual. Mande um e-mail para o suporte da revolution que acredito que o >>> problema foi sua implementação. Eles respondem rápido. Já recebei uma >>> resposta no domingo depois de 20 minutos! >>> >>> >>>- para cada procedimento que se quer fazer, pode ser necessário >>>instalar um pacote diferente, eu tentei pedir para o administrador do >>>meu trabalho que uma pasta minha fiique livre para eu instalar >>>pacotes, mas não sei se isso vai ser feito. >>> >>> >>>Acho que já foi falado aqui, mas não custa lembrar: não é necessário >>> privilégio de administrador para instalar pacotes. Caso você não seja >>> administrador basta dar um install.packages('pacote') que ele vai criar >>> um biblioteca local de pacotes. Esse não é um problema. >>> >>> >>>No mais, é normal para aprender qualquer software ter que gastar >>>tempo, mas as questões que citei acima se quiserem podem considerar >>>como gosto pessoal mesmo. Mas acho importante eu citar essas questões >>>para os empregadores saberem, pelo menos. Pois creio que essas >>>questões não permitem ainda dizer para todos que o R é uma opção fácil >>>de programa livre, porque não é tão fácil assim. Mas ele é ótimo, não >>>que não seja ótimo. Apenas acho que ainda há argumentos para que uma >>>empresa compre tanto o SPSS, quando achar que precisa, quanto o SAS, >>>quando for necessário. >>> >>> >>>Na minha opinião existem casos e casos para migrar do SAS ou SPSS para o >>> R. Uma questão importante é o treinamento: se a empresa já tem >>> profissionais bem treinados usando SAS, muito código pronto que seria >>> difícil migrar a melhor solução é continuar a usar essas ferramentas. >>> Agora no caso de um consultor, uma startup ou uma empresa pública que vai >>> optar por um pacote estatístico acredito que o R é sim uma melhor opção. >>> >>> >>>O R é muito utilizado no Google por exemplo. Se é bom para o Google dá >>> para usar em qualquer empresa. >>> >>> >>>Não advogo todo mundo usar o R, mas a linguagem tem muitas características >>> interessantes, é bem intuitiva para quem programa em outras linguagens >>> (python, C, Java por ex) e tem custo muito menor que o SAS ou SPSS. E é >>> open source. >>> >>> >>>No seu caso Richard que acredito que deve ter uma boa experiência em SAS, >>> trocar pode não ser a melhor opção para você. Se você tem acesso a uma >>> licença do SAS ou do SPSS e tem experiência com essas ferramentas não >>> vejo porque trocar. No entanto é interessante ir aprendendo o R em >>> paralelo. >>> >>> >>>Flavio Barros >>>Contact me: flaviomargarito(a)gmail.com >>> >>>“ "Sing like no one's listening, love like you've never been hurt, dance >>> like nobody's watching, and live like its heaven on earth." - Mark Twain >>> ” Get this email app! >>>Create your free signature: CLICK HERE! >>> >>> >>> >>>2012/11/22 Richard Santos <jamesrichardsantos(a)gmail.com> >>> >>> >>>>Caros, >>>> >>>>Apenas gostaria de > informar que eu, como > estatístico, > particularmente, >>>>não vou divulgar para as > empresas e instituições > que o R é uma opção >>>>com os mesmos níveis de > facilidade do que o SAS e > o SPSS. Eu não tenho >>>>dúvidas de que alguém que > fique expert no R consiga > fazer tudo e ainda >>>>mais rápido do que o SAS. >>>> >>>>Apenas que eu pessoalmente > não vou optar pelo R > agora, pelo menos por >>>>enquanto, ainda que eu > tenha que trazer notebook > para o trabalho, por >>>>alguns motivos que não > tive adesão por questão de > gosto pessoal: >>>> >>>>- sempre prefiro a > programação, mas a > programação do R é mais > extensa >>>>do que o SAS e o SPSS. Por > exemplo, para um comando > de igualdade no >>>>SAS ou no SPSS se digita > X4=1 e no R se digita > dados$x4 != 1; >>>> >>>>- para grandes bases, > tem-se algumas opções como > Revolution, mas mesmo >>>>o Revolution parece que > apresenta um problema que > encontrei similar ao >>>>Minitab, que as variáveis > envolvidas tem que ter o > mesmo tamanho ou >>>>formato. Eu não queria > encontrar nenhuma > limitação quanto a ler >>>>grandes bases; >>>> >>>>- para cada procedimento > que se quer fazer, pode > ser necessário >>>>instalar um pacote > diferente, eu tentei pedir > para o administrador do >>>>meu trabalho que uma pasta > minha fiique livre para eu > instalar >>>>pacotes, mas não sei se > isso vai ser feito. >>>> >>>>No mais, é normal para > aprender qualquer software > ter que gastar >>>>tempo, mas as questões que > citei acima se quiserem > podem considerar >>>>como gosto pessoal mesmo. > Mas acho importante eu > citar essas questões >>>>para os empregadores > saberem, pelo menos. Pois > creio que essas >>>>questões não permitem > ainda dizer para todos que > o R é uma opção fácil >>>>de programa livre, porque > não é tão fácil assim. Mas > ele é ótimo, não >>>>que não seja ótimo. Apenas > acho que ainda há > argumentos para que uma >>>>empresa compre tanto o > SPSS, quando achar que > precisa, quanto o SAS, >>>>quando for necessário. >>>> >>>>Obrigado. >>>> >>>>Abraços, >>>>Att. >>>>Richard >>>> >>> >>> >> > >

1 0

ajuda do R
by alanarocha＠sapo.pt 23 Nov '12

23 Nov '12

Bom dia queria a vossa ajuda. tenho um trabalho onde tenho uma variavel condicional com dois niveis "H" e "c" para relacionar com duas variaveis. e fiz a analise de regressão e anova. não sei se preciso fazer anova neste caso. O meu orientador queria que eu fizesse um grafico de dispersão a mostrar so os H e outra so a mostrar os C. Como eu sou cega nao consigo saber se o grafico mostra isso e queria saber se tenho de fazer mais algum grafico. mostro aqui os meus dados e codigo. Doente Grupo Pthi Dose_Cinacalcet 1 54020637 H 221.60 1080 2 54020698 H 317.00 2790 3 54020681 H 412.60 930 4 54131227 H 462.00 900 5 54030338 H 474.90 930 6 54220370 H 547.00 930 7 54110496 H 547.50 930 8 54110160 H 589.00 930 9 54140258 H 709.00 930 10 54110468 H 804.90 930 11 54130380 H 842.00 930 12 54300217 H 855.40 930 13 54110542 H 911.50 930 14 54300021 H 979.30 930 15 54020732 H 1000.00 930 16 54300160 H 1054.00 840 17 54300040 H 1178.00 900 18 54131239 H 1409.00 930 19 54140376 H 1579.00 420 20 54020303 H 2794.00 1860 21 54030135 H 204.00 420 22 54140295 H 236.00 2790 23 54140011 H 253.00 930 24 54210252 H 334.00 930 25 54130307 H 358.00 390 26 54300160 H 360.80 930 27 54210085 H 426.00 930 28 54030451 H 494.50 1860 29 54020681 H 535.90 900 30 54030344 H 594.60 1860 31 54110542 H 615.70 900 32 54110255 H 659.00 1860 33 54110468 H 666.10 900 34 54020637 H 705.70 480 35 54300040 H 765.20 930 36 54300217 H 875.80 930 37 54110496 H 953.20 900 38 54131239 H 982.00 930 39 54300021 H 1275.00 930 40 54110160 H 1297.00 900 41 54110070 H 1312.00 1440 42 54130380 H 1330.00 930 43 54020732 H 1388.00 900 44 54020303 H 2060.00 2790 45 54110386 H 2432.00 900 46 54210217 H 63.50 1860 47 54110255 H 86.19 1710 48 54210104 H 236.00 2040 49 54210333 H 249.00 1860 50 54300040 H 256.00 390 51 54030135 H 265.90 60 52 54140295 H 322.70 2790 53 54030451 H 354.60 1860 54 54110542 H 436.00 930 55 54060515 H 458.80 510 56 54110468 H 460.40 930 57 54210252 H 480.10 930 58 54300160 H 480.90 900 59 54210085 H 627.80 930 60 54110496 H 700.50 930 61 54110160 H 700.90 930 62 54020732 H 739.90 870 63 54030344 H 760.70 1860 64 54020681 H 808.60 870 65 54300021 H 830.10 900 66 54210233 H 883.00 930 67 54140326 H 1033.00 390 68 54300217 H 1074.00 900 69 54210181 H 1138.00 2820 70 54130380 H 1145.00 900 71 54110070 H 1174.00 1350 72 54020637 H 1308.00 870 73 54220660 H 1401.00 1860 74 54300248 H 1462.00 1860 75 54020303 H 1536.00 2700 76 54131031 C 76.50 780 77 54210217 C 159.30 390 78 54210252 C 161.10 780 79 54210181 C 268.00 2070 80 54300021 C 280.90 390 81 54030451 C 364.10 720 82 54300160 C 574.30 390 83 54110468 C 603.90 390 84 54030135 C 651.80 390 85 54300255 C 727.10 360 86 54110542 C 738.60 390 87 54110255 C 745.00 390 88 54210233 C 764.50 1170 89 54020681 C 771.20 420 90 54300217 C 772.30 840 91 54140011 C 788.70 540 92 54030344 C 807.50 780 93 54210085 C 813.60 900 94 54110160 C 856.30 780 95 54300268 C 865.10 840 96 54110070 C 882.70 810 97 54140326 C 883.90 540 98 54140295 C 916.30 600 99 54110496 C 923.70 780 100 54210333 C 984.70 780 101 54300068 C 1044.00 420 102 54210104 C 1084.00 810 103 54020637 C 1287.00 2970 104 54020732 C 1332.00 390 105 54130380 C 1334.00 600 106 54210180 C 1401.00 420 107 54300248 C 1452.00 390 108 54110386 C 2664.00 780 109 54020303 C 2980.00 1620 110 54300068 C 176.70 390 111 54210217 C 217.00 390 112 54210104 C 221.60 1080 113 54110160 C 369.70 720 114 54300160 C 373.00 420 115 54030451 C 484.40 840 116 54300156 C 512.10 390 117 54030135 C 574.50 390 118 54210252 C 613.00 780 119 54210181 C 613.30 1080 120 54130380 C 626.30 840 121 54060515 C 706.10 840 122 54140326 C 777.30 720 123 54030344 C 793.20 780 124 54110468 C 819.00 360 125 54300255 C 822.60 390 126 54110255 C 855.20 390 127 54140011 C 861.50 960 128 54300021 C 948.70 960 129 54110496 C 984.50 660 130 54300268 C 1069.00 780 131 54140295 C 1130.00 720 132 54110070 C 1177.00 780 133 54210180 C 1250.00 1170 134 54300217 C 1318.00 780 135 54210333 C 1382.00 870 136 54300248 C 1404.00 570 137 54110542 C 1446.00 360 138 54210085 C 2230.00 1170 139 54110386 C 2577.00 720 140 54020303 C 2924.00 1440 141 54210180 C 116.80 1080 142 54210181 C 229.70 960 143 54210217 C 267.90 270 144 54210252 C 285.80 690 145 54300068 C 311.80 390 146 54300156 C 470.20 420 147 54110255 C 485.90 390 148 54030451 C 526.40 780 149 54110374 C 628.10 720 150 54110070 C 687.10 780 151 54210333 C 702.70 1170 152 54030344 C 722.20 840 153 54140326 C 765.30 840 154 54130380 C 811.80 780 155 54300217 C 845.20 720 156 54030135 C 853.80 420 157 54300255 C 949.30 570 158 54210104 C 978.10 780 159 54140295 C 1076.00 780 160 54300268 C 1240.00 780 161 54210085 C 1391.00 1260 162 54300248 C 1738.00 780 grupo<-dados$Grupo pthi<-dados$Pthi dose<-dados$Dose_Cinacalcet grupo pthi dose #'pthi ~ dose+grupo' indica: modele pthi como função estatística das variáveis dose e grupo (efeito aditivo dos modelos lineares); #a variável grupo entrou alterando o intercepto da regressão entre (pthi) e (dose). class(grupo) [1] "factor" class(pthi) [1] "numeric" class(dose) [1] "integer" modelo <- lm(pthi~dose+grupo) modelo Call: lm(formula = pthi ~ dose + grupo) Coefficients: (Intercept) dose grupoH 760.1821 0.1729 -159.5969 summary(modelo3) Call: lm(formula = pthi ~ dose + grupo) Residuals: Min 1Q Median 3Q Max -858.70 -322.92 -64.76 219.93 1939.70 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 760.18206 84.59910 8.986 7.13e-16 *** dose 0.17291 0.08286 2.087 0.0385 * grupoH -159.59687 92.60311 -1.723 0.0868 . --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Residual standard error: 541.9 on 159 degrees of freedom Multiple R-squared: 0.0325, Adjusted R-squared: 0.02033 F-statistic: 2.671 on 2 and 159 DF, p-value: 0.07231 anova(modelo3) Analysis of Variance Table Response: pthi Df Sum Sq Mean Sq F value Pr(>F) dose 1 696390 696390 2.3711 0.12559 grupo 1 872380 872380 2.9703 0.08675 . Residuals 159 46698728 293703 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 x11() par(col="green") plot(pthi~dose+grupo,x1lab="dose",x2lab="grupo",ylab="pthi", main="grafico de dispersao") #Adicionar a reta da regressão no gráfico: abline(modelo3) #recta da regressão text("y=760.1821+0.1729dose-159.5969grupo,p=0.07231") gráfico: n cnsigo colar aqui o grafico apara que mil posso enviar em anexo?o Ana

1 0

Parcela Subdividida com dados desbalanceados
by ROSIANA RODRIGUES ALVES 23 Nov '12

23 Nov '12

Olá Pessoal, Estou com problemas em analisar dados de um experimento com parcelas subdivididas no tempo desbalanceado. Geralmente eu faço a analise o seguindo modelo modelo <- aov(Bloco + FatorA *Tempo + Error(FatorA:Bloco)), mas quando eu faço com os dados com parcela perdida a somas de quadrados não estão saindo correntamente eles saem assim : Error: FatorA:Bloco Df Sum Sq Mean Sq F value Pr(>F) Bloco 3 0.1158 0.0386 0.278 0.8400 FatorA 3 2.1162 0.7054 5.073 0.0355 * Tempo 1 0.0124 0.0124 0.089 0.7743 Cultivar:Tempo 1 0.0152 0.0152 0.110 0.7504 Residuals 7 0.9733 0.1390 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Tempo 1 1.7152 1.7152 24.943 0.000542 *** Cultivar:Tempo 3 0.5318 0.1773 2.578 0.112073 Residuals 10 0.6877 0.0688 preciso de sugestões para fazer essa analise. Obrigada desde já.

3 2

Parcela Subdividida
by tiago souza marçal 22 Nov '12

22 Nov '12

Pessoal, gostaria de saber se em parcela subdividida eu devo realizar o teste normalidade dos resíduos (shapiro.test()) para a parcela e subparcela separadamente. Agradeço desde já pela ajuda. Att. Tiago.

2 1

Normalidade dos resíduos em parcela subdividida
by tiago souza marçal 22 Nov '12

22 Nov '12

Pessoal, gostaria de saber se em parcela subdividida eu tenho que verificar a normalidade dos resíduos (shapiro.test()) da parcela e da subparcela? Agradeço desde já por qualquer ajuda. Att. Tiago.

1 0

rect.hclust
by Mauricio Figueira 22 Nov '12

22 Nov '12

Oi, estou tentando delimitar os grupos no Cluster analysis, porém a função rect.hclust extrapola os limites dos retângulos, deixando-os abertos. O script que estou usando segue abaixo. #dissimilaridade EUCLIDIANA dis<-vegdist (dados, method = "euclidean") dis range(dis) clust<- hclust(dis, "average") plot(clust, main="Análise de Agrupamento", sub="Ligação pela média", xlab="Locais e orgãos", ylab="Dissimilaridade (Distância Euclidiana)") cor(dis, cophenetic(clust)) rect.hclust(clust, k=2, border="red")#desenhar dois retângulos x <- rect.hclust(clust, h=40)#cortar na altura de 40 Se alguém puder me ajudar, desde já agradeço. Maurício Figueira Mestrando em Engenharia Florestal (UFSM) Herbário do Departamento de Ciências Florestais (HDCF) Currículo: http://lattes.cnpq.br/1022600405818850 Tel.:(55) 9663-4920

2 2

Remover linhas do banco de dados, sujeito a algum critério da variável
by Abel Brasil Ramos da Silva 22 Nov '12

22 Nov '12

Olá a todos, Preciso excluir as linhas do banco de dados segundo algum critério da variável, por exemplo: x1=c(seq(1,10),NA,seq(1,10)) x2=c(seq(1,9),NA,seq(1,11)) x3=c(seq(1,8),NA,seq(1,12)) x4=c(seq(1,7),NA,seq(1,13)) x5=c(seq(1,6),NA,seq(1,14)) dados=data.frame(x1,x2,x3,x4,x5) >dados x1 x2 x3 x4 x5 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 5 5 5 5 5 5 6 6 6 6 6 6 7 7 7 7 7 NA 8 8 8 8 NA 1 9 9 9 NA 1 2 10 10 NA 1 2 3 11 NA 1 2 3 4 12 1 2 3 4 5 13 2 3 4 5 6 14 3 4 5 6 7 15 4 5 6 7 8 16 5 6 7 8 9 17 6 7 8 9 10 18 7 8 9 10 11 19 8 9 10 11 12 20 9 10 11 12 13 21 10 11 12 13 14 se eu utilizar > na.exclude(dados) x1 x2 x3 x4 x5 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 5 5 5 5 5 5 6 6 6 6 6 6 12 1 2 3 4 5 13 2 3 4 5 6 14 3 4 5 6 7 15 4 5 6 7 8 16 5 6 7 8 9 17 6 7 8 9 10 18 7 8 9 10 11 19 8 9 10 11 12 20 9 10 11 12 13 21 10 11 12 13 14 as linhas 7 a 11 foram excluídas, porém como eu trabalho com banco de dados de 50 variáveis, gostaria de excluir , por exemplo, as linhas cuja variáveis x4 e x5 sejam igual a NA (não somente NA, mas outro critério numérico ou caracteres), produzindo o seguinte resultado: x1 x2 x3 x4 x5 1 1 1 1 1 1 2 2 2 2 2 2 3 3 3 3 3 3 4 4 4 4 4 4 5 5 5 5 5 5 6 6 6 6 6 6 9 9 9 NA 1 2 10 10 NA 1 2 3 11 NA 1 2 3 4 12 1 2 3 4 5 13 2 3 4 5 6 14 3 4 5 6 7 15 4 5 6 7 8 16 5 6 7 8 9 17 6 7 8 9 10 18 7 8 9 10 11 19 8 9 10 11 12 20 9 10 11 12 13 21 10 11 12 13 14 Obrigado.

3 8

Ajuda gráfico
by Sérgio Henrique almeida da silva ju 22 Nov '12

22 Nov '12

Olá Amigos Estou rodando um gráfico através da função mosaicplot: exemplo: CAT=c(rep('CI',60),rep('CII',30),rep('CIII',10)) SEXO=c(rep('F',28),rep('M',32),rep('F',20),rep('M',10),rep('F',7),rep('M',3)) MOSAICO = cbind(as.data.frame(CAT),as.data.frame(SEXO)) edit(MOSAICO) tabela =table(MOSAICO$CAT,MOSAICO$SEXO) mosaicplot(tabela) Porém gostaríamos de que os grupos (CI, CII e CIII) tivessem os mesmos tamanhos (Como o gráfico anexo). O que podemos fazer na função ou no comando para que os gráficos fiquem da forma anexa? Obrigado mais uma vez! http://www.datafilehost.com/download-fe09c339.html -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 94429486/78101651 id: 123*20942

1 1

Modelo Misto + LSD de Fisher
by Juliana Passos 22 Nov '12

22 Nov '12

Olá, eu gostaria de usar o teste de comparações múltiplas LSD de Fisher após rodar um modelo misto, porém não encontrei nenhuma função para tal. Os comandos que tentei são: mod1 <- aov(resposta ~ fator + Error(sujeito/fator), data=dados) require(nlme) mod2 <- lme(resposta ~ fator, random = ~1|sujeito/fator,data=dados) Usando o mod2 consigo rodar o teste de Tukey com a seguinte função: glht(mod2,linfct=mcp(fator="Tukey")) Porém eu preciso do LSD de Fisher. Encontrei a função LSD.test, mas ela não funciona para modelo misto require(agricolae) LSD.test(mod1,"fator", group=TRUE) Alguém sabe me indicar alguma função que faça o teste LSD de Fisher considerando o mod1, mod2, ou outro modelo com mesmos resultados? Desde já agradecida, Juliana

8 9

dúvidas em outputs
by alanarocha＠sapo.pt 22 Nov '12

22 Nov '12

Olá, eu tenho um output da função summary de 4 variaveis ou seja tenho 4 valores de medias diferentes medianas diferentes desvios padroes diferentes. eu ando a procura de algo que me esplique por exemplo qual dos 4 vaores de medias é imais sgnificativo, e o mesmo relativamente, a mediana quantis desvio padrão variancia. obrigada desde já. vai aqui meu exemplo para ser mais facil. clínicas X1s09 X1s10 X1s11 X1s12 1 Ponte da Barca 70.5 73.9 75.7 81.0 2 Braga 72.5 79.9 77.3 73.4 3 Fafe 74.7 74.7 71.1 74.4 4 Maia 76.7 72.7 77.3 78.1 5 VNGaia 84.6 84.3 83.4 81.8 6 SMFeira 73.0 73.8 76.7 79.4 7 Viseu 67.4 75.7 75.6 74.5 8 Guarda 70.6 65.4 68.7 78.9 9 CovilhÆ 76.7 76.0 76.4 81.3 10 Coimbra 77.5 77.9 77.0 77.7 11 Abrantes 84.0 81.2 82.9 79.1 12 Entroncamento 74.4 73.5 74.1 72.0 13 Santar‚m 75.2 69.0 80.6 83.3 14 TVedras 70.6 74.1 81.1 78.5 15 VFXira 70.7 76.7 78.3 79.7 16 Alverca 74.3 80.9 84.4 80.2 17 Amadora 69.5 70.5 74.7 75.5 18 Lumiar 71.4 72.3 73.2 76.9 19 Restelo 71.2 69.7 72.5 74.0 20 CUF 64.2 65.5 71.6 66.5 21 SAMS 78.8 82.9 84.0 82.3 22 APDP 0.0 0.0 0.0 57.6 23 Montijo 77.3 71.4 69.8 84.3 24 Barreiro 78.6 78.6 77.8 77.6 25 Almada 76.5 73.1 74.8 78.2 26 Seixal 78.6 75.4 80.9 80.0 27 Set£bal 76.3 79.8 78.2 74.8 28 Portalegre 77.0 80.0 80.1 77.7 29 vora 81.2 79.6 80.9 81.8 30 Grƒndola 75.9 78.4 80.2 78.8 31 Beja 73.9 76.4 75.0 67.8 32 Tavira 73.0 66.5 66.9 75.1 33 Faro 70.6 78.5 79.1 77.6 34 PortimÆo 71.9 80.2 79.0 80.6 dim(dados) [1] 34 5 #estatísticas descritivas sapply(dados[,2:4],mean) média: X1s09 X1s10 X1s11 X1s121s09 72.332 73.19 1 74.685 74.685 sapply(dados[,2:5],median) medianas: X1s09 X1s10 X1s11 X1s12 74.35 75.55 77.15 78.15 sapply(dados[,2:5],min) mínimo: X1s09 X1s10 X1s11 X1s12 0.0 0.0 0.0 57.6 sapply(dados[,2:5],max) máximo: X1s09 X1s10 X1s11 X1s12 84.6 84.3 84.4 84.3 sapply(dados[,2:5],sd) desvio padrão: X1s09 X1s10 X1s11 X1s12 13.501958 13.796999 13.890047 5.241185 sapply(dados[,2:5],var) variância: X1s09 X1s10 X1s11 X1s12 182.30286 190.35719 192.93341 27.47002 sapply(dados[,2:5],quantile) quantis: X1s09 X1s10 X1s11 X1s12 0% 0.000 0.00 0.000 57.600 25% 70.825 72.40 74.250 74.875 50% 74.350 75.55 77.150 78.150 75% 76.925 79.35 80.175 80.150 100% 84.600 84.30 84.400 84.300 diff(range(dados$X1s09, na.rm=T)) # amplitude amplitude total X1s09 X1s10 X1s11 X1s12 84.6 diff(range(dados$X1s10, na.rm=T)) # amplitude [1] 84.3 diff(range(dados$X1s11, na.rm=T)) # amplitude [1] 84.4 diff(range(dados$X1s12, na.rm=T)) # amplitude [1] 26.7 Ana Rocha

1 0

Criando Rotina
by Giovane Pretto Azambuja 22 Nov '12

22 Nov '12

Olá, pessoal Sou um pouco novato em R, mas vamos tentar. Tenho de fazer a análise estatística, pela função summary(), de uma pancada de objetos. Gostaria de saber se existe uma maneira em que o R faça essas análises, e as armazene em um objeto diferente para cada análise, de maneira automática. Tipo uma rotina. Creio que consigo fazer isso com o comando for, porém como são vários objetos não me veio em mente como organiza-lo. OBS: Não precisa formular a rotina inteira, apenas dê uma luz. Agradeço quem puder ajudar Giovane P. Azambuja Enviado do Email do Windows

3 2

Enc: Problemas com rgl!
by Ivan Bezerra Allaman 22 Nov '12

22 Nov '12

Boa tarde Daniel! Tenho a seguinte mensagem como camando: Lendo listas de pacotes... Pronto Construindo árvore de dependências Lendo informação de estado... Pronto E: Você deve colocar algumas URIs 'source' em seu sources.list O que isso quer dizer eu não sei, só sei que tenho a seguinte conteúdo no meu sources.list: deb http://mint-packages.c3sl.ufpr.br/ maya main upstream import deb http://archive.ubuntu.com/ubuntu/ precise main restricted universe multiverse deb http://archive.ubuntu.com/ubuntu/ precise-updates main restricted universe multiverse deb http://security.ubuntu.com/ubuntu/ precise-security main restricted universe multiverse deb http://archive.canonical.com/ubuntu/ precise partner deb http://packages.medibuntu.org/ precise free non-free # deb http://archive.getdeb.net/ubuntu precise-getdeb apps # deb http://archive.getdeb.net/ubuntu precise-getdeb games #Repositório do R deb http://cran-r.c3sl.ufpr.br/bin/linux/ubuntu precise/ deb http://cran-r.c3sl.ufpr.br/ precise-backports main restricted universe Inté, \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature} ________________________________ De: Daniel C Bezerra <danielcbezerra(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br; Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> Enviadas: Sábado, 8 de Setembro de 2012 15:57 Assunto: Re: [R-br] Problemas com rgl! Tente rodar no terminal: sudo apt-get build-dep r-cran-rgl On Sat, Sep 8, 2012 at 3:30 PM, Ivan Bezerra Allaman <ivanalaman(a)yahoo.com.br> wrote: Boa tarde senhores! > > >Ao rodar um gráfico em 3D com a rgl percebi que estava muito lento, travando. Então descobri que não estava instalado o driver de vídeo da minha máquina. Ao instalar e tentar usar a rgl obtive erros. Então removi o pacote que estava dentro do usr blablab e tentei reinstalar o pacote sem sucesso com o seguinte erro: > > >* installing *source* package ‘rgl’ ... >** package ‘rgl’ successfully unpacked and MD5 sums checked >checking for gcc... gcc -std=gnu99 >checking whether the C compiler works... yes >checking for C compiler default output file name... a.out >checking for suffix of executables... >checking whether we are cross compiling... no >checking for suffix of object files... o >checking whether we are using the GNU C compiler... yes >checking whether gcc -std=gnu99 accepts -g... yes >checking for gcc -std=gnu99 option to accept ISO C89... none needed >checking how to run the C preprocessor... gcc -std=gnu99 -E >checking for gcc... (cached) gcc -std=gnu99 >checking whether we are using the GNU C compiler... (cached) yes >checking whether gcc -std=gnu99 accepts -g... (cached) yes >checking for gcc -std=gnu99 option to accept ISO C89... (cached) none needed >checking for libpng-config... yes >configure: using libpng-config >configure: using libpng dynamic linkage >checking for X... libraries , headers >checking GL/gl.h usability... yes >checking GL/gl.h presence... yes >checking for GL/gl.h... yes >checking GL/glu.h usability... yes >checking GL/glu.h presence... yes >checking for GL/glu.h... yes >checking for glEnd in -lGL... yes >checking for gluProject in -lGLU... yes >checking for freetype-config... yes >configure: using Freetype and FTGL >configure: creating ./config.status >config.status: creating src/Makevars > >. >. >. > >Error : .onLoad failed in loadNamespace() for 'rgl', details: > call: dyn.load(file, DLLpath = DLLpath, ...) > error: impossível carregar objeto compartilhado '/usr/local/lib/R/site-library/rgl/libs/rgl.so': > libGL.so.1: cannot open shared object file: No such file or directory >Erro: loading failed >Execução interrompida >ERROR: loading failed >* removing ‘/usr/local/lib/R/site-library/rgl’ > >The downloaded source packages are in > ‘/tmp/RtmpS6wZHV/downloaded_packages’ >Mensagens de aviso perdidas: >In install.packages("rgl") : > installation of package ‘rgl’ had non-zero exit status > > >Já pesquisei em vários lugares, instalar várias coisas e o erro persiste. Grato se alguém puder ajudar. Utilizo o mint 13 maya ubuntu. > > >Abraços! > > > >\begin{signatu ><<>>= >Prof. Dr. Ivan Bezerra Allaman >Universidade Estadual de Santa Cruz >Departamento de Ciências Exatas e Tecnológicas >Ilhéus/BA - Brasil >Fone: +55 73 3680-5596 >E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com >@ >\end{signature} >_______________________________________________ >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. >

3 3

Regressão linear
by Narede Golaco 22 Nov '12

22 Nov '12

Amigos gostaria de saber como se analisa uma regressão linear apartir do test T ? Agradeço a ajuda e qualquer exemplo, assim como a indicação de material, pois ja coloquei no google " teste T e regressão linear '' e não aparece nada de relevante que avalie a regressão com o teste T. Agradeço desde ja a ajuda.

4 6

Comparar graficamente densidades estimadas
by Matheus Monteiro 21 Nov '12

21 Nov '12

Obtive as densidades aproximadas da 't' e da 'normal', gostaria de compara-las através de graficos dessas densidades estimadas. Qual seria o melhor comando para isso? Uso o que está entre parênteses ou o outro? Obrigado. *t * > fitdistr(x[,1],"t") m s df -0.001874414 0.009122738 3.522064232 ( 0.001726835) ( 0.001862692) ( 2.234067399)* **normal* fitdistr(x[,1],"normal") mean sd 0.0001040816 0.0128201843 (0.0018314549) (0.0012950342) x são dados simulados.

2 4

duvida
by Hélio Gallo Rocha 21 Nov '12

21 Nov '12

Boa noite, Espero não estar enviando este e-mail para todos, mas sim para algum admin da lista. Não estou conseguindo configurar para receber todas as msg separadamente e não resumos. Grato pela atenção Hélio

2 1

Re: [R-br] INLA Geoestatística
by Elias T. Krainski 21 Nov '12

21 Nov '12

Oi Taynãna, DIC negativo pode ocorrer com distribuições contínuas sem problemas. O mesmo vale para log-Verossimilhança. Veja sobre DIC em http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/dicpage.shtml ### exemplo de log-verossimilhança negativa y1 <- rnorm(100, 0, 100) ### log-verossimilhança em 0, 100 sum(dnorm(y1, 0, 100, log=TRUE)) ### exemplo de log-verossimilhança positiva y2 <- rnorm(100, 0, 0.01) ### log-verossimilhança em 0, 0.01 sum(dnorm(y2, 0, 0.01, log=TRUE)) ### para distribuições discretas, o ### log-verossimilhança será sempre negativa y3 <- rpois(100, 3) ### log lik sum(dpois(y3, 3, log=TRUE)) att. Elias T. Krainski > >Ei, Elias. >Agora sim, consegui rodar, mas o DIC está dando negativo, com qualquer das estratégias utilizadas (gaussiana, laplace ou laplace simplificada)... >Percebeu isso? > >Abraço, Taynãna > > >Taynãna, fiz > inla.upgrade(testing=TRUE) >atualizou para a versao de 5 de novembro >e consegui executar todo o script SwissRainfall.R sem problemas... > > > > Att. >Elias T. Krainski > >>Oi, Elias. >> >> >>O Havard Rue respondeu dizendo p tentar o comando inla.upgrade(testing=TRUE), mas não deu certo. >>A versão que estou usando é de 5 de novembro. >> >> >>Abraço, Taynãna >> >>Oi Taynãna, >> >>qual a versão do seu INLA. >> inla.version() >>Eu consegui executar sem problemas na versão de 18 de agosto >> >> >> >>Elias T. Krainski >> >>> >>> >>>Pessoal, >>> >>>Estou tentando reproduzir um exemplo de modelagem espacial no INLA com dados de geoestatística, do artigo >>> >>>de Blangiardo et al (2012). O código em R está disponível e o exemplo 3 (SwissRainfall.R). >>> >>>Ao rodar o modelo aparece a seguinte mensagem: >>> >>>> mod = inla(formula, >>>+ data=inla.stack.data(stack, spde=spde),family="gaussian", >>>+ control.predictor=list(A=inla.stack.A(stack), compute=TRUE), >>>+ >>control.compute=list(cpo=TRUE, dic=TRUE), >>>+ keep=FALSE, verbose=TRUE) >>>Erro em .local(x, i, j, ..., value) : >>> not-yet-implemented 'Matrix[<-' method >>>Calls: inla -> inla.predictor.section -> [<- -> [<- -> .local >>> >>>Alguém tem idéia do quê possa estar acontecendo? >>> >>>Att, >>>Taynãna >>>Pesquisadora em Saúde Pública - ENSP/FIOCRUZ >

2 1

Geoestatística: pacotes para o R [iniciantes]
by Ligia Souza 21 Nov '12

21 Nov '12

Boa tarde, Estou procurando, novamente, indicações. Desta vez de pacotes do R para usar geoestatística, de preferência para iniciantes hehehe Abraços, -- *LIGIA PEREIRA DE SOUZA* *Instituto Brasileiro para Educação e Negócios - IBEN * * * *Graduanda em Ecologia* * * *UNESP - Rio Claro/SP* *Lattes - http://lattes.cnpq.br/3635101224851131<https://wwws.cnpq.br/cvlattesweb/PKG_MENU.menu?f_cod=729600D125BF53DC3FC4A0…> * * *

3 2

problemas da correlação resolvidos
by alanarocha＠sapo.pt 21 Nov '12

21 Nov '12

muito obrigado Daniel Ana

2 1

matriz de correlação
by alanarocha＠sapo.pt 21 Nov '12

21 Nov '12

Olá Daniel, auentão quer dizer q vou ter de fazer varios cor.test pra relacionar X1s09 primeiro com X1s10 depois com Xs11 e depois com X1s12 e assim sucessivamente? Ana

2 1

Re: [R-br] AKAIKE
by Eder Comunello 21 Nov '12

21 Nov '12

Helio, bom dia! Rodei seu script, mas só atentei pra parte do AIC e BIC. O cálculo me parece estar correto. Acredito que o procedimento também, sendo aplicado AIC e BIC em um objeto resultante de 'lm' (data~predicted oriundos da validação cruzada). O BIC é uma derivação do AIC. Tanto que: BIC(x) == AIC(x, k = log(nobs(object))) ou BIC(reglin_ols) == AIC(reglin_ols, k = log(31)) adaptando para um de seus procedimentos (31 observações). Fiz o teste com seus dados e está OK... A diferença entre AIC e BIC é notada com a diferenciação quanto ao número de parâmetros do modelo. No teu caso, você compara um mesmo modelo linear simples sucessivas vezes, não variando o número de parâmetros. Fora isso, por ser um modelo simples não deve haver diferença na penalização por AIC ou BIC, diferindo apenas o cálculo inicial dos índices. Isso é só um 'palpite' e precisaria dar uma estudada na teoria dos índices pra entender melhor. Talvez algum dos colegas da lista possa contribuir mais nesse sentido... ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00

1 0

matriz de correlação
by alanarocha＠sapo.pt 21 Nov '12

21 Nov '12

O problema dos pacotes esta resolvido. agora gostava de saber porque é que me dá erro na instrução return e como o resolvo. obrigada correl<-cor(dados[,2:5],method=c("pearson")) correl result=correl return(list(sumario=summary(dados[,2:5]),resultado=result)) > result=correl > return(list(sumario=summary(dados[,2:5]),resultado=result)) Error: no function to return from, jumping to top level e este outro erro: cor_with_p_test(dados[2:5]) ct <- cor.test(dados[,2:5]) ct$p ct$v cov(dados[,2:5], na.rm=TRUE) Error in cov(dados[, 2:5], na.rm = TRUE) : unused argument(s) (na.rm = TRUE) Ana

2 1

matriz de correlação
by alanarocha＠sapo.pt 21 Nov '12

21 Nov '12

olá, como indico ao R que quero usar este package? cor.ps package=unknown R Documentation **Description: Calculo dos coeficientes de correlação de Pearson/Spearman de um data.frame com duas ou mais variáveis, mostra matriz de espalhagrama e análise exploratória de dados omitindo os valores de NA quando presentes e mostrando um boxplot dos dados iniciais. **Usage: cor.ps(x, Pearson=TRUE,NAs=TRUE) **Arguments: x: data.frame dos dados que se deseja calcular o coeficiente de correlação de Pearson/Spearman Pearson: lógico. Se TRUE, calcula o coeficiente de correlação de Pearson. Se FALSE, calcula o coeficiente de correlação de Spearman. NAs: lógico. Se TRUE, indica presença de NAs nos dados e retira para calcular o coeficiente de correlação. Se FALSE, indica ausência de NAs nos dados e ele calcula o coeficiente de correlação sem retirar os valores de NAs. **Details: O input dessa função deve ser um data.frame contendo duas ou mais variáveis, todas com mesmo comprimento. x - no primeiro argumento da função, quando deseja calcular o coeficiente de correlação de algumas variáveis do data.frame, deve-se indexar o valor correspondente ao das colunas que se queira calcular a correlação. (Ex: dados[,3:6] - dos dados desejo calcular a correlação das colunas 3 a 6). Como análise exploratória dos dados iniciais, além do sumário, um boxplot dos dados é gerado com o intuito de melhorar a visualização dos dados, bem como identificar presença de outliers que possam interferir nos resultados e uma matriz de espalhagrama dos dados originais. Matriz de espalhagrama dos resultados - representação gráfica dos resultados da matriz de correlação de Pearson/Spearman **Value: Return da função em formato de lista: sumario: summary dos dados do data.frame original resultado: result do coeficiente de correlação de Pearson/Spearman **Warning: Para calcular o coeficiente de correlação de algumas das variáveis do data.frame, utilizar indexação, como mostrado nos detalhes do help. Para calcular o coeficiente de correlação, os dados das variáveis devem ter mesmo comprimento. **Author(s): Mariana Rodrigues dos Santos santos.mariana(a)usp.br

3 2

matriz de correlação
by alanarocha＠sapo.pt 21 Nov '12

21 Nov '12

Olá Natália, o problema é esse o R não esta a reconhecer a package cor.ps Ana

2 1

Dúvida em summary em regressão logística
by Alexandre Santos 21 Nov '12

21 Nov '12

Boa tarde Pessoal, Estou ajustando uma regressão logística e me deparei com a seguinte dúvida: ################################################################################################### #Regressão logística entre a proporção de herbívoros predados e a densidade de herbívoros oferecidos # pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2 ,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26 ,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros predados dens<-sort(rep(2^(2:6),10))#### Densidade de herbivoros oferecidos ## ##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ exp(P0+P1*N+P2*N^2+P3*N3) p.model1<-glm(pred/par~dens+I(dens^2)+I(dens^3),family="binomial") summary(p.model1) # Minha dúvida é se o coeficiente linear que aparece no summary esta transformado em exp(x)/1+exp(x) ou trata-se do valor sem transformação? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso - Campus Cáceres Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

4 6

Re: [R-br] excluir mar no ETOPO
by Eder Comunello 21 Nov '12

21 Nov '12

Carlos, boa noite! Você chegou a ler as imagens com raster ou rgdal? Em que ponto exatamente está com problemas? Na leitura ou no fatiamento? Dá pra fatiar o objeto criado por raster ao invés daquele criado por rgdal... O link que segue tem algumas considerações sobre o uso de rgdal() e raster para leitura de imagens do ETOPO1. Comenta sobre causas de erro e dá algumas opções. http://r-sig-geo.2731867.n2.nabble.com/Import-Part-of-a-GeoTiff-using-readG… Espero que ajude, ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00 Em 20 de novembro de 2012 16:51, <r-br-request(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br-request@listas.c3sl…> > escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br-request@listas.c3sl…> > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br-owner@listas.c3sl.u…> > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: excluir mar no ETOPO (Thiago Veloso) > 2. Re: excluir mar no ETOPO (D'Apolito) > 3. Graficos de decomposição (kaue veras) > 4. Re: Gráfico de distribuiçoes estimadas (Matheus Monteiro) > 5. Re: (sem assunto) (Edson Lira) > 6. Re: Usando o R em grandes bases de dados (Leonard de Assis) > 7. matrizes de correlação (alanarocha(a)sapo.pt<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alanarocha@sapo.pt> > ) > 8. ExpDes - Comparações múltiplas (Andre Oliveira) > 9. Re: ExpDes - Comparações múltiplas (Walmes Zeviani) > 10. Re: Dúvida sobre IC simultâneos para modelagem de variáveis > binarias no R (multcomp) (Walmes Zeviani) > 11. Re: Problema para reinstalar pacote (ExpDes) no Linux Ubuntu > 10.04 (Andre Oliveira) > 12. Re: Problema para reinstalar pacote (ExpDes) no Linux Ubuntu > 10.04 (Andre Oliveira) > 13. Dúvida em summary em regressão logística (Alexandre Santos) > 14. Dúvida em summary em regressão logística (Alexandre Santos) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 20 Nov 2012 07:22:57 -0800 (PST) > From: Thiago Veloso <thi_veloso(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=thi_veloso@yahoo.com.br> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: Re: [R-br] excluir mar no ETOPO > Message-ID: > <1353424977.92989.YahooMailNeo(a)web121903.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353424977.92989.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Carlos, > > Eu tentaria algo assim: > > library (raster) #carrega pacote exigido > r <- raster('ETOPO1_Ice_g_geotiff.tif') #carrega arquivo de trabalho > r [ r<0 ] <- NA #atribui NA a todos os valores negativos no raster. > > E caso você queira atribuir NA a um escopo de valores, pode ser mais fácil > criar uma função: > > > vals <- c(0,24,32,56) # estes são os valores que você quer eliminar da > tabela de dados do raster > fun <- function (x) { > x <- ifelse (x %in% vals, 1, NA) # aqui você substitui os valores acima > por "1" e atribui NA para todo o resto. Você pode substituir o "1" por vals > } > > E então aplicá-la sobre o raster: > r.new <- calc (r, fun) > > > Uma vantagem de usar funções é pode realizar uma tarefa para diversos > arquivos sem precisar repetir comandos para cada um deles. > > Espero que ajude, > Thiago. > > > ________________________________ > From: D'Apolito <carlosdapolito(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=carlosdapolito@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Sent: Monday, November 19, 2012 8:15 PM > Subject: [R-br] excluir mar no ETOPO > > > Prezados, > > Estou graficando mapas com os dados do ETOPO01 de elevação e gostaria de > excluir o mar, para que os dados iniciem em zero ao invés de números > negativos (profundidades dos oceanos). > Qualquer ajuda é bem vinda, desde já agradeço > > Att > Carlos > > #código: > > library(raster) > library(rgdal) > > etopo<-"ETOPO1_Ice_g_geotiff.tif" > r.etopo<-raster(etopo) > > ext <- extent(-95,-35,-31,31) > > zoom(r.etopo, ext=ext,col=rev(gray.colors(12))) > > > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/675bda4e/atta… > > > > ------------------------------ > > Message: 2 > Date: Tue, 20 Nov 2012 15:22:24 +0000 > From: "D'Apolito" <carlosdapolito(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=carlosdapolito@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Subject: Re: [R-br] excluir mar no ETOPO > Message-ID: > <CAObczOaGmp3ug9u_ktRu7ZfSgoQ3XUFD9V8yg= > 440j-wZhVkuQ(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=440j-wZhVkuQ@mail.gmail.…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Ola Eder, > > O arquivo de origem dos dados (ETOPO1) `e muito grande, de maneira que > tenho obtido erro constantemente, Tentei aumentar a capacidade da memoria > com memory.limit() mas tambem nao funcionou. > > > att > Carlos > > > Em 20 de novembro de 2012 11:30, Eder Comunello <ecomunel(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ecomunel@gmail.com> > >escreveu: > > > Carlos e colegas, bom dia! > > > > Vou aproveitar um script anterior, fazendo algumas alterações. Como da > > outra vez, creio que dá pra adaptar facilmente pro seu caso. > > > > Vou trabalhar extraindo o vetor de interesse, substituindo este vetor no > > objeto origem ao final... > > > > ##### <begin> ##### > > require(sp); require(rgdal) > > > > ### Download do exemplo (399 Kb) > > download.file('https://dl.dropbox.com/u/117618178/dourados/DDOS.tif', > > file.path(getwd(),'DDOS.tif'), mode='wb') > > > > alt = readGDAL("DDOS.tif"); gridded(alt) ### 'importa' o geotiff > > image(alt, col=terrain.colors(21)) ## visualização > > slotNames(alt) ### nome dos slots (S4) > > proj4string(alt) ### projeçãoo definida no geotiff > > head(alt@data) ### observe os valores dos dados > > head(alt@data$band1) ### observe os valores dos dados > > > > ### Opção para fatiar (0 ou 1) > > > > alt2 <- alt ### para manipular uma cópia do objeto > > z <- alt2@data$band1 ### para manipular apenas o vetor de interesse > > n <- length(z); n > > > > for(i in 1:n){if (z[i]<400){z[i]=0} else{z[i]=1}}; head(z) ### 0 ou 1 > > alt2@data$band1 <- z ### substitui o vetor manipulado no objeto > > image(alt2) ## visualização > > > > ### Opção para mascarar valores abaixo de critério (com 0 ou NA) > > > > alt3 <- alt ### para manipular outra cópia do objeto > > z <- alt3@data$band1 ### para manipular apenas o vetor de interesse > > n <- length(z); n > > > > for(i in 1:n){if (z[i]<400){z[i]=0}}; head(z) ### 0 ou 1 > > #for(i in 1:n){if (z[i]<400){z[i]=NA}}; head(z) ### 0 ou 1 > > alt3@data$band1 <- z ### substitui o vetor manipulado no objeto > > image(alt3) ## visualização > > ##### <end> ##### > > > > > > ================================================ > > Éder Comunello > > > > Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) > > Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] > > > > Researcher at Embrapa Western Region Agriculture > > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > > ================================================ > > UTC-03:00 > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > > > > > -- > _____________________________________ > *Carlos D'Apolito - PhD candidate - ID:1239480* > Earth Sciences > School of Geography & Environmental Sciences > Aston Webb Building > University of Birmingham > Birmingham > B15 2TT > UK > > phone: 44 - 07 919 673 173 > CXD280(a)bham.ac.uk<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=CXD280@bham.ac.uk> > CV: http://lattes.cnpq.br/4557754746424026 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/615af15d/atta… > > > > ------------------------------ > > Message: 3 > Date: Tue, 20 Nov 2012 14:14:30 -0200 > From: kaue veras <kaueveras(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kaueveras@gmail.com> > > > To: R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@listas.c3sl.ufpr.br> > Subject: [R-br] Graficos de decomposição > Message-ID: > <CAENvtmjrMnC52i4-JYn2NdhrwWMNw= > Qu0HJLeTOqz2s5q2i6uQ(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=Qu0HJLeTOqz2s5q2i6uQ@mai…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde prezados, > > Estou realizando um trabalho em que utilizo a decomposição, análises de > previsão usando o modelo ARIMA e Holt-Winters e goataria de mudar os > títulos do gráfico e titulos dos eixos. Com os comandos main, xlab e ylab > não esta dando certo, poderiam me dar uma luz ou me ajudar como faço para > gerar essas modificações? > > Desde já agradeço. > > Atenciosamente, > > Kauê P. Veras da Cunha > Estatística - UERJ > MSN: kaueveras(a)hotmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kaueveras@hotmail.com> > E-mail: kaueveras(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kaueveras@gmail.com> > Cel: (21) 8254-9601 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/1154d3e8/atta… > > > > ------------------------------ > > Message: 4 > Date: Tue, 20 Nov 2012 14:34:00 -0200 > From: Matheus Monteiro <matheusmont(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=matheusmont@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Subject: Re: [R-br] Gráfico de distribuiçoes estimadas > Message-ID: > < > CAAG8985KZaYb24ckqaGFxeon85H3g18ny9gEtBc79DXK7p8mcg(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=CAAG8985KZaYb24ckqaGFxeo…> > > > Content-Type: text/plain; charset="iso-8859-1" > > mais alguem pode me ajudar? > > Em 10 de novembro de 2012 13:21, Edson Lira <edinhoestat(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=edinhoestat@yahoo.com.br> > >escreveu: > > > Mateus, no pacote MASS acho que tem alguns gráficos empírica x teórica. > > > > [. ]'s. > > Edson Lira > > Estatístico > > Ma-Am > > > > Em 10/11/2012, às 02:20, Matheus Monteiro <matheusmont(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=matheusmont@gmail.com> > > > > escreveu: > > > > Boa noite, obtive as densidades aproximadas da 't' e da 'normal', > gostaria > > de compara-las através de graficos dessas densidades estimadas. > > Alguma sugestão, como faria isso? > > > > *t * > > > > > fitdistr(x[,1],"t") > > > > m s df > > > > -0.001874414 0.009122738 3.522064232 > > > > ( 0.001726835) ( 0.001862692) ( 2.234067399)* > > > > **normal* > > > > fitdistr(x[,1],"normal") > > > > mean sd > > > > 0.0001040816 0.0128201843 > > > > (0.0018314549) (0.0012950342) > > > > > > x são dados simulados. > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/1c03edac/atta… > > > > ------------------------------ > > Message: 5 > Date: Tue, 20 Nov 2012 09:06:00 -0800 (PST) > From: Edson Lira <edinhoestat(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=edinhoestat@yahoo.com.br> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: Re: [R-br] (sem assunto) > Message-ID: > <1353431160.32913.YahooMailNeo(a)web161205.mail.bf1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353431160.32913.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Mas Eder a idéia do NA é não considerar estas respostas na contagem e/ou > no percentual. > > Mas valeu e muito obrigado. > > Peço desculpas a todos, esqueci de colocar o assunto na mensagem > > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: Eder Comunello <ecomunel(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ecomunel@gmail.com> > > > Para: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Enviadas: Terça-feira, 20 de Novembro de 2012 7:39 > Assunto: Re: [R-br] (sem assunto) > > > Edson e colegas, bom dia! > > Ao invés de deixar o campo 'motivo' em branco (NA), experimente atribuir > um valor, p. ex. 'sem resposta', e rode novamente. Assim estes casos serão > considerados e aparecerão na saída. Vale pros outros campos. > > ================================================ > Éder Comunello > > Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) > Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] > > Researcher at Embrapa Western Region Agriculture > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > ================================================ > UTC-03:00 > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/a9ea6c58/atta… > > > > ------------------------------ > > Message: 6 > Date: Tue, 20 Nov 2012 15:30:17 -0200 > From: Leonard de Assis <assis.leonard(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=assis.leonard@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Subject: Re: [R-br] Usando o R em grandes bases de dados > Message-ID: <50ABBE29.4030805(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=50ABBE29.4030805@gmail.c…> > > > Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" > > Bigmemory eu acredito que lide somente com matrizes, não com data > frames, resumo: todos os campos devem ter o mesmo tipo. > > []s > Leonard de Assis > http://about.me/ldeassis > > Em 20/11/2012 03:06, Vinicius Brito Rocha escreveu: > > > > fiquei brincando com o bigmemory e carreguei uma tabela mas os campos > > character vieram como NA. > > > > Quando usava a opção type="char" piorava e alterava os valores dos > > campos numéricos. O melhor ( ou menos pior resultado) foi carregar a > > tabela com a opção type="double" > > > > Alguem sabe o que pode ser? > > > > > > Abs. > > > > Vinicius Brito Rocha > > > > Em 19 de novembro de 2012 15:40, Richard Santos > > <jamesrichardsantos(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=jamesrichardsantos@gmail…><mailto: > jamesrichardsantos(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=jamesrichardsantos@gmail…> > >> > > escreveu: > > > > Valeu!! Toda e qq ajuda nesse sentido é bem vinda! > > Abç. > > Richard > > > > Em 19/11/12, Daniel C Bezerra<danielcbezerra(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=danielcbezerra@gmail.com> > > <mailto:danielcbezerra@gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=danielcbezerra@gmail.com>>> > escreveu: > > > Eu encontrei estas apresentações e achei legal dividir com vcs: > > > > > > http://goo.gl/l6fuC > > > > > > Tem link para os slides, código do R e tudo mais por lá. > > > > > > Abs, > > > > > > Daniel > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@listas.c3sl.ufpr.br><mailto: > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > > > > > > > > > -- > > /Vinicius Brito Rocha./ > > /Estatístico e Atuário (IM / UFRJ)// > > Mestre em Pesquisa Operacional (COPPE / UFRJ)/ > > > > www.aplicademic.blogspot.com <http://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. > > > > > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/b8d49acf/atta… > > > > ------------------------------ > > Message: 7 > Date: Tue, 20 Nov 2012 17:51:25 +0000 > From: alanarocha(a)sapo.pt<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alanarocha@sapo.pt> > To: R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@listas.c3sl.ufpr.br> > Subject: [R-br] matrizes de correlação > Message-ID: > <20121120175125.Horde.X_LLT3CXF3VQq8MdQ0ASvYA(a)mail.sapo.pt<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=20121120175125.Horde.X_L…> > > > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > Olá, a actualização do R resultou para o meu problema. > agora preciso de ajuda pois tou a tentar fazer uma matriz de > correlação e outra de covariancias de 4 variaveis de um dataframe e > dame estes erros: > ##matriz de correlação > > cor.p<-cor.ps(dados[,2:4], Pearson=TRUE,NAs=TRUE) > cor.p > dados=na.omit(dados[,2:4]) > > cat("valores NA retirados\n") > x11() > pairs(dados) > x11() > boxplot(dados) > dados=dados[,2:4] > x11() > boxplot(dados) > correl=cor(dados[2:4]) > result=correl > return(list(sumario=summary(dados),resultado=result)) > cor(cbind(X1s09, X1s10, X1s11, X1s12)) > cor.test(dados[,2:4], method='pearson')$p.value > cor_with_p_test(cbind(X1s09, X1s10, X1s11, X1s12)) > ct <- cor.test(X1s10, X1s12) > ct$p > ct$v > > #matriz de covariancias > cov(dados[,2:4], na.rm=TRUE) > > > > cor.p<-cor.ps(dados[,2:4], Pearson=TRUE,NAs=TRUE) > Error: could not find function "cor.ps" > > cor.p > Error: object 'cor.p' not found > > dados=na.omit(dados[,2:4]) > > > > cat("valores NA retirados\n") > valores NA retirados > > > > > > > > > > correl=cor(dados[2:4]) > Error in `[.data.frame`(dados, 2:4) : undefined columns selected > > result=correl > Error: object 'correl' not found > > return(list(sumario=summary(dados),resultado=result)) > Error: object 'result' not found > > cor(cbind(X1s09, X1s10, X1s11, X1s12)) > Error in cbind(X1s09, X1s10, X1s11, X1s12) : object 'X1s09' not found > > cor.test(dados[,2:4], method='pearson')$p.value > Error in `[.data.frame`(dados, , 2:4) : undefined columns selected > > cor_with_p_test(cbind(X1s09, X1s10, X1s11, X1s12)) > Error: could not find function "cor_with_p_test" > > ct <- cor.test(X1s10, X1s12) > Error in cor.test(X1s10, X1s12) : object 'X1s10' not found > > ct$p > Error: object 'ct' not found > > ct$v > Error: object 'ct' not found > > > > #matriz de covariancias > > cov(dados[,2:4], na.rm=TRUE) > Error in cov(dados[, 2:4], na.rm = TRUE) : > unused argument(s) (na.rm = TRUE) > > > cumprimentos e obrigada Ana Rocha > > > > > > ------------------------------ > > Message: 8 > Date: Tue, 20 Nov 2012 09:57:15 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: [R-br] ExpDes - Comparações múltiplas > Message-ID: > <1353434235.56527.YahooMailNeo(a)web125902.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353434235.56527.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde, > estou aqui tentando reproduzir algumas análises do ExpDes. A analise > abaixo não me deu as comparações múltiplas como faço para tais comparações? > > require(ExpDes) > Carregando pacotes exigidos: ExpDes > data(ex1) > attach(ex1) > dic(trat, ig, quali = FALSE, sigF = 0.05) > > > > André Oliveira Souza > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/696d59ab/atta… > > > > ------------------------------ > > Message: 9 > Date: Tue, 20 Nov 2012 16:24:10 -0200 > From: Walmes Zeviani <walmeszeviani(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmeszeviani@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>, > Andre Oliveira > <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Subject: Re: [R-br] ExpDes - Comparações múltiplas > Message-ID: > <CAFU= > EkZQNQHMTUM2nUKMnJhKTcCoDJ7ngqJ+_69-XZ6mf3FJDQ(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=EkZQNQHMTUM2nUKMnJhKTcCo…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Você usou quali=FALSE, ou seja, está declarando que os níveis são não > qualitativos. Use quali=TRUE e leia o manual. > > À disposiçã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<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br> > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/6b5bc51b/atta… > > > > ------------------------------ > > Message: 10 > Date: Tue, 20 Nov 2012 16:28:49 -0200 > From: Walmes Zeviani <walmeszeviani(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmeszeviani@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Subject: Re: [R-br] Dúvida sobre IC simultâneos para modelagem de > variáveis binarias no R (multcomp) > Message-ID: > <CAFU=EkYKvcE65FEpjv=qb0xnJ=mME7ec_8c7DD28wQNPF= > 1_yg(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1_yg@mail.gmail.com> > > > Content-Type: text/plain; charset="iso-8859-1" > > Seu exemplo não é reproduzível. O guia de postagem não recomenda envio de > dados como anexo das mensagens. Hospede os dados na núvem e acrescente > comandos para leitura do arquivo. Assim rapidamente podemos reproduzir sua > análise e ajudar você. > > À disposiçã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<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br> > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/7d4a7d34/atta… > > > > ------------------------------ > > Message: 11 > Date: Tue, 20 Nov 2012 10:30:37 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > To: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > >, > "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > Message-ID: > <1353436237.16303.YahooMailNeo(a)web125901.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353436237.16303.YahooMa…> > > > Content-Type: text/plain; charset="utf-8" > > Versão ExpDes_1.1.1.tar.gz agora sim consegui instalar no linux. > > > > André Oliveira Souza > > > ________________________________ > De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Para: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Enviadas: Quinta-feira, 5 de Janeiro de 2012 11:05 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > O Walmes me disse, aqui mesmo na lista, ha algum tempo que realmente o > pacote para linux em português ele não consegui colocar para funcionar. E > no caso do windows eu não sabia que estava em inglês! Obrigado > > > > André Oliveira Souza > > > ________________________________ > De: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Para: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Enviadas: Quinta-feira, 5 de Janeiro de 2012 10:55 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > Olá Andre, > > Que bom que não desistiu. Vamos por partes que solucionaremos todos os > seus problemas. > Quanto ao Windows, > não se esqueça que a versão do CRAN está em ingles. Portanto, após > carregar, voce deve pedir o help da função cdr (completely randomized > design) e não dic, ou seja > ?crd > Veja se funciona. > > > Quanto ao linux, talvez voce deva escrever para o Walmes, que faz parte > dessa lista. Ele usa o ExpDes no Linux e diz que funciona perfeitamente. > > abs > > Eric. > > > Em 5 de janeiro de 2012 10:42, Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br>> > escreveu: > > Eric, > > > >eu acompanhei o pacote no seu site e nunca consegui coloca-lo para > funcionar no Linux a versão em português. Agora com ele no CRAN voltei a > tentar colocar o pacote para funcionar no Linux, mas sem sucesso. E para > minha surpresa no windows não está funcionado também. Olha a minha > tentativa no windows! > >require(ExpDes) > >Carregando pacotes exigidos: ExpDes > >> ?dic > >No documentation for ?dic? in specified packages and libraries: > >you could try ???dic? > > > > > > Estou ansioso para colocar este pacote para funcionar por aqui. Já li o > TCC e o manual sobre o mesmo e vi a importância do mesmo. Visto que antes > tinha que usar muitos pacotes para a mesma análise e com o ExpDes é muito > mais prático. > > > > > > > > > > > > > > > > > > > > > >André Oliveira Souza > > > > > > > > > > > > > >________________________________ > > De: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > >Para: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>; > Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > >Enviadas: Quarta-feira, 4 de Janeiro de 2012 18:25 > >Assunto: Re: [R-br] Problema para reinstalar pacote (EspDes) no Linux > Ubuntu 10.04 > > > > > >Olá André, > > > >Se quiser trocar uma ideia em particular, posso tentar te ajudar. Sou o > criador do ExpDes, mas não usu Linux. Não sei se serei útil. > > > >A princípio, posso te dizer que existe o pacote no CRAN e na minha página > pessoal: https://sites.google.com/site/ericbferreira/unifal/downloads-1 > > > >Tente pegar o pacote de lá para ver se o erro se repete. > >Peço gentilmente aos usuárioas de (R no) Linux que me ajudem a exclarecer > o André. > > > >Atenciosamente, > > > >Eric. > > > > > >Em 4 de janeiro de 2012 16:39, Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br>> > escreveu: > > > >Opa vacilo meu erro no nome do pacote > >> > >> > >> install.packages("ExpDes", dep=T) > >> > >> > >> > >> > >> > >> > >>André Oliveira Souza > >> > >> > >>________________________________ > >> De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > >>Para: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > >>Enviadas: Quarta-feira, 4 de Janeiro de 2012 16:36 > >>Assunto: [R-br] Problema para reinstalar pacote (EspDes) no Linux Ubuntu > 10.04 > >> > >> > >> > >>Pessoal estou com problemas para reinstalar o EspDes !! Alguém saberia > dizer o que pode ser? > >> > >> > >> > >> > >>install.packages("EspDes", dep=T) > >>Installing package(s) into ?/home/andre/R/i486-pc-linux-gnu-library/2.14? > >>(as ?lib? is unspecified) > >>--- Please select a CRAN mirror for use in this session --- > >>Carregando a interface Tcl/Tk ... pronto > >>Mensagens de aviso perdidas: > >>In getDependencies(pkgs, dependencies, available, lib) : > >> package ?EspDes? is not available (for R version 2.14.1) > >> > >> > >> > >> > >> > >> > >> > >>André Oliveira Souza > >>obrigado > >> > >> > >>_______________________________________________ > >>R-br mailing list > >>R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > >> > >> > >>_______________________________________________ > >>R-br mailing list > >>R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > >> > > > > > >-- > >Dr Eric B Ferreira > >Exact Sciences Department > >Federal University of Alfenas > >Brazil > > > > > > > > > -- > Dr Eric B Ferreira > Exact Sciences Department > Federal University of Alfenas > Brazil > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/20c67bb9/atta… > > > > ------------------------------ > > Message: 12 > Date: Tue, 20 Nov 2012 10:37:54 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>>, > Eric > Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Subject: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > Message-ID: > <1353436674.31607.YahooMailNeo(a)web125904.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353436674.31607.YahooMa…> > > > Content-Type: text/plain; charset="utf-8" > > Em português no linux! > > André Oliveira Souza > > > ________________________________ > De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Para: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com>>; > "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Enviadas: Terça-feira, 20 de Novembro de 2012 16:30 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > Versão ExpDes_1.1.1.tar.gz agora sim consegui instalar no linux. > > > > André Oliveira Souza > > > ________________________________ > De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Para: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Enviadas: Quinta-feira, 5 de Janeiro de 2012 11:05 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > O Walmes me disse, aqui mesmo na lista, ha algum tempo que realmente o > pacote para linux em português ele não consegui colocar para funcionar. E > no caso do windows eu não sabia que estava em inglês! Obrigado > > > > André Oliveira Souza > > > ________________________________ > De: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Para: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Enviadas: Quinta-feira, 5 de Janeiro de 2012 10:55 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > Olá Andre, > > Que bom que não desistiu. Vamos por partes que solucionaremos todos os > seus problemas. > Quanto ao Windows, > não se esqueça que a versão do CRAN está em ingles. Portanto, após > carregar, voce deve pedir o help da função cdr (completely randomized > design) e não dic, ou seja > ?crd > Veja se funciona. > > > Quanto ao linux, talvez voce deva escrever para o Walmes, que faz parte > dessa lista. Ele usa o ExpDes no Linux e diz que funciona perfeitamente. > > abs > > Eric. > > > Em 5 de janeiro de 2012 10:42, Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br>> > escreveu: > > Eric, > > > >eu acompanhei o pacote no seu site e nunca consegui coloca-lo para > funcionar no Linux a versão em português. Agora com ele no CRAN voltei a > tentar colocar o pacote para funcionar no Linux, mas sem sucesso. E para > minha surpresa no windows não está funcionado também. Olha a minha > tentativa no windows! > >require(ExpDes) > >Carregando pacotes exigidos: ExpDes > >> ?dic > >No documentation for ?dic? in specified packages and libraries: > >you could try ???dic? > > > > > > Estou ansioso para colocar este pacote para funcionar por aqui. Já li o > TCC e o manual sobre o mesmo e vi a importância do mesmo. Visto que antes > tinha que usar muitos pacotes para a mesma análise e com o ExpDes é muito > mais prático. > > > > > > > > > > > > > > > > > > > > > >André Oliveira Souza > > > > > > > > > > > > > >________________________________ > > De: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > >Para: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>; > Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > >Enviadas: Quarta-feira, 4 de Janeiro de 2012 18:25 > >Assunto: Re: [R-br] Problema para reinstalar pacote (EspDes) no Linux > Ubuntu 10.04 > > > > > >Olá André, > > > >Se quiser trocar uma ideia em particular, posso tentar te ajudar. Sou o > criador do ExpDes, mas não usu Linux. Não sei se serei útil. > > > >A princípio, posso te dizer que existe o pacote no CRAN e na minha página > pessoal: https://sites.google.com/site/ericbferreira/unifal/downloads-1 > > > >Tente pegar o pacote de lá para ver se o erro se repete. > >Peço gentilmente aos usuárioas de (R no) Linux que me ajudem a exclarecer > o André. > > > >Atenciosamente, > > > >Eric. > > > > > >Em 4 de janeiro de 2012 16:39, Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br>> > escreveu: > > > >Opa vacilo meu erro no nome do pacote > >> > >> > >> install.packages("ExpDes", dep=T) > >> > >> > >> > >> > >> > >> > >>André Oliveira Souza > >> > >> > >>________________________________ > >> De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > >>Para: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > >>Enviadas: Quarta-feira, 4 de Janeiro de 2012 16:36 > >>Assunto: [R-br] Problema para reinstalar pacote (EspDes) no Linux Ubuntu > 10.04 > >> > >> > >> > >>Pessoal estou com problemas para reinstalar o EspDes !! Alguém saberia > dizer o que pode ser? > >> > >> > >> > >> > >>install.packages("EspDes", dep=T) > >>Installing package(s) into ?/home/andre/R/i486-pc-linux-gnu-library/2.14? > >>(as ?lib? is unspecified) > >>--- Please select a CRAN mirror for use in this session --- > >>Carregando a interface Tcl/Tk ... pronto > >>Mensagens de aviso perdidas: > >>In getDependencies(pkgs, dependencies, available, lib) : > >> package ?EspDes? is not available (for R version 2.14.1) > >> > >> > >> > >> > >> > >> > >> > >>André Oliveira Souza > >>obrigado > >> > >> > >>_______________________________________________ > >>R-br mailing list > >>R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > >> > >> > >>_______________________________________________ > >>R-br mailing list > >>R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > >> > > > > > >-- > >Dr Eric B Ferreira > >Exact Sciences Department > >Federal University of Alfenas > >Brazil > > > > > > > > > -- > Dr Eric B Ferreira > Exact Sciences Department > Federal University of Alfenas > Brazil > > > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/6a5b8ac7/atta… > > > > ------------------------------ > > Message: 13 > Date: Tue, 20 Nov 2012 10:49:58 -0800 (PST) > From: Alexandre Santos <alexandresantosbr(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alexandresantosbr@yahoo.…> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: [R-br] Dúvida em summary em regressão logística > Message-ID: > <1353437398.81186.YahooMailNeo(a)web121206.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353437398.81186.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde Pessoal, > Estou ajustando uma regressão logística e me deparei com a seguinte > dúvida: > > > ################################################################################################### > #Regressão logística entre a proporção de herbívoros predados e a > densidade de herbívoros oferecidos > # > pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2 > ,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26 > ,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros > predados > dens<-sort(rep(2^(2:6),10))#### Densidade de herbivoros oferecidos > ## > ##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ > exp(P0+P1*N+P2*N^2+P3*N3) > p.model1<-glm(pred/par~dens+I(dens^2)+I(dens^3),family="binomial") > summary(p.model1) > # > > > Minha dúvida é se o coeficiente linear que aparece no summary esta > transformado em exp(x)/1+exp(x) > ou trata-se do valor sem transformação? > > Obrigado, > > -- > ====================================================================== > Alexandre dos Santos > Proteção Florestal > IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso > - Campus Cáceres > Avenida dos Ramires, s/n > Bairro: Distrito Industrial > Cáceres - MT CEP: 78.200-000 > Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) > e-mails:alexandresantosbr@yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=e-mails%3Aalexandresanto…> > alexandre.santos(a)cas.ifmt.edu.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alexandre.santos@cas.ifm…>====================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/0b764a8d/atta… > > > > ------------------------------ > > Message: 14 > Date: Tue, 20 Nov 2012 10:51:20 -0800 (PST) > From: Alexandre Santos <alexandresantosbr(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alexandresantosbr@yahoo.…> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: [R-br] Dúvida em summary em regressão logística > Message-ID: > <1353437480.93196.YahooMailNeo(a)web121202.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353437480.93196.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde Pessoal, > Estou ajustando uma regressão logística e me deparei com a seguinte > dúvida: > > > ################################################################################################### > #Regressão logística entre a proporção de herbívoros predados e a > densidade de herbívoros oferecidos > # > pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2 > ,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26 > ,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros > predados > dens<-sort(rep(2^(2:6),10))#### Densidade de herbivoros oferecidos > ## > ##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ > exp(P0+P1*N+P2*N^2+P3*N3) > p.model1<-glm(pred/par~dens+I(dens^2)+I(dens^3),family="binomial") > summary(p.model1) > # > > > Minha dúvida é se o coeficiente linear que aparece no summary esta > transformado em exp(x)/1+exp(x) > ou trata-se do valor sem transformação? > > Obrigado, > > -- > ====================================================================== > Alexandre dos Santos > Proteção Florestal > IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso > - Campus Cáceres > Avenida dos Ramires, s/n > Bairro: Distrito Industrial > Cáceres - MT CEP: 78.200-000 > Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) > e-mails:alexandresantosbr@yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=e-mails%3Aalexandresanto…> > alexandre.santos(a)cas.ifmt.edu.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alexandre.santos@cas.ifm…>====================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/e2cfa427/atta… > > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@listas.c3sl.ufpr.br> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 21, assunto 47 > ***************************************** >

2 1

Preencher um serie temporal
by Daniel Marcelino 21 Nov '12

21 Nov '12

Caros, Eu gostaria de expandir os dados abaixo para realizar uma análise temporal de base diária. O meu período inicial é '2012-07-01' até '2012-10-06'. É possível eu fazer isso com meu data frame ou eu terei que criar um novo data frame com a sequencia de dias e juntar as duas tabelas? Meu data frame final deveria repetir os valores apresentados nas colunas até que uma nova medida surgisse já que os pontos de medida que tenho são exporádicos. Alguém já expandiu um data frame assim pra me dizer como fazer? structure(list(Date.beg = structure(c(15635, 15617, 15615, 15610, 15609, 15605, 15604, 15601, 15593, 15593, 15586, 15581, 15580, 15577, 15572, 15565, 15552, 15540, 15530, 15516), class = "Date"), Date.end = structure(c(15619, 15619, 15616, 15615, 15610, 15607, 15604, 15602, 15595, 15594, 15587, 15582, 15581, 15579, 15572, 15567, 15554, 15541, 15533, 15517), class = "Date"), Pollster = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 4L, 1L, 2L, 1L, 1L, 2L, 1L, 4L, 1L, 2L, 2L, 1L, 3L, 1L), .Label = c("Datafolha", "Ibope", "Veritá", "Vox Populi"), class = "factor"), Serra.PSDB = c(24, 22, 23, 19, 22, 17, 17, 21, 19, 20, 21, 20, 22, 22, 27, 26, 26, 30, 31.4, 31), Russomanno.PRB = c(23, 22, 25, 27, 30, 34, 34, 35, 35, 32, 35, 31, 31, 31, 31, 26, 25, 26, 17.7, 24), Haddad.PT = c(20, 22, 19, 18, 18, 18, 17, 15, 15, 17, 16, 16, 14, 14, 8, 9, 6, 7, 9.5, 6), Chalita.PMDB = c(11, 11, 11, 10, 9, 7, 5, 8, 6, 8, 7, 5, 7, 5, 6, 5, 5, 6, 4.3, 6), Others = c(8, 7, 7, 4, 7, 6, 4, 5, 6, 7, 9, 5, 8, 6, 10, 12, 15, 15, 8.6, 16), Null = c(8, 8, 8, NA, 8, 10, 10, NA, 13, 9, 8, 12, 10, NA, 10, 12, 14, 11, 13.1, 11), Undecided = c(6, 8, 6, NA, 6, 8, 13, NA, 6, 7, 4, 9, 7, 13, 6, 10, 9, 6, 13.3, 5), Round = structure(c(3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L), .Label = c("year before", "off campaign", "first", "second"), class = "factor"), Stage = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("0", "1", "2"), class = "factor"), Serra.Haddad = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("", "1", "0"), class = "factor"), N = c(3959, 1204, 2099, 1204, 1799, 1204, 2000, 1802, 1001, 1221, 1078, 1001, 1069, 1200, 1077, 805, 805, 1075, 1331, 1081), Err = c(2, 3, 2, 3, 2, 3, 2.2, 2, 3, 3, 3, 3, 3, 2.8, 3, 3, 3, 3, 2.7, 3)), .Names = c("Date.beg", "Date.end", "Pollster", "Serra.PSDB", "Russomanno.PRB", "Haddad.PT", "Chalita.PMDB", "Others", "Null", "Undecided", "Round", "Stage", "Serra.Haddad", "N", "Err"), row.names = c("6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25"), class = "data.frame") -- "Small steps toward a much better world"

2 1

Usando o R em grandes bases de dados
by Daniel C Bezerra 20 Nov '12

20 Nov '12

Eu encontrei estas apresentações e achei legal dividir com vcs: http://goo.gl/l6fuC Tem link para os slides, código do R e tudo mais por lá. Abs, Daniel

5 6

Re: [R-br] (sem assunto)
by Eder Comunello 20 Nov '12

20 Nov '12

Edson, boa noite! Desculpe a falta de atenção, mas na verdade por padrão o NA é desconsiderado. Você deve definir se quer computá-lo, valendo a dica do colega Fernando Mayer. Para demonstrar montei o script abaixo com os dados que vc forneceu... ### criando um objeto com seus dados h <- structure(list(id = 1:31, expect_exc = structure(c(1L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 1L, 2L, 2L, 1L), .Label = c("não", "sim"), class = "factor"), motivo = structure(c(NA, 1L, 1L, 4L, 5L, 7L, NA, 2L, NA, 1L, 1L, 1L, 6L, 4L, 5L, 3L, 3L, 7L, 5L, 7L, 4L, 7L, 4L, 5L, 5L, 4L, 7L, 1L, 5L, 5L, NA), .Label = c("atendimento", "despertar", "espera por apto", "mensageiro", "recepção", "reserva", "todos"), class = "factor"), visitaria.novamente = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, NA, 1L, 1L, 1L, 1L, 1L, 1L, NA, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = "sim", class = "factor"), nao_pq = structure(c(NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, NA, 2L, NA, NA, NA, NA, NA, 1L, 3L), .Label = c("agilidade na troca do apto", "falta de preparo", "hospitalidade e simpatia"), class = "factor")), .Names = c("id", "expect_exc", "motivo", "visitaria.novamente", "nao_pq"), class = "data.frame", row.names = c(NA, -31L)) ### término do objeto head(h) length(h$motivo) ### retorna 31 observações length(h$expect_exc) ### idem sum(table(h$motivo,h$expect_exc)) ### retorna 27 (pois os 4 NA ficaram de fora) table(h$motivo,h$expect_exc) ### tabela de dados sem NA sum(table(h$motivo,h$expect_exc,useNA='always')) ### inclui o NA e portanto retorna 31 observações table(h$motivo,h$expect_exc,useNA='always') ### tabela de dados com NA 100*prop.table(table(h$motivo,h$expect_exc)) ### logo estes valores não consideram o NA (27 obs) 100*prop.table(table(h$motivo,h$expect_exc, useNA='always')) ### como fica se considerar o NA (31 obs) Roda aí pra confirmar! ================================================ Éder Comunello Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] Researcher at Embrapa Western Region Agriculture Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] ================================================ UTC-03:00 Em 20 de novembro de 2012 16:51, <r-br-request(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br-request@listas.c3sl…> > escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br-request@listas.c3sl…> > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br-owner@listas.c3sl.u…> > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Re: excluir mar no ETOPO (Thiago Veloso) > 2. Re: excluir mar no ETOPO (D'Apolito) > 3. Graficos de decomposição (kaue veras) > 4. Re: Gráfico de distribuiçoes estimadas (Matheus Monteiro) > 5. Re: (sem assunto) (Edson Lira) > 6. Re: Usando o R em grandes bases de dados (Leonard de Assis) > 7. matrizes de correlação (alanarocha(a)sapo.pt<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alanarocha@sapo.pt> > ) > 8. ExpDes - Comparações múltiplas (Andre Oliveira) > 9. Re: ExpDes - Comparações múltiplas (Walmes Zeviani) > 10. Re: Dúvida sobre IC simultâneos para modelagem de variáveis > binarias no R (multcomp) (Walmes Zeviani) > 11. Re: Problema para reinstalar pacote (ExpDes) no Linux Ubuntu > 10.04 (Andre Oliveira) > 12. Re: Problema para reinstalar pacote (ExpDes) no Linux Ubuntu > 10.04 (Andre Oliveira) > 13. Dúvida em summary em regressão logística (Alexandre Santos) > 14. Dúvida em summary em regressão logística (Alexandre Santos) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Tue, 20 Nov 2012 07:22:57 -0800 (PST) > From: Thiago Veloso <thi_veloso(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=thi_veloso@yahoo.com.br> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: Re: [R-br] excluir mar no ETOPO > Message-ID: > <1353424977.92989.YahooMailNeo(a)web121903.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353424977.92989.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Carlos, > > Eu tentaria algo assim: > > library (raster) #carrega pacote exigido > r <- raster('ETOPO1_Ice_g_geotiff.tif') #carrega arquivo de trabalho > r [ r<0 ] <- NA #atribui NA a todos os valores negativos no raster. > > E caso você queira atribuir NA a um escopo de valores, pode ser mais fácil > criar uma função: > > > vals <- c(0,24,32,56) # estes são os valores que você quer eliminar da > tabela de dados do raster > fun <- function (x) { > x <- ifelse (x %in% vals, 1, NA) # aqui você substitui os valores acima > por "1" e atribui NA para todo o resto. Você pode substituir o "1" por vals > } > > E então aplicá-la sobre o raster: > r.new <- calc (r, fun) > > > Uma vantagem de usar funções é pode realizar uma tarefa para diversos > arquivos sem precisar repetir comandos para cada um deles. > > Espero que ajude, > Thiago. > > > ________________________________ > From: D'Apolito <carlosdapolito(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=carlosdapolito@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Sent: Monday, November 19, 2012 8:15 PM > Subject: [R-br] excluir mar no ETOPO > > > Prezados, > > Estou graficando mapas com os dados do ETOPO01 de elevação e gostaria de > excluir o mar, para que os dados iniciem em zero ao invés de números > negativos (profundidades dos oceanos). > Qualquer ajuda é bem vinda, desde já agradeço > > Att > Carlos > > #código: > > library(raster) > library(rgdal) > > etopo<-"ETOPO1_Ice_g_geotiff.tif" > r.etopo<-raster(etopo) > > ext <- extent(-95,-35,-31,31) > > zoom(r.etopo, ext=ext,col=rev(gray.colors(12))) > > > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/675bda4e/atta… > > > > ------------------------------ > > Message: 2 > Date: Tue, 20 Nov 2012 15:22:24 +0000 > From: "D'Apolito" <carlosdapolito(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=carlosdapolito@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Subject: Re: [R-br] excluir mar no ETOPO > Message-ID: > <CAObczOaGmp3ug9u_ktRu7ZfSgoQ3XUFD9V8yg= > 440j-wZhVkuQ(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=440j-wZhVkuQ@mail.gmail.…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Ola Eder, > > O arquivo de origem dos dados (ETOPO1) `e muito grande, de maneira que > tenho obtido erro constantemente, Tentei aumentar a capacidade da memoria > com memory.limit() mas tambem nao funcionou. > > > att > Carlos > > > Em 20 de novembro de 2012 11:30, Eder Comunello <ecomunel(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ecomunel@gmail.com> > >escreveu: > > > Carlos e colegas, bom dia! > > > > Vou aproveitar um script anterior, fazendo algumas alterações. Como da > > outra vez, creio que dá pra adaptar facilmente pro seu caso. > > > > Vou trabalhar extraindo o vetor de interesse, substituindo este vetor no > > objeto origem ao final... > > > > ##### <begin> ##### > > require(sp); require(rgdal) > > > > ### Download do exemplo (399 Kb) > > download.file('https://dl.dropbox.com/u/117618178/dourados/DDOS.tif', > > file.path(getwd(),'DDOS.tif'), mode='wb') > > > > alt = readGDAL("DDOS.tif"); gridded(alt) ### 'importa' o geotiff > > image(alt, col=terrain.colors(21)) ## visualização > > slotNames(alt) ### nome dos slots (S4) > > proj4string(alt) ### projeçãoo definida no geotiff > > head(alt@data) ### observe os valores dos dados > > head(alt@data$band1) ### observe os valores dos dados > > > > ### Opção para fatiar (0 ou 1) > > > > alt2 <- alt ### para manipular uma cópia do objeto > > z <- alt2@data$band1 ### para manipular apenas o vetor de interesse > > n <- length(z); n > > > > for(i in 1:n){if (z[i]<400){z[i]=0} else{z[i]=1}}; head(z) ### 0 ou 1 > > alt2@data$band1 <- z ### substitui o vetor manipulado no objeto > > image(alt2) ## visualização > > > > ### Opção para mascarar valores abaixo de critério (com 0 ou NA) > > > > alt3 <- alt ### para manipular outra cópia do objeto > > z <- alt3@data$band1 ### para manipular apenas o vetor de interesse > > n <- length(z); n > > > > for(i in 1:n){if (z[i]<400){z[i]=0}}; head(z) ### 0 ou 1 > > #for(i in 1:n){if (z[i]<400){z[i]=NA}}; head(z) ### 0 ou 1 > > alt3@data$band1 <- z ### substitui o vetor manipulado no objeto > > image(alt3) ## visualização > > ##### <end> ##### > > > > > > ================================================ > > Éder Comunello > > > > Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) > > Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] > > > > Researcher at Embrapa Western Region Agriculture > > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > > ================================================ > > UTC-03:00 > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > > > > > -- > _____________________________________ > *Carlos D'Apolito - PhD candidate - ID:1239480* > Earth Sciences > School of Geography & Environmental Sciences > Aston Webb Building > University of Birmingham > Birmingham > B15 2TT > UK > > phone: 44 - 07 919 673 173 > CXD280(a)bham.ac.uk<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=CXD280@bham.ac.uk> > CV: http://lattes.cnpq.br/4557754746424026 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/615af15d/atta… > > > > ------------------------------ > > Message: 3 > Date: Tue, 20 Nov 2012 14:14:30 -0200 > From: kaue veras <kaueveras(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kaueveras@gmail.com> > > > To: R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@listas.c3sl.ufpr.br> > Subject: [R-br] Graficos de decomposição > Message-ID: > <CAENvtmjrMnC52i4-JYn2NdhrwWMNw= > Qu0HJLeTOqz2s5q2i6uQ(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=Qu0HJLeTOqz2s5q2i6uQ@mai…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde prezados, > > Estou realizando um trabalho em que utilizo a decomposição, análises de > previsão usando o modelo ARIMA e Holt-Winters e goataria de mudar os > títulos do gráfico e titulos dos eixos. Com os comandos main, xlab e ylab > não esta dando certo, poderiam me dar uma luz ou me ajudar como faço para > gerar essas modificações? > > Desde já agradeço. > > Atenciosamente, > > Kauê P. Veras da Cunha > Estatística - UERJ > MSN: kaueveras(a)hotmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kaueveras@hotmail.com> > E-mail: kaueveras(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kaueveras@gmail.com> > Cel: (21) 8254-9601 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/1154d3e8/atta… > > > > ------------------------------ > > Message: 4 > Date: Tue, 20 Nov 2012 14:34:00 -0200 > From: Matheus Monteiro <matheusmont(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=matheusmont@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Subject: Re: [R-br] Gráfico de distribuiçoes estimadas > Message-ID: > < > CAAG8985KZaYb24ckqaGFxeon85H3g18ny9gEtBc79DXK7p8mcg(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=CAAG8985KZaYb24ckqaGFxeo…> > > > Content-Type: text/plain; charset="iso-8859-1" > > mais alguem pode me ajudar? > > Em 10 de novembro de 2012 13:21, Edson Lira <edinhoestat(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=edinhoestat@yahoo.com.br> > >escreveu: > > > Mateus, no pacote MASS acho que tem alguns gráficos empírica x teórica. > > > > [. ]'s. > > Edson Lira > > Estatístico > > Ma-Am > > > > Em 10/11/2012, às 02:20, Matheus Monteiro <matheusmont(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=matheusmont@gmail.com> > > > > escreveu: > > > > Boa noite, obtive as densidades aproximadas da 't' e da 'normal', > gostaria > > de compara-las através de graficos dessas densidades estimadas. > > Alguma sugestão, como faria isso? > > > > *t * > > > > > fitdistr(x[,1],"t") > > > > m s df > > > > -0.001874414 0.009122738 3.522064232 > > > > ( 0.001726835) ( 0.001862692) ( 2.234067399)* > > > > **normal* > > > > fitdistr(x[,1],"normal") > > > > mean sd > > > > 0.0001040816 0.0128201843 > > > > (0.0018314549) (0.0012950342) > > > > > > x são dados simulados. > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/1c03edac/atta… > > > > ------------------------------ > > Message: 5 > Date: Tue, 20 Nov 2012 09:06:00 -0800 (PST) > From: Edson Lira <edinhoestat(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=edinhoestat@yahoo.com.br> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: Re: [R-br] (sem assunto) > Message-ID: > <1353431160.32913.YahooMailNeo(a)web161205.mail.bf1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353431160.32913.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Mas Eder a idéia do NA é não considerar estas respostas na contagem e/ou > no percentual. > > Mas valeu e muito obrigado. > > Peço desculpas a todos, esqueci de colocar o assunto na mensagem > > Edson Lira > Estatístico > Manaus-Amazonas > > > ________________________________ > De: Eder Comunello <ecomunel(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ecomunel@gmail.com> > > > Para: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Enviadas: Terça-feira, 20 de Novembro de 2012 7:39 > Assunto: Re: [R-br] (sem assunto) > > > Edson e colegas, bom dia! > > Ao invés de deixar o campo 'motivo' em branco (NA), experimente atribuir > um valor, p. ex. 'sem resposta', e rode novamente. Assim estes casos serão > considerados e aparecerão na saída. Vale pros outros campos. > > ================================================ > Éder Comunello > > Ph.D. Student in Agricultural Systems Engineering (USP/ESALQ) > Piracicaba, SP, Brazil [22 42.7'S, 47 37.8'W] > > Researcher at Embrapa Western Region Agriculture > Dourados, MS, Brazil [22 16.5'S, 54 49.0'W] > ================================================ > UTC-03:00 > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/a9ea6c58/atta… > > > > ------------------------------ > > Message: 6 > Date: Tue, 20 Nov 2012 15:30:17 -0200 > From: Leonard de Assis <assis.leonard(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=assis.leonard@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Subject: Re: [R-br] Usando o R em grandes bases de dados > Message-ID: <50ABBE29.4030805(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=50ABBE29.4030805@gmail.c…> > > > Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" > > Bigmemory eu acredito que lide somente com matrizes, não com data > frames, resumo: todos os campos devem ter o mesmo tipo. > > []s > Leonard de Assis > http://about.me/ldeassis > > Em 20/11/2012 03:06, Vinicius Brito Rocha escreveu: > > > > fiquei brincando com o bigmemory e carreguei uma tabela mas os campos > > character vieram como NA. > > > > Quando usava a opção type="char" piorava e alterava os valores dos > > campos numéricos. O melhor ( ou menos pior resultado) foi carregar a > > tabela com a opção type="double" > > > > Alguem sabe o que pode ser? > > > > > > Abs. > > > > Vinicius Brito Rocha > > > > Em 19 de novembro de 2012 15:40, Richard Santos > > <jamesrichardsantos(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=jamesrichardsantos@gmail…><mailto: > jamesrichardsantos(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=jamesrichardsantos@gmail…> > >> > > escreveu: > > > > Valeu!! Toda e qq ajuda nesse sentido é bem vinda! > > Abç. > > Richard > > > > Em 19/11/12, Daniel C Bezerra<danielcbezerra(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=danielcbezerra@gmail.com> > > <mailto:danielcbezerra@gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=danielcbezerra@gmail.com>>> > escreveu: > > > Eu encontrei estas apresentações e achei legal dividir com vcs: > > > > > > http://goo.gl/l6fuC > > > > > > Tem link para os slides, código do R e tudo mais por lá. > > > > > > Abs, > > > > > > Daniel > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@listas.c3sl.ufpr.br><mailto: > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > > > > > > > > > -- > > /Vinicius Brito Rocha./ > > /Estatístico e Atuário (IM / UFRJ)// > > Mestre em Pesquisa Operacional (COPPE / UFRJ)/ > > > > www.aplicademic.blogspot.com <http://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. > > > > > > > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/b8d49acf/atta… > > > > ------------------------------ > > Message: 7 > Date: Tue, 20 Nov 2012 17:51:25 +0000 > From: alanarocha(a)sapo.pt<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alanarocha@sapo.pt> > To: R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@listas.c3sl.ufpr.br> > Subject: [R-br] matrizes de correlação > Message-ID: > <20121120175125.Horde.X_LLT3CXF3VQq8MdQ0ASvYA(a)mail.sapo.pt<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=20121120175125.Horde.X_L…> > > > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > Olá, a actualização do R resultou para o meu problema. > agora preciso de ajuda pois tou a tentar fazer uma matriz de > correlação e outra de covariancias de 4 variaveis de um dataframe e > dame estes erros: > ##matriz de correlação > > cor.p<-cor.ps(dados[,2:4], Pearson=TRUE,NAs=TRUE) > cor.p > dados=na.omit(dados[,2:4]) > > cat("valores NA retirados\n") > x11() > pairs(dados) > x11() > boxplot(dados) > dados=dados[,2:4] > x11() > boxplot(dados) > correl=cor(dados[2:4]) > result=correl > return(list(sumario=summary(dados),resultado=result)) > cor(cbind(X1s09, X1s10, X1s11, X1s12)) > cor.test(dados[,2:4], method='pearson')$p.value > cor_with_p_test(cbind(X1s09, X1s10, X1s11, X1s12)) > ct <- cor.test(X1s10, X1s12) > ct$p > ct$v > > #matriz de covariancias > cov(dados[,2:4], na.rm=TRUE) > > > > cor.p<-cor.ps(dados[,2:4], Pearson=TRUE,NAs=TRUE) > Error: could not find function "cor.ps" > > cor.p > Error: object 'cor.p' not found > > dados=na.omit(dados[,2:4]) > > > > cat("valores NA retirados\n") > valores NA retirados > > > > > > > > > > correl=cor(dados[2:4]) > Error in `[.data.frame`(dados, 2:4) : undefined columns selected > > result=correl > Error: object 'correl' not found > > return(list(sumario=summary(dados),resultado=result)) > Error: object 'result' not found > > cor(cbind(X1s09, X1s10, X1s11, X1s12)) > Error in cbind(X1s09, X1s10, X1s11, X1s12) : object 'X1s09' not found > > cor.test(dados[,2:4], method='pearson')$p.value > Error in `[.data.frame`(dados, , 2:4) : undefined columns selected > > cor_with_p_test(cbind(X1s09, X1s10, X1s11, X1s12)) > Error: could not find function "cor_with_p_test" > > ct <- cor.test(X1s10, X1s12) > Error in cor.test(X1s10, X1s12) : object 'X1s10' not found > > ct$p > Error: object 'ct' not found > > ct$v > Error: object 'ct' not found > > > > #matriz de covariancias > > cov(dados[,2:4], na.rm=TRUE) > Error in cov(dados[, 2:4], na.rm = TRUE) : > unused argument(s) (na.rm = TRUE) > > > cumprimentos e obrigada Ana Rocha > > > > > > ------------------------------ > > Message: 8 > Date: Tue, 20 Nov 2012 09:57:15 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: [R-br] ExpDes - Comparações múltiplas > Message-ID: > <1353434235.56527.YahooMailNeo(a)web125902.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353434235.56527.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde, > estou aqui tentando reproduzir algumas análises do ExpDes. A analise > abaixo não me deu as comparações múltiplas como faço para tais comparações? > > require(ExpDes) > Carregando pacotes exigidos: ExpDes > data(ex1) > attach(ex1) > dic(trat, ig, quali = FALSE, sigF = 0.05) > > > > André Oliveira Souza > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/696d59ab/atta… > > > > ------------------------------ > > Message: 9 > Date: Tue, 20 Nov 2012 16:24:10 -0200 > From: Walmes Zeviani <walmeszeviani(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmeszeviani@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>, > Andre Oliveira > <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Subject: Re: [R-br] ExpDes - Comparações múltiplas > Message-ID: > <CAFU= > EkZQNQHMTUM2nUKMnJhKTcCoDJ7ngqJ+_69-XZ6mf3FJDQ(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=EkZQNQHMTUM2nUKMnJhKTcCo…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Você usou quali=FALSE, ou seja, está declarando que os níveis são não > qualitativos. Use quali=TRUE e leia o manual. > > À disposiçã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<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br> > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/6b5bc51b/atta… > > > > ------------------------------ > > Message: 10 > Date: Tue, 20 Nov 2012 16:28:49 -0200 > From: Walmes Zeviani <walmeszeviani(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmeszeviani@gmail.com> > > > To: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > Subject: Re: [R-br] Dúvida sobre IC simultâneos para modelagem de > variáveis binarias no R (multcomp) > Message-ID: > <CAFU=EkYKvcE65FEpjv=qb0xnJ=mME7ec_8c7DD28wQNPF= > 1_yg(a)mail.gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1_yg@mail.gmail.com> > > > Content-Type: text/plain; charset="iso-8859-1" > > Seu exemplo não é reproduzível. O guia de postagem não recomenda envio de > dados como anexo das mensagens. Hospede os dados na núvem e acrescente > comandos para leitura do arquivo. Assim rapidamente podemos reproduzir sua > análise e ajudar você. > > À disposiçã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<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br> > skype: walmeszeviani > twitter: @walmeszeviani > homepage: http://www.leg.ufpr.br/~walmes > linux user number: 531218 > ========================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/7d4a7d34/atta… > > > > ------------------------------ > > Message: 11 > Date: Tue, 20 Nov 2012 10:30:37 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > To: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > >, > "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > Message-ID: > <1353436237.16303.YahooMailNeo(a)web125901.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353436237.16303.YahooMa…> > > > Content-Type: text/plain; charset="utf-8" > > Versão ExpDes_1.1.1.tar.gz agora sim consegui instalar no linux. > > > > André Oliveira Souza > > > ________________________________ > De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Para: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Enviadas: Quinta-feira, 5 de Janeiro de 2012 11:05 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > O Walmes me disse, aqui mesmo na lista, ha algum tempo que realmente o > pacote para linux em português ele não consegui colocar para funcionar. E > no caso do windows eu não sabia que estava em inglês! Obrigado > > > > André Oliveira Souza > > > ________________________________ > De: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Para: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Enviadas: Quinta-feira, 5 de Janeiro de 2012 10:55 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > Olá Andre, > > Que bom que não desistiu. Vamos por partes que solucionaremos todos os > seus problemas. > Quanto ao Windows, > não se esqueça que a versão do CRAN está em ingles. Portanto, após > carregar, voce deve pedir o help da função cdr (completely randomized > design) e não dic, ou seja > ?crd > Veja se funciona. > > > Quanto ao linux, talvez voce deva escrever para o Walmes, que faz parte > dessa lista. Ele usa o ExpDes no Linux e diz que funciona perfeitamente. > > abs > > Eric. > > > Em 5 de janeiro de 2012 10:42, Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br>> > escreveu: > > Eric, > > > >eu acompanhei o pacote no seu site e nunca consegui coloca-lo para > funcionar no Linux a versão em português. Agora com ele no CRAN voltei a > tentar colocar o pacote para funcionar no Linux, mas sem sucesso. E para > minha surpresa no windows não está funcionado também. Olha a minha > tentativa no windows! > >require(ExpDes) > >Carregando pacotes exigidos: ExpDes > >> ?dic > >No documentation for ?dic? in specified packages and libraries: > >you could try ???dic? > > > > > > Estou ansioso para colocar este pacote para funcionar por aqui. Já li o > TCC e o manual sobre o mesmo e vi a importância do mesmo. Visto que antes > tinha que usar muitos pacotes para a mesma análise e com o ExpDes é muito > mais prático. > > > > > > > > > > > > > > > > > > > > > >André Oliveira Souza > > > > > > > > > > > > > >________________________________ > > De: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > >Para: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>; > Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > >Enviadas: Quarta-feira, 4 de Janeiro de 2012 18:25 > >Assunto: Re: [R-br] Problema para reinstalar pacote (EspDes) no Linux > Ubuntu 10.04 > > > > > >Olá André, > > > >Se quiser trocar uma ideia em particular, posso tentar te ajudar. Sou o > criador do ExpDes, mas não usu Linux. Não sei se serei útil. > > > >A princípio, posso te dizer que existe o pacote no CRAN e na minha página > pessoal: https://sites.google.com/site/ericbferreira/unifal/downloads-1 > > > >Tente pegar o pacote de lá para ver se o erro se repete. > >Peço gentilmente aos usuárioas de (R no) Linux que me ajudem a exclarecer > o André. > > > >Atenciosamente, > > > >Eric. > > > > > >Em 4 de janeiro de 2012 16:39, Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br>> > escreveu: > > > >Opa vacilo meu erro no nome do pacote > >> > >> > >> install.packages("ExpDes", dep=T) > >> > >> > >> > >> > >> > >> > >>André Oliveira Souza > >> > >> > >>________________________________ > >> De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > >>Para: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > >>Enviadas: Quarta-feira, 4 de Janeiro de 2012 16:36 > >>Assunto: [R-br] Problema para reinstalar pacote (EspDes) no Linux Ubuntu > 10.04 > >> > >> > >> > >>Pessoal estou com problemas para reinstalar o EspDes !! Alguém saberia > dizer o que pode ser? > >> > >> > >> > >> > >>install.packages("EspDes", dep=T) > >>Installing package(s) into ?/home/andre/R/i486-pc-linux-gnu-library/2.14? > >>(as ?lib? is unspecified) > >>--- Please select a CRAN mirror for use in this session --- > >>Carregando a interface Tcl/Tk ... pronto > >>Mensagens de aviso perdidas: > >>In getDependencies(pkgs, dependencies, available, lib) : > >> package ?EspDes? is not available (for R version 2.14.1) > >> > >> > >> > >> > >> > >> > >> > >>André Oliveira Souza > >>obrigado > >> > >> > >>_______________________________________________ > >>R-br mailing list > >>R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > >> > >> > >>_______________________________________________ > >>R-br mailing list > >>R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > >> > > > > > >-- > >Dr Eric B Ferreira > >Exact Sciences Department > >Federal University of Alfenas > >Brazil > > > > > > > > > -- > Dr Eric B Ferreira > Exact Sciences Department > Federal University of Alfenas > Brazil > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/20c67bb9/atta… > > > > ------------------------------ > > Message: 12 > Date: Tue, 20 Nov 2012 10:37:54 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>>, > Eric > Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Subject: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > Message-ID: > <1353436674.31607.YahooMailNeo(a)web125904.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353436674.31607.YahooMa…> > > > Content-Type: text/plain; charset="utf-8" > > Em português no linux! > > André Oliveira Souza > > > ________________________________ > De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Para: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com>>; > "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Enviadas: Terça-feira, 20 de Novembro de 2012 16:30 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > Versão ExpDes_1.1.1.tar.gz agora sim consegui instalar no linux. > > > > André Oliveira Souza > > > ________________________________ > De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Para: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Enviadas: Quinta-feira, 5 de Janeiro de 2012 11:05 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > O Walmes me disse, aqui mesmo na lista, ha algum tempo que realmente o > pacote para linux em português ele não consegui colocar para funcionar. E > no caso do windows eu não sabia que estava em inglês! Obrigado > > > > André Oliveira Souza > > > ________________________________ > De: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > Para: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > Enviadas: Quinta-feira, 5 de Janeiro de 2012 10:55 > Assunto: Re: [R-br] Problema para reinstalar pacote (ExpDes) no Linux > Ubuntu 10.04 > > > Olá Andre, > > Que bom que não desistiu. Vamos por partes que solucionaremos todos os > seus problemas. > Quanto ao Windows, > não se esqueça que a versão do CRAN está em ingles. Portanto, após > carregar, voce deve pedir o help da função cdr (completely randomized > design) e não dic, ou seja > ?crd > Veja se funciona. > > > Quanto ao linux, talvez voce deva escrever para o Walmes, que faz parte > dessa lista. Ele usa o ExpDes no Linux e diz que funciona perfeitamente. > > abs > > Eric. > > > Em 5 de janeiro de 2012 10:42, Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br>> > escreveu: > > Eric, > > > >eu acompanhei o pacote no seu site e nunca consegui coloca-lo para > funcionar no Linux a versão em português. Agora com ele no CRAN voltei a > tentar colocar o pacote para funcionar no Linux, mas sem sucesso. E para > minha surpresa no windows não está funcionado também. Olha a minha > tentativa no windows! > >require(ExpDes) > >Carregando pacotes exigidos: ExpDes > >> ?dic > >No documentation for ?dic? in specified packages and libraries: > >you could try ???dic? > > > > > > Estou ansioso para colocar este pacote para funcionar por aqui. Já li o > TCC e o manual sobre o mesmo e vi a importância do mesmo. Visto que antes > tinha que usar muitos pacotes para a mesma análise e com o ExpDes é muito > mais prático. > > > > > > > > > > > > > > > > > > > > > >André Oliveira Souza > > > > > > > > > > > > > >________________________________ > > De: Eric Ferreira <ericbferreira(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=ericbferreira@gmail.com> > > > >Para: r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>; > Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > >Enviadas: Quarta-feira, 4 de Janeiro de 2012 18:25 > >Assunto: Re: [R-br] Problema para reinstalar pacote (EspDes) no Linux > Ubuntu 10.04 > > > > > >Olá André, > > > >Se quiser trocar uma ideia em particular, posso tentar te ajudar. Sou o > criador do ExpDes, mas não usu Linux. Não sei se serei útil. > > > >A princípio, posso te dizer que existe o pacote no CRAN e na minha página > pessoal: https://sites.google.com/site/ericbferreira/unifal/downloads-1 > > > >Tente pegar o pacote de lá para ver se o erro se repete. > >Peço gentilmente aos usuárioas de (R no) Linux que me ajudem a exclarecer > o André. > > > >Atenciosamente, > > > >Eric. > > > > > >Em 4 de janeiro de 2012 16:39, Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br>> > escreveu: > > > >Opa vacilo meu erro no nome do pacote > >> > >> > >> install.packages("ExpDes", dep=T) > >> > >> > >> > >> > >> > >> > >>André Oliveira Souza > >> > >> > >>________________________________ > >> De: Andre Oliveira <andreolsouza(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=andreolsouza@yahoo.com.br> > > > >>Para: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > >>Enviadas: Quarta-feira, 4 de Janeiro de 2012 16:36 > >>Assunto: [R-br] Problema para reinstalar pacote (EspDes) no Linux Ubuntu > 10.04 > >> > >> > >> > >>Pessoal estou com problemas para reinstalar o EspDes !! Alguém saberia > dizer o que pode ser? > >> > >> > >> > >> > >>install.packages("EspDes", dep=T) > >>Installing package(s) into ?/home/andre/R/i486-pc-linux-gnu-library/2.14? > >>(as ?lib? is unspecified) > >>--- Please select a CRAN mirror for use in this session --- > >>Carregando a interface Tcl/Tk ... pronto > >>Mensagens de aviso perdidas: > >>In getDependencies(pkgs, dependencies, available, lib) : > >> package ?EspDes? is not available (for R version 2.14.1) > >> > >> > >> > >> > >> > >> > >> > >>André Oliveira Souza > >>obrigado > >> > >> > >>_______________________________________________ > >>R-br mailing list > >>R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > >> > >> > >>_______________________________________________ > >>R-br mailing list > >>R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > >> > > > > > >-- > >Dr Eric B Ferreira > >Exact Sciences Department > >Federal University of Alfenas > >Brazil > > > > > > > > > -- > Dr Eric B Ferreira > Exact Sciences Department > Federal University of Alfenas > Brazil > > > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/6a5b8ac7/atta… > > > > ------------------------------ > > Message: 13 > Date: Tue, 20 Nov 2012 10:49:58 -0800 (PST) > From: Alexandre Santos <alexandresantosbr(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alexandresantosbr@yahoo.…> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: [R-br] Dúvida em summary em regressão logística > Message-ID: > <1353437398.81186.YahooMailNeo(a)web121206.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353437398.81186.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde Pessoal, > Estou ajustando uma regressão logística e me deparei com a seguinte > dúvida: > > > ################################################################################################### > #Regressão logística entre a proporção de herbívoros predados e a > densidade de herbívoros oferecidos > # > pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2 > ,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26 > ,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros > predados > dens<-sort(rep(2^(2:6),10))#### Densidade de herbivoros oferecidos > ## > ##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ > exp(P0+P1*N+P2*N^2+P3*N3) > p.model1<-glm(pred/par~dens+I(dens^2)+I(dens^3),family="binomial") > summary(p.model1) > # > > > Minha dúvida é se o coeficiente linear que aparece no summary esta > transformado em exp(x)/1+exp(x) > ou trata-se do valor sem transformação? > > Obrigado, > > -- > ====================================================================== > Alexandre dos Santos > Proteção Florestal > IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso > - Campus Cáceres > Avenida dos Ramires, s/n > Bairro: Distrito Industrial > Cáceres - MT CEP: 78.200-000 > Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) > e-mails:alexandresantosbr@yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=e-mails%3Aalexandresanto…> > alexandre.santos(a)cas.ifmt.edu.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alexandre.santos@cas.ifm…>====================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/0b764a8d/atta… > > > > ------------------------------ > > Message: 14 > Date: Tue, 20 Nov 2012 10:51:20 -0800 (PST) > From: Alexandre Santos <alexandresantosbr(a)yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alexandresantosbr@yahoo.…> > > > To: "r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br>" > <r-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=r-br@listas.c3sl.ufpr.br> > > > Subject: [R-br] Dúvida em summary em regressão logística > Message-ID: > <1353437480.93196.YahooMailNeo(a)web121202.mail.ne1.yahoo.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=1353437480.93196.YahooMa…> > > > Content-Type: text/plain; charset="iso-8859-1" > > Boa tarde Pessoal, > Estou ajustando uma regressão logística e me deparei com a seguinte > dúvida: > > > ################################################################################################### > #Regressão logística entre a proporção de herbívoros predados e a > densidade de herbívoros oferecidos > # > pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2 > ,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26 > ,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros > predados > dens<-sort(rep(2^(2:6),10))#### Densidade de herbivoros oferecidos > ## > ##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ > exp(P0+P1*N+P2*N^2+P3*N3) > p.model1<-glm(pred/par~dens+I(dens^2)+I(dens^3),family="binomial") > summary(p.model1) > # > > > Minha dúvida é se o coeficiente linear que aparece no summary esta > transformado em exp(x)/1+exp(x) > ou trata-se do valor sem transformação? > > Obrigado, > > -- > ====================================================================== > Alexandre dos Santos > Proteção Florestal > IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso > - Campus Cáceres > Avenida dos Ramires, s/n > Bairro: Distrito Industrial > Cáceres - MT CEP: 78.200-000 > Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) > e-mails:alexandresantosbr@yahoo.com.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=e-mails%3Aalexandresanto…> > alexandre.santos(a)cas.ifmt.edu.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=alexandre.santos@cas.ifm…>====================================================================== > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/e2cfa427/atta… > > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=R-br@listas.c3sl.ufpr.br> > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 21, assunto 47 > ***************************************** >

1 0

Dúvida em summary em regressão logística
by Alexandre Santos 20 Nov '12

20 Nov '12

Boa tarde Pessoal, Estou ajustando uma regressão logística e me deparei com a seguinte dúvida: ################################################################################################### #Regressão logística entre a proporção de herbívoros predados e a densidade de herbívoros oferecidos # pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2 ,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26 ,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros predados dens<-sort(rep(2^(2:6),10))#### Densidade de herbivoros oferecidos ## ##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ exp(P0+P1*N+P2*N^2+P3*N3) p.model1<-glm(pred/par~dens+I(dens^2)+I(dens^3),family="binomial") summary(p.model1) # Minha dúvida é se o coeficiente linear que aparece no summary esta transformado em exp(x)/1+exp(x) ou trata-se do valor sem transformação? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso - Campus Cáceres Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

1 0