Outubro 2014 - R-br

Transformação reversa
by Felinto COSTA 29 Out '14

29 Out '14

Prezados. Alguém já implementou uma rotina para realizar a transformação reversa dos resultados de uma krigagem ordinária cujos dados originais foram transformados, cf. o Prof. Jorge Yamamoto propôs?

2 2

Submissão ao CRAN
by Waldelene Moura 29 Out '14

29 Out '14

Olá Pessoal Gostaria de submeter ao CRAN, um pacote que implementei no meu mestrado. Alguém que já domine o assunto pode me dar umas dicas de como faze-lo para diminuir a possibilidade de erros? Tenho algumas dúvidas, por exemplo: 1. O manual de documentação necessariamente tem que ser em ingles? 2. Dentre as funções do pacote existe uma que foi implementada no Matlab por outra pessoa e adaptada por mim no R. Preciso pedir autorização ao autor para submeter ou só referencio ele na documentação?. Agradeço a quem puder me ajudar. Waldelene Moura Analista de Sistemas Universidade Federal do Maranhão

1 0

mesclar dataframes
by Leonard de Assis 29 Out '14

29 Out '14

O melhor jeito de juntar dois dataframes ao bom e velho estilo SQL (SELECT ... INNER JOIN ...) é com merge mesmo né?

7 7

R - matlab
by Fernanda De Bastiani 29 Out '14

29 Out '14

Caros, alguém tem algum exemplo prático de como usar o pacote R.matlab para poder ler e rodar funções do matlab no R? Ou conhece algum outro pacote? Obrigada. Fernanda

2 1

Fazer rasterbricks a partir de rasterstacks que estão no ambiente do R e na pasta de trabalho.
by Maraisa Braga 28 Out '14

28 Out '14

