- R-br - DINF UFPR

Tinn-R 2.3.7.4 pre-release
by Jose Claudio Faria 01 Nov '12

01 Nov '12

Prezados usuários do Tinn-R, Aproveitando que ando me recuperando de uma cirurgia de hérnia inguinal dupla convencional em casa, com calma e bastante tempo (não estou indo trabalhar na Universidade), andei fazendo umas alterações no Tinn-R e creio que ficaram bem boas. Contudo, gostaria de escutar as opiniões do usuários da casa antes de liberar para todos os usuários no SourceForge. Estou usando a nova versão já a alguns dias e já testei em dois sabores do WIndows 7, tudo parece estável. Mas seria bom testes em outras plataformas e máquinas. Os uauários que quiserem testar podem entrar em contato direto comigo: joseclaudio.faria(a)gmail.com que mando o link de onde pode ser baixado. As principais alterações estão abaixo (em Inglês): 2.3.7.4 (Nov/01/2012) - The graphical interface was updated with some improvements, mainly the Application options. - A new toolbar "Format" was added to the main toolbar. - A new resource allowing reformat R code (selection or whole file) using formatR package was added. The icon resource was placed in the Format toolbar, so that from this version on, the formatR package will be necessary together with the already traditional TinnR and svSocket packages. - Due to new resource related to reformat code the variable .trPaths was changed. As a result, it will be necessary to run again the R/Configure/Permanent (Rprofile.site). In this case, do not forget to remove any prior script generated by Tinn-R in the Rprofile.site file. - The R/Hotkeys interface was deeply reworked and it now has two tabs: Default and Custom. - Default: Contains the already traditional instructions of Tinn-R; - Custom: Allows the user to customize any instructions to be send to R interpreter (thanks to Philemon Lenherr for the suggestion). The instructions must be as follows: - Simple: search(). The R interpreter will receive > search(); - Replace word or small selection: View(%s, title='Test'). If the editor cursor is over the word iris or it is selected, the R interpreter will receive > View(iris, title='Test') - Replace whole file: source(%f, echo=TRUE, verbose=TRUE). The R interpreter will receive > source(.trPaths[4], echo=TRUE, verbose=TRUE). All rules related to send file will be preserved. - Sorry, due to deeply changes in R/Hotkeys, all hotkeys configured prior to this version will be lost. It will be necessary to reconfigure all. - The chapter Some secrets for an efficient use is being revised by Ricardo Pietrobon and soon will be completed. Many thanks for his hard work! Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

Problemas na integração R e Tinn-R!
by Danilo Cuzzuol Pedrini 01 Nov '12

01 Nov '12

Prezados, Instalei as últimas versões do R e do Tinn-R, ao enviar uma linha de comando do Tinn-R ao R, aparece a seguinte mensagem de erro: "source(.trPaths[4]) Erro em source(.trPaths[4]) : objeto '.trPaths' não encontrado" Atenciosamente, Danilo Cuzzuol Pedrini www.danilocp.blogspot.com

2 1

[Dúvida] Função Gama Incompleta.
by Pedro Rafael 01 Nov '12

01 Nov '12

Pessoal qual o comando certo ou pacote que trabalha com a função gamma incompleta? Estou obtendo valores muito diferentes da gamma incompleta no R e no mathematica. Por exemplo: No R: pgamma(.10,1.5) = 0.0224107 No Mathematica: N[Gamma[.10, 1.5], 10] = 0.107921 Qual seria o problema que não estou enxergando? [ ], Pedro Rafael Diniz Marinho.

2 1

Problema ao ler arquivo UTF-16 com nulos
by Roberto de Pinho 01 Nov '12

01 Nov '12

