Outubro 2012 - R-br

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 30 Out '12

30 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 29 Out '12

29 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 28 Out '12

28 Out '12

4 5