- R-br - DINF UFPR

Gráfico no R como eixo secundário.
by Alisson Lucrecio 22 Jan '13

22 Jan '13

Caros Colegas r-br, Boa tarde. Eu preciso fazer um gráfico como eixo secundário no R como o da figura em anexo, alguém poderia me dar umas dicas? Obrigado. Alisson Lucrécio da Costa

3 3

nível de confiança pgamma
by Cynthia Lara 22 Jan '13

22 Jan '13

Prezados, A função pgamma retorna a probabilidade acumulada à esquerda de uma Distribuição Gamma, certo? Ou seja, se utilizamos a função pgamma (x, alfa, beta), o resultado é a probabilidade de se obter um valor menor que x. Como posso obter o nível de confiança desse resultado? É algo em torno de 60,82% ? Obrigada pela atenção, Cynthia Oliveira Lara

3 2

quantreg: predição e erro-padrão
by Luna Hidalgo 21 Jan '13

21 Jan '13

Boa tarde. Ajustei uma regressão quantílica usando a função rq do pacote quantreg para tau=0.1 e desejo saber os valores preditos por ela, bem como os erros-padrão associados a tais valores preditos. Porém parece que o argumento se.fit não vale para a função predict quando é de um objeto rq, isto é, predict(regression, se.fit=TRUE) me gera um vetor idêntico a predict(regression, se.fit=FALSE). Alguém sabe alguma forma de obter os erros-padrões dos valores ajustados de uma regressão quantílica? --- library(quantreg) data(engel) attach(engel) regression <- rq(foodexp~income,tau=0.10, data=engel) predito <- predict(regression, se.fit=TRUE) --- Muito obrigada. Atenciosamente, Luna Hidalgo Mestranda em Estudos Populacionais e Pesquisas Sociais Escola Nacional de Ciências Estatísticas - ENCE

1 0

Construindo data frames com dados de indivíduos de outro data frame
by Fátima Lima Paula 21 Jan '13

21 Jan '13

Acho que já pedi ajuda para isso, mas vou tentar de novo. Tenho uma lista de nomes (lista) da qual quero pegar dados que pertencem a outro data frame (fati) para construir um novo data frame, apenas com aqueles indivíduos. Exemplo: lista=c("ana","basilio","clovis","duda","ellen","guto") nome=c("ana","basilio","aurora","carol","duda","helio","angela","clovis","ellen","fada","guto") idade=c("12","24","31","92","68","59","81","12","40","62","11") sexo=c("1","2","1","1","1","2","1","2","1","1","2") fati=cbind(nome,idade,sexo) fati=as.data.frame(fati) fati Quero construir um data frame com as informações contidas, que ficaria desta forma: nome idade sexo 1 ana 12 1 2 basilio 24 2 3 duda 68 1 4 clovis 12 2 5 ellen 40 1 6 guto 11 2 Alguém conseguiria me ajudar? Obrigada Fátima

3 3

Escrevendo em mais de uma linha
by Fátima Lima Paula 20 Jan '13

20 Jan '13

Já consegui obter como se faz para escrever linhas usei a dica do Rodrigo e fiz exatamente assim: for(i in naoread) { writeLines(strwrap(paste("cada=subset(dados2,nome==\"",i,'\")',sep=''))) } Agora quero acrescentar duas linhas. Fiz da seguinte forma: for(i in naoread) { writeLines(strwrap(paste("cada=subset(dados2,nome==\"",i,'\")',sep='',"\n","read=cada","\n","read=rbind(read,cada)"))) } Não sei por que, mas o texto "read=cada" está sendo escrito na mesma linha que o cada=subset, como abaixo: cada=subset(dados2,nome=="ZUMILRA ") read=cada read=rbind(read,cada) Alguém saberia dizer o porquê? Obrigada.

1 0

Plot (x,y) com desvios erro padrão da média como linha vertical
by ASANTOS 20 Jan '13

20 Jan '13

Boa noite pessoal, Gostaria de plotar um gráfico com pontos que representam as médias e que o erro padrão da média previamente fornecido aparecesse sobre os pontos, sendo meu CRM: ##Gráfico casal10<-c(16.03636,16.14,16.62,23.05)##Média 10 casais sd.casal10<-(0.8402676,0.69461420,0.5656461,1.479279)#SD para 10 casais casal1<-c(14.14000,14.21,17.37,17.26667)## Média 1casal sd.casal1<-(0.1166190,0.07520343,1.1059486,0.768739)#SD para 1 casal ger<-c(3,6,9,12)## Geração plot(ger,casal10,ylim=c(0,25),pch=16) lines(ger,casal10) points(ger,casal1, col='gray',pch=2) lines(ger,casal1, col='gray') # Eu pergunto: tem como sd.casal10 e sd.casal1 aparecerem como uma barra vertical ("arrows") sobre os pontos que representam as média? Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

4 6

R e C++
by ctrucios 19 Jan '13

19 Jan '13

Olá pessoal, Estou tentando de passar parte de um código que tenho no R para C++ (quero utilizar o Rcpp para que meu programa rode mais rápido) embora to meio enferrujado com C, alguém poderia me ajudar com esta linha do R e me dizer como teria que escrever ela em C++? (pensei em utilizar for mas fiquei pensando se existe uma forma mais direta) Codigo R: coeff[1] + sum(coeff[2:(p+1)]*y_boot[m - (1:p)]^2) + sum(coeff[(p+2):(aux3[2]+1)]*s2_boot[m - (1:q)]) Desde já agradeço. *Carlos Trucíos Maza ** : ctrucios(a)gmail.com

4 10

R + Windows 7
by Leonard de Assis 18 Jan '13

18 Jan '13

Instalei Windows em um Notebook, cou colocar o R. Me desejem sorte, vou precisar

3 5

recorte
by Hélio Gallo Rocha 18 Jan '13

18 Jan '13

Boa tarde a todos. Tendo de preparar um mapa com duas áreas, como consigo separar as duas áreas sem que fiquem ligadas por uma linha limite=read.table("rec.txt",header=T) limite plot(limite,type="l") anexei o arquivo rec.txt A finalidade é fazer um recorte em um mapa de krigagem grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 3

geoestatistica e script 100
by Hélio Gallo Rocha 18 Jan '13

18 Jan '13

Prezados. Propus um grupo de estudo de geoestatistica com um script, tá "meio" parado. Alguém conseguiu rodar? Alguma ideia para trocar? []´s -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 1

Otimização de código
by Leandro Marino 17 Jan '13

17 Jan '13

Caros, tudo bem? Estou com um código que gostaria de otimizar. Notem que faço um for para fazer tabulações simples. Como posso otimizar este código? Tenho que realizar 3k~4k tabelas... temp <- list() for(i in 1:20){ temp[[i]] <- data.frame(matrix(ncol=3,nrow=100)) colnames(temp[[i]]) <- c('ID','Tipo','Valor') temp[[i]]$ID <- paste(i,sprintf("%03d",1:100),sep='') temp[[i]]$Tipo <- factor(sample(0:10,100,replace=T),levels=c(0:10)) temp[[i]]$Valor <- rnorm(100,100,20) } dado <- data.frame(matrix(ncol=13,nrow=20)) colnames(dado) <- c('ID',paste('Tipo',sprintf('%02d',0:10),sep=''),'TipoNA') for(i in 1:20){ dado[i,] <- c(i,table(temp[[i]]$Tipo,exclude=NULL)) } Pensei algo com do.call ou sapply, mas não sei como definir apenas uma coluna... Abraços, Leandro

4 5

Pacote no R
by Ari Clecius 17 Jan '13

17 Jan '13

Boa tarde, alguém conhece algum pacote no R para modelagem de colunas de adsorção por elementos finitos? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)87687786 (085)33669042

1 0

substituir caracteres especiais de strings
by Paulo Nogueira 17 Jan '13

17 Jan '13

Olá a todos. Estou com um problema simples que não consigo resolver de maneira prática. Dado um vetor com nomes de bairos: 1- transformar os caracteres especiais ç, ã, õ, ü, etc, em caracteres normais c, a, o, u, etc; 2- substituir eventuais espaços entre nomes compostos por hífen; 3- transformar todas as letras para caracteres minúsculos Eu sei resolver 2 e 3, mas empaquei no item número 1. Como resolver de forma direta, sem utilizar gsub várias vezes? Como CMR, tenho o seguinte vetor de exemplo: c("Santa Cecília", "Anhangüera", "Freguesia Do Ó") E preciso que o resultado seja: santa-cecilia anhanguera freguesia-do-o Desde já agradeço a ajuda, Abraços, Paulo Nogueira Starzynski

3 3

Como calcular o p-valor...
by andrebvs＠bol.com.br 17 Jan '13

17 Jan '13

4 6

Fwd:
by bruno vidigal 17 Jan '13

17 Jan '13

http://www.arsnota.it/xqvrzo.php

3 2

Leitura e Gravação de Imagens usando URL
by Henrique Leme Felizatti 16 Jan '13

16 Jan '13

Caros, bom dia. Alguém poderia me ajudar a carregar uma imagem jpg/png/gif usando url e depois salvá-la em disco ? Também gostaria de sugestões de pacotes ou funções para processar essas imagens. Um exemplo de URL seria url <- "http://img0.cptec.inpe.br/~rclima/historicos/mensal/brasil/jan/abrchuvat011…" e o arquivo de saida poderia ser arquivo <- "C:\\Projetos\\TesteR\\pluviometria.gif" Grato desde já. Henrique Leme Felizatti

2 2

Re: [R-br] Junção de dados no mclapply
by Junior Beleti 15 Jan '13

15 Jan '13

Desculpa Benilton, mas seria esse o código completo: grid <- expand.grid(i=1:ncol(m2), j=1:nrow(m1)) ## versao 1 f <- function(idx) m1[grid[idx, 2],] %*% m2[, grid[idx, 1]] matrix(unlist(mclapply(1:nrow(grid), f)), nc=ncol(m2), byrow=T) ## versao 2 g <- function(idx){ aux = 0 X = m1[grid[idx, 2], ] Y = m2[, grid[idx, 1]] for (i in 1:length(X)) aux = aux + X[i]*Y[i] aux } matrix(unlist(mclapply(1:nrow(grid), g)), nc=ncol(m2), byrow=TRUE) Obrigado, Junior.

2 3

Escrevendo linhas
by Fátima Lima Paula 14 Jan '13

14 Jan '13