Caros, Estou tentando ler o arquivo de servidores disponível no portal da transparência. ( http://arquivos.portaldatransparencia.gov.br/PortalTransparenciaEscolheTipo… - mais de 200MB) É um CSV codificado em UTF-16. No entanto, ele reconhece prematuramente fins de linhas em algumas linhas quando encontra campos com nulo [00 00]. Há valores não nulos não sendo lidos após campos nulos que provocam o reconhecimento de fim de linha. Alguém saberia como poderia fazer com que o R possa ler este arquivo? Estou lendo o arquivo com: servidores = read.table("teste.csv", header = T, sep="\t",quote="",stringsAsFactors=T,fileEncoding="UTF-16", fill=T) o fill=T estou sando para evitar que dê erro de diferente numero de campos. Envio, em anexo, um extrato das primeiras linhas do arquivo. Obrigado pela atenção. ==== Sequencia de bytes que provocam o reconhecimento de nova linha: 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 00 00 09 00 -- Roberto de Pinho robertodepinho(a)gmail.com http://www.ascoisas.com http://lattes.cnpq.br/4816166073408660

7 20

LER CSV N0 MAC OS
by Edson Lira 31 Out '12

31 Out '12

Caros amigos da lista, meu pc é MacBookPro, OS X vefrsão 10.6.8, processador I7, com 2,7 HZ, ram de 4gb. Ao ler arquivos no formato csv, tenho tido dificuldades, no último que estou tentando ler, inclusive já postei em um email anterior uma amostra do arquivo. Quais os cuidados mínimos que devo tomar para não ter problemas? Com a rotina abaixo, h<-read.csv("/Users/edsonlira/Edinho/Edinho/Consulta/Ciesa/BSC_2307.csv",row.names = 1) Está me dando o erro abaixo: Erro em read.table(file = file, header = header, sep = sep, quote = quote, : more columns than column names Edson Lira Estatístico Manaus-Amazonas

3 2

Re: [R-br] em plot.xy(xy.coords(x, y), type = type, ...) : objeto 'membros' não encontrado
by Eder Comunello 31 Out '12

31 Out '12

Kin, boa tarde! O erro listado é por atribuir ao parâmetro pch o objeto 'membros' que não existe ou não está acessível (veja seu arquivo de dados e declarações). Acredito que a ideia seria diferenciar os pontos (tipos) de acordo com o valor em membros... Acredito que o exemplo que segue faça algo parecido com o que você deseja, embora usando somente a biblioteca sp(). Inclui os pontos que você forneceu. ########### require(sp) br.df <- read.table("http://dl.dropbox.com/u/117618178/brlite.dat", sep=";", head=T) class(br.df); head(br.df); str(br.df) br.lin <- Lines(list(Line(br.df)), ID = "Brasil") br.sl <- SpatialLines(list(br.lin)) proj4string(br.sl) <- CRS('+proj=longlat +datum=WGS84') summary(br.sl); class(br.sl) est.x <- c(-53.35,-55.60,-51.80,-51.05,-50.13,-52.41,-52.10,-52.51) est.y <- c(-33.51,-30.83,-29.45,-29.68,-29.05,-31.78,-32.03,-30.53) est.id <- c(1:8) est.df <- data.frame(est.x,est.y,est.id) est.sp <- est.df coordinates(est.sp) <- c('est.x','est.y') proj4string(br.sl) <- CRS('+proj=longlat +datum=WGS84') est.sp #png(filename="c:/temp/mapa.png",width=600,height=800) plot(br.sl, axes=T) plot(est.sp, pch=21, bg=3, add=T) #dev.off() ########### -- ================================================ É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 > Date: Tue, 30 Oct 2012 21:09:26 -0200 > From: kin Huts <kim-met(a)hotmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=kim-met@hotmail.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] em plot.xy(xy.coords(x, y), type = type, ...) : > objeto 'membros' não encontrado > Message-ID: <SNT141-W57AFF1FD679C8C33FDDE18F1620(a)phx.gbl> > Content-Type: text/plain; charset="iso-8859-1" > > > > Olá, > Gostaria de uma ajuda para entender o erro encontrado e como corrigir esse > script,Estou tentando gerar um mapa para o RS com a localização das > estações que estou utilizando. > > > require(MASS)require(maps) > > # GERAÇÃO DO MAPA COM AS ESTAÇÕES > # carregando divisão política brasileiramapaBrasil <- > read.table("BASE_BR.txt",sep="") > # carregando informações de localização das estações > meteorológicaslatlonRS <- read.csv2("coordenadas.txt",h=F,sep="") > # gerando mapa em PNGpng(filename="mapa_grupos.png",width=600,height=800) > # criando grafico em branco p/ o estado do > RSmap(xlim=c(-80,-35),ylim=c(-45,10), type="n")map.axes() > # plotando mapalines(mapaBrasil[,2], mapaBrasil[,1]) > # plotando estaçõespoints( latlonRS[,2], latlonRS[,1], pch=membros) > # fechando PNGdev.off() > > dados que estou usando são: V1 V21 -33.51 -53.352 -30.83 > -55.63 -29.45 -51.84 -29.68 -51.055 -29.05 -50.136 -31.78 -52.417 -32.03 > -52.18 -30.53 -52.51 > e o erro que aparece é: > > Carregando pacotes exigidos: MASSCarregando pacotes exigidos: mapsErro em > plot.xy(xy.coords(x, y), type = type, ...) : objeto 'membros' não > encontrado > >

1 0

Re: [R-br] coordenadas geograficas
by Eder Comunello 31 Out '12

31 Out '12

Natalia e colegas, bom dia! Ainda não estou muito familiarizado com o R, mas acho que este exemplo que rodei irá ajudar... Não utilizei o project(), mas sim o spTransform()... ############################################################# #install.packages("rgdal", dep=T) #install.packages("sp", dep=T) require(sp); require(rgdal) ### long coord.x = c(-54.607845, -54.599831, -54.430241, -54.258459, -54.344408, -54.333258, -53.607359, -53.589245, -53.546762, -53.105123, -52.972362, -52.107579, -52.089572, -51.940453, -51.718151, -51.263728, -50.741531, -50.734130, -50.667779, -50.662230, -49.986102, -49.970018, -49.901474, -49.911543, -49.678260, -49.654858, -49.610676, -49.599658, -49.511570, -49.200709, -49.222508, -49.317621, -49.311839, -49.275882, -48.581791, -48.578971, -48.582604, -48.568999, -48.528527, -48.411012, -48.251622, -48.027824, -48.029886, -48.023501, -48.590796, -48.594522, -48.643601, -49.175230, -49.555033, -49.567207, -49.689784, -49.941462, -49.942167, -50.571750, -50.591636, -50.553035, -50.718453, -51.208345, -51.216170, -51.292329, -51.239319, -51.258597, -51.283372, -51.411480, -53.551638, -53.552753, -53.643791, -53.835953, -53.842011, -53.858677, -53.955544, -54.099576, -54.108952, -54.429627, -54.438295, -54.593660, -54.593683, -54.619090, -54.617907, -54.607845) ### lat coord.y = c(-25.446263, -25.434436, -25.159459, -24.357597, -24.140630, -24.116248, -22.951327, -22.936230, -22.890090, -22.621923, -22.570479, -22.516685, -22.521237, -22.567335, -22.669007, -22.669047, -22.961387, -22.962502, -22.897375, -22.895686, -22.897501, -22.910601, -22.999147, -23.051110, -23.165054, -23.210043, -23.851684, -23.857080, -23.925858, -24.344464, -24.362727, -24.566630, -24.663615, -24.693608, -24.671609, -24.685526, -25.051554, -25.051240, -25.100784, -24.980119, -24.978739, -25.221284, -25.226189, -25.230964, -25.976982, -25.977087, -25.955943, -26.001844, -26.237326, -26.233284, -26.188456, -26.008879, -26.008876, -26.003790, -26.022134, -26.055436, -26.244779, -26.297107, -26.313987, -26.404261, -26.628383, -26.625825, -26.653990, -26.717283, -26.292744, -26.292051, -26.253072, -25.971930, -25.687925, -25.680270, -25.646984, -25.495677, -25.495742, -25.695430, -25.691293, -25.592352, -25.544262, -25.457555, -25.452072, -25.446263) pr.df <- data.frame(coord.x,coord.y); head(pr.df) ### criar um objeto "SpatialPoints" pr.spgeo <- pr.df; head (pr.spgeo) coordinates(pr.spgeo) <- c('coord.x','coord.y') ### definir projeção e datum de entrada proj4string(pr.spgeo) <- CRS('+proj=longlat +datum=WGS84') class(pr.spgeo) head(pr.spgeo) plot(pr.spgeo) ### lista parâmetros disponíveis projInfo("proj") projInfo("datum") projInfo("ellps") ### transforma para utm ### atenção para o uso do parâmetro '+south' ### ter claro que o Paraná utiliza mais de uma Zona UTM e maior parte do estado em 22S pr.sputm <- spTransform(pr.spgeo, CRS('+proj=utm +zone=22 +south')); head(pr.sputm) plot(pr.sputm) ############################################################# Esse exemplo assumiu o WGS84 como datum padrão de entrada e saída. Para datums (data) diferentes tem que dar uma estudada mais a fundo, porque tem poucas opções pré-configuradas no pacote que usei (como pode ser visto usando <projInfo("datum")>). -- ================================================ É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 > ---------------------------------------------------------------------- > > > 2012/10/29 Natalia Martins <nsmbarreto(a)gmail.com<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=nsmbarreto@gmail.com> > > > > > Walmes, bom dia. > > Estou tentando transformar minhas coordenadas geograficas (das estações > > metereologicas do PR) > > em utm, procurando no help do R encontrei: > > > > data(state) > > s <- project(state.center, "+proj=merc") > > > > mas eu não estou conseguindo entender qual proj eu devo usar, no link > > > > http://www.remotesensing.org/geotiff/proj_list/transverse_mercator.html > > > > até encontri mas quando tento fazer aparece o seguinte erro > > > > s <- project(state.center, "+proj=utm") > > Erro em project(state.center, "+proj=utm") : elliptical usage required > > > > > > se vc puder me ajudar... > > > > Muito obrigada!! > > > > Em 28 de outubro de 2012 11:11, Walmes Zeviani <walmes(a)ufpr.br<https://mail.google.com/mail/?view=cm&fs=1&tf=1&to=walmes@ufpr.br>> > escreveu: > > > > Sim, para calculo de distâncias você tem que ter coordenadas decimais e > >> não exagesimais (º ' ''). > >> > >> > ========================================================================== > >> 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 > >> > ========================================================================== > >> > > > > > > > > -- > > > > 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 > > > > > > Mensagem Verificada Pelo AntiSpam UFPR > > > > >

2 1

Carregar dados pelo Rcmdr_ ubuntu 12.04
by Andre Oliveira 31 Out '12

31 Out '12

Boa noite, eu sempre uso R via terminal mesmo. Mas pensando em dar uma aula para o ADM pensei no Rcmdr. Porém estou com um problema não consigo carregar arquivos .txt, .xls e .csv. Vou em carregar arquivo seleciono o arquivo e tem o aviso. 1] NOTA: Versão do R Commander 1.8-4: Fri Oct 26 20:09:58 2012 [2] ERRO: número mágico do arquivo de restauração inválido (arquivo pode estar corrompido) -- nenhum dado carregado Alguma dica? André Oliveira Souza

4 10

Modelos econométricos no R
by Gilbert Queiroz 31 Out '12

31 Out '12

Olá, amigos e amigas, Alguém poderia me indicar alguma leitura sobre Modelos Econométricos no R??? Saudações. Gilbert.

4 4

Ajuda para plotar dados (usando o ggplot2)
by Thiago Veloso 30 Out '12

30 Out '12

Pessoal, Alguém domina bem a produção de gráficos usando o pacote ggplot2? Meus dados contém a frequência de dados em 13 classes diferentes ao longo de vários anos (simulação climática). Segue o link para baixar o arquivo contendo os dados: https://www.dropbox.com/s/50loxmtd2z2eqco/rcp26.csv (a primeira coluna pode ser ignorada). Gostaria de fazer algo parecido com esse gráfico: https://www.dropbox.com/s/8p3g42c54w3dnof/Untitled.png. Meus dados contém todas as informações para montar o gráfico: as classes (eixo x) e suas frequências (ao invés de hectares - eixo y) e os anos (três barras). Por enquanto, o código mínimo reprodutível que eu consigo enviar é apenas a abertura dos arquivos: a <- read.csv ('~/Dropbox/web/rcp26.csv') b <- read.csv ('~/Dropbox/web/rcp85.csv') Preciso de ajuda para reproduzir a imagem acima usando meus dados, com os anos 1700, 2005 e 2100 no lugar de 1985, 1993 e 2005. Agradeço desde já, Thiago.

3 5

em plot.xy(xy.coords(x, y), type = type, ...) : objeto 'membros' não encontrado
by kin Huts 30 Out '12

30 Out '12

Olá, Gostaria de uma ajuda para entender o erro encontrado e como corrigir esse script,Estou tentando gerar um mapa para o RS com a localização das estações que estou utilizando. require(MASS)require(maps) # GERAÇÃO DO MAPA COM AS ESTAÇÕES # carregando divisão política brasileiramapaBrasil <- read.table("BASE_BR.txt",sep="") # carregando informações de localização das estações meteorológicaslatlonRS <- read.csv2("coordenadas.txt",h=F,sep="") # gerando mapa em PNGpng(filename="mapa_grupos.png",width=600,height=800) # criando grafico em branco p/ o estado do RSmap(xlim=c(-80,-35),ylim=c(-45,10), type="n")map.axes() # plotando mapalines(mapaBrasil[,2], mapaBrasil[,1]) # plotando estaçõespoints( latlonRS[,2], latlonRS[,1], pch=membros) # fechando PNGdev.off() dados que estou usando são: V1 V21 -33.51 -53.352 -30.83 -55.63 -29.45 -51.84 -29.68 -51.055 -29.05 -50.136 -31.78 -52.417 -32.03 -52.18 -30.53 -52.51 e o erro que aparece é: Carregando pacotes exigidos: MASSCarregando pacotes exigidos: mapsErro em plot.xy(xy.coords(x, y), type = type, ...) : objeto 'membros' não encontrado

2 1

Comando para os Graficos
by T Branquinho 30 Out '12

30 Out '12

Bom dia, Para uma regressão multipla linear qualquer: > set.seed(1) > B1 <- data.frame(x1=rpois(16,1),x2=rpois(16,2),x3=rpois(16,4),y=rpois(16,10)) > lm.B1=lm(y~x1+x2+x3) > lm.B1=lm(y~x1+x2+x3,data=B1) > plot(lm.B1) Aperte <Enter> para ver o próximo gráfico: Resisuals versus Fitted Aperte <Enter> para ver o próximo gráfico: Normal Q-Q Aperte <Enter> para ver o próximo gráfico: Scale Location Aperte <Enter> para ver o próximo gráfico: Residuals versus Leverage Problema: Possuo um modelo não linear e gostaria de obter os mesmos gráficos para uma apresentação e para comparação dos dados, desta forma quais os ->comandos<- e dados necessários para plotar esses gráficos para qualquer modelo de predição? No mais agradeço, Tales

2 3

Re: [R-br] [R_STAT] Funções 'sort' e 'order' no R....
by Ary Casado 30 Out '12

30 Out '12

Olá A lista mudou para o endereço r-br(a)listas.c3sl.ufpr.br Algumas considerações: 1) sort sort faz uma classificação (crescente ou decrescente) simples e tem a sintaxe: sort(x, decreasing = FALSE, ...) onde x é o vetor a ser classificado, decreasing = false define que a classificação será em ordem crescente e ... significa que outros argumentos podem ser incrementados. Por exemplo: na.last = NA significa que valores que não foram definidos não vão ser listados. Outras opções seriam na.last = TRUE (os valores faltantes ficariam após o último valor classificado) e na.last=FALSE (os valores faltantes ficariam antes do primeiro valor classificado) Tente os exemplos abaixo > amostra <- c(4, 10, 3, -2, NA, 0, NA, 30, -5) > sort(amostra) #É a classificação simples com os argumentos com seus valores padrões [1] -5 -2 0 3 4 10 30 > sort(amostra, decreasing=TRUE, na.last=TRUE) [1] 30 10 4 3 0 -2 -5 NA NA 2) order A função order faz uma permutação que rearranja os valores em ordem crescente ou decrescente. A sintaxe da função order é order(..., na.last = TRUE, decreasing = FALSE) onde ... é qualquer sequência válida. Talvez fique mais fácil você colocar um exemplo com caracteres de texto: > letras1 <-c('g', 'a', 'z', 'c') > order(letras1) [1] 2 4 1 3 Observe que o vetor tem as letras fora da ordem, então a função order as ordena mostrando em qual posição está cada uma: a está na posição 2 c está na posição 4 g está na posição 1 z está na posição 3 Se você utiliza o operador [ ] você seleciona um elemento através da posição dele na sequência ou no vetor. então fica claro se você fizer > letras1[order(letras1)] o resultado será: [1] "a" "c" "g" "z" pois os colchetes determinam qual a posição do elemento que será retornado da lista. quem está na posição 2? a quem está na posição 4? c quem está na posição 1? g quem está na posição 3? z NÃO SEI SE DÁ PRA FAZER A FUNÇÃO SORT SER MODIFICADA PARA FAZER O PAPEL DE ORDER. Espero ter ajudado. Ary Nascimento Engenheiro Químico Consultor Especialista em Informática Mestrando em Ciências Nucleares - Proten - UFPE ________________________________ De: relatividade <alissonpereira_14(a)hotmail.com> Para: R_STAT(a)yahoogrupos.com.br Enviadas: Terça-feira, 30 de Outubro de 2012 8:00 Assunto: [R_STAT] Funções 'sort' e 'order' no R.... Olá amigos tudo bem ? Espero que sim. Sou novo aqui e gostaria muito da ajuda dos senhores para resolver um pequeno probleminha no R. Então, estou com um probleminha com uma atividade que relaciona duas funções no R, é o seguinte: Primeiro através da função sample é criado um vetor aleatótrio: x <- sample(1:100, 20) Depois aplico a função sort que tem como objetivo ordenar os elementos em ordem crescente sort(x) Ai o problema aplica a função order em x: order(x) ---------------------------------------------------------- Ai o problema cita duas questões: 1) Qual é o objetivo da função order ? 2) Como emular, ou seja, como fazer que a função order realize o mesmo objetivo que a função sort mas sem utilizar a função sort. E vice-versa ? Para a questão 2 obtive a ajuda de um colega onde este me deu à dica de fazer a primeira parte da questão da seguinte maneira: Para fazer com que a função order execute o mesmo objetivo que a função sort, faz-se o seguinte: x[order(x)] *Se plotarmos este comando no R realmente o objetivo da primeira pate da questão é cumprido, porém não entendi o pq desses dois colchetes ? ** E como fazer com que a função sort realize o objetivo da função order mas sem utilizar o order ?? Agradeço muito mesmo se puderem me ajudar na questão 1 e nas questões acima (*) e (**) __._,_.___ | através de email | Responder através da web | Adicionar um novo tópico Mensagens neste tópico (1) Atividade nos últimos dias: * Novos usuários 3 Visite seu Grupo Não passe mais um Dia dos Namorados sozinho(a). Conheça o Y! Encontros.<img width="1" height="1" alt=""> ________________________________ O Dia dos Namorados está chegando. Encontre alguém especial até lá!<img width="1" height="1" alt=""> ________________________________ Super Descontos HP <img width="1" height="1" alt=""> Trocar para: SÃ³ Texto, Resenha DiÃ¡ria • Sair do grupo • Termos de uso . __,_._,___

1 0

Soundex para português
by Fátima Lima Paula 29 Out '12

29 Out '12

Prezados, alguém sabe como conseguir um algoritmo Soundex para Português para utilizar no RecordLinkage? Obrigada, Fátima

2 1

ERRO ao instalar pacote rJava
by Emerson Cotta Bodevan 29 Out '12

29 Out '12