Olá a todos! Sou nova por aqui, então se cometer alguma gafe, por favor me avisem. Preciso fazer rasterbricks de objetos que são rasterstacks e estão no meu ambiente do R com arquivos que tb são rasterstacks e estão no meu working diretory. Então, tenho isso: JJA<-list.files(pattern="JJA.grd") seasonal<-apropos(what="seasonal_") anual<-apropos(what="anual_") >JJA [1] "CCSM_0k_JJA.grd" "CCSM_26_JJA.grd" "CCSM_85_JJA.grd" > seasonal [1] "seasonal_CCSM_0k" "seasonal_CCSM_26" "seasonal_CCSM_85" > anual [1] "anual_CCSM_0k" "anual_CCSM_26" "anual_CCSM_85" Assim, os objetos "anual" e "seasonal" estão no ambiente do R e o objeto "JJA" está no working diretory. Tentei fazer isso, mas não dá certo: >brick(stack(JJA[1]), anual[1], seasonal[1]) Error in round(nl) : non-numeric argument to mathematical function > brick(stack(JJA[1]), stack(anual[1]), stack(seasonal[1])) Error in .local(.Object, ...) : `D:\...\seasonal_CCSM_0k' does not exist in the file system, and is not recognised as a supported dataset name. Error in .rasterObjectFromFile(x, objecttype = "RasterBrick", ...) : Cannot create a RasterLayer object from this file. (file does not exist) Alguma sugestão? Maraísa

4 6

GeoR <likfit>
by Felinto COSTA 28 Out '14

28 Out '14

Caros colegas. Ao se usar a função <likfit>, os valores por ela estimados da anisotropia geométrica (anisotropy angle e anisotropy ratio) já estão considerados no variograma ajustado ou devemos usá-los como base para transformar as coordenadas via <coords.aniso>, corrigir a anisotropia detectada e então estimar novamente os parâmetros do variograma teórico. FCosta

2 2

geoR: variog4, posição da legenda
by Helder Gramacho 28 Out '14

28 Out '14

Bom dia pessoal, Sempre que utilizo a função variog4 do pacote geoR as legendas estão ficando sobrepostas às curvas dos semivariogramas direcionais. Tentei posiciona-la adequadamente mas não consegui. É possível posicionar a legenda em outro local? [image: Imagem inline 1] Desde já agradeço, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 2

análise de agrupamento
by Aline Lipsky 28 Out '14

28 Out '14

Prezados colegas, Estou fazendo uma análise de agrupamento hierárquico pelo método de Ward com o seguinte script: safra2012 <- read.csv2("safra2012.csv") tiff("safra2012.tiff", width = 2000, height = 1100) par( mfrow = c (1,1), cex.main = 1.6, cex.lab = 1.2, cex.axis = 1.2, cex = 1.4) d <- vegdist(safra2012, method = "jaccard", binary = TRUE) dissimilaridade <- hclust(d, method="ward") plot(dissimilaridade, hang=-1, main = "Safra 2012") dev.off( ) Eu gostaria de plotar o dendograma sem o nome de cada objeto da minha matriz. Grata pela atenção. Atenciosamente, Aline Lipsky

5 5

[OFF-topic] Processar uma base de dados grande - ENEM 2012 - 6 milhões de registros - só 40 segundos
by walmes . 27 Out '14

27 Out '14

Saudações, Essa mensagem é apenas para compartilhar uma experiência que tive e considero que o problema aqui discutido não seja tão particular. Fui procurado por alguém que precisava processar um arquivo com resultados da prova do ENEM 2012 para aplicar teoria de resposta ao item. Precisava-se processar o arquivo original para chegar em um que tivesse apenas 5 das 80 colunas (notas para as competências sobre a redação) e somente às linhas que tivessem o conteúdo "P" (de presente) para a coluna "IN_STATUS_REDACAO". O arquivo é um csv, delimitador vírgula, aspas nas strings, 5791066 de linhas, 3.8 GB. Fiz o que me veio na cabeça de primeira: "dividir para conquistar". Lia porções de 30 mil linhas do arquivo (usando na read.table() o skip= e nrow=), separava as 5 colunas, mantinha só as linhas com registro P e escrevia para um arquivo com write.table(..., append=TRUE). Isso dentro de um for(). Na minha implementação, um tanto ingenua e corrida porque fiz enquanto era assistido a programar, quanto apliquei ao arquivo alvo levou 3h40. Esperava que fosse demorar mas não tanto. Minha máquina tem 16 GB de RAM. Sabendo que o bash do linux é muito eficiente para tarefas como operações em arquivos texto, eu fui buscar na internet como fazer o mesmo. Resultado é que tudo se resolveu com 37 segundos!!! É por isso que achei interessante compartilhar com a lista, inclusive para saber se alguem tem meios alternativos para solucionar o problema. Seguem os comandos que apliquei no terminal do Linux. Os dados fornecidos no exemplo são apenas as 10 mil primeiras linhas do arquivo original. ##----------------------------------------------------------------------------- ## 0. Aquisição dos dados. Apenas 10000 registros dos 5791066 do arquivo ## original. wget http://www.leg.ufpr.br/~walmes/data/DADOS_ENEM_2012_10millinhas.csv -O DADOS_ENEM_2012.csv ls wc -l DADOS_ENEM_2012.csv ##----------------------------------------------------------------------------- ## 1. Filtrar só às colunas de interesse. Da 74 à 79. Na 74 tem-se os ## valores P, F e outros. cut -d , -f 74-79 DADOS_ENEM_2012.csv > file.csv head -10 file.csv wc -l file.csv ##----------------------------------------------------------------------------- ## 2. Manter só as linhas com ocorrência do P. grep "P" file.csv > fileP.csv wc -l fileP.csv head -10 fileP.csv ##----------------------------------------------------------------------------- ## 3. Remover à colunas com P, ou seja, manter da 2 à 6. cut -d , -f 2-6 fileP.csv > file.csv head -10 file.csv wc -l file.csv ##----------------------------------------------------------------------------- ## 4. Remover as aspas para salvar espaço em disco. sed 's/\"//g' file.csv > fileP.csv head -10 fileP.csv ##----------------------------------------------------------------------------- ## 5. Como contar o número de linhas repetidas. Gerar tabela de frequência dos ## padrões de resposta. Economiza espaço e é a informação mínima ## necessária para ajustar TRI. Remover o cabeçalho. sed 1d fileP.csv > file.csv sort file.csv | uniq --count > fileP.csv head -10 fileP.csv tail -10 fileP.csv ##----------------------------------------------------------------------------- ## 6. Eliminar espaços no ínicio e colocar uma vírgula entre a ## frequência e o primeiro registro para que todos os campos sejam ## separados por `,`. sed -e 's/^ *//;s/ /,/' fileP.csv > fileFreq.csv head -10 fileFreq.csv wc -l fileFreq.csv ##----------------------------------------------------------------------------- À 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 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

8 9

Tabela 2 x 2
by Taynãna César Simões 27 Out '14

27 Out '14

Bom dia, pessoal. Gostaria de identificar os labels em uma tabela 2x2, através de uma condição. Em uma tabela comum, consigo usar: names(tabela)[condicao] Por exemplo: names(tabela)[sum(tabela)>1] Não consigo usar a mesma ideia com uma tabela 2x2: names(table(variável1,variável2))[condicao] Já tentei colnames e rownames, sem sucesso. Alguma ideia? Grata, Taynãna FIOCRUZ-Minas

3 5