Prezados, preciso escrever 1353 linhas do formato abaixo, colocando entre aspas os nomes de uma lista. cadaum=subset(dados1,nome=="AIDA PEREIRA") Consegui fazer apenas uma parte, dessa forma: for(i in lista) { writeLines(strwrap(paste("cadaum=subset(dados1,nome==",i))) } Exemplo lista=c("jose","pedro","maria","teresa","Rosa") for(i in lista) { writeLines(strwrap(paste("cadaum=subset(dados1,nome==",i))) } cadaum=subset(dados1,nome== jose cadaum=subset(dados1,nome== pedro cadaum=subset(dados1,nome== maria cadaum=subset(dados1,nome== teresa cadaum=subset(dados1,nome== Rosa Como podem ver, está faltando as aspas e o parêntese do final. Só consegui chegar até aí porque vi um exemplo parecido ao que eu quero no google, é claro. Alguém saberia como faço isso? Caso contrário terei que completar com aspas e parêntese as 1353 linhas. Obrigada Fátima

2 2

Suprimir labels do eixo de X.
by Alessandro Corrêa 13 Jan '13

13 Jan '13

Prezados, Gostaria de gerar uma curva normal sem os valores do eixo de X. Consegui um resultado satisfatório da seguinte maneira (col.axis="white",tck=0) : ---------------------------- x <- seq(-3, 3, by=0.1) fx<-dnorm(x,mean=0,sd=1) plot(x, fx, col.axis="white",tck=0, main="",xlab="",ylab="",type="l",yaxt="n", col="blue",lwd=2,frame.plot=FALSE, ) ------------------------------ Haveria alguma maneira mais eficiente e elegante de suprimir os valores do eixo de X, que não tornando a suas cores brancas? Já tentei exis(labels=FALSE) , mas não está dando certo. Obrigado, Alessandro

2 1

Exclusão de linha e reordenar as posições
by Fátima Lima Paula 13 Jan '13

13 Jan '13

Caros, quando eu excluo uma linha o número dela some, como no exemplo abaixo. No banco 2 ficam as linhas 1,3,5, e 6. Como fazer para isso não acontecer? Para que no banco2 passe a ser 1,2,3,4 nome=c("Ricardo","Maria","Tereza","Pedro","José", "Germano") sobrenome=c("Cunha","Andrade","Silva","Soares","Silva","Lima") idade=c(68,67,89,78,65,67) ano=c(1944,1945,1923,1934,1947,1945) dataint=c(20121004,20120907,20120703,20120301,20120805,20121209) banco1=cbind(nome,sobrenome,idade,ano,dataint) banco1=data.frame(banco1) banco1 nome sobrenome idade ano dataint 1 Ricardo Cunha 68 1944 20121004 2 Maria Andrade 67 1945 20120907 3 Tereza Silva 89 1923 20120703 4 Pedro Soares 78 1934 20120301 5 José Silva 65 1947 20120805 6 Germano Lima 67 1945 20121209 > banco2=banco1[-c(2,4),] > banco2 nome sobrenome idade ano dataint 1 Ricardo Cunha 68 1944 20121004 3 Tereza Silva 89 1923 20120703 5 José Silva 65 1947 20120805 6 Germano Lima 67 1945 20121209 Obrigada Fátima

2 2

Looping
by Fátima Lima Paula 11 Jan '13

11 Jan '13

Pessoal, tenho uma lista de nomes num banco, que pertencem a outro banco. Exemplo. lista=c("a","b","c","d","e","f","g") nome=c("a","b","a","c","d","h","a","g","e","f","g") inf=c("1","2","3","9","6","5","8","12","4","6","11") banco=cbind(nome,inf) banco=as.data.frame(banco) Quando faço, por exemplo which(banco=="a") Retorna 1 3 7 Preciso criar uma lista com os ids desses nomes no banco. Assim, a lista para o meu exemplo ficaria: lista2= 1 3 7 2 4 5 9 10 8 11 Na verdade, pensei num looping, mas não sei fazer. A minha ideia seria lista2=which(banco=="n") sendo n=cada um dos nomes da 'lista' Teria que fazer isto com todos os n da 'lista' na ordem e ir guardando na lista2 Isso é possível? Alguém poderia ajudar? Obrigada Fátima

3 5

Adicionar um intervalo (quebra - axis.break) no eixo y em gráficos no R
by Wenceslau Teixeira 11 Jan '13

11 Jan '13

*Caros colegas da lista R,* * * *Envio este post para saber se alguém pode me auxiliar ou dar algumas dicas de bibliografia sobre como colocar uma quebra (intervalo) no eixo y em gráficos.* * * *Segue o código reproduzível abaixo cujo intervalo ultrapassa o eixo y e não gostaria de aumentar os limites de Maximo no eixo y pois a maioria dos dados acaba ficando restrita num pequeno intervalo que dificulta a visualização. Eu encontrei o comando (**axis.break) mas não consegui implementar no grafico. * * * *Sds estatísticas e antecipadamente agradeço alguma orientação.* * * *Wenceslau Teixeira * *Embrapa Solos* * * *#X SÃO AS DATAS* * * *x <- c("30/04/2010","13/05/2010","28/05/2010",* * "15/06/2010","30/06/2010","15/07/2010",* * "30/07/2010","15/08/2010","30/08/2010",* * "15/09/2010","30/09/2010","15/10/2010",* * "30/10/2010","15/11/2010","30/11/2010",* * "15/12/2010","30/12/2010","15/01/2011",* * "30/01/2011","15/02/2011","28/02/2011",* * "15/03/2011","30/03/2011","15/04/2011",* * "30/04/2011","15/05/2011","30/05/2011",* * "15/06/2011","30/06/2011","15/07/2011")* * * *#Y* *y <- c(0.067,0.094,0.039,0.048,0.123,0.080,0.036,0.179,0.069,0.601,NA,NA,0.057,NA,NA,0.028,0.273,0.052,0.045,0.476,0.154,0.314,0.120,0.100,0.079,1.208,8.388,3.307,0.774,0.244) * * * *#SDXX E O DESVIO PADRAO DO ELEMENTO FE A 10 CM QUE SERA FEITO O GRAFICO* * * *SDFE10 <-c(0.018,0.028,0.013,0.047,0.132,0.056,0.011,0.177,0.008,0.000,NA,NA,0.028,NA,NA,0.000,0.387,0.030,0.000,0.622,0.037,0.157,0.054,0.030,0.031,1.907,9.203,4.629,1.041,0.189) * * * *# colocar num data.frame, eliminar NA, converter para data* *da <- data.frame(x=as.POSIXlt(x, format="%d/%m/%Y"), y=y,cex.axis=0.5, sd=SDFE10)* *#da <- na.omit(da)* *da$x* * * *#fazer o grafico * * * *plot(y~x,tck=0.01, data=da,ylab=expression(Concentração~de~Ferro~na~solução~do~solo~(mg~L^{-1})),tck=0.01,xlab="", pch=24,cex.axis=0.5, type="p", xaxt="n",ylim=c(-0.2,0.8))* *with(da, arrows(x, y-sd, x, y+sd,cex.axis=0.5, code=3, angle=90,length=0.03))* * * *axis.POSIXct(side=1, at=da$x, las=2, format="%d/%m/%Y",cex.axis=0.5)* * * *#adicionado os outros dados da prof50 cm * * * *par(new=TRUE)* * * *#Y SÃO OS DADOS DO FE50CM* * * *y1 <- c(0.097,0.029,0.026,0.053,0.056,0.014,0.015,0.122,NA,0.140,0.160,NA,0.088,NA,NA,NA,NA,NA,0.017,0.070,0.052,0.173,0.079,0.037,0.050,0.420,0.280,0.030,0.057,0.037) * * * *#SDXX E O DESVIO PADRAO FERRO A 50 CM FEITO O GRAFICO * * * *SDFE50 <- c(0.078,0.000,0.012,0.057,0.011,0.000,0.000,0.000,NA,0.000,0.000,,0.000,NA,NA,NA,NA,NA,0.000,0.027,0.029,0.062,0.019,0.006,0.000,0.114,0.227,0.000,0.000,0.000) * *# break axis* *#axis.break(2, 22, style="zigzag")* * * * * *#fazer o segundo grafico * * * *plot(y1~x, data=da,cex.axis=0.5,* * ylab=expression(Concentração~de~Ferro~na~solução~do~solo~(mg~L^{-1})),xlab="",tck=0.01, pch=19, type="p", xaxt="n",ylim=c(-0.2,0.8))* *with(da, arrows(x, y1-sd, x, y1+sd, code=3, angle=90,length = 0.03))* * * *#Adicionando uma legenda* * * *legend("topright", c("10 cm", "50 cm"), pch = c(24,19))* * * *#informação sobre O tratamento * *legend("topleft", c("Tratamento - TENS Dendê"))* * * -- Wenceslau Geraldes Teixeira Pesquisador Embrapa Embrapa Solos (www.cnps.embrapa.br) Rua Jardim Botânico, 1024 CEP 22460-000 Rio de Janeiro - RJ - Brasil Tel: 55 21 21794583 Skype: wgt007 E-mail: wenceslau.teixeira(a)embrapa.br wenceslau.teixeira(a)pq.cnpq.br CV - http://lattes.cnpq.br/4136381865367906 ________________________________________________________ II ENCONTRO BRASILEIRO DE FÍSICA DO SOLO - 2013 2nd Brazilian Soil Physics Meeting Rio de Janeiro, 06 a 10 de Maio de 2013 http://www.bspm2013.org/ -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

2 1

encontrar num banco a linha de um nome
by Fátima Lima Paula 11 Jan '13

11 Jan '13

Acho que estou surtando. Não consigo achar o indivíduo pelo nome em um data frame. Tenho o nome dele e quero achar a linha em que ele está. Alguém poderia me ajudar, por faor? Obrigada Fátima

3 4

data.frame para planilha de dados
by Nicolay Cunha 11 Jan '13

11 Jan '13

Prezados, Gostaria de saber se existe alguma maneira automatizada de eu reverter um data.frame para uma planilha de dados como no exemplo abaixo. Grato! Nicolay #################################################################################################################### # Data.frame # a partir disso! Ind Tipo L1 L2 L3 M1 M2 M2 C1 C2 C3 S 1 LONGO NA NA NA 20 19 18 15 14 13 30 2 MÉDIO 32 31 30 NA NA NA 16 15 13 22 3 CURTO 33 32.5 31 23 22 21.5 NA NA NA 14 #################################################################################################################### # Planilha de dados # obter isso! Ind Tipo Ant Pol S 1 LONGO M1 20 30 1 LONGO M2 19 30 1 LONGO M3 18 30 1 LONGO C1 15 30 1 LONGO C2 14 30 1 LONGO C3 13 30 2 MÉDIO L1 32 22 2 MÉDIO L2 31 22 2 MÉDIO L3 30 22 2 MÉDIO C1 16 22 2 MÉDIO C2 15 22 2 MÉDIO C3 13 22 3 CURTO L1 33 14 3 CURTO L2 32.5 14 3 CURTO L3 31 14 3 CURTO M1 23 14 3 CURTO M2 22 14 3 CURTO M3 21.5 14 -- 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

2 2

Re: [R-br] Junção de dados no mclapply
by Junior Beleti 10 Jan '13

10 Jan '13

Entendido, Benilton. Apenas na instrução: X = m1[grid[idx, 2], ] Está criando um objeto X que seria um grid (malha) com os elementos de m1, certo? Mas "idx, 2" teria qual significado? Muito obrigado, Junior.

2 1

Excluir nomes sem sobrenome
by Fátima Lima Paula 10 Jan '13

10 Jan '13

Descobri que não consigo separar os nomes no bancão poque têm vários nomes que não têm sobrenome nenhum. Só que são muitos. Existe uma forma de excluir essas linhas? Aí vai o exemplo com a Maria. Não roda. Quando peço um summary no "separado", retorna que Maria tem length=1. Seria retirar todos que têm length=1. nome=c("Teresa dos Santos","Paulo da Silva","Denis Henrique de Paula","Pedro dos Santos","Pedro Maia", "Maria") idade=c("80","89","78","78","78","85") dataint=c("26","20","25","25","25","20") dados1=data.frame(nome,idade,dataint) dados1 str(dados1) nome1=dados1$nome # remover da,das,do,dos,de nome1=gsub(" d(a|e|o)s? ", " ", nome1) separado1=strsplit(as.character(nome1), split=" ") head(separado1) fa1=sapply(separado1, function(i){ n=length(i) if(n>2) return(c(i[1], paste(i[2:(n-1)],collapse=" "), i[n])) if(n==2) return(c(i[1], NA, i[n])) }) nomenovo1=t(fa1) nomenovo1=as.data.frame(nomenovo1) dim(nomenovo1) colnames(nomenovo1)=c("nome1","nome2","nome3") dados2=cbind(dados1,nomenovo1) head(dados2) Por favor, alguém pode ajudar?

3 3

Diferenças nas saídas
by Fátima Lima Paula 10 Jan '13

10 Jan '13

Pessoal, desculpem, mas não consigo rodar um script que o Walmes me passou. No exemplo, fica tudo bem, mas na real dá zebra. Acho que os objetos estão sendo considerados diferentemente. No exemplo tudo estava com factor. Passei tudo para factor. Nos exemplos os dois bancos eram data.frames. Mas quando faço a função em um a saída vem diferente da do outro. Por favor, alguém poderia me ajudar. Exemplo que apresentei nome=c("Teresa dos Santos","Paulo da Silva","Denis Henrique de Paula","Pedro dos Santos","Pedro Maia", "Maria do Carmo") idade=c("80","89","78","78","78","85") dataint=c("26","20","25","25","25","20") dados1=data.frame(nome,idade,dataint) dados1 str(dados1) nome1=dados1$nome # remover da,das,do,dos,de nome1=gsub(" d(a|e|o)s? ", " ", nome1) separado1=strsplit(as.character(nome1), split=" ") head(separado1) fa1=sapply(separado1, function(i){ n=length(i) if(n>2) return(c(i[1], paste(i[2:(n-1)],collapse=" "), i[n])) if(n==2) return(c(i[1], NA, i[n])) }) nomenovo1=t(fa1) nomenovo1=as.data.frame(nomenovo1) dim(nomenovo1) colnames(nomenovo1)=c("nome1","nome2","nome3") dados2=cbind(dados1,nomenovo1) head(dados2) Quando faço head para o fa1 no meu exemplo, retorna assim, que pelo pouco que sei de R, parece uma matriz > head(fa1) [,1] [,2] [,3] [,4] [,5] [,6] [1,] "Teresa" "Paulo" "Denis" "Pedro" "Pedro" "Maria" [2,] NA NA "Henrique" NA NA NA [3,] "Santos" "Silva" "Paula" "Santos" "Maia" "Carmo" Já quando peço o head do fa do meu bancão, retorna assim e zebra tudo a partir daí. > head(fa) [[1]] [1] "HERALDO" "TAVARES" "SILVA" [[2]] [1] "JUPITER" NA "SILVA" [[3]] [1] "JACI" NA "SANTOS" [[4]] [1] "MARCOS" "SILVA" "FILHO" [[5]] [1] "PAULO" "SILVA" "PORTO" [[6]] [1] "INACIO" "AVELINO" "BELO" Alguém tem ideia do que seja? Obrigada Fátima

3 3

Busca textual
by Clécio da Silva Ferreira 10 Jan '13

10 Jan '13

Olá. Tem uma mestranda de Linguística interessada em fazer busca textual de ocorrências de associações de uma palavra com outra na mesma sentença. Essa busca é feita em textos (livros, transcrições de falas, etc.). Existe alguma ferramenta no R para isso? Atte, xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx x Prof. Clécio da Silva Ferreira x x Depto Estatística - UFJF x xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx

4 3

Separando nomes e retirando de da dos das
by Fátima Lima Paula 09 Jan '13

09 Jan '13

Prezados, há um tempo Walmes me passou um script para separar nomes. Como vcs podem verificar o nome do meio ficou com um problema (Pedro Maia), mas eu só ia utilizar o primeiro e último nome, por isso não me importei. Mas agora preciso de tudo. Assim, eu gostaria de saber se alguém poderia me ajudar desta forma: acertando o problema que retornou para quem só tem dois nomes e retirando os "de", "das", "da", etc dos nomes. Não sei programar e nem imagino como se faz para acrescentar isso no script. Você poderia me ajudar? nome=c("Teresa dos Santos","Paulo da Silva","Denis Henrique de Paula","Pedro dos Santos","Pedro Maia", "Maria do Carmo") idade=c("80","89","78","78","78","85") dataint=c("26","20","25","25","25","20") dados=data.frame(nome,idade,dataint) dados Ulitizei seu script assim: nome=dados$nome separado=strsplit(as.character(nome),split=" ") fa=sapply(separado, function(i){ n=length(i) c(i[1],paste(i[2:(n-1)],collapse=" "), i[n]) }) nomenovo=t(fa) nomenovo=as.data.frame(nomenovo) dim(nomenovo) colnames(nomenovo)=c("nome1","nome2","nome3") dados1=cbind(dados,nomenovo) head(dados1) Obrigada Fátima

2 4

Excluindo "da" "de" ao mesmo tempo
by Fátima Lima Paula 09 Jan '13

09 Jan '13

Olá, pessoal. O Daniel Marcelino me passou esse script com a solução para extrair os "da" "de" dos nomes. Mas como sou um horror em R não consigo eliminar tudo ao mesmo tempo. Se rodo o script dele, funciona para o "de". Mas quando repito para "da" retorna um problema. Como não manjo nada, tentei usar list(c("de","da")). Não funcionou. Tentei usar ("de" | "da"). Nada. Alguém pode me ajudar, por favor. Obrigada Fátima mydata <- data.frame( name = a<-c("Ricardo de Almeida","Maria José da Cruz","Maria Tereza dos Santos","Paulo da Rocha","José Maria", "Germano das Neves"), age = b<-c(67,89,78,65,68,67), birth = c <-c(1945,1923,1934,1947,1944,1945), date = d<-c(20120907,20120703,20120301,20120805,20121004,20121209) ) mydata$name2 <-strsplit(as.character(mydata[,1]),split=c(" de "))

2 1

Re: [R-br] Não estou conseguindo instalar o pacote ggplot2
by Cynthia Lara 09 Jan '13

09 Jan '13

Obrigada Marcelo! Consegui resolver o problema utilizando a função setInternet2(use = TRUE). Atenciosamente, Cynthia Oliveira Lara

1 0

[off-topic] Criar uma página para armazenar material
by Andre Oliveira 09 Jan '13

09 Jan '13

Boa tarde, estou querendo criar uma página simples para eu armazenar lista de exercícios avisos para minhas turmas de estatística entre outras coisas que um professor precisa. Algo bem simples e que possa trabalhar com linux com facilidade. Alguém sugere algum serviço gratuito? Conheço pouco sobre criação de páginas de preferência algo bem fácil de se trabalhar. obrigado André Oliveira Souza

4 3

Soundex
by Fátima Lima Paula 09 Jan '13

09 Jan '13

Olá Marcelino, quando você fez o Soundex, não retirou os "de", "da", "do", e plurais dos nomes, não é? Observei que usei o nome do meio para fazer o linkage e ele não captou um par porque num dos bancos o nome estava sem o "de". Exemplo: Fátima de Lima e Fátima Lima Abs Fátima

2 2

DESDOBRAMENTO DE UM FATORIAL TRIPLO
by Tiago Souza Marçal 09 Jan '13

09 Jan '13

Boa tarde pessoal, Estou tentando desdobrar a interação tripla do exemplo abaixo, entretanto, as somas de quadrado e os graus de liberdade estão dando 0. Ex: anos<-c(rep(1,144),rep(2,144)) ambientes<-c(rep(1,48),rep(2,48),rep(3,48),rep(1,48),rep(2,48),rep(3,48)) genotipos<-c(rep(c(rep(1,6),rep(2,6),rep(3,6),rep(4,6),rep(5,6),rep(6,6),rep(7,6),rep(8,6)),6)) blocos<-c(rep(c(1,2,3,4,5,6),48)) resp<-c(195,196,192,208,180,213,216,194,198,209,238,236,307,341,305,316,335,340,435,403,499,422,487,483,218,192,222,221,218,219,227,422,429,448,463,461,134,178,188,285,294,295,707,803,828,714,848,808,141,153,168,270,146,250,244,134,142,253,261,234,342,359,349,359,369,330,448,465,431,449,459,440,243,163,250,244,210,231,246,440,444,454,462,464,141,154,172,269,249,259,759,840,838,741,845,806,569,460,433,373,480,384,490,191,406,500,424,314,842,703,881,804,912,929,739,749,629,630,848,721,931,849,988,880,893,993,290,204,380,311,215,295,190,109,239,137,200,146,200,218,132,249,235,202,294,324,386,382,319,395,625,620,719,619,627,620,626,540,668,662,639,575,882,984,994,796,802,709,525,515,540,409,589,467,180,256,260,198,292,109,365,372,364,274,396,302,929,783,894,866,772,872,130,394,399,203,387,213,211,396,192,208,234,238,303,442,302,316,332,449,436,608,492,429,487,489,218,393,225,220,214,217,226,424,426,541,462,467,139,279,389,288,296,298,701,801,823,913,849,805,195,390,199,239,217,180,306,380,488,372,309,375,502,534,514,726,682,459,273,324,278,296,302,293,304,315,316,337,291,312,273,355,286,283,331,403,305,296,283,318,309,481,283,255,306,313,274,484) dados<-data.frame(anos,ambientes,genotipos,blocos,resp) #Estou usando os seguintes comandos: anos<-factor(anos) ambientes<-factor(ambientes) genotipos<-factor(genotipos) blocos<-factor(blocos) av<-aov(resp~anos*ambientes/genotipos+(anos:ambientes:blocos)) des.tab<-sapply(paste("anos:ambientes", levels(anos:ambientes), sep=""), simplify=F, grep, x=names(coef(av)[av$assign==4])) des.tab1<-summary(av, split=list("anos:ambientes:genotipos"=des.tab)) Agradeço desde já por qualquer ajuda. Att. Tiago.

2 4

Não estou conseguindo instalar o pacote ggplot2
by Cynthia Lara 08 Jan '13

08 Jan '13

Prezados, bom dia! Não estou conseguindo instalar o pacote ggplot2. A seguinte mensagem é apresentada: Aviso: unable to access index for repository http://cran.fiocruz.br/bin/windows/contrib/2.13 Aviso: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.13 Mensagens de aviso perdidas: 1: In open.connection(con, "r") : não foi possível conectar a 'cran.r-project.org' na porta 80. 2: In getDependencies(pkgs, dependencies, available, lib) : package ggplot2 is not available (for R version 2.13.2) Tentei instalar outras versões, como a 2.15.2 e a 2.14.0, mas não obtive sucesso. Alguém pode me ajudar? Obrigada! Cynthia Oliveira Lara

2 1

Pegar variáveis de um banco para outro com restrições
by Fátima Lima Paula 08 Jan '13

08 Jan '13

Pessoal, desculpem, mas preciso urgente de ajuda. Quando eliminei as linhas repetidas com unique, usei apenas duas variáveis. Agora fiquei com um banco muito menor, mas sem as variáveis que eu quero utilizar. Como faço para trazer as variáveis do primeiro banco para o segundo. Aí vai o script: nome=c("Teresa","Paulo","Pedro","Pedro","Pedro", "Maria","Teresa","Pedro") idade=c("80","89","78","78","78","85","80","78") dataint=c("26","20","25","25","25","20","26","30") banco=data.frame(nome,idade,dataint) banco banco1=banco[,c(1,3)] banco2=unique(banco1) banco2 Como faço agora para pegar as idades do banco para utilizar no banco2? O que tenho é um banco muito grande com muitas variáveis e se faço o unique com ele completo não elimino linha nenhuma já que as outras variáveis são diferentes. Então, utilizei apenas as duas variáveis que me interessavam para fazer o unique. Agora, no banco sem duplicatas, preciso das outras variáveis. Obrigada

2 1

Trazendo de volta variáveis depois do unique
by Fátima Lima Paula 08 Jan '13

08 Jan '13

Pessoal, desculpem, mas aí vai mais uma dúvida elementar. Quando eliminei as linhas repetidas com unique, usei apenas duas variáveis. Agora fiquei com um banco muito menor, mas sem as variáveis que eu quero utilizar. Como faço para trazer as variáveis do primeiro banco para o segundo. Aí vai o script: nome=c("Teresa","Paulo","Pedro","Pedro","Pedro", "Maria","Teresa","Pedro") idade=c("80","89","78","78","78","85","80","78") dataint=c("26","20","25","25","25","20","26","30") banco=data.frame(nome,idade,dataint) banco banco1=banco[,c(1,3)] banco2=unique(banco1) banco2 Como faço agora para pegar as idades do banco para utilizar no banco2? Obrigada Fátima

2 4

excluir linhas com duas condições
by Fátima Lima Paula 08 Jan '13

08 Jan '13

Prezados, preciso de uma ajuda crucial. Tenho um banco enorme de indivíduos com suas datas de internação. Preciso excluir as datas repetidas do mesmo indivíduo. Exemplo: nome=c("Teresa","Paulo","Pedro","Pedro","Pedro", "Maria","Teresa","Pedro") idade=c("80","89","78","78","78","85","80","78") dataint=c("26","20","25","25","25","20","26","30") banco=data.frame(nome,idade,dataint) banco Preciso excluir os nomes iguais que tenham dataint iguais. Por favor, alguém poderia me ajudar? Obrigada Fátima

2 2

Fwd:
by Wagner Pinheiro 07 Jan '13

07 Jan '13

http://englisheasy4u.com/lpvuolf.php

1 0

salvar gráficos
by Luiz André Zardo 06 Jan '13

06 Jan '13

Prezados, Poderiam me indicar qual comando, pode ser utilizado para salvar os gráficos construídos. Obrigado.

3 2

Criar colunas com id
by Fátima Lima Paula 06 Jan '13

06 Jan '13

Prezados, tenho um banco e preciso criar uma coluna com o número dos respectivos id de cada um, da seguinte forma nome=c("Maria","Teresa","Pedro") idade=c("89","76","67") id2=c("3","1","6",) banco1=cbind(nome,idade,id,read) nome=c("Teresa","Ana","Maria","Joyce","Isabela","Pedro") idade=c("76","60","89","90","69","67") banco2=cbind(nome,idade) Quero criar uma coluna no banco2 chamada reinter, em que os nomes que aparecem no banco 1 tenham "sim" e o restante tenham "não". Para isso preciso dizer que os ids do banco 2 que tenham valor do id2 do banco 1 levarão "sim" e o restante levará "não". Deu para entender? O resultado ficaria assim: nome idade reinter 1 Teresa 76 s 2 Ana 60 n 3 Maria 89 s 4 Joyce 90 n 5 Isabela 69 n 6 Pedro 67 s Alguém poderia me ajudar? Obrigada Fátima

2 6

Revista Sigmae
by Eric Ferreira 05 Jan '13

05 Jan '13

Caros colegas useRs, É com satisfação que comunico o lançamento do primeiro volume da revista científica Sigmae: http://publicacoes.unifal-mg.edu.br/revistas/index.php/sigmae/issue/current Dedicada a artigos em Matemática e Estatística (puras, aplicadas e educações), trabalhos feitos em R não constantes. A revista é exclusivamente eletrônica, totalmente gratuita e admite artigos em Português e Inglês. Aos pesquisadores, peço que considerem a revista Sigmae quando forem publicar seus artigos. Abraços, -- Dr Eric B Ferreira Exact Sciences Department Federal University of Alfenas Brazil

1 0

TESTE DE SHAPIRO-WILK
by Tiago Souza Marçal 05 Jan '13

05 Jan '13

Boa tarde pessoal, Gostaria de saber se é possivel fixar um nível alpha para a realização do teste de Shapiro-Wilk. Agradeço desde já. Att. Tiago.

3 6

Fwd: R version 3.0.0
by Leandro Marino 04 Jan '13

04 Jan '13

Enviado de um dispositivo móvel. Sent from a mobile device. ---------- Mensagem encaminhada ---------- De: "Peter Dalgaard" <pd.mes(a)cbs.dk> Data: 01/01/2013 18:18 Assunto: R version 3.0.0 Para: "r-announce(a)lists.r-project.org" <r-announce(a)r-project.org> Cc: This is no secret to those who read the NEWS file of the development version regularly, as the following has been in place since December 12th: \section{\Rlogo CHANGES IN R-devel}{ \subsection{SIGNIFICANT USER-VISIBLE CHANGES}{ \itemize{ \item It is intended that this version will be released as \R 3.0.0. .... However, it seems reasonable to supplement this with a more direct public announcement. The intended timing is to follow the annual release schedule and have R 3.0.0 around April 1 and a finalizing 2.15.3 a month earlier. Major R releases have not previously marked great landslides in terms of new features. Rather, they represent that the codebase has developed to a new level of maturity. This is not going to be an exception to the rule. Version 1.0.0 was released at a point in time when we felt that we had reached a level of completeness and stability high enough to characterize a full statistical system, which could be put to production use. Version 2.0.0 came out after strong enhancements of the memory management subsystem as well as several major features, including Sweave. Version 3.0.0, as of this writing, contains only really major new feature: The inclusion of long vectors (containing more than 2^31-1 elements!). More changes are likely to make it into the final release, but the main reason for having it as a new major release is that R over the last 8.5 years has reached a new level: we now have 64 bit support on all platforms, support for parallel processing, the Matrix package, and much more. On behalf of the R Core Team, Peter D. Happy New Year! -- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd.mes(a)cbs.dk Priv: PDalgd(a)gmail.com _______________________________________________ R-announce(a)r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-announce

1 0

R em Nuvem
by Tiago Vieira 04 Jan '13

04 Jan '13

Olá colegas usuários do R, Desta vez meu problema é o seguinte: Estou escrevendo uma dissertação e preciso rodar um script com uma simulação com várias iterações. O problema é que em meus computadores está sendo inviável executar esta tarefa, devido a exigência de memória e tempo, visto que, dependendo o configuração da máquina, estimo que a execução deve demorar alguns meses para ser concluída. Então pergunto: Vocês conhecem alguma funcionalidade que me permita executar esta tarefa em uma nuvem e não em uma máquina local ou até mesmo alguma empresa que preste este tipo de serviço? Pois só encontro opções de armazenamento de dados na Rede, mas não de execução de programas (que o caso é o R). Desde já, agradeço pela atenção e ajuda! Att. Tiago Vieira.

7 8

Ajuda para escrever um arquivo texto.txt usando a função sink()
by Augusto Ribas 03 Jan '13

03 Jan '13

Ola a todos e feliz ano novo. Eu estou com um problema para exportar um resultado do R. Eu gostaria de escrever um arquivo de texto .txt assim: ATGCAACT A: 5 1 0 0 5 5 0 0 C: 0 0 1 4 2 0 6 1 G: 1 1 6 3 0 1 0 0 T: 1 5 0 0 0 1 1 6 Eu tenho um string chamado consensus que são a sequência de bases ali em cima: ATGCAACT E tenho uma matriz que é aquela parte ali de baixo: A: 5 1 0 0 5 5 0 0 C: 0 0 1 4 2 0 6 1 G: 1 1 6 3 0 1 0 0 T: 1 5 0 0 0 1 1 6 O problema é que quando a matriz tem muitas colunas, eu não consigo produzir o resultado que queria. Eu tentei usar o comando sink() pra escrever o arquivo. Dai eu fiz o seguinte: ########################## sink("resultado.txt") cat(paste(consensus,collapse="")) print(resul) sink() ############################# E tive como resultado: ################################################################################################# CCAAGATCGTCACCAGTTCACGCCGTACATTCCCAAAGTTCTAGCCATCAAGTGTGCAACAAGTCACAGACTATAAAGGACTT... 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 A 1 3 3 3 3 3 3 2 2 3 2 3 1 1 4 2 3 2 2 6 1 1 1 0 2 1 4 3 5 2 3 2 2 2 3 3 5 2 2 C 3 4 2 2 2 1 1 4 3 3 3 2 5 5 2 3 1 3 4 1 5 3 3 4 3 1 4 5 3 1 2 3 4 3 2 1 0 2 1 G 3 1 3 2 4 3 2 2 4 0 3 2 3 1 2 5 2 1 1 2 4 5 3 3 5 3 1 0 1 2 1 3 1 2 3 3 3 4 3 T 3 2 2 3 1 3 4 2 1 4 2 3 1 3 2 0 4 4 3 1 0 1 3 3 0 5 1 2 1 5 4 2 3 3 2 3 2 2 4 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 A 2 1 2 3 2 0 2 4 2 2 3 3 2 2 2 1 1 1 4 4 0 3 3 2 1 0 5 3 3 1 3 2 2 3 3 4 C 2 4 3 2 2 4 5 2 2 4 3 3 2 2 1 1 3 4 1 2 5 3 3 0 3 7 1 4 3 2 2 3 2 2 0 3 G 1 2 1 3 4 4 2 2 0 4 1 3 4 1 4 2 4 2 1 1 3 1 2 5 2 0 3 1 1 5 2 2 2 2 2 1 T 5 3 4 2 2 2 1 2 6 0 3 1 2 5 3 6 2 3 4 3 2 3 2 3 4 3 1 2 3 2 3 3 4 3 5 2 . . . ###################################################################################################3 Mas o que eu queria era o seguinte (... é que a linha continua): ###################################################################################################3 CCAAGATCGTCACCAGTTCACGCCGTACATTCCCAAAGTTCTAGCCATCAAGTGTGCAACAAGTCACAGACTATAAAGGACTT... A 1 3 3 3 3 3 3 2 2 3 2 3 1 1 4 2 3 2 2 6 1 1 1 0 2 1 4 3 5 2 3 2 2 2 3 3 5 2 2 2 1 2 3 2 0 2 4 2 2... C 3 4 2 2 2 1 1 4 3 3 3 2 5 5 2 3 1 3 4 1 5 3 3 4 3 1 4 5 3 1 2 3 4 3 2 1 0 2 1 2 1 2 3 2 0 2 4 2 2 ... G 3 1 3 2 4 3 2 2 4 0 3 2 3 1 2 5 2 1 1 2 4 5 3 3 5 3 1 0 1 2 1 3 1 2 3 3 3 4 3 1 2 1 3 4 4 2 2 0 4... T 3 2 2 3 1 3 4 2 1 4 2 3 1 3 2 0 4 4 3 1 0 1 3 3 0 5 1 2 1 5 4 2 3 3 2 3 2 2 4 5 3 4 2 2 2 1 2 6 0 ... ###################################################################################################3 Resumindo eu tenho duas duvidas: 01. Como escrever a matriz num arquivo txt sem quebrar linhas, independente do tamanho dela. Eu tentei mudar o argumento width no comando print e também no options() e ainda não consegui fazer a matrix ficar sem quebras no arquivo .txt. 02. Como excluir a linha de nome de colunas na hora de escrever o .txt, para ficar como no exemplo acima. Essa eu ainda nem tenho ideia de como proceder. Eu estou usando o sink() para escrever o arquivo txt pois pareceu o mais simples, não sei também se existem outras formar para escrever varias coisas juntas. Os dados estão no fim da msg, é uma matriz chamada resul e um vetor chamado consensus. Obrigado pela atenção: consensus<-c("C", "C", "A", "A", "G", "A", "T", "C", "G", "T", "C", "A", "C", "C", "A", "G", "T", "T", "C", "A", "C", "G", "C", "C", "G", "T", "A", "C", "A", "T", "T", "C", "C", "C", "A", "A", "A", "G", "T", "T", "C", "T", "A", "G", "C", "C", "A", "T", "C", "A", "A", "G", "T", "G", "T", "G", "C", "A", "A", "C", "A", "A", "G", "T", "C", "A", "C", "A", "G", "A", "C", "T", "A", "T", "A", "A", "A", "G", "G", "A", "C", "T", "T", "A", "T", "A", "T", "A", "A", "A", "G", "T", "A", "A", "A", "C", "C", "A", "G", "A", "T", "A", "A", "A", "C", "T", "A", "A", "C", "A", "G", "T", "T", "G", "C", "G", "G", "A", "G", "T", "C", "T", "C", "A", "C", "G", "A", "A", "C", "A", "G", "G", "A", "G", "C", "A", "A", "C", "G", "C", "C", "G", "A", "C", "A", "C", "A", "G", "A", "G", "A", "A", "C", "G", "C", "T", "G", "T", "G", "G", "T", "C", "G", "C", "G", "G", "A", "G", "C", "A", "C", "T", "C", "C", "G", "G", "G", "T", "C", "G", "C", "G", "A", "C", "A", "T", "C", "G", "A", "G", "G", "A", "T", "A", "C", "C", "G", "T", "A", "G", "C", "C", "A", "G", "C", "A", "A", "G", "C", "A", "C", "G", "C", "A", "T", "C", "A", "C", "G", "A", "C", "C", "A", "G", "G", "A", "A", "C", "T", "T", "A", "T", "C", "G", "G", "C", "A", "T", "A", "T", "T", "C", "A", "A", "A", "A", "G", "A", "C", "A", "T", "A", "A", "T", "G", "G", "G", "T", "G", "A", "G", "G", "G", "A", "C", "T", "T", "A", "T", "C", "G", "A", "A", "C", "A", "A", "C", "C", "G", "C", "G", "A", "C", "A", "T", "T", "G", "A", "G", "G", "A", "C", "G", "C", "G", "A", "T", "T", "A", "A", "A", "T", "A", "A", "C", "T", "A", "G", "G", "G", "T", "A", "A", "T", "G", "C", "C", "A", "G", "A", "T", "A", "T", "T", "A", "T", "G", "A", "A", "A", "G", "G", "C", "T", "T", "G", "T", "G", "C", "G", "A", "A", "T", "G", "T", "T", "C", "G", "A", "C", "G", "C", "T", "A", "G", "T", "G", "G", "T", "C", "A", "A", "A", "C", "C", "A", "C", "T", "C", "T", "T", "C", "C", "T", "T", "A", "A", "T", "A", "G", "T", "T", "G", "T", "G", "C", "C", "C", "G", "A", "G", "C", "T", "A", "C", "G", "A", "T", "C", "C", "A", "T", "G", "C", "T", "T", "C", "C", "G", "T", "C", "A", "C", "A", "C", "A", "G", "A", "A", "G", "T", "C", "A", "A", "A", "A", "A", "C", "T", "T", "T", "A", "G", "G", "C", "A", "G", "A", "G", "A", "G", "T", "G", "A", "A", "G", "G", "A", "A", "C", "G", "C", "C", "G", "C", "G", "C", "A", "T", "C", "C", "C", "C", "T", "T", "T", "A", "A", "A", "G", "C", "A", "G", "A", "G", "C", "A", "C", "G", "C", "G", "A", "C", "A", "A", "T", "C", "G", "G", "T", "C", "C", "C", "C", "A", "C", "T", "G", "C", "A", "T", "C", "C", "A", "T", "C", "G", "A", "C", "C", "A", "A", "C", "T", "C", "C", "A", "G", "C", "G", "G", "C", "C", "C", "A", "A", "C", "A", "C", "G", "G", "T", "C", "T", "A", "C", "G", "G", "C", "A", "A", "G", "G", "C", "A", "A", "T", "A", "G", "A", "A", "T", "G", "A", "T", "G", "G", "C", "A", "T", "T", "A", "A", "T", "G", "T", "G", "T", "A", "G", "G", "C", "C", "T", "C", "C", "T", "C", "C", "A", "C", "A", "G", "A", "C", "A", "A", "T", "C", "A", "G", "A", "G", "C", "A", "A", "A", "T", "A", "C", "A", "T", "C", "A", "C", "T", "T", "A", "A", "A", "C", "G", "C", "C", "T", "G", "A", "C", "G", "G", "G", "A", "C", "A", "G", "T", "G", "C", "C", "A", "A", "A", "A", "A", "A", "C", "G", "G", "G", "C", "C", "C", "T", "T", "G", "C", "A", "C", "G", "G", "C", "A", "A", "C", "C", "A", "A", "G", "T", "A", "A", "A", "C", "A", "C", "T", "G", "A", "G", "A", "T", "G", "C", "A", "T", "G", "A", "T", "A", "G", "A", "C", "C", "T", "C", "A", "T", "G", "A", "A", "A", "G", "G", "G", "A", "A", "G", "G", "G", "A", "C", "G", "G", "G", "T", "A", "A", "T", "T", "A", "T", "C", "A", "G", "A", "T", "G", "A", "G", "G", "G", "G", "C", "G", "C", "A", "A", "A", "C", "C", "C", "G", "G", "A", "T", "A", "A", "C", "T", "A", "A", "A", "G", "C", "G", "A", "A", "T", "G", "A", "C", "G", "A", "A", "T", "A", "T", "A", "A", "A", "G", "C", "C", "C", "C", "A", "A", "G", "C", "A", "C", "C", "G", "A", "G", "T", "A", "G", "T", "G", "A", "A", "T", "G", "C", "A", "G", "C", "C", "A", "A", "T", "T", "C", "T", "T", "G", "C", "A", "G", "A", "G", "T", "G", "G", "C", "C", "G", "C", "T", "C", "A", "G", "C", "A", "C", "C", "C", "A", "C", "G", "A", "T", "C", "C", "T", "T", "G", "G", "G", "T", "T", "A", "C", "T", "G", "A", "G", "C", "A", "T", "A", "G", "A", "G", "G", "A", "T", "G", "C", "G", "G", "C", "T", "T", "C", "A", "A", "T", "A", "T", "C", "C", "C", "A", "A", "G", "T", "A", "A", "A", "A", "A", "T", "A", "C", "T", "A", "G", "A", "A", "A", "A", "A", "T", "A", "A", "G", "C", "C", "G", "C", "A", "C", "A", "G", "A", "G", "A", "T", "C", "G", "T", "T", "T", "G", "G", "C", "G", "A", "G", "T", "G", "A", "C", "T", "G", "G", "T", "T", "A", "C", "G", "C", "A", "C", "C", "A", "T", "T", "G", "A", "A", "C", "A", "A", "A", "T", "C", "C", "C", "T", "A", "A", "G", "A", "A", "A", "C", "C", "C", "T", "A", "A", "A", "T", "C") resul<-structure(c(1, 3, 3, 3, 3, 4, 1, 2, 3, 2, 3, 2, 3, 2, 2, 3, 3, 2, 4, 1, 3, 1, 3, 3, 3, 1, 2, 4, 2, 4, 2, 2, 2, 3, 4, 1, 3, 3, 0, 4, 2, 3, 3, 2, 3, 2, 2, 3, 1, 5, 3, 1, 1, 5, 1, 3, 4, 2, 2, 2, 2, 3, 5, 0, 3, 1, 2, 4, 2, 3, 1, 4, 2, 4, 1, 3, 6, 1, 2, 1, 1, 5, 4, 0, 1, 3, 5, 1, 1, 3, 3, 3, 0, 4, 3, 3, 2, 3, 5, 0, 1, 1, 3, 5, 4, 4, 1, 1, 3, 5, 0, 2, 5, 3, 1, 1, 2, 1, 2, 5, 3, 2, 1, 4, 2, 3, 3, 2, 2, 4, 1, 3, 2, 3, 2, 3, 3, 2, 3, 2, 3, 1, 3, 3, 5, 0, 3, 2, 2, 2, 4, 2, 2, 1, 3, 4, 2, 2, 1, 5, 1, 4, 2, 3, 2, 3, 1, 4, 3, 2, 3, 2, 2, 2, 4, 2, 0, 4, 4, 2, 2, 5, 2, 1, 4, 2, 2, 2, 2, 2, 0, 6, 2, 4, 4, 0, 3, 3, 1, 3, 3, 3, 3, 1, 2, 2, 4, 2, 2, 2, 1, 5, 2, 1, 4, 3, 1, 1, 2, 6, 1, 3, 4, 2, 1, 4, 2, 3, 4, 1, 1, 4, 4, 2, 1, 3, 0, 5, 3, 2, 3, 3, 1, 3, 3, 3, 2, 2, 2, 0, 5, 3, 1, 3, 2, 4, 0, 7, 0, 3, 5, 1, 3, 1, 3, 4, 1, 2, 3, 3, 1, 3, 1, 2, 5, 2, 3, 2, 2, 3, 2, 3, 2, 3, 2, 2, 2, 4, 3, 2, 2, 3, 3, 0, 2, 5, 4, 3, 1, 2, 3, 3, 3, 1, 6, 0, 0, 4, 2, 1, 5, 2, 1, 2, 5, 2, 3, 2, 3, 2, 2, 4, 3, 1, 2, 2, 2, 4, 1, 2, 3, 4, 4, 1, 2, 3, 1, 2, 3, 4, 4, 1, 3, 2, 0, 1, 3, 6, 3, 2, 3, 2, 5, 3, 0, 2, 3, 2, 3, 2, 3, 2, 4, 1, 2, 1, 3, 4, 4, 1, 3, 2, 3, 3, 3, 1, 4, 3, 2, 1, 1, 3, 3, 3, 2, 6, 2, 0, 3, 2, 2, 3, 3, 2, 4, 1, 4, 3, 1, 2, 0, 2, 3, 5, 4, 3, 1, 2, 3, 1, 3, 3, 3, 2, 2, 3, 1, 3, 3, 3, 2, 1, 2, 5, 4, 3, 1, 2, 3, 3, 3, 1, 3, 4, 1, 2, 3, 2, 3, 2, 3, 0, 5, 2, 2, 2, 2, 4, 2, 2, 1, 5, 3, 2, 4, 1, 2, 3, 2, 3, 2, 2, 3, 3, 1, 1, 4, 4, 5, 1, 1, 3, 0, 2, 4, 4, 2, 2, 1, 5, 1, 4, 4, 1, 1, 3, 2, 4, 3, 5, 1, 1, 4, 3, 3, 0, 1, 6, 1, 2, 3, 0, 4, 3, 4, 2, 4, 0, 3, 3, 1, 3, 2, 3, 2, 3, 5, 2, 0, 3, 2, 2, 5, 1, 2, 2, 3, 3, 3, 2, 2, 3, 2, 3, 4, 1, 2, 3, 3, 2, 4, 3, 2, 1, 4, 1, 1, 4, 1, 4, 3, 2, 2, 3, 5, 0, 3, 5, 2, 0, 1, 5, 2, 2, 2, 2, 4, 2, 3, 1, 3, 3, 1, 4, 3, 2, 4, 3, 2, 1, 2, 3, 3, 2, 3, 2, 3, 2, 3, 3, 4, 0, 4, 2, 2, 2, 2, 2, 3, 3, 3, 1, 3, 3, 4, 3, 3, 0, 2, 3, 2, 3, 3, 1, 6, 0, 2, 4, 3, 1, 2, 3, 1, 4, 3, 1, 4, 2, 2, 2, 2, 4, 3, 1, 4, 2, 3, 3, 4, 0, 2, 1, 1, 6, 2, 4, 4, 0, 1, 3, 5, 1, 2, 4, 2, 2, 0, 2, 7, 1, 0, 2, 6, 2, 3, 3, 3, 1, 2, 3, 4, 1, 2, 6, 0, 2, 4, 1, 3, 2, 4, 5, 0, 1, 1, 2, 3, 4, 0, 4, 4, 2, 2, 4, 2, 2, 2, 2, 4, 2, 1, 3, 5, 1, 2, 2, 3, 3, 2, 1, 3, 4, 2, 5, 3, 0, 1, 2, 4, 3, 2, 5, 2, 1, 2, 2, 3, 3, 3, 2, 3, 2, 1, 4, 1, 4, 4, 2, 3, 1, 2, 3, 1, 4, 2, 4, 3, 1, 2, 0, 4, 4, 3, 3, 1, 3, 3, 3, 4, 0, 2, 2, 3, 3, 4, 3, 2, 1, 3, 1, 2, 4, 3, 2, 3, 2, 1, 4, 3, 2, 1, 5, 3, 1, 1, 1, 7, 1, 2, 1, 3, 4, 5, 2, 0, 3, 1, 2, 5, 2, 1, 5, 2, 2, 2, 4, 3, 1, 3, 2, 2, 3, 1, 1, 5, 3, 2, 3, 2, 3, 3, 2, 2, 3, 4, 4, 0, 2, 2, 2, 4, 2, 1, 4, 3, 2, 4, 0, 2, 4, 1, 3, 3, 3, 2, 2, 3, 3, 2, 3, 2, 3, 4, 1, 2, 3, 0, 3, 3, 4, 2, 5, 1, 2, 4, 2, 2, 2, 2, 4, 3, 1, 2, 2, 4, 2, 4, 1, 3, 2, 2, 4, 3, 1, 2, 3, 3, 2, 4, 0, 4, 2, 1, 3, 5, 1, 2, 2, 4, 2, 6, 1, 0, 3, 4, 2, 2, 2, 2, 3, 2, 3, 2, 2, 1, 5, 2, 3, 1, 4, 4, 1, 3, 2, 2, 2, 2, 4, 3, 4, 2, 1, 2, 2, 4, 2, 3, 3, 4, 0, 1, 4, 1, 4, 4, 2, 2, 2, 2, 1, 2, 5, 3, 1, 3, 3, 2, 2, 1, 5, 2, 2, 2, 4, 0, 6, 1, 3, 3, 2, 3, 2, 4, 3, 1, 2, 5, 1, 0, 4, 3, 3, 2, 2, 3, 2, 4, 1, 4, 2, 2, 2, 1, 4, 3, 2, 3, 3, 1, 3, 1, 3, 2, 4, 4, 1, 4, 1, 4, 2, 2, 2, 2, 3, 0, 5, 2, 2, 3, 3, 3, 3, 4, 0, 1, 3, 5, 1, 3, 3, 0, 4, 0, 2, 4, 4, 4, 4, 2, 0, 2, 1, 4, 3, 1, 1, 5, 3, 3, 2, 5, 0, 5, 3, 1, 1, 1, 4, 1, 4, 1, 3, 2, 4, 2, 1, 3, 4, 4, 2, 1, 3, 3, 1, 2, 4, 2, 4, 1, 3, 2, 2, 4, 2, 4, 2, 1, 3, 4, 2, 1, 3, 1, 4, 3, 2, 4, 4, 0, 2, 5, 1, 2, 2, 3, 5, 0, 2, 1, 6, 1, 2, 1, 2, 4, 3, 1, 4, 2, 3, 0, 3, 5, 2, 4, 2, 1, 3, 2, 3, 2, 3, 3, 3, 2, 2, 3, 1, 2, 4, 0, 2, 3, 5, 2, 2, 4, 2, 5, 2, 1, 2, 1, 2, 5, 2, 1, 4, 5, 0, 4, 2, 2, 2, 2, 4, 1, 3, 3, 2, 4, 1, 1, 6, 1, 2, 1, 2, 4, 3, 3, 2, 3, 2, 3, 0, 3, 4, 2, 3, 1, 4, 4, 1, 4, 1, 4, 2, 4, 0, 5, 1, 1, 3, 2, 3, 0, 5, 6, 2, 1, 1, 3, 3, 1, 3, 1, 4, 1, 4, 2, 3, 1, 4, 5, 0, 3, 2, 2, 1, 5, 2, 0, 2, 6, 2, 2, 3, 4, 1, 3, 1, 1, 5, 6, 2, 1, 1, 4, 1, 4, 1, 0, 2, 3, 5, 2, 1, 4, 3, 2, 3, 2, 3, 1, 4, 2, 3, 3, 1, 3, 3, 2, 0, 6, 2, 4, 1, 3, 2, 3, 0, 2, 5, 3, 3, 3, 1, 2, 3, 1, 4, 1, 3, 1, 5, 5, 1, 2, 2, 0, 1, 3, 6, 1, 1, 8, 0, 5, 1, 1, 3, 3, 2, 3, 2, 3, 2, 3, 2, 1, 3, 5, 1, 1, 3, 4, 2, 0, 4, 3, 3, 3, 1, 2, 4, 0, 3, 2, 5, 2, 2, 4, 2, 1, 3, 0, 6, 1, 0, 5, 4, 3, 4, 2, 1, 3, 2, 4, 1, 3, 2, 3, 2, 3, 3, 2, 2, 2, 3, 1, 4, 2, 3, 4, 1, 2, 3, 1, 4, 3, 1, 1, 5, 3, 4, 2, 1, 2, 1, 4, 3, 3, 2, 3, 2, 2, 3, 2, 3, 2, 2, 4, 2, 3, 4, 2, 1, 2, 1, 3, 4, 5, 2, 3, 0, 2, 2, 4, 2, 2, 1, 3, 4, 2, 2, 3, 3, 4, 0, 5, 1, 1, 2, 2, 5, 2, 4, 2, 2, 4, 0, 3, 3, 4, 1, 4, 1, 5, 0, 1, 4, 1, 5, 2, 2, 0, 6, 2, 2, 3, 3, 2, 2, 1, 6, 2, 1, 1, 1, 3, 5, 3, 5, 1, 1, 3, 3, 0, 4, 2, 3, 1, 4, 3, 4, 1, 2, 1, 4, 2, 3, 2, 2, 2, 4, 1, 3, 2, 4, 4, 3, 1, 2, 3, 3, 2, 2, 2, 2, 2, 4, 3, 2, 3, 2, 1, 1, 4, 4, 0, 4, 1, 5, 3, 2, 0, 5, 2, 1, 5, 2, 1, 2, 2, 5, 2, 2, 3, 3, 1, 4, 2, 3, 2, 4, 3, 1, 2, 4, 1, 3, 0, 2, 6, 2, 5, 3, 1, 1, 2, 2, 3, 3, 2, 4, 2, 2, 0, 1, 2, 7, 4, 3, 2, 1, 4, 5, 1, 0, 2, 1, 4, 3, 7, 2, 1, 0, 1, 1, 2, 6, 1, 4, 3, 2, 2, 5, 3, 0, 4, 2, 1, 3, 2, 2, 2, 4, 3, 1, 6, 0, 2, 3, 2, 3, 1, 3, 2, 4, 2, 1, 2, 5, 1, 4, 2, 3, 2, 5, 0, 3, 1, 1, 4, 4, 1, 3, 2, 4, 2, 3, 2, 3, 3, 2, 2, 3, 2, 3, 2, 3, 4, 1, 4, 1, 2, 3, 2, 3, 4, 3, 2, 1, 2, 1, 5, 2, 3, 3, 1, 3, 3, 2, 2, 3, 2, 3, 4, 1, 2, 3, 1, 4, 3, 4, 1, 2, 4, 1, 3, 2, 3, 3, 1, 3, 5, 0, 3, 2, 3, 3, 1, 3, 4, 1, 2, 3, 2, 3, 3, 2, 2, 2, 2, 4, 4, 1, 0, 5, 3, 0, 3, 4, 5, 1, 2, 2, 1, 1, 5, 3, 3, 1, 5, 1, 2, 3, 2, 3, 4, 2, 2, 2, 0, 3, 5, 2, 3, 2, 3, 2, 2, 3, 5, 0, 3, 3, 2, 2, 2, 2, 3, 3, 2, 1, 1, 6, 3, 2, 4, 1, 7, 1, 1, 1, 5, 2, 1, 2, 2, 2, 3, 3, 1, 2, 5, 2, 4, 1, 1, 4, 4, 1, 3, 2, 1, 5, 2, 2, 2, 2, 3, 3, 2, 4, 3, 1, 1, 5, 2, 2, 3, 1, 5, 1, 3, 5, 2, 0, 3, 1, 6, 0, 0, 5, 3, 2, 4, 3, 1, 2, 0, 4, 1, 5, 0, 4, 3, 3, 0, 5, 2, 3, 1, 5, 3, 1, 1, 4, 2, 3, 3, 1, 2, 4, 1, 2, 0, 7, 1, 3, 2, 4, 4, 1, 3, 2, 3, 3, 2, 2, 5, 2, 2, 1, 3, 1, 6, 0, 2, 4, 1, 3, 3, 2, 2, 3, 2, 1, 5, 2, 3, 1, 3, 3, 2, 1, 4, 3, 2, 3, 2, 3, 4, 2, 1, 3, 3, 4, 2, 1, 2, 1, 5, 2, 2, 4, 3, 1, 2, 1, 4, 3, 4, 1, 3, 2, 2, 3, 3, 2, 3, 2, 3, 2, 3, 3, 1, 3, 3, 2, 1, 4, 2, 5, 1, 2, 0, 2, 5, 3, 2, 1, 5, 2, 1, 1, 2, 6, 1, 3, 3, 3, 1, 5, 1, 3, 2, 3, 3, 2, 3, 4, 0, 3, 6, 3, 1, 0, 2, 5, 1, 2, 2, 2, 1, 5, 2, 0, 7, 1, 1, 4, 2, 3, 4, 1, 2, 3, 1, 2, 3, 4, 3, 4, 1, 2, 1, 3, 3, 3, 3, 2, 3, 2, 1, 2, 3, 4, 2, 4, 4, 0, 2, 2, 3, 3, 3, 2, 3, 2, 2, 3, 2, 3, 3, 4, 1, 2, 4, 3, 2, 1, 3, 1, 3, 3, 2, 3, 2, 3, 1, 2, 3, 4, 2, 3, 3, 2, 2, 3, 3, 2, 3, 1, 3, 3, 3, 2, 5, 0, 1, 3, 3, 3, 3, 1, 4, 2, 2, 2, 3, 3, 2, 5, 2, 1, 3, 4, 1, 2, 1, 4, 2, 3, 5, 0, 2, 3, 4, 3, 1, 2, 2, 4, 3, 1, 3, 2, 2, 3, 0, 4, 3, 3, 2, 0, 4, 4, 3, 2, 5, 0, 3, 1, 2, 4, 1, 5, 3, 1, 0, 3, 3, 4, 3, 2, 3, 2, 2, 3, 3, 2, 2, 2, 4, 2, 2, 2, 6, 0, 2, 3, 2, 3, 4, 3, 1, 2, 3, 2, 3, 2, 1, 3, 5, 1, 1, 3, 4, 2, 0, 4, 2, 4, 4, 2, 3, 1, 3, 3, 1, 3, 0, 2, 3, 5, 4, 1, 2, 3, 1, 3, 4, 2, 4, 1, 1, 4, 4, 3, 1, 2, 1, 3, 2, 4, 3, 2, 4, 1, 3, 3, 1, 3, 1, 2, 3, 4, 2, 3, 5, 0, 2, 0, 4, 4, 1, 3, 3, 3, 3, 1, 3, 3, 1, 2, 2, 5, 2, 0, 1, 7, 5, 3, 2, 0, 6, 0, 1, 3, 1, 3, 2, 4, 1, 1, 4, 4, 2, 2, 2, 4, 2, 2, 4, 2, 3, 0, 3, 4, 3, 3, 3, 1, 1, 1, 4, 4, 2, 1, 7, 0, 3, 4, 1, 2, 2, 5, 2, 1, 1, 3, 2, 4, 4, 5, 1, 0, 2, 4, 1, 3, 3, 1, 1, 5, 2, 5, 1, 2, 2, 4, 3, 1, 3, 3, 3, 1, 2, 3, 3, 2, 3, 2, 2, 3, 1, 2, 4, 3, 4, 2, 2, 2, 2, 3, 3, 2, 5, 2, 2, 1, 6, 0, 4, 0, 2, 3, 1, 4, 2, 3, 3, 2, 5, 2, 2, 1, 2, 2, 5, 1, 3, 2, 3, 2, 2, 3, 5, 0, 3, 4, 2, 1, 4, 3, 3, 0, 4, 2, 1, 3, 4, 2, 4, 0, 2, 2, 2, 4, 4, 1, 2, 3, 2, 4, 2, 2, 4, 0, 2, 4, 2, 2, 2, 4, 2, 4, 2, 2, 3, 2, 2, 3, 2, 3, 3, 2, 2, 2, 2, 4, 1, 3, 2, 4, 4, 1, 4, 1, 4, 4, 1, 1, 4, 3, 1, 2, 2, 4, 2, 2, 2, 2, 3, 3, 1, 5, 2, 2, 1, 4, 4, 1, 1, 2, 3, 4, 0, 1, 7, 2, 4, 2, 1, 3, 2, 3, 2, 3, 2, 2, 3, 3, 2, 2, 3, 3, 2, 1, 4, 3, 4, 1, 3, 2, 3, 4, 1, 2, 3, 3, 1, 3, 2, 2, 4, 2, 3, 1, 2, 4, 3, 3, 4, 0, 2, 6, 1, 1, 1, 4, 2, 3, 4, 0, 4, 2, 6, 1, 2, 1, 3, 2, 2, 3, 5, 1, 2, 2, 4, 4, 2, 0, 3, 1, 3, 3, 2, 3, 2, 3, 2, 2, 6, 0, 3, 2, 4, 1, 0, 2, 4, 4, 1, 5, 2, 2, 3, 4, 3, 0, 2, 4, 2, 2, 1, 2, 3, 4, 3, 2, 1, 4, 1, 3, 4, 2, 2, 4, 2, 2, 3, 3, 2, 2, 2, 4, 0, 4, 2, 2, 3, 3, 2, 1, 6, 1, 2, 4, 2, 2, 4, 3, 3, 0, 4, 0, 2, 4, 1, 3, 3, 3, 1, 7, 1, 1, 4, 2, 2, 2, 5, 2, 3, 0, 3, 2, 5, 0, 2, 3, 1, 4, 4, 3, 1, 2, 4, 3, 1, 2, 3, 2, 2, 3, 2, 5, 3, 0, 5, 4, 0, 1, 4, 5, 1, 0, 1, 2, 2, 5, 2, 3, 4, 1, 4, 3, 0, 3, 1, 3, 4, 2, 5, 2, 3, 0, 1, 3, 1, 5, 1, 2, 4, 3, 2, 4, 3, 1, 3, 3, 3, 1, 1, 2, 3, 4, 2, 1, 5, 2, 3, 3, 1, 3, 0, 3, 3, 4, 4, 3, 1, 2, 2, 2, 4, 2, 4, 2, 2, 2, 2, 5, 1, 2, 3, 4, 1, 2, 1, 2, 3, 4, 3, 5, 1, 1, 4, 1, 3, 2, 1, 2, 3, 4, 0, 4, 5, 1, 3, 3, 2, 2, 3, 3, 2, 2, 5, 2, 2, 1, 1, 2, 5, 2, 3, 1, 5, 1, 2, 3, 4, 1, 3, 2, 2, 3, 4, 3, 1, 2, 3, 2, 4, 1, 2, 2, 3, 3, 2, 0, 5, 3, 4, 1, 3, 2, 1, 3, 3, 3, 0, 3, 4, 3, 2, 1, 4, 3, 3, 1, 4, 2, 1, 1, 2, 6, 4, 3, 3, 0, 5, 1, 0, 4, 2, 0, 3, 5, 3, 1, 2, 4, 4, 2, 2, 2, 2, 2, 1, 5, 2, 4, 1, 3, 3, 3, 2, 2, 2, 2, 4, 2, 4, 1, 2, 3, 3, 0, 2, 5, 2, 2, 3, 3, 5, 0, 4, 1, 2, 2, 5, 1, 1, 3, 5, 1, 2, 2, 4, 2, 2, 2, 3, 3, 2, 3, 2, 3, 2, 2, 3, 3, 1, 4, 3, 2, 5, 0, 3, 2, 3, 3, 2, 2, 5, 3, 1, 1, 3, 4, 2, 1, 1, 3, 3, 3, 1, 4, 3, 2, 2, 2, 4, 2, 2, 2, 4, 2, 4, 2, 4, 0, 2, 2, 1, 5, 4, 3, 1, 2, 6, 3, 0, 1, 1, 5, 1, 3, 3, 3, 0, 4, 4, 1, 1, 4, 5, 1, 3, 1, 3, 2, 3, 2, 2, 2, 3, 3, 0, 5, 4, 1, 1, 1, 6, 2, 4, 2, 2, 2, 4, 3, 1, 2, 2, 1, 2, 5, 1, 2, 5, 2, 5, 2, 0, 3, 1, 6, 2, 1, 3, 1, 4, 2, 4, 1, 1, 4, 5, 2, 2, 1, 3, 2, 1, 4, 4, 2, 3, 1, 2, 2, 2, 4, 3, 3, 2, 2, 5, 0, 3, 2, 3, 2, 2, 3, 3, 2, 4, 1, 2, 5, 0, 3, 1, 3, 3, 3, 0, 7, 2, 1, 1, 4, 3, 2, 3, 3, 1, 3, 4, 2, 3, 1, 2, 2, 3, 3, 1, 5, 3, 1, 4, 2, 3, 1, 2, 4, 1, 3, 0, 6, 2, 2, 0, 3, 5, 2, 5, 0, 4, 1, 2, 2, 4, 2, 3, 2, 1, 4, 3, 2, 3, 2, 1, 3, 4, 2, 3, 2, 1, 4, 2, 1, 4, 3, 4, 1, 1, 4, 3, 3, 2, 2, 1, 2, 3, 4, 2, 2, 3, 3, 1, 5, 1, 3, 5, 2, 2, 1, 3, 2, 4, 1, 2, 3, 3, 2, 3, 5, 1, 1, 4, 3, 2, 1, 5, 1, 1, 3, 2, 2, 2, 4, 2, 2, 2, 4, 2, 5, 3, 0, 2, 1, 3, 4, 2, 2, 1, 5, 1, 4, 5, 0, 2, 3, 3, 2, 4, 3, 2, 1, 2, 3, 4, 1, 6, 2, 2, 0, 0, 2, 6, 2, 1, 2, 2, 5, 3, 0, 4, 3, 1, 2, 4, 3, 1, 5, 3, 1, 2, 4, 2, 2, 0, 3, 5, 2, 2, 3, 3, 2, 3, 2, 1, 4, 0, 5, 2, 3, 3, 3, 2, 2, 2, 1, 4, 3, 2, 4, 2, 2, 3, 3, 3, 1, 2, 4, 3, 1, 1, 3, 3, 3, 1, 5, 1, 3, 4, 1, 4, 1, 2, 3, 3, 2, 2, 2, 5, 1, 3, 3, 1, 3, 3, 2, 1, 4, 1, 4, 2, 3, 3, 4, 3, 0, 1, 3, 2, 4, 2, 3, 0, 5, 4, 1, 5, 0, 3, 2, 4, 1, 3, 1, 4, 2, 1, 2, 3, 4, 3, 2, 1, 4, 3, 3, 1, 3, 2, 3, 2, 3, 2, 2, 2, 4, 2, 1, 4, 3, 4, 2, 2, 2, 2, 2, 3, 3, 2, 4, 1, 3, 4, 1, 2, 3, 1, 3, 2, 4, 4, 3, 2, 1, 2, 2, 3, 3, 3, 2, 2, 3, 0, 3, 5, 2, 1, 2, 4, 3, 5, 3, 1, 1, 0, 1, 3, 6, 2, 2, 3, 3, 2, 3, 3, 2, 2, 2, 3, 3, 1, 2, 4, 3, 2, 4, 2, 2, 3, 2, 1, 4, 2, 1, 0, 7, 2, 4, 0, 4, 3, 3, 2, 2, 4, 1, 1, 4, 2, 2, 2, 4, 5, 3, 0, 2, 2, 1, 3, 4, 2, 3, 2, 3, 2, 4, 1, 3, 3, 5, 2, 0, 5, 0, 3, 2, 3, 2, 3, 2, 3, 2, 4, 1, 2, 1, 2, 5, 4, 0, 4, 2, 5, 2, 2, 1, 4, 2, 2, 2, 3, 3, 1, 3, 4, 1, 2, 3, 0, 1, 4, 5, 4, 1, 2, 3, 0, 6, 2, 2, 0, 3, 3, 4, 5, 0, 2, 3, 3, 0, 5, 2, 3, 3, 3, 1, 3, 2, 3, 2, 4, 1, 2, 3, 4, 2, 2, 2, 3, 3, 2, 2, 1, 1, 2, 6, 3, 1, 3, 3, 4, 1, 4, 1, 1, 3, 4, 2, 1, 4, 2, 3, 1, 4, 1, 4, 1, 2, 4, 3, 2, 3, 2, 3, 4, 4, 2, 0, 2, 4, 1, 3, 5, 3, 0, 2, 2, 2, 4, 2, 3, 3, 3, 1, 2, 2, 3, 3, 3, 3, 2, 2, 3, 0, 2, 5, 0, 4, 2, 4, 3, 1, 5, 1, 2, 3, 1, 4, 0, 4, 1, 5, 2, 1, 3, 4, 2, 2, 5, 1, 1, 1, 4, 4, 1, 4, 1, 4, 2, 2, 5, 1, 3, 3, 1, 3, 1, 2, 4, 3, 3, 2, 1, 4, 2, 2, 4, 2, 5, 2, 2, 1, 2, 3, 3, 2, 3, 1, 2, 4, 2, 1, 5, 2, 3, 0, 5, 2, 3, 1, 1, 5, 3, 2, 1, 4, 4, 1, 3, 2, 2, 4, 3, 1, 2, 1, 5, 2, 2, 3, 3, 2, 4, 1, 1, 4, 2, 6, 2, 0, 1, 5, 2, 2, 5, 1, 1, 3, 3, 1, 1, 5, 2, 3, 1, 4, 1, 1, 5, 3, 4, 2, 1, 3, 3, 1, 3, 3, 1, 3, 3, 3, 5, 2, 1, 2, 4, 1, 4, 1, 4, 2, 2, 2, 2, 3, 1, 4, 1, 4, 1, 4, 3, 4, 2, 1, 1, 4, 3, 2, 2, 1, 3, 4, 4, 3, 0, 3, 5, 1, 0, 4, 3, 2, 5, 0, 4, 0, 2, 4, 3, 2, 3, 2, 3, 3, 3, 1, 2, 3, 2, 3, 4, 5, 1, 0, 1, 4, 3, 2, 2, 1, 0, 7, 3, 3, 1, 3, 3, 3, 3, 1, 4, 3, 1, 2, 3, 1, 2, 4, 2, 5, 3, 0), .Dim = c(4L, 957L), .Dimnames = list( c("A", "C", "G", "T"), c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", "26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", "37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", "48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", "59", "60", "61", "62", "63", "64", "65", "66", "67", "68", "69", "70", "71", "72", "73", "74", "75", "76", "77", "78", "79", "80", "81", "82", "83", "84", "85", "86", "87", "88", "89", "90", "91", "92", "93", "94", "95", "96", "97", "98", "99", "100", "101", "102", "103", "104", "105", "106", "107", "108", "109", "110", "111", "112", "113", "114", "115", "116", "117", "118", "119", "120", "121", "122", "123", "124", "125", "126", "127", "128", "129", "130", "131", "132", "133", "134", "135", "136", "137", "138", "139", "140", "141", "142", "143", "144", "145", "146", "147", "148", "149", "150", "151", "152", "153", "154", "155", "156", "157", "158", "159", "160", "161", "162", "163", "164", "165", "166", "167", "168", "169", "170", "171", "172", "173", "174", "175", "176", "177", "178", "179", "180", "181", "182", "183", "184", "185", "186", "187", "188", "189", "190", "191", "192", "193", "194", "195", "196", "197", "198", "199", "200", "201", "202", "203", "204", "205", "206", "207", "208", "209", "210", "211", "212", "213", "214", "215", "216", "217", "218", "219", "220", "221", "222", "223", "224", "225", "226", "227", "228", "229", "230", "231", "232", "233", "234", "235", "236", "237", "238", "239", "240", "241", "242", "243", "244", "245", "246", "247", "248", "249", "250", "251", "252", "253", "254", "255", "256", "257", "258", "259", "260", "261", "262", "263", "264", "265", "266", "267", "268", "269", "270", "271", "272", "273", "274", "275", "276", "277", "278", "279", "280", "281", "282", "283", "284", "285", "286", "287", "288", "289", "290", "291", "292", "293", "294", "295", "296", "297", "298", "299", "300", "301", "302", "303", "304", "305", "306", "307", "308", "309", "310", "311", "312", "313", "314", "315", "316", "317", "318", "319", "320", "321", "322", "323", "324", "325", "326", "327", "328", "329", "330", "331", "332", "333", "334", "335", "336", "337", "338", "339", "340", "341", "342", "343", "344", "345", "346", "347", "348", "349", "350", "351", "352", "353", "354", "355", "356", "357", "358", "359", "360", "361", "362", "363", "364", "365", "366", "367", "368", "369", "370", "371", "372", "373", "374", "375", "376", "377", "378", "379", "380", "381", "382", "383", "384", "385", "386", "387", "388", "389", "390", "391", "392", "393", "394", "395", "396", "397", "398", "399", "400", "401", "402", "403", "404", "405", "406", "407", "408", "409", "410", "411", "412", "413", "414", "415", "416", "417", "418", "419", "420", "421", "422", "423", "424", "425", "426", "427", "428", "429", "430", "431", "432", "433", "434", "435", "436", "437", "438", "439", "440", "441", "442", "443", "444", "445", "446", "447", "448", "449", "450", "451", "452", "453", "454", "455", "456", "457", "458", "459", "460", "461", "462", "463", "464", "465", "466", "467", "468", "469", "470", "471", "472", "473", "474", "475", "476", "477", "478", "479", "480", "481", "482", "483", "484", "485", "486", "487", "488", "489", "490", "491", "492", "493", "494", "495", "496", "497", "498", "499", "500", "501", "502", "503", "504", "505", "506", "507", "508", "509", "510", "511", "512", "513", "514", "515", "516", "517", "518", "519", "520", "521", "522", "523", "524", "525", "526", "527", "528", "529", "530", "531", "532", "533", "534", "535", "536", "537", "538", "539", "540", "541", "542", "543", "544", "545", "546", "547", "548", "549", "550", "551", "552", "553", "554", "555", "556", "557", "558", "559", "560", "561", "562", "563", "564", "565", "566", "567", "568", "569", "570", "571", "572", "573", "574", "575", "576", "577", "578", "579", "580", "581", "582", "583", "584", "585", "586", "587", "588", "589", "590", "591", "592", "593", "594", "595", "596", "597", "598", "599", "600", "601", "602", "603", "604", "605", "606", "607", "608", "609", "610", "611", "612", "613", "614", "615", "616", "617", "618", "619", "620", "621", "622", "623", "624", "625", "626", "627", "628", "629", "630", "631", "632", "633", "634", "635", "636", "637", "638", "639", "640", "641", "642", "643", "644", "645", "646", "647", "648", "649", "650", "651", "652", "653", "654", "655", "656", "657", "658", "659", "660", "661", "662", "663", "664", "665", "666", "667", "668", "669", "670", "671", "672", "673", "674", "675", "676", "677", "678", "679", "680", "681", "682", "683", "684", "685", "686", "687", "688", "689", "690", "691", "692", "693", "694", "695", "696", "697", "698", "699", "700", "701", "702", "703", "704", "705", "706", "707", "708", "709", "710", "711", "712", "713", "714", "715", "716", "717", "718", "719", "720", "721", "722", "723", "724", "725", "726", "727", "728", "729", "730", "731", "732", "733", "734", "735", "736", "737", "738", "739", "740", "741", "742", "743", "744", "745", "746", "747", "748", "749", "750", "751", "752", "753", "754", "755", "756", "757", "758", "759", "760", "761", "762", "763", "764", "765", "766", "767", "768", "769", "770", "771", "772", "773", "774", "775", "776", "777", "778", "779", "780", "781", "782", "783", "784", "785", "786", "787", "788", "789", "790", "791", "792", "793", "794", "795", "796", "797", "798", "799", "800", "801", "802", "803", "804", "805", "806", "807", "808", "809", "810", "811", "812", "813", "814", "815", "816", "817", "818", "819", "820", "821", "822", "823", "824", "825", "826", "827", "828", "829", "830", "831", "832", "833", "834", "835", "836", "837", "838", "839", "840", "841", "842", "843", "844", "845", "846", "847", "848", "849", "850", "851", "852", "853", "854", "855", "856", "857", "858", "859", "860", "861", "862", "863", "864", "865", "866", "867", "868", "869", "870", "871", "872", "873", "874", "875", "876", "877", "878", "879", "880", "881", "882", "883", "884", "885", "886", "887", "888", "889", "890", "891", "892", "893", "894", "895", "896", "897", "898", "899", "900", "901", "902", "903", "904", "905", "906", "907", "908", "909", "910", "911", "912", "913", "914", "915", "916", "917", "918", "919", "920", "921", "922", "923", "924", "925", "926", "927", "928", "929", "930", "931", "932", "933", "934", "935", "936", "937", "938", "939", "940", "941", "942", "943", "944", "945", "946", "947", "948", "949", "950", "951", "952", "953", "954", "955", "956", "957"))) -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

2 1

Estimativa da Variância Global - Krigagem (krige.conv)
by Rodrigo Sant'Ana 03 Jan '13

03 Jan '13

Boa noite Pessoal, minha dúvida é mais teórica do que prática. Por isso não apresento aqui nenhum CMR. Mas acredito que vocês possam me ajudar. Estou com uma dúvida na estimativa da variância global resultante do processo de krigagem por intermédio da função krige.conv() prevista no pacote geoR. Poderia considerar que a variância total da minha predição (krige) seria o simples somatório das variâncias preditas para cada ponto do meu grid de predição? Desculpem-me pela pergunta teórica, mas é uma dúvida consideravelmente importante para o meu trabalho, e estou tendo dificuldades em esclarecer essa questão. Muito obrigado pela atenção de todos. Cordialmente, _______________ Rodrigo Sant'Ana Oceanógrafo

3 3

Re: [R-br] eliminando linhas específicas
by Humberto Hazin 03 Jan '13

03 Jan '13

Olá Fátima basta apenas você digitar o comando banco[!duplicated(banco), ] ou unique(banco) Humberto Enviado do Email do Windows De: Fátima Lima Paula Enviado: ‎3‎ de ‎janeiro‎ de ‎2013 ‎13‎:‎48 Para: r-br(a)listas.c3sl.ufpr.br Assunto: [R-br] eliminando linhas específicas Prezados, aproveito para desejar um ótimo ano a todos e agradecer a ajuda do ano que passou. Já começo o ano precisando de vocês. Estou com uma tabela que possui dados repetidos de indivíduos, como o exemplo: id=c("1","2","2","2","3","4","4") nome=c("Maria","Jose","Jose","Jose","Paulo","Malu","Malu") idade=c("65","70","70","70","68","80","80") banco=data.frame(id,nome,idade) banco Quero eliminar as linhas repetidas e deixar apenas uma de cada tipo, tendo como referência "id", de forma que fique assim: id nome idade 1 1 Maria 65 2 2 Jose 70 3 3 Paulo 68 4 4 Malu 80 Alguém pode me ajudar, por favor. Obrigada Fátima

2 1

eliminando linhas específicas
by Fátima Lima Paula 03 Jan '13

03 Jan '13

Prezados, aproveito para desejar um ótimo ano a todos e agradecer a ajuda do ano que passou. Já começo o ano precisando de vocês. Estou com uma tabela que possui dados repetidos de indivíduos, como o exemplo: id=c("1","2","2","2","3","4","4") nome=c("Maria","Jose","Jose","Jose","Paulo","Malu","Malu") idade=c("65","70","70","70","68","80","80") banco=data.frame(id,nome,idade) banco Quero eliminar as linhas repetidas e deixar apenas uma de cada tipo, tendo como referência "id", de forma que fique assim: id nome idade 1 1 Maria 65 2 2 Jose 70 3 3 Paulo 68 4 4 Malu 80 Alguém pode me ajudar, por favor. Obrigada Fátima

1 0

Mudar cor da variável em plot
by Sérgio Henrique almeida da silva ju 03 Jan '13

03 Jan '13

Olá Tem como eu mudar a cor das variáveis plotadas em um gráfico? exemplo: exemplo: CAT=c(*rep*('Categoria I',60),*rep*('Categoria II',30),rep('Categoria III',10)) SEXO=c(rep('Feminino',28),rep('Masculino',32),rep('Feminino',20),rep('Masculino',10),rep('Feminino',7),rep('Masculino',3)) dados = cbind(as.data.frame(CAT),as.data.frame(SEXO)) grafico =table(dados$CAT,dados$SEXO) plot(grafico) Nesse caso tem como por exemplo minha variável sexo ser azul e a variável categoria ser vermelha? Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383

4 5