Prezados, minha intenção é instalar o pacote speedR. No entanto, ao tentar instalá-lo, aparece a mensagem de que está faltando o rJava. Bom... ao tentar instalar o rJava aparece a mensagem abaixo: > install.packages() Carregando a interface Tcl/Tk ... pronto --- Please select a CRAN mirror for use in this session --- Installing package(s) into ‘/home/ecbodevan/R/i686-pc-linux-gnu-library/2.14’ (as ‘lib’ is unspecified) tentando a URL 'http://cran-r.c3sl.ufpr.br/src/contrib/rJava_0.9-3.tar.gz' Content type 'application/x-gzip' length 537153 bytes (524 Kb) URL aberta ================================================== downloaded 524 Kb * installing *source* package ‘rJava’ ... ** package ‘rJava’ successfully unpacked and MD5 sums checked checking for gcc... gcc -std=gnu99 checking for C compiler default output file name... a.out checking whether the C compiler works... yes checking whether we are cross compiling... no checking for suffix of executables... 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 grep that handles long lines and -e... /bin/grep checking for egrep... /bin/grep -E checking for ANSI C header files... yes checking for sys/wait.h that is POSIX.1 compatible... yes checking for sys/types.h... yes checking for sys/stat.h... yes checking for stdlib.h... yes checking for string.h... yes checking for memory.h... yes checking for strings.h... yes checking for inttypes.h... yes checking for stdint.h... yes checking for unistd.h... yes checking for string.h... (cached) yes checking sys/time.h usability... yes checking sys/time.h presence... yes checking for sys/time.h... yes checking for unistd.h... (cached) yes checking for an ANSI C-conforming const... yes checking whether time.h and sys/time.h may both be included... yes configure: checking whether gcc -std=gnu99 supports static inline... yes checking whether setjmp.h is POSIX.1 compatible... yes checking whether sigsetjmp is declared... yes checking whether siglongjmp is declared... yes checking Java support in R... present: interpreter : '/usr/bin/java' archiver : '' compiler : '' header prep.: '' cpp flags : '' java libs : '-L/usr/lib/jvm/java-6-openjdk-i386/jre/lib/i386/server -L/usr/lib/jvm/java-6-openjdk-i386/jre/lib/i386 -L/usr/lib/jvm/java-6-openjdk-i386/jre/../lib/i386 -L/usr/java/packages/lib/i386 -L/usr/lib/i386-linux-gnu/jni -L/lib/i386-linux-gnu -L/usr/lib/i386-linux-gnu -L/usr/lib/jni -L/lib -L/usr/lib -ljvm' *configure: error: Java Development Kit (JDK) is missing or not registered in R Make sure R is configured with full Java support (including JDK). Run R CMD javareconf as root to add Java support to R.* If you don't have root privileges, run R CMD javareconf -e to set all Java-related variables and then install rJava. *ERROR: configuration failed for package ‘rJava’ * removing ‘/home/ecbodevan/R/i686-pc-linux-gnu-library/2.14/rJava’ * The downloaded packages are in ‘/tmp/RtmpTC1ljU/downloaded_packages’ Mensagens de aviso perdidas: In install.packages() : installation of package ‘rJava’ had non-zero exit status Alguém sabe como resolver? Agradeço qualquer ajuda. Atenciosamente, Emerson. OBS.: ao tentar instalar o speedR e o rJava, eu rodava o R como "sudo". Uso o ubuntu.

4 4

Re: [R-br] coordenadas geograficas
by Walmes Zeviani 29 Out '12

29 Out '12

Natália, Seu caso não é minha área de competência. Estou certo que na R-br têm membros com tais competências e disposição para te ajudar. À 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 skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== 2012/10/29 Natalia Martins <nsmbarreto(a)gmail.com> > Walmes, bom dia. > Estou tentando transformar minhas coordenadas geograficas (das estações > metereologicas do PR) > em utm, procurando no help do R encontrei: > > data(state) > s <- project(state.center, "+proj=merc") > > mas eu não estou conseguindo entender qual proj eu devo usar, no link > > http://www.remotesensing.org/geotiff/proj_list/transverse_mercator.html > > até encontri mas quando tento fazer aparece o seguinte erro > > s <- project(state.center, "+proj=utm") > Erro em project(state.center, "+proj=utm") : elliptical usage required > > > se vc puder me ajudar... > > Muito obrigada!! > > Em 28 de outubro de 2012 11:11, Walmes Zeviani <walmes(a)ufpr.br> escreveu: > > Sim, para calculo de distâncias você tem que ter coordenadas decimais e >> não exagesimais (º ' ''). >> >> ========================================================================== >> Walmes Marques Zeviani >> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) >> Departamento de Estatística - Universidade Federal do Paraná >> fone: (+55) 41 3361 3573 >> VoIP: (3361 3600) 1053 1173 >> e-mail: walmes(a)ufpr.br >> skype: walmeszeviani >> twitter: @walmeszeviani >> homepage: http://www.leg.ufpr.br/~walmes >> linux user number: 531218 >> ========================================================================== >> > > > > -- > > 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 > > > Mensagem Verificada Pelo AntiSpam UFPR > >

2 1

Separa primeiro e último nome
by Fátima Lima Paula 29 Out '12

29 Out '12

Prezados, consegui o que eu queria, graças a vocês. Estava com dificuldade porque o que eu tinha depois de fazer o t(fa) não era um data.frame. Depois que transformei em data.frame foi só juntar o antigo e criado por mim e deu certo. Muito obrigada a todos. Fátima Segue o passo a passo: idoso200801=read.csv("idoso200801.csv",header=T) teste=idoso200801 nome=teste$AH_PACIENT library(stringr) separado=strsplit(as.character(nome),split=" ") head(separado) fa=sapply(separado, function(i){ n=length(i) c(i[1],paste(i[2:(n-1)],collapse=" "), i[n]) }) head(fa) nomenovo=t(fa) is.data.frame(nomenovo) as.data.frame(nomenovo) dim(nomenovo) colnames(nomenovo)=c("nome1","nome2","nome3") head(nomenovo) testenovo=cbind(teste,nomenovo) head(testenovo)

1 0

Plotar vetor com cores diferentes para grupo de pontos.
by Carlos Pereira 29 Out '12

29 Out '12

Pessoal estou com uma duvida eu tenho um vetor e quero plotar pontos em cores diferentes, consegui plotar o vetor de forma que ele agrupa os dados em subgrupos estou com duvida apenas como diferencia esses grupos por cores! Carlos Pereira da Silva Bacharel em Estatística – UEPB Mestrando em Estatística e Experimentação Agropecuária - UFLA Lavras - MG Cel: 035 - 9119 - 2245 Email: ccpsilva81(a)hotmail.com __________________________________________ Tenha Um Dia Abençoado Por DEUS!!! Sl – 34, 8

2 1

Re: [R-br] Transformação de dados no Dataframe
by fabrizzioconde＠gmail.com 28 Out '12

28 Out '12

Caros Augusto e Rodrigo, Muito obrigado pela dica de usar o reshape. Funcionou perfeitamente. []s, Fabrízzio -------Original Message------- From: r-br-request(a)listas.c3sl.ufpr.br Date: 10/26/12 17:18:01 To: r-br(a)listas.c3sl.ufpr.br Subject: Digest R-br, volume 20, assunto 40 Enviar submissões para a lista de discussão R-br para r-br(a)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 Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br 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: Transformação de dados no dataframe (Rodrigo Coster) 2. Separa primeiro e último nome (Fátima Lima Paula) 3. Re: Separa primeiro e último nome (Walmes Zeviani) 4. Re: Separa primeiro e último nome (Alisson Lucrecio) ---------------------------------------------------------------------- Message: 1 Date: Fri, 26 Oct 2012 16:16:18 -0200 From: Rodrigo Coster <rcoster(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Transformação de dados no dataframe Message-ID: <CAKU4wotutucRnLJMOaocbUL_L=RGY+tYKgMGiyAERLZE03QHYw(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" Tenta com o reshape() > dados V1 V2 V3 1 Característica A Indivíduo 1 0 2 Característica B Indivíduo 1 1 3 Característica C Indivíduo 1 2 4 Característica A Indivíduo 2 2 5 Característica B Indivíduo 2 0 6 Característica C Indivíduo 2 1 > reshape(dados, idvar='V2', timevar='V1',direction='wide') V2 V3.Característica A V3.Característica B V3.Característica C 1 Indivíduo 1 0 1 2 4 Indivíduo 2 2 0 1 []'s 2012/10/26 fabrizzioconde(a)gmail.com <fabrizzioconde(a)gmail.com> > Caro Augusto, > > > Na realidade, cada característica na tabela abaixo é um marcador genético > e estou analisando a variação desse marcador dentro da amostra de > indivíduos. Os números na terceira coluna da tabela abaixo são os alelos ( > 0 = AA, 1 = BB e 2 = AB). Os dados que tenham ESTão no formato abaixo: > > > Característica A****** > > Indivíduo 1**** > > 0**** > > Característica B**** > > Indivíduo 1**** > > 1**** > > Característica C**** > > Indivíduo 1**** > > 2**** > > Característica A**** > > Indivíduo 2**** > > 2**** > > Característica B**** > > Indivíduo 2**** > > 0**** > > Característica C**** > > Indivíduo 2**** > > 1**** > > Eu preciso mudar o posicionamento da tabela acima para o formato da tabela > abaixo: (para aplicar o modelo de regressão) > > > ** ** > > Característica A**** > > Característica B**** > > Característica C**** > > Indivíduo 1**** > > 0**** > > 1**** > > 2**** > > Indivíduo 2**** > > 2**** > > 0**** > > 1**** > > > Desta forma, a segunda tabela é de certa forma uma "transposição" da > primeira, mas não uma contagem do número de vezes que aparece o 0, 1 ou 2. > > Preciso colocar as características nas colunas do dataframe, pois irei > aplicar uma regressão linear logística, sendo as características A, B e C > as variáveis explanatórias. Tenho os dados da variável dicotomica y em > outro arquivo como um vetor coluna. > > Muito obrigado. > > Abraços, > Fabrízzio > > > > > > > > *-------Original Message-------* > > *From:* r-br-request(a)listas.c3sl.ufpr.br > *Date:* 10/26/12 14:43:29 > *To:* r-br(a)listas.c3sl.ufpr.br > *Subject:* Digest R-br, volume 20, assunto 38 > > Enviar submissões para a lista de discussão R-br para > r-br(a)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 > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > 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: Coeficiente de Correlação Cofenética (T Branquinho) > 2. Re: Transformação de dados no dataframe > (fabrizzioconde(a)gmail.com) > 3. Re: Transformação de dados no dataframe (Augusto Ribas) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 26 Oct 2012 13:36:07 -0200 > From: T Branquinho <tbranquinho1(a)hotmail.com> > To: R <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Coeficiente de Correlação Cofenética > Message-ID: <BAY167-W1202591253BD2194857FEC6917E0(a)phx.gbl> > Content-Type: text/plain; charset="utf-8" > > > Esse material pode lhe ajudarestá no pacote stat que já vem no R ao > instalá-lo. > http://stat.ethz.ch/R-manual/R-devel/library/stats/html/cophenetic.html > > Date: Fri, 26 Oct 2012 10:34:25 -0200 > From: caetano.vidigal(a)gmail.com > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Coeficiente de Correlação Cofenética > > Bom dia, > alguém sabe se existe algum pacote no R que calcula o coeficiente de > correlação cofenética ao se fazer um agrupamento?Obrigado > -- > Bruno C. Vidigal > > > > > _______________________________________________ > 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. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/ecdb53b6/attachment-0001.html > > > > ------------------------------ > > Message: 2 > Date: Fri, 26 Oct 2012 13:38:39 -0200 (Horário brasileiro de verão) > From: "fabrizzioconde(a)gmail.com" <fabrizzioconde(a)gmail.com> > To: <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Transformação de dados no dataframe > Message-ID: <508AAE7F.000009.01192@FABRIZZIO-PC> > Content-Type: text/plain; charset="iso-8859-1" > > > Caro Augusto, > > O comando table do R faz uma contagem em fatores, assim, esse comando não > solucionou o meu problema, pois o que preciso fazer é reposicionar os dados > na tabela abaixo: > > Característica AIndivíduo 10 > Característica BIndivíduo 11 > Característica CIndivíduo 12 > Característica AIndivíduo 22 > Característica BIndivíduo 20 > Característica CIndivíduo 21 > > > Tentei usar o comando reshape com o seguinte código: > > Wide <- reshape(dados, v.names = "col3", idvar = "col2", timevar = "col 1 ) > > Mas esse comando está travando o R. > > A tabela final que preciso é: > > Característica ACaracterística BCaracterística C > Indivíduo 1012 > Indivíduo 2201 > > > Li o help do comando table, mas não consegui adaptá-lo para resolver esse > problema. > > Será que terei que implementar uma função para realizar essa tarefa? > > Abraços e muito obrigado pela ajuda e pela paciencia. > > Fabrízzio Condé de Oliveira > > -------Original Message------- > > From: r-br-request(a)listas.c3sl.ufpr.br > Date: 10/24/12 12:00:14 > To: r-br(a)listas.c3sl.ufpr.br > Subject: Digest R-br, volume 20, assunto 30 > > Enviar submissões para a lista de discussão R-br para > r-br(a)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 > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > 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: Transformação de dados no dataframe (Augusto Ribas) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 24 Oct 2012 09:50:12 -0300 > From: Augusto Ribas <ribas.aca(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Transformação de dados no dataframe > Message-ID: > <CACMkfRyM6o+=mvtJ7GuQXn2xMiNYgGdo1Za0m0erZ6aw6Xc09g(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Vamos supor um exemplo de dados: > > #dados de exemplo > dados<-data.frame(col1=round(runif(20,1,4)),col2=sample(c("A","B"),20 > replace=T),col3=sample(c("A","B"),20,replace=T)) > dados > > > #primeiro com o comando paste, vamos juntar a coluna 2 com a 3, criando os > AA AB e BB > dados$col4<-factor(paste(dados$col2,dados$col3,sep="")) > > #Depois como os walmes nos ensinou aqui na lista, a gente usa o match() > #para recodificar nos valores que queremos, no caso como vc falou > > levels(dados$col4) > recod<-c(0,1,1,2) > dados$col5<-recod[match(dados$col4,levels(dados$col4))] > dados > > > #agora que vc te todos os ingredientes, é so usar o comando table() para > fazer essa tabela > #eu fiz um exemplo de tabela com a coluna 5 e a coluna 5, mas pelo que vi > vc quer a coluna > #5 e os nomes das linhas, rownames(sua planilha), entao use os nomes da > linha na tabela > table(dados$col5,dados$col4) > > Espero que ajude. > Abraços > > Em 24 de outubro de 2012 09:16, fabrizzioconde(a)gmail.com < > fabrizzioconde(a)gmail.com> escreveu: > > > Bom dia, > > > > Existe algum script que faça o seguinte: > > > > 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na > quarta > > coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. > > > > AA1 1 A B AA2 1 A A AA3 1 B B AA4 1 B B AA1 2 A A AA2 2 A B AA3 2 A B > > AA4 2 A A AA1 3 B B AA2 3 B B AA3 3 A A AA4 3 A A AA1 4 B B AA2 4 A B AA3 > > 4 A A AA4 4 A B > > > > Assim, a matriz se transformaria em: > > > > AA1 1 A B 2 AA2 1 A A 0 AA3 1 B B 1 AA4 1 B B 1 AA1 2 A A 0 AA2 2 A B > 2 > > AA3 2 A B 2 AA4 2 A A 0 AA1 3 B B 1 AA2 3 B B 1 AA3 3 A A 0 AA4 3 A A 0 > > AA1 4 B B 1 AA2 4 A B 2 AA3 4 A A 1 AA4 4 A B 2 > > > > 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima > > (talvez em um novo dataframe) e redimensioná-la para o seguinte formato: > > > > AA1 AA2 AA3 AA4 1 2 0 1 1 2 0 2 2 0 3 1 1 0 0 4 1 2 1 2 > > > > Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R > > como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de > > regressão logística penalizada stepPlr e considerar essas variáveis > > preditoras. > > > > Muito obrigado. > > Fabrízzio > > [image: FREE Christmas Animations for your email ? by IncrediMail! > > Click Here!]<http://www.incredimail.com/?id=620188&did=10500&ppd=2747 > 201111071523,9,[TypeID],[IM_UPN2]&rui=139771782&sd=20121024> > > > > _______________________________________________ > > 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. > > > > > > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://augustoribas.heliohost.org > Lattes: http://lattes.cnpq.br/7355685961127056 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr > br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.html> > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: xmas_girl_en.gif > Tipo: image/gif > Tamanho: 61041 bytes > Descrição: não disponível > URL: <http://listas.inf.ufpr > br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.gif> > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 20, assunto 30 > ***************************************** > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/e3a78863/attachment-0001.html > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: xmas_girl_en.gif > Tipo: image/gif > Tamanho: 61041 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/e3a78863/attachment-0001.gif > > > > ------------------------------ > > Message: 3 > Date: Fri, 26 Oct 2012 13:43:13 -0300 > From: Augusto Ribas <ribas.aca(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Transformação de dados no dataframe > Message-ID: > <CACMkfRxYhrSakx+OjZeB-_Drp44H=MpuYsnKRa-OaF8C5ZaDyA(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Eu não entendi perfeitamente o que é o que, desculpe. > Vc não consegue fazer um exemplo de dados semelhante ao que eu fiz, e > colocar o nome das colunas o que é Caracteristica e o que é Individuo? > > Em 26 de outubro de 2012 12:38, fabrizzioconde(a)gmail.com < > fabrizzioconde(a)gmail.com> escreveu: > > > > > Caro Augusto, > > > > O comando *table* do R faz uma contagem em fatores, assim, esse comando > > não solucionou o meu problema, pois o que preciso fazer é reposicionar os > > dados na tabela abaixo: > > > > > > Característica A****** > > > > Indivíduo 1**** > > > > 0**** > > > > Característica B**** > > > > Indivíduo 1**** > > > > 1**** > > > > Característica C**** > > > > Indivíduo 1**** > > > > 2**** > > > > Característica A**** > > > > Indivíduo 2**** > > > > 2**** > > > > Característica B**** > > > > Indivíduo 2**** > > > > 0**** > > > > Característica C**** > > > > Indivíduo 2**** > > > > 1**** > > > > Tentei usar o comando *reshape* com o seguinte código: > > > > Wide <- reshape(dados, v.names = "col3", idvar = "col2", timevar = "col > 1") > > > > Mas esse comando está travando o R. > > > > A tabela final que preciso é: > > > > > > ** ** > > > > Característica A**** > > > > Característica B**** > > > > Característica C**** > > > > Indivíduo 1**** > > > > 0**** > > > > 1**** > > > > 2**** > > > > Indivíduo 2**** > > > > 2**** > > > > 0**** > > > > 1**** > > > > Li o help do comando table, mas não consegui adaptá-lo para resolver esse > > problema. > > > > Será que terei que implementar uma função para realizar essa tarefa? > > > > Abraços e muito obrigado pela ajuda e pela paciencia. > > > > Fabrízzio Condé de Oliveira > > > > *-------Original Message-------* > > > > *From:* r-br-request(a)listas.c3sl.ufpr.br > > *Date:* 10/24/12 12:00:14 > > *To:* r-br(a)listas.c3sl.ufpr.br > > *Subject:* Digest R-br, volume 20, assunto 30 > > > > Enviar submissões para a lista de discussão R-br para > > r-br(a)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 > > > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > > endereço > > r-br-owner(a)listas.c3sl.ufpr.br > > > > 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: Transformação de dados no dataframe (Augusto Ribas) > > > > > > ---------------------------------------------------------------------- > > > > Message: 1 > > Date: Wed, 24 Oct 2012 09:50:12 -0300 > > From: Augusto Ribas <ribas.aca(a)gmail.com> > > To: r-br(a)listas.c3sl.ufpr.br > > Subject: Re: [R-br] Transformação de dados no dataframe > > Message-ID: > > <CACMkfRyM6o+=mvtJ7GuQXn2xMiNYgGdo1Za0m0erZ6aw6Xc09g(a)mail.gmail.com> > > Content-Type: text/plain; charset="windows-1252" > > > > Vamos supor um exemplo de dados: > > > > #dados de exemplo > > > > > dados<-data.frame(col1=round(runif(20,1,4)),col2=sample(c("A","B"),20 replace=T),col3=sample(c("A","B"),20,replace=T)) > > dados > > > > > > #primeiro com o comando paste, vamos juntar a coluna 2 com a 3, criando > os > > AA AB e BB > > dados$col4<-factor(paste(dados$col2,dados$col3,sep="")) > > > > #Depois como os walmes nos ensinou aqui na lista, a gente usa o match() > > #para recodificar nos valores que queremos, no caso como vc falou > > > > levels(dados$col4) > > recod<-c(0,1,1,2) > > dados$col5<-recod[match(dados$col4,levels(dados$col4))] > > dados > > > > > > #agora que vc te todos os ingredientes, é so usar o comando table() para > > fazer essa tabela > > #eu fiz um exemplo de tabela com a coluna 5 e a coluna 5, mas pelo que vi > > vc quer a coluna > > #5 e os nomes das linhas, rownames(sua planilha), entao use os nomes da > > linha na tabela > > table(dados$col5,dados$col4) > > > > Espero que ajude. > > Abraços > > > > Em 24 de outubro de 2012 09:16, fabrizzioconde(a)gmail.com < > > fabrizzioconde(a)gmail.com> escreveu: > > > > > Bom dia, > > > > > > Existe algum script que faça o seguinte: > > > > > > 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na > > quarta > > > coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. > > > > > > AA1 1 A B AA2 1 A A AA3 1 B B AA4 1 B B AA1 2 A A AA2 2 A B AA3 2 A > B > > > AA4 2 A A AA1 3 B B AA2 3 B B AA3 3 A A AA4 3 A A AA1 4 B B AA2 4 A B > AA3 > > > 4 A A AA4 4 A B > > > > > > Assim, a matriz se transformaria em: > > > > > > AA1 1 A B 2 AA2 1 A A 0 AA3 1 B B 1 AA4 1 B B 1 AA1 2 A A 0 AA2 2 A > B > > 2 > > > AA3 2 A B 2 AA4 2 A A 0 AA1 3 B B 1 AA2 3 B B 1 AA3 3 A A 0 AA4 3 A A 0 > > > AA1 4 B B 1 AA2 4 A B 2 AA3 4 A A 1 AA4 4 A B 2 > > > > > > 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima > > > (talvez em um novo dataframe) e redimensioná-la para o seguinte > formato: > > > > > > AA1 AA2 AA3 AA4 1 2 0 1 1 2 0 2 2 0 3 1 1 0 0 4 1 2 1 2 > > > > > > Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R > > > como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de > > > regressão logística penalizada stepPlr e considerar essas variáveis > > > preditoras. > > > > > > Muito obrigado. > > > Fabrízzio > > > [image: FREE Christmas Animations for your email ? by IncrediMail! > > > Click Here!]< > > http://www.incredimail.com/?id=620188&did=10500&ppd=2747,201111071523,9 > > ,[TypeID],[IM_UPN2]&rui=139771782&sd=20121024> > > > > > > _______________________________________________ > > > 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. > > > > > > > > > > > -- > > Grato > > Augusto C. A. Ribas > > > > Site Pessoal: http://augustoribas.heliohost.org > > Lattes: http://lattes.cnpq.br/7355685961127056 > > -------------- Próxima Parte ---------- > > Um anexo em HTML foi limpo... > > URL: < > > > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.html > > > > > -------------- Próxima Parte ---------- > > Um anexo não-texto foi limpo... > > Nome: xmas_girl_en.gif > > Tipo: image/gif > > Tamanho: 61041 bytes > > Descrição: não disponível > > URL: < > > > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.gif > > > > > > > ------------------------------ > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > > > > Fim da Digest R-br, volume 20, assunto 30 > > ***************************************** > > [image: FREE Christmas Animations for your email ? by IncrediMail! > > Click Here!]< > http://www.incredimail.com/?id=620188&did=10500&ppd=2747,201111071523,9 > ,[TypeID],[IM_UPN2]&rui=139771782&sd=20121026> > > > > _______________________________________________ > > 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. > > > > > > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://recologia.blogspot.com.br/< > http://augustoribas.heliohost.org> > Lattes: http://lattes.cnpq.br/7355685961127056 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/7ff2d7cb/attachment.html > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: xmas_girl_en.gif > Tipo: image/gif > Tamanho: 61041 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/7ff2d7cb/attachment.gif > > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 20, assunto 38 > ***************************************** > [image: FREE Christmas Animations for your email ? by IncrediMail! > Click Here!]<http://www.incredimail.com/?id=620189&did=10500&ppd=2748 201111071527,9,[TypeID],[IM_UPN2]&rui=139771782&sd=20121026> > > _______________________________________________ > 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. > -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/41350124/attachment-0001.html> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: xmas_tree_en.gif Tipo: image/gif Tamanho: 48428 bytes Descrição: não disponível URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/41350124/attachment-0001.gif> ------------------------------ Message: 2 Date: Fri, 26 Oct 2012 16:28:09 -0200 From: Fátima Lima Paula <fatima.lima.paula(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Separa primeiro e último nome Message-ID: <CAJzNC0vRiJ-E-+npGr+bYwsoDKOijg5jeEbzxQUAhqRinw0h7w(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Prezados, tenho uma coluna com nomes completos. Existe uma forma de separar o nome em 3 colunas: uma com primeiro nome, uma com nome do meio e uma com último nome? Obrigada Fátima -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/b6aa797b/attachment-0001.html> ------------------------------ Message: 3 Date: Fri, 26 Oct 2012 17:15:01 -0200 From: Walmes Zeviani <walmeszeviani(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Separa primeiro e último nome Message-ID: <CAFU=EkYppn=gqoKQ6yhkrSzVeXywXTSFkWCzvyQnhSNUswwR_w(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Tem sim, tudo depende do padrão apresentado. Todos os nomes são compostos por 3 partes? Bem, á pra usar a função strsplit() e quebrar no espaço em branco, o resultado é uma lista. Um probleminha surge com nomes que tem "da", "do", "das", "dos", "de" (da Silva, dos Santos, de Oliveira, do Nascimento, das Flores), mas pode ser revolvido. Veja nomes <- c("Antonio Teixeira Batista", "Suellen Cristina Albuquerque", "Paulo Matias dos Santos") # separa pelo espaço x1 <- strsplit(nomes, split=" ") x1 # junta o meio sapply(x1, function(i){ n <- length(i) c(i[1], paste(i[2:(n-1)], collapse=" "), i[n]) }) À 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 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/20121026/da795aec/attachment-0001.html> ------------------------------ Message: 4 Date: Fri, 26 Oct 2012 12:17:44 -0700 (PDT) From: Alisson Lucrecio <alissonluc(a)yahoo.com.br> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Separa primeiro e último nome Message-ID: <1351279064.24241.YahooMailNeo(a)web160705.mail.bf1.yahoo.com> Content-Type: text/plain; charset="iso-8859-1" Olá, vc nao colocou muito bem como quer fazer essa divisão, mas no caso se os nomes forem separados por uma "/" use esse script, mas se for outro separador é só fazer uma adaptação. Saudações Alisson Lucrécio da Costa library(stringr) df names3 <- str_split(df$names3, pattern="/") names3 <- do.call(rbind, names3) names3 <- data.frame(names3) names(names3)[which(names(names3) == "X1")] <- "names1" names(names3)[which(names(names3) == "X2")] <- "names2" names(names3)[which(names(names3) == "X3")] <- "names3" names3$id <- with(names, paste(names1, names2, names3, sep=".")) df2 <- cbind(names3, df[,-1]) ________________________________ From: Fátima Lima Paula <fatima.lima.paula(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Sent: Friday, October 26, 2012 4:28 PM Subject: [R-br] Separa primeiro e último nome Prezados, tenho uma coluna com nomes completos. Existe uma forma de separar o nome em 3 colunas: uma com primeiro nome, uma com nome do meio e uma com último nome? Obrigada Fátima _______________________________________________ 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. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/4334d389/attachment.html> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 20, assunto 40 *****************************************

1 0

levelplot e contourplot
by D'Apolito 28 Out '12

28 Out '12

Olá todos, Estou produzindo mapas de nível topográfico com o Lattice, usando as funções levelplot() e contourplot(). Depois de pronto, tento colocar pontos no mapa, usando points() e não funciona. Se alguém souber onde está o erro, agradeço qualquer ajudar. Att Carlos o código é o seguinte: elevation.loess = loess(z ~ x*y, data=data.s, degree = 2, span = .3, family="gaussian") elevation.fit = expand.grid(list(x = seq(65.86,71.6,0.01), y = seq(1.81,6.845,0.01))) z = predict(elevation.loess, newdata = elevation.fit) z1<-ifelse(z<0,0,z) height<-(elevation.fit$Height = as.numeric(z1)) library(lattice) levelplot(height~-x*-y, data = elevation.fit, xlab = "Latitude", ylab = "Longitude", main = "Surface elevation data", col.regions = terrain.colors(100),cuts=40,contour=T, plot.new=T, points=T) contourplot(height~-x*-y, data = elevation.fit, cuts = 20,labels = TRUE,col.regions = terrain.colors(100), contour = TRUE, pretty = TRUE,region = T) points(-70,-4,pch=20,col="blue") # aqui é onde não funciona

3 4

ANOVA para mais de 3 fatores para dados não normais....
by andrebvs 27 Out '12

27 Out '12

4 5

Dúvida em Modelo linear
by Luciano F. Sgarbi 27 Out '12

27 Out '12

Boa noite, R's Estou com uma dúvida, Como eu faço para saber o tamanho da diferença entre dois tratamentos (no caso "Alto" e "Baixo") em um modelo linear. No exemplo abaixo, não existe efeito do "x" sobre o "y", entretanto existe efeito de "t" (diferença entre os tratamentos "Alto" e "Baixo"). E é nisso que eu estou interessado. Qual seria a forma mais adequada para encontrar essa diferença? Eu usei "lm(y~x+factor(t))" mas não consegui chegar a nada. Alguém poderia me ajudar, por favor? O script é o mesmo que eu estou usando para gerar diferentes padrões e tentar entender os resultados. Abraços Luciano #-------------------------------------------------------------------------------- SD=0.1 n=50 #----- aA=0.8#intercepto Alto bA=-0#inclinação Alta aB=0.3#intercepto Baixo bB=-0#inclinação Baixo xA<-runif(n) yA<-(aA+(bA*xA)) yA<-yA+rnorm(n,0,sd=SD) xA<-xA*100 xB<-runif(n) yB<-(aB+(bB*xB)) yB<-yB+rnorm(n,0,sd=SD) xB<-xB*100 x<-c(xB,xA) y<-c(yB,yA) v<-rep(c("BAIXO","ALTO"),c(length(xB),length(xA))) plot(y~x,pch="",xlab="Eixo X",ylab="Eixo Y") points(xA,yA,col="pink",pch=16) points(xB,yB,col="grey",,pch=16) abline(lm(yB~xB),lwd=4,col="black") abline(lm(yA~xA),lwd=4,col="red") legend("topright", c("ALTO","BAIXO"), col = c("red","black"),pch=16,merge = TRUE, bg = 'gray90',lty=1) model<-lm(y~x+factor(v)) summary(model) #-------------------------------------------------------------------------------- -- Luciano F. Sgarbi Mestrando em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos Cel. (62)8174-2262 Lab. (62)3521-1732

2 4

Separa primeiro e último nome
by Fátima Lima Paula 27 Out '12

27 Out '12

Prezados, tenho uma coluna com nomes completos. Existe uma forma de separar o nome em 3 colunas: uma com primeiro nome, uma com nome do meio e uma com último nome? Obrigada Fátima

4 9

[Dúvida] Pacote para trabalhar com função de verossimilhança.
by Pedro Rafael 26 Out '12

26 Out '12

Pessoal qual o pacote que vocês me indicam para trabalhar com maximização de função log-varossimilhança. Em geral eu escrevo a função e uso a função optim a depender da função utilizo uma método específico para maximizar. Gostaria de saber qual é um bom pacote confiável que apresenta objetos e métodos para esse propósito. [ ], Pedro Rafael Diniz Marinho.

4 5

integral de dados discretos
by Luigi Cavalcanti 26 Out '12

26 Out '12

Desculpem-me pela mancada, segue email anterior, saída do dput() e CMR... Prezados colegas, tenho alguns dados de produção de metano ao longo do tempo em intervalos de 12 minutos. Gostaria da opinião de vocês sobre a melhor forma de integrar a área sob a curva formada. Já pesquisei algumas coisas e acredito que usar a função integrate() sobre o vetor gerado com as funções approxfun() ou splinefun() podem resolver o problema ao ligarem os pontos, embora redundem em resultados diferentes e não contemplem exatamente todos os pontos. Segue o banco de dados e as funções usadas por mim... ### Banco de dados data<-structure(list(Sample = 1:100, CH4A = c(0.02492, 0.0272, 0.03284, 0.04048, 0.04256, 0.0474, 0.05128, 0.05348, 0.05292, 0.05156, 0.05532, 0.05608, 0.05152, 0.05234, 0.05222, 0.05266, 0.056, 0.05008, 0.05038, 0.05148, 0.0485, 0.04748, 0.04272, 0.04294, 0.04132, 0.03824, 0.03672, 0.04056, 0.03592, 0.03492, 0.03818, 0.03628, 0.03492, 0.03588, 0.03196, 0.0338, 0.03456, 0.03284, 0.03684, 0.03282, 0.03116, 0.03082, 0.03072, 0.03488, 0.0349, 0.035, 0.03516, 0.03256, 0.02988, 0.03408, 0.02854, 0.03392, 0.02844, 0.0327, 0.031, 0.03074, 0.0308, 0.02808, 0.02796, 0.02822, 0.04016, 0.03704, 0.04814, 0.04636, 0.04824, 0.0514, 0.0506, 0.05602, 0.0623, 0.06436, 0.06132, 0.06396, 0.0612, 0.05466, 0.05832, 0.06704, 0.06776, 0.06076, 0.05698, 0.06096, 0.06566, 0.06282, 0.05626, 0.05456, 0.05716, 0.05634, 0.05258, 0.05564, 0.05608, 0.05284, 0.05228, 0.05008, 0.05322, 0.05018, 0.04752, 0.04478, 0.04732, 0.04524, 0.03788, 0.04192)), .Names = c("Sample", "CH4A"), row.names = c(NA, 100L), class = "data.frame") #### Plotagem ----- plot(CH4A~Sample,data,type='l') #### integrando com approxfun() f2 <- approxfun(data$Sample, data$CH4A) integrate(f2, min(data$Sample), max(data$Sample),subdivisions=1000) #### integrando com splinefun() f3 <- splinefun(data$Sample, data$CH4A) integrate(f3, min(data$Sample), max(data$Sample),subdivisions=1000) #### plotagem final -------- plot(data$CH4A~data$Sample, type = "n") curve(f2(x), add = TRUE, col=1) curve(f3(x), add = TRUE, col=3) f4 <- splinefun(data$Sample, data$CH4A, method ="natural") curve(f4(x), add = TRUE, col=4) ### Percebam que algumas variações bruscas não são contempladas pelas curvas ajustadas, o que provavelmente afeta a área.

2 1

integral de dados discretos
by Luigi Cavalcanti 26 Out '12

26 Out '12

Prezados colegas, tenho alguns dados de produção de metano ao longo do tempo em intervalos de 12 minutos. Gostaria da opinião de vocês sobre a melhor forma de integrar a área sob a curva formada. Já pesquisei algumas coisas e acredito que usar a função integrate() sobre o vetor gerado com as funções approxfun() ou splinefun() podem resolver o problema ao ligarem os pontos, embora redundem em resultados diferentes e não contemplem exatamente todos os pontos. Segue o banco de dados e as funções usadas por mim... Banco de dados: data Sample CH4A 1 1 0.02492 2 2 0.02720 3 3 0.03284 4 4 0.04048 5 5 0.04256 6 6 0.04740 7 7 0.05128 8 8 0.05348 9 9 0.05292 10 10 0.05156 11 11 0.05532 12 12 0.05608 13 13 0.05152 14 14 0.05234 15 15 0.05222 16 16 0.05266 17 17 0.05600 18 18 0.05008 19 19 0.05038 20 20 0.05148 21 21 0.04850 22 22 0.04748 23 23 0.04272 24 24 0.04294 25 25 0.04132 26 26 0.03824 27 27 0.03672 28 28 0.04056 29 29 0.03592 30 30 0.03492 31 31 0.03818 32 32 0.03628 33 33 0.03492 34 34 0.03588 35 35 0.03196 36 36 0.03380 37 37 0.03456 38 38 0.03284 39 39 0.03684 40 40 0.03282 41 41 0.03116 42 42 0.03082 43 43 0.03072 44 44 0.03488 45 45 0.03490 46 46 0.03500 47 47 0.03516 48 48 0.03256 49 49 0.02988 50 50 0.03408 51 51 0.02854 52 52 0.03392 53 53 0.02844 54 54 0.03270 55 55 0.03100 56 56 0.03074 57 57 0.03080 58 58 0.02808 59 59 0.02796 60 60 0.02822 61 61 0.04016 62 62 0.03704 63 63 0.04814 64 64 0.04636 65 65 0.04824 66 66 0.05140 67 67 0.05060 68 68 0.05602 69 69 0.06230 70 70 0.06436 71 71 0.06132 72 72 0.06396 73 73 0.06120 74 74 0.05466 75 75 0.05832 76 76 0.06704 77 77 0.06776 78 78 0.06076 79 79 0.05698 80 80 0.06096 81 81 0.06566 82 82 0.06282 83 83 0.05626 84 84 0.05456 85 85 0.05716 86 86 0.05634 87 87 0.05258 88 88 0.05564 89 89 0.05608 90 90 0.05284 91 91 0.05228 92 92 0.05008 93 93 0.05322 94 94 0.05018 95 95 0.04752 96 96 0.04478 97 97 0.04732 98 98 0.04524 99 99 0.03788 100 100 0.04192 #### Plotagem ----- plot(CH4A~Sample,data,type='l') #### integrando com approxfun() f2 <- approxfun(data$Sample, data$CH4A) integrate(f2, min(data$Sample), max(data$Sample),subdivisions=1000) #### integrando com splinefun() f3 <- splinefun(data$Sample, data$CH4A) integrate(f3, min(data$Sample), max(data$Sample),subdivisions=1000) #### plotagem final -------- plot(data$CH4A~data$Sample, type = "n") curve(f2(x), add = TRUE, col=1) curve(f3(x), add = TRUE, col=3) f4 <- splinefun(data$Sample, data$CH4A, method ="natural") curve(f4(x), add = TRUE, col=4) ### Percebam que algumas variações bruscas não são contempladas pelas curvas ajustadas, o que provavelmente afeta a área. Grato pela atenção, Luigi Cavalcanti Doutorando em Zootecnia UFMG/Texas A&M University

2 1

Re: [R-br] Transformação de dados no dataframe
by fabrizzioconde＠gmail.com 26 Out '12

26 Out '12

Caro Augusto, Na realidade, cada característica na tabela abaixo é um marcador genético e estou analisando a variação desse marcador dentro da amostra de indivíduos. Os números na terceira coluna da tabela abaixo são os alelos ( 0 = AA, 1 = BB e 2 = AB). Os dados que tenham ESTão no formato abaixo: Característica AIndivíduo 10 Característica BIndivíduo 11 Característica CIndivíduo 12 Característica AIndivíduo 22 Característica BIndivíduo 20 Característica CIndivíduo 21 Eu preciso mudar o posicionamento da tabela acima para o formato da tabela abaixo: (para aplicar o modelo de regressão) Característica ACaracterística BCaracterística C Indivíduo 1012 Indivíduo 2201 Desta forma, a segunda tabela é de certa forma uma "transposição" da primeira, mas não uma contagem do número de vezes que aparece o 0, 1 ou 2. Preciso colocar as características nas colunas do dataframe, pois irei aplicar uma regressão linear logística, sendo as características A, B e C as variáveis explanatórias. Tenho os dados da variável dicotomica y em outro arquivo como um vetor coluna. Muito obrigado. Abraços, Fabrízzio -------Original Message------- From: r-br-request(a)listas.c3sl.ufpr.br Date: 10/26/12 14:43:29 To: r-br(a)listas.c3sl.ufpr.br Subject: Digest R-br, volume 20, assunto 38 Enviar submissões para a lista de discussão R-br para r-br(a)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 Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br 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: Coeficiente de Correlação Cofenética (T Branquinho) 2. Re: Transformação de dados no dataframe (fabrizzioconde(a)gmail.com) 3. Re: Transformação de dados no dataframe (Augusto Ribas) ---------------------------------------------------------------------- Message: 1 Date: Fri, 26 Oct 2012 13:36:07 -0200 From: T Branquinho <tbranquinho1(a)hotmail.com> To: R <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Coeficiente de Correlação Cofenética Message-ID: <BAY167-W1202591253BD2194857FEC6917E0(a)phx.gbl> Content-Type: text/plain; charset="utf-8" Esse material pode lhe ajudarestá no pacote stat que já vem no R ao instalá-lo. http://stat.ethz.ch/R-manual/R-devel/library/stats/html/cophenetic.html Date: Fri, 26 Oct 2012 10:34:25 -0200 From: caetano.vidigal(a)gmail.com To: r-br(a)listas.c3sl.ufpr.br Subject: [R-br] Coeficiente de Correlação Cofenética Bom dia, alguém sabe se existe algum pacote no R que calcula o coeficiente de correlação cofenética ao se fazer um agrupamento?Obrigado -- Bruno C. Vidigal _______________________________________________ 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. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/ecdb53b6/attachment-0001.html> ------------------------------ Message: 2 Date: Fri, 26 Oct 2012 13:38:39 -0200 (Horário brasileiro de verão) From: "fabrizzioconde(a)gmail.com" <fabrizzioconde(a)gmail.com> To: <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Transformação de dados no dataframe Message-ID: <508AAE7F.000009.01192@FABRIZZIO-PC> Content-Type: text/plain; charset="iso-8859-1" Caro Augusto, O comando table do R faz uma contagem em fatores, assim, esse comando não solucionou o meu problema, pois o que preciso fazer é reposicionar os dados na tabela abaixo: Característica AIndivíduo 10 Característica BIndivíduo 11 Característica CIndivíduo 12 Característica AIndivíduo 22 Característica BIndivíduo 20 Característica CIndivíduo 21 Tentei usar o comando reshape com o seguinte código: Wide <- reshape(dados, v.names = "col3", idvar = "col2", timevar = "col 1") Mas esse comando está travando o R. A tabela final que preciso é: Característica ACaracterística BCaracterística C Indivíduo 1012 Indivíduo 2201 Li o help do comando table, mas não consegui adaptá-lo para resolver esse problema. Será que terei que implementar uma função para realizar essa tarefa? Abraços e muito obrigado pela ajuda e pela paciencia. Fabrízzio Condé de Oliveira -------Original Message------- From: r-br-request(a)listas.c3sl.ufpr.br Date: 10/24/12 12:00:14 To: r-br(a)listas.c3sl.ufpr.br Subject: Digest R-br, volume 20, assunto 30 Enviar submissões para a lista de discussão R-br para r-br(a)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 Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br 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: Transformação de dados no dataframe (Augusto Ribas) ---------------------------------------------------------------------- Message: 1 Date: Wed, 24 Oct 2012 09:50:12 -0300 From: Augusto Ribas <ribas.aca(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Transformação de dados no dataframe Message-ID: <CACMkfRyM6o+=mvtJ7GuQXn2xMiNYgGdo1Za0m0erZ6aw6Xc09g(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" Vamos supor um exemplo de dados: #dados de exemplo dados<-data.frame(col1=round(runif(20,1,4)),col2=sample(c("A","B"),20 replace=T),col3=sample(c("A","B"),20,replace=T)) dados #primeiro com o comando paste, vamos juntar a coluna 2 com a 3, criando os AA AB e BB dados$col4<-factor(paste(dados$col2,dados$col3,sep="")) #Depois como os walmes nos ensinou aqui na lista, a gente usa o match() #para recodificar nos valores que queremos, no caso como vc falou levels(dados$col4) recod<-c(0,1,1,2) dados$col5<-recod[match(dados$col4,levels(dados$col4))] dados #agora que vc te todos os ingredientes, é so usar o comando table() para fazer essa tabela #eu fiz um exemplo de tabela com a coluna 5 e a coluna 5, mas pelo que vi vc quer a coluna #5 e os nomes das linhas, rownames(sua planilha), entao use os nomes da linha na tabela table(dados$col5,dados$col4) Espero que ajude. Abraços Em 24 de outubro de 2012 09:16, fabrizzioconde(a)gmail.com < fabrizzioconde(a)gmail.com> escreveu: > Bom dia, > > Existe algum script que faça o seguinte: > > 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na quarta > coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. > > AA1 1 A B AA2 1 A A AA3 1 B B AA4 1 B B AA1 2 A A AA2 2 A B AA3 2 A B > AA4 2 A A AA1 3 B B AA2 3 B B AA3 3 A A AA4 3 A A AA1 4 B B AA2 4 A B AA3 > 4 A A AA4 4 A B > > Assim, a matriz se transformaria em: > > AA1 1 A B 2 AA2 1 A A 0 AA3 1 B B 1 AA4 1 B B 1 AA1 2 A A 0 AA2 2 A B 2 > AA3 2 A B 2 AA4 2 A A 0 AA1 3 B B 1 AA2 3 B B 1 AA3 3 A A 0 AA4 3 A A 0 > AA1 4 B B 1 AA2 4 A B 2 AA3 4 A A 1 AA4 4 A B 2 > > 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima > (talvez em um novo dataframe) e redimensioná-la para o seguinte formato: > > AA1 AA2 AA3 AA4 1 2 0 1 1 2 0 2 2 0 3 1 1 0 0 4 1 2 1 2 > > Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R > como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de > regressão logística penalizada stepPlr e considerar essas variáveis > preditoras. > > Muito obrigado. > Fabrízzio > [image: FREE Christmas Animations for your email ? by IncrediMail! > Click Here!]<http://www.incredimail.com/?id=620188&did=10500&ppd=2747 201111071523,9,[TypeID],[IM_UPN2]&rui=139771782&sd=20121024> > > _______________________________________________ > 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. > -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.html> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: xmas_girl_en.gif Tipo: image/gif Tamanho: 61041 bytes Descrição: não disponível URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.gif> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 20, assunto 30 ***************************************** -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/e3a78863/attachment-0001.html> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: xmas_girl_en.gif Tipo: image/gif Tamanho: 61041 bytes Descrição: não disponível URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/e3a78863/attachment-0001.gif> ------------------------------ Message: 3 Date: Fri, 26 Oct 2012 13:43:13 -0300 From: Augusto Ribas <ribas.aca(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Transformação de dados no dataframe Message-ID: <CACMkfRxYhrSakx+OjZeB-_Drp44H=MpuYsnKRa-OaF8C5ZaDyA(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" Eu não entendi perfeitamente o que é o que, desculpe. Vc não consegue fazer um exemplo de dados semelhante ao que eu fiz, e colocar o nome das colunas o que é Caracteristica e o que é Individuo? Em 26 de outubro de 2012 12:38, fabrizzioconde(a)gmail.com < fabrizzioconde(a)gmail.com> escreveu: > > Caro Augusto, > > O comando *table* do R faz uma contagem em fatores, assim, esse comando > não solucionou o meu problema, pois o que preciso fazer é reposicionar os > dados na tabela abaixo: > > > Característica A****** > > Indivíduo 1**** > > 0**** > > Característica B**** > > Indivíduo 1**** > > 1**** > > Característica C**** > > Indivíduo 1**** > > 2**** > > Característica A**** > > Indivíduo 2**** > > 2**** > > Característica B**** > > Indivíduo 2**** > > 0**** > > Característica C**** > > Indivíduo 2**** > > 1**** > > Tentei usar o comando *reshape* com o seguinte código: > > Wide <- reshape(dados, v.names = "col3", idvar = "col2", timevar = "col 1 ) > > Mas esse comando está travando o R. > > A tabela final que preciso é: > > > ** ** > > Característica A**** > > Característica B**** > > Característica C**** > > Indivíduo 1**** > > 0**** > > 1**** > > 2**** > > Indivíduo 2**** > > 2**** > > 0**** > > 1**** > > Li o help do comando table, mas não consegui adaptá-lo para resolver esse > problema. > > Será que terei que implementar uma função para realizar essa tarefa? > > Abraços e muito obrigado pela ajuda e pela paciencia. > > Fabrízzio Condé de Oliveira > > *-------Original Message-------* > > *From:* r-br-request(a)listas.c3sl.ufpr.br > *Date:* 10/24/12 12:00:14 > *To:* r-br(a)listas.c3sl.ufpr.br > *Subject:* Digest R-br, volume 20, assunto 30 > > Enviar submissões para a lista de discussão R-br para > r-br(a)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 > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > 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: Transformação de dados no dataframe (Augusto Ribas) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 24 Oct 2012 09:50:12 -0300 > From: Augusto Ribas <ribas.aca(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Transformação de dados no dataframe > Message-ID: > <CACMkfRyM6o+=mvtJ7GuQXn2xMiNYgGdo1Za0m0erZ6aw6Xc09g(a)mail.gmail.com> > Content-Type: text/plain; charset="windows-1252" > > Vamos supor um exemplo de dados: > > #dados de exemplo > > dados<-data.frame(col1=round(runif(20,1,4)),col2=sample(c("A","B"),20 replace=T),col3=sample(c("A","B"),20,replace=T)) > dados > > > #primeiro com o comando paste, vamos juntar a coluna 2 com a 3, criando os > AA AB e BB > dados$col4<-factor(paste(dados$col2,dados$col3,sep="")) > > #Depois como os walmes nos ensinou aqui na lista, a gente usa o match() > #para recodificar nos valores que queremos, no caso como vc falou > > levels(dados$col4) > recod<-c(0,1,1,2) > dados$col5<-recod[match(dados$col4,levels(dados$col4))] > dados > > > #agora que vc te todos os ingredientes, é so usar o comando table() para > fazer essa tabela > #eu fiz um exemplo de tabela com a coluna 5 e a coluna 5, mas pelo que vi > vc quer a coluna > #5 e os nomes das linhas, rownames(sua planilha), entao use os nomes da > linha na tabela > table(dados$col5,dados$col4) > > Espero que ajude. > Abraços > > Em 24 de outubro de 2012 09:16, fabrizzioconde(a)gmail.com < > fabrizzioconde(a)gmail.com> escreveu: > > > Bom dia, > > > > Existe algum script que faça o seguinte: > > > > 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na > quarta > > coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. > > > > AA1 1 A B AA2 1 A A AA3 1 B B AA4 1 B B AA1 2 A A AA2 2 A B AA3 2 A B > > AA4 2 A A AA1 3 B B AA2 3 B B AA3 3 A A AA4 3 A A AA1 4 B B AA2 4 A B AA3 > > 4 A A AA4 4 A B > > > > Assim, a matriz se transformaria em: > > > > AA1 1 A B 2 AA2 1 A A 0 AA3 1 B B 1 AA4 1 B B 1 AA1 2 A A 0 AA2 2 A B > 2 > > AA3 2 A B 2 AA4 2 A A 0 AA1 3 B B 1 AA2 3 B B 1 AA3 3 A A 0 AA4 3 A A 0 > > AA1 4 B B 1 AA2 4 A B 2 AA3 4 A A 1 AA4 4 A B 2 > > > > 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima > > (talvez em um novo dataframe) e redimensioná-la para o seguinte formato: > > > > AA1 AA2 AA3 AA4 1 2 0 1 1 2 0 2 2 0 3 1 1 0 0 4 1 2 1 2 > > > > Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R > > como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de > > regressão logística penalizada stepPlr e considerar essas variáveis > > preditoras. > > > > Muito obrigado. > > Fabrízzio > > [image: FREE Christmas Animations for your email ? by IncrediMail! > > Click Here!]< > http://www.incredimail.com/?id=620188&did=10500&ppd=2747,201111071523,9 > ,[TypeID],[IM_UPN2]&rui=139771782&sd=20121024> > > > > _______________________________________________ > > 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. > > > > > > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://augustoribas.heliohost.org > Lattes: http://lattes.cnpq.br/7355685961127056 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.html > > > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: xmas_girl_en.gif > Tipo: image/gif > Tamanho: 61041 bytes > Descrição: não disponível > URL: < > http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.gif > > > > ------------------------------ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > Fim da Digest R-br, volume 20, assunto 30 > ***************************************** > [image: FREE Christmas Animations for your email ? by IncrediMail! > Click Here!]<http://www.incredimail.com/?id=620188&did=10500&ppd=2747 201111071523,9,[TypeID],[IM_UPN2]&rui=139771782&sd=20121026> > > _______________________________________________ > 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. > -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.blogspot.com.br/<http://augustoribas heliohost.org> Lattes: http://lattes.cnpq.br/7355685961127056 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/7ff2d7cb/attachment.html> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: xmas_girl_en.gif Tipo: image/gif Tamanho: 61041 bytes Descrição: não disponível URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121026/7ff2d7cb/attachment.gif> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 20, assunto 38 *****************************************

2 1

Coeficiente de Correlação Cofenética
by bruno vidigal 26 Out '12

26 Out '12

Bom dia, alguém sabe se existe algum pacote no R que calcula o coeficiente de correlação cofenética ao se fazer um agrupamento? Obrigado -- Bruno C. Vidigal

2 2

Re: [R-br] Transformação de dados no dataframe
by fabrizzioconde＠gmail.com 26 Out '12

26 Out '12

Caro Augusto, O comando table do R faz uma contagem em fatores, assim, esse comando não solucionou o meu problema, pois o que preciso fazer é reposicionar os dados na tabela abaixo: Característica AIndivíduo 10 Característica BIndivíduo 11 Característica CIndivíduo 12 Característica AIndivíduo 22 Característica BIndivíduo 20 Característica CIndivíduo 21 Tentei usar o comando reshape com o seguinte código: Wide <- reshape(dados, v.names = "col3", idvar = "col2", timevar = "col 1") Mas esse comando está travando o R. A tabela final que preciso é: Característica ACaracterística BCaracterística C Indivíduo 1012 Indivíduo 2201 Li o help do comando table, mas não consegui adaptá-lo para resolver esse problema. Será que terei que implementar uma função para realizar essa tarefa? Abraços e muito obrigado pela ajuda e pela paciencia. Fabrízzio Condé de Oliveira -------Original Message------- From: r-br-request(a)listas.c3sl.ufpr.br Date: 10/24/12 12:00:14 To: r-br(a)listas.c3sl.ufpr.br Subject: Digest R-br, volume 20, assunto 30 Enviar submissões para a lista de discussão R-br para r-br(a)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 Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br 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: Transformação de dados no dataframe (Augusto Ribas) ---------------------------------------------------------------------- Message: 1 Date: Wed, 24 Oct 2012 09:50:12 -0300 From: Augusto Ribas <ribas.aca(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Transformação de dados no dataframe Message-ID: <CACMkfRyM6o+=mvtJ7GuQXn2xMiNYgGdo1Za0m0erZ6aw6Xc09g(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" Vamos supor um exemplo de dados: #dados de exemplo dados<-data.frame(col1=round(runif(20,1,4)),col2=sample(c("A","B"),20 replace=T),col3=sample(c("A","B"),20,replace=T)) dados #primeiro com o comando paste, vamos juntar a coluna 2 com a 3, criando os AA AB e BB dados$col4<-factor(paste(dados$col2,dados$col3,sep="")) #Depois como os walmes nos ensinou aqui na lista, a gente usa o match() #para recodificar nos valores que queremos, no caso como vc falou levels(dados$col4) recod<-c(0,1,1,2) dados$col5<-recod[match(dados$col4,levels(dados$col4))] dados #agora que vc te todos os ingredientes, é so usar o comando table() para fazer essa tabela #eu fiz um exemplo de tabela com a coluna 5 e a coluna 5, mas pelo que vi vc quer a coluna #5 e os nomes das linhas, rownames(sua planilha), entao use os nomes da linha na tabela table(dados$col5,dados$col4) Espero que ajude. Abraços Em 24 de outubro de 2012 09:16, fabrizzioconde(a)gmail.com < fabrizzioconde(a)gmail.com> escreveu: > Bom dia, > > Existe algum script que faça o seguinte: > > 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na quarta > coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. > > AA1 1 A B AA2 1 A A AA3 1 B B AA4 1 B B AA1 2 A A AA2 2 A B AA3 2 A B > AA4 2 A A AA1 3 B B AA2 3 B B AA3 3 A A AA4 3 A A AA1 4 B B AA2 4 A B AA3 > 4 A A AA4 4 A B > > Assim, a matriz se transformaria em: > > AA1 1 A B 2 AA2 1 A A 0 AA3 1 B B 1 AA4 1 B B 1 AA1 2 A A 0 AA2 2 A B 2 > AA3 2 A B 2 AA4 2 A A 0 AA1 3 B B 1 AA2 3 B B 1 AA3 3 A A 0 AA4 3 A A 0 > AA1 4 B B 1 AA2 4 A B 2 AA3 4 A A 1 AA4 4 A B 2 > > 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima > (talvez em um novo dataframe) e redimensioná-la para o seguinte formato: > > AA1 AA2 AA3 AA4 1 2 0 1 1 2 0 2 2 0 3 1 1 0 0 4 1 2 1 2 > > Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R > como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de > regressão logística penalizada stepPlr e considerar essas variáveis > preditoras. > > Muito obrigado. > Fabrízzio > [image: FREE Christmas Animations for your email ? by IncrediMail! > Click Here!]<http://www.incredimail.com/?id=620188&did=10500&ppd=2747 201111071523,9,[TypeID],[IM_UPN2]&rui=139771782&sd=20121024> > > _______________________________________________ > 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. > -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.html> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: xmas_girl_en.gif Tipo: image/gif Tamanho: 61041 bytes Descrição: não disponível URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.gif> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 20, assunto 30 *****************************************

2 1

Re: [R-br] Ajuda para criar um grafico
by Ivan Bezerra Allaman 25 Out '12

25 Out '12

Utilize as funções plot, segments e par para obter êxito! O link abaixo será útil também. http://gallery.r-enthusiasts.com/ (S,f,P) Allaman \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}

5 13

Erro ao carregar Rcmdr
by Andre Oliveira 25 Out '12

25 Out '12

Pessoas, estou aqui com ubuntu 12.04. Instalei o Rcmdr e os pacotes dependentes, mas o rgl não instala de jeito nenhum e ele é requerido pela Rcmdr. ======================================================================================================== install.packages("rgl", dep=T) Installing package(s) into ‘/home/andre/R/i686-pc-linux-gnu-library/2.15’ (as ‘lib’ is unspecified) --- Please select a CRAN mirror for use in this session --- Carregando a interface Tcl/Tk ... pronto tentando a URL 'http://www.vps.fmvz.usp.br/CRAN/src/contrib/rgl_0.92.894.tar.gz' Content type 'application/x-gzip' length 1713982 bytes (1.6 Mb) URL aberta ================================================== downloaded 1.6 Mb * 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... no configure: error: X11 not found but required, configure aborted. ERROR: configuration failed for package ‘rgl’ * removing ‘/home/andre/R/i686-pc-linux-gnu-library/2.15/rgl’ The downloaded source packages are in ‘/tmp/RtmpfivhTH/downloaded_packages’ Mensagens de aviso perdidas: In install.packages("rgl", dep = T) : installation of package ‘rgl’ had non-zero exit status ========================================================= Obrigado André Oliveira Souza

3 2

Modelos segmentados de regressão não linear...
by andrebvs 25 Out '12

25 Out '12

3 11

simulacao gama
by Natalia Martins 25 Out '12

25 Out '12

Prezados redistas, boa tarde estou tentando simular uma gama, usei o comando abaixo no R, no entanto eu gostaria de por um valor maximo e um valor minimo (truncá-la), alguem pode me ajudar? Muito obrigada. radiacao <- rgamma(n, shape=3, rate=2) -- 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

4 3

Re: [R-br] Transformação de dados no dataframe
by fabrizzioconde＠gmail.com 25 Out '12

25 Out '12

Caro Augusto, Muito obrigado pelas dicas de codificação, resolveram muito bem o problema, mas ainda estou com problema na montagem da tabela final deste e-mail, pois o comando table que sugeriu conta o número de casos de AA, BB e AB na tabela final e o que necessito é posicionar de maneira diferente o dado e não contar o número de vezes que o mesmo aparece nos fatores. AA11AB2 AA21AA0 AA31BB1 AA41BB1 AA12AA0 AA22AB2 AA32AB2 AA42AA0 AA13BB1 AA23BB1 AA33AA0 AA43AA0 AA14BB1 AA24AB2 AA34AA1 AA44AB2 Desta forma, note que o elemento (1,5) na tabela acima foi colocado na posição (1,1) da tabela abaixo. O elemento (5,5) foi colocado na posição (2 1). O elemento (9,5) foi inserido na posição (3,1) e o elemento (13,5) está na posição (4,1). Assim, consigo capturar a variável AA1 de cada indivíduo (1, 2, 3 e 4). Da mesma forma, o elemento (2,5) foi colocado na posição (1 2), o elemento (6,5) foi inserido na posição (2,2) e assim sucessivamente. AA1AA2AA3AA4 12011 20220 31100 41212 Não consigo fazer esse redimensionamento da tabela antiga para a nova. Desde já agradeço imensamente sua ajuda. Abraços, Fabrízzio -------Original Message------- From: r-br-request(a)listas.c3sl.ufpr.br Date: 24/10/2012 12:00:14 To: r-br(a)listas.c3sl.ufpr.br Subject: Digest R-br, volume 20, assunto 30 Enviar submissões para a lista de discussão R-br para r-br(a)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 Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br 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: Transformação de dados no dataframe (Augusto Ribas) ---------------------------------------------------------------------- Message: 1 Date: Wed, 24 Oct 2012 09:50:12 -0300 From: Augusto Ribas <ribas.aca(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Transformação de dados no dataframe Message-ID: <CACMkfRyM6o+=mvtJ7GuQXn2xMiNYgGdo1Za0m0erZ6aw6Xc09g(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" Vamos supor um exemplo de dados: #dados de exemplo dados<-data.frame(col1=round(runif(20,1,4)),col2=sample(c("A","B"),20 replace=T),col3=sample(c("A","B"),20,replace=T)) dados #primeiro com o comando paste, vamos juntar a coluna 2 com a 3, criando os AA AB e BB dados$col4<-factor(paste(dados$col2,dados$col3,sep="")) #Depois como os walmes nos ensinou aqui na lista, a gente usa o match() #para recodificar nos valores que queremos, no caso como vc falou levels(dados$col4) recod<-c(0,1,1,2) dados$col5<-recod[match(dados$col4,levels(dados$col4))] dados #agora que vc te todos os ingredientes, é so usar o comando table() para fazer essa tabela #eu fiz um exemplo de tabela com a coluna 5 e a coluna 5, mas pelo que vi vc quer a coluna #5 e os nomes das linhas, rownames(sua planilha), entao use os nomes da linha na tabela table(dados$col5,dados$col4) Espero que ajude. Abraços Em 24 de outubro de 2012 09:16, fabrizzioconde(a)gmail.com < fabrizzioconde(a)gmail.com> escreveu: > Bom dia, > > Existe algum script que faça o seguinte: > > 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na quarta > coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. > > AA1 1 A B AA2 1 A A AA3 1 B B AA4 1 B B AA1 2 A A AA2 2 A B AA3 2 A B > AA4 2 A A AA1 3 B B AA2 3 B B AA3 3 A A AA4 3 A A AA1 4 B B AA2 4 A B AA3 > 4 A A AA4 4 A B > > Assim, a matriz se transformaria em: > > AA1 1 A B 2 AA2 1 A A 0 AA3 1 B B 1 AA4 1 B B 1 AA1 2 A A 0 AA2 2 A B 2 > AA3 2 A B 2 AA4 2 A A 0 AA1 3 B B 1 AA2 3 B B 1 AA3 3 A A 0 AA4 3 A A 0 > AA1 4 B B 1 AA2 4 A B 2 AA3 4 A A 1 AA4 4 A B 2 > > 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima > (talvez em um novo dataframe) e redimensioná-la para o seguinte formato: > > AA1 AA2 AA3 AA4 1 2 0 1 1 2 0 2 2 0 3 1 1 0 0 4 1 2 1 2 > > Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R > como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de > regressão logística penalizada stepPlr e considerar essas variáveis > preditoras. > > Muito obrigado. > Fabrízzio > [image: FREE Christmas Animations for your email ? by IncrediMail! > Click Here!]<http://www.incredimail.com/?id=620188&did=10500&ppd=2747 201111071523,9,[TypeID],[IM_UPN2]&rui=139771782&sd=20121024> > > _______________________________________________ > 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. > -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.html> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: xmas_girl_en.gif Tipo: image/gif Tamanho: 61041 bytes Descrição: não disponível URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.gif> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 20, assunto 30 *****************************************

4 4

Ajuda para criar um grafico
by Humberto 24 Out '12

24 Out '12

Ola pessoal, gostaria de saber se alguem poderia me ajudar na criacao desse grafico em anexo (link para baixar) a partir da tabela abaixo: Species lci or uci 1 Bigeye tuna 0.0952 0.1293 0.176 2 Yellowfin tuna 0.0916 0.1255 0.172 3 Albacore 0.0913 0.1562 0.267 4 Swordfish 0.1311 0.1565 0.187 5 Sailfish 0.0368 0.0923 0.231 6 White marlim 0.1300 0.2161 0.359 7 Blue marlim 0.0738 0.1471 0.293 8 Blue shark 0.1303 0.1653 0.210 9 Crocodile shark 0.1123 0.1792 0.286 10 Pelagic stingray 0.1074 0.2052 0.392 11 Turtle oliva 0.1275 0.3802 1.134 12 Oceanic whitetip shark 0.1658 0.3788 0.866 desde ja agradeco Humberto Eu inseri um link para 1 arquivo nesta mensagem: figure.docx <http://www.yousendit.com/download/WUJhQ3Q1Qk50NjlMWE1UQw>(172KB)YouSendIt <https://www.yousendit.com>http://www.yousendit.com/download/WUJhQ3Q1Qk50NjlMWE1UQw O Mozilla Thunderbird <http://www.getthunderbird.com> torna fácil compartilhar grandes arquivos pelo e-mail.

1 0

Transformação de dados no dataframe
by fabrizzioconde＠gmail.com 24 Out '12

24 Out '12

Bom dia, Existe algum script que faça o seguinte: 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na quarta coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. AA11AB AA21AA AA31BB AA41BB AA12AA AA22AB AA32AB AA42AA AA13BB AA23BB AA33AA AA43AA AA14BB AA24AB AA34AA AA44AB Assim, a matriz se transformaria em: AA11AB2 AA21AA0 AA31BB1 AA41BB1 AA12AA0 AA22AB2 AA32AB2 AA42AA0 AA13BB1 AA23BB1 AA33AA0 AA43AA0 AA14BB1 AA24AB2 AA34AA1 AA44AB2 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima (talvez em um novo dataframe) e redimensioná-la para o seguinte formato: AA1AA2AA3AA4 12011 20220 31100 41212 Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de regressão logística penalizada stepPlr e considerar essas variáveis preditoras. Muito obrigado. Fabrízzio

2 1

Gráficos animado no R
by Jodavid Ferreira 23 Out '12

23 Out '12

Gostaria de saber como fazer um gráfico animado no R tipo, um gráfico, que vai se transformando em outro! -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

3 2

Regressão Logística
by Edson Lira 23 Out '12

23 Out '12

Estou usando as rotinas abaixo, PARA LER OS DADOS h<-read.csv2("C:\\Consulta\\CIESA\\bd.csv") h[h==""]<-NA names(h) table(h$quem_elab) str(h) h$tempo<-2012-h$ano_ini h$melh_orig<-h$niv_melhora[drop=T] h$niv_melhora<-h$niv_melhora[drop=T] h$bsc<-h$bsc[drop=T] levels(h$niv_melhora) levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("elevado","muito elevado")] <- "Melhora" levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("nenhum","pouco", "intermediário")] <- "Sem melhora" #levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("intermediário")] <-"" levels(h$quem_elab) h$part<-h$quem_elab levels(h$part)[levels(h$part)%in% c("alta direção")]<- "Direção" levels(h$part)[levels(h$part)%in% c("alta direção\ngerentes", "alta direção\ngerentes\nfuncionários","alta direção\ngerentes\noutros","gerentes", "funcionários")] <- "Gerentes" table(h$part) #PARA SELECIONAR SOMENTE AS VARIÁVEIS QUE PENSO USAR NO MODELO nb<-data.frame(melhora=h$melh_orig,bsc=h$bsc,champion=h$champion,ap_dir=h$apoio_bsc, niv_tre=h$niv_treina,cons_ext=h$cext_bsc,reu_f_est=h$reu_func_estr, com_est=h$comu_estr,ajuste=h$ajustes,ind_apr=h$ind_aprop,ind_desemp=h$ind_desd_dep, mapa_est=h$mapa_estr,tempo=h$tempo) levels(nb$melhora) levels(nb$melhora)[levels(nb$melhora)%in% c("elevado","muito elevado")] <- 1 levels(nb$melhora)[levels(nb$melhora)%in% c("nenhum","pouco", "intermediário")] <- 0 nb$desemp<-as.numeric(as.character(nb$melhora)) #PARA ELIMINAR OS NA'S cd <- subset(nb, subset=!is.na(desemp) & !is.na(bsc)& !is.na(champion)& !is.na(ap_dir)& !is.na(niv_tre)& !is.na(cons_ext) & !is.na(reu_f_est)& !is.na(com_est)& !is.na(ajuste)& !is.na(ind_apr)& !is.na(ind_desemp) & !is.na(mapa_est)& !is.na(tempo)& !is.na(desemp)) dim(cd) #MODELO DE REGRESSÃO glm3<-glm(desemp~bsc+champion+ap_dir+niv_tre+cons_ext+reu_f_est+com_est+ ajuste+ind_apr+ind_desemp+mapa_est+tempo+desemp,family=binomial,data=cd) modelstep <- step(glm3, direction = "both") logistic.display(modelstep) summary(modelstep) Estou tentando aplicar uma regressão logística, quero estimar um modelo aonde a variável resposta é o desempenho de determinada empresa, aonde tenho os níveis melhora=1, sem melhora=0 com outras variáveis. No link abaixo tem uma amostra do banco de dados. http://www.datafilehost.com/download-996309b6.html Alguém poderia dar uma avaliada no que estou fazendo? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 1

Fwd: Re: Teste no R
by Leonard de Assis 22 Out '12

22 Out '12

Esqueci de postar na lista > median.test<-function(y1,y2, na.rm = TRUE){ + z <- c(y1,y2) + g <- rep(1:2, c(length(y1), length(y2))) + m <- median(z, na.rm) + fisher.test(z<m,g)$p.value + } > turma1=c(72,79,93,91,70,95,82,80,74,86, NA) > turma2=c(78,66,65,84,69,73,71,75,68,90,76) > median.test(turma1,turma2) [1] 0.1983806 []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com Em 22/10/2012 17:40, Zé Augusto escreveu: > Luiz Henrique > > Seu script retorna o erro > "Erro em median.default(dados, na.rm = TRUE) : need numeric data" > no R2.15.0 64bits. > > ##################################################### > > Leonardo > > Desculpe a minha ignorância, mas seu script, para > y1 = turma1 > y2 = turma2 > não retorna nada. > > Grato por qualquer informação. > José Augusto > > > Em 22 de outubro de 2012 17:05, Leonard de Assis > <assis.leonard(a)gmail.com> escreveu: > median.test<-function(y1,y2){ > z<-c(y1,y2) > g <- rep(1:2, c(length(y1),length(y2))) > m<-median(z) > fisher.test(z<m,g)$p.value > } >> Fonte: Lista R >> >> []s >> Leonard de Assis >> assis <dot> leonard <at> gmail <dot> com >> >> Em 22/10/2012 16:36, Luiz Henrique da Conceicao Leal escreveu: >> >> # Dados >> turma1=c(72,79,93,91,70,95,82,80,74,86,NA) >> turma2=c(78,66,65,84,69,73,71,75,68,90,76) >> Notas=data.frame(turma1,turma2) >> >> >> # Teste da mediana >> teste.da.mediana=function(dados,alfa=0.05){ >> M=median(dados,na.rm=TRUE) >> Tabela=table(dados[,1]<=M) >> for(j in 2:ncol(dados)) Tabela=rbind(Tabela,table(dados[,j]<=M)) >> colnames(Tabela)=c("Acima da mediana","Abaixo da mediana") >> rownames(Tabela)=colnames(dados) >> if(length(table(Tabela<=5))==1) {Teste=chisq.test(Tabela,correct = FALSE)} >> else Teste=fisher.test(Tabela) >> Conclusão=if(Teste$p.value<alfa) {"Rejeita-se a hipótese de igualdade das >> medianas"} else "Não se pode rejeitar a hipótese de igualdade das medianas" >> Medianas=apply(Notas,2,median,na.rm=TRUE) >> list(Tabela=Tabela,Teste=Teste,Conclusão=Conclusão,Medianas=Medianas) >> } >> >> teste.da.mediana(Notas,alfa=0.05) >> >> >> Em 22 de outubro de 2012 15:54, Márcia Varela da Silva >> <marcinhavarela(a)hotmail.com> escreveu: >>> Como fazer o "Teste da Mediana" no R? >>> >>> >>> >>> _______________________________________________ >>> 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. >> >> >> >> _______________________________________________ >> 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

Teste no R
by Márcia Varela da Silva 22 Out '12

22 Out '12

Como fazer o "Teste da Mediana" no R?

5 4

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

22 Out '12

Saudações, Encaminho link do edital da UTFPR sobre divulgação de abertura de 03 vagas para Professor de Estatística dedicação exclusiva. O concurso tem vagas para outras áreas do conhecimento também. http://www.utfpr.edu.br/concursos/campi/ct/arquivos-gerais/Edital_068_2012_… Aos interessados, boa sorte. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

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

22 Out '12

Saudações, Encaminho link do edital da UTFPR sobre divulgação de abertura de 03 vagas para Professor de Estatística dedicação exclusiva. O concurso tem vagas para outras áreas do conhecimento também. http://www.utfpr.edu.br/concursos/campi/ct/arquivos-gerais/Edital_068_2012_… Aos interessados, boa sorte. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0

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

22 Out '12

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

4 8

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

22 Out '12

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

4 6

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

21 Out '12

1 0

read.csv
by Edson Lira 20 Out '12

20 Out '12

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

2 1

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

20 Out '12

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

4 5

ggplot2 ou Lattice
by Leonard de Assis 20 Out '12

20 Out '12

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

5 6

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

19 Out '12

Prezados, Ha tempos venho tendo problemas com a digitação/edição/uso de dados digitados em planilhas do Excel 2010 (formato . xlsx). Quando salvo e uso os dados em .txt ou .csv, o R usa os dados em problemas. No entanto, quando abro esses mesmos arquivos novamente no Excel 2010, o programa bagunça todos os valores que originalmente eram números fracionados Por exemplo, se tenho em uma coluna qualquer uma valore de, digamos, 5.677899 na planilha original, quando abro a mesma planilha em formato .txt ou . csv esse mesmo valor se transforma em 5.677.899. E esse arquivo, se usado posteriormente no R, gera resultados errôneos. Alguém já passou por isso e teria alguma dica para compartilhar? Obrigado desde logo. R. -- =8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8-->=8--> Rodrigo Silvestre Martins, PhD Bolsista FAPESP Pós-Doutorado Laboratório de Ecossistemas Pesqueiros (LabPesq) Universidade de São Paulo, Instituto Oceanográfico Praça do Oceanográfico, 191. Cidade Universitaria (sala 107-A/B) Butantã - São Paulo/SP, Brasil 05508-900 Tel: +55 11 3091 6549 Email: rodrigo.plei(a)gmail.com ; ocersm(a)lycos.com; rsmartins(a)usp.br CV Lattes: http://lattes.cnpq.br/5350064124902777 Google Scholar: http://scholar.google.com.br/citations?user=zX_EezEAAAAJ&hl=pt-BR

5 5

mclapply
by clara luz 18 Out '12

18 Out '12

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

3 3