Agosto 2015 - R-br

tiff() pixels, dpi e megabytes
by Pedro Emmanuel Alvarenga Americano do Brasil 07 Ago '15

07 Ago '15

Amigos de R, Peço paciencia nesse email porque eu não sei reproduzir o meu problema. Estou tentando submeter um paper de pesquisa a um periodico. No sitio deles a orientação é que as figuras devem ser tiff, ter no mínimo 300 dpi, e no maximo 40 megapixels. O tamanho de todos os arquivos não deve exceder 100 MB. O tamanho do arquivo não é problema porque dá pra comprimir. Usualmente eu uso algo do tipo (480 de altura e largura é o padrao e res 72 é o padrao.) tiff("meugrafico.tiff,width=480*10,heigh=480*10,res=72*10,compress="lzw") Entendendo que manter essas proporções de 480 x 480 com res = 72 é importante para manter o gráfico sem distorção. Só que um dos meus graficos é grande na horizontal e na vertical o que significa que pra não distorcer o grafico, mater tudo dentro das margens e mantendo a mesma resulução e devo fazer algo do tipo tiff("meugrafico.tiff,width=480*10,heigh=480*17.5,res=72,compression="lzw") Então o grafico já possui um montão de megapixels e poucos dpi. Quando eu tento aumentar a resolução para 300 dpi com algo como tiff("meugrafico.tiff,width=480*10*4.2,heigh=480*17.5*4.2,res=72*4.2,compression="lzw") o grafico excede a quantidade de megapixels que o arquivo pode ter. Eu descobri que o paint pode reduzir os megapixels e manter o dpi, mas essa figura tem tantos megapixels que o paint dá erro na hora que abre. Então até agora, eu estou tendo que submeter uma figura com dpi aquem do desejado e megapixels dentro do esperado. A pergunta é: há alguma forma de aumentar a resolução, sem ditorcer a imagem e sem aumentar a quantidade de megapixles, preferencialmente com o R? Abraço forte a todos, Pedro Brasil

3 5

Erro
by Pedro Mestre 06 Ago '15

06 Ago '15

Prezad@s, boa tarde. Estou tentando atualizar o R pelo "installr" seguindo os passos disponíveis no link ( http://www.r-bloggers.com/a-step-by-step-screenshots-tutorial-for-upgrading…) porém, está dando o seguinte erro: Erro em file(con, "r") : não é possível abrir a conexão As linhas de comando que executei (inclusive com as explicações) foram essas: # installing/loading the latest installr package: install.packages("installr"); library(installr) # install+load installr updateR() # updating R. Alguém sabe o que pode ser? PS. Desculpem se a dúvida é básica, mas não domino linguagem de programação. Abraços e obrigado, -- Pedro Mestre Ferreira Alves Assistente Técnico de Pesquisa Instituto de Pesca Av. Bartolomeu de Gusmão, 192 Ponta da Praia - Santos - SP CEP: 11030-906 Fone/Fax: (13) 3261 2342 www.pesca.sp.gov.br/pmestre@pesca.sp.gov.br Sent with MailTrack <https://mailtrack.io/install?source=signature&lang=en&referral=pedromestref…>

3 5

[OffTopic] Termos de mesma hierarquia em DoE (Delineamento de Experimentos)
by Cleber N.Borges 04 Ago '15

04 Ago '15

Prezados, Gostaria de uma indicação de bibliografia que dê detalhes (teóricos) sobre o porquê quando uma variável participa de uma *Interação* que se mostra significativa, mantém se no modelo o *Efeito Principal* mesmo que este não for. Agradeço pelas informação desde já. Cleber --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

2 1

geoR parana
by Hélio Gallo Rocha 04 Ago '15

04 Ago '15

Boa tarde listeiros, Estudando os dados de precipitação parana, nativos no geoR, observa-se tendencia no semivariogramas . Nos mapas exploratórios,observa-se tendencia na direção sudoeste, o que, a meu modo de ver, é natural, pois os dados climáticos apresentam variação mais ampla do que apenas área de um estado. Utilizando da Verossimilhança para estimar os parâmetros ( efeito pepita, patamar e alcance) com os dados sem retirada de tendencia e com três modos de retirada de tendencia, observou-se grande variação entre os parâmetros estimadas por vero entre as três retiradas de tendencia. Os parâmetros estimados para os dados sem retirada de tendencia foram mais discrepantes ainda, quando comparados com os com retirada de tendencia (acréscimos em torno de 300%). Interessantes foram os quatro mapas estimados por krigagem, que apresentaram pouca variação. Algum listeiro já estudou estes dados? Alguma explicação para parâmetros tão diferentes, produzirem mapas com pouca variação? -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Qualidade de vida
by Edson Lira 04 Ago '15

04 Ago '15

Estou tentando analisar um banco que mede qualidade de vida. Abaixo um script usado no SPSS (http://www.ufrgs.br/psiquiatria/psiq/whoqol86.html)Instrumentos de Avaliação de Qualidade de vida - OMS - versão em português | | | | | | | | | | | Instrumentos de Avaliação de Qualidade de vida - OMS - v...Versão em português dos Instrumentos de Avaliação de Qualidade de Vida - WHOQOL-100 e WHOQOL-bref - da Organização Mundial de Saúde. Projeto de... | | | | Visualizar em www.ufrgs.br | Visualizado por Yahoo | | | | | STEPS FOR CHECKING AND CLEANING DATA AND COMPUTING FACET AND DOMAIN SCORES Recode q3 q4 q26 (1=5) (2=4) (3=3) (4=2) (5=1). (scores based on a 4-20 scale) compute dom 1= (mean.6 (q3, q4, q10, q15, q16, q17, q18)) * 4. compute dom 2= (mean.5 (q5, q6, q7, q11, q19, q26)) * 4 compute dom 3= (mean.2 (q20,q21,q22)) * 4 compute dom 4= (mean.6 (q8, q9, q12, q13, q14, q23, q24, q25)) * 4 compute overall= (mean.2 (q1, q2)) * 4 (scores transformed to a 0-100 scale) compute dom1b= (dom1 - 4) * (100/16). compute dom2b= (dom2 - 4) * (100/16). compute dom3b= (dom3 - 4) * (100/16). compute dom4b= (dom4 - 4) * (100/16). compute q1b= (q1 - 1) * (100/16). compute q2b= (q2 - 1) * (100/16). Gostaria de replicar este script, alguém poderia me ajudar? O link do banco: Download qv.csv | | | | | | | | | | | Download qv.csvDownload qv.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehos... | Visualizado por Yahoo | | | | | [ ]'s Edson LiraEstatístico Manaus-Amazonas

2 1

Exportar resultado de ctree para txt
by Paulo Nogueira Starzynski 03 Ago '15

03 Ago '15

Olá colegas de lista. Estou criando algumas árvores de decisão usando party::ctree e gostaria de gravar o resultado em um arquivo texto, para documentar o resultado rapidamente. Um simples txt, sem formatação, já bastaria. Tentei: library(party) airq <- subset(airquality, !is.na(Ozone)) airct <- ctree(Ozone ~ ., data = airq, controls = ctree_control(maxsurrogate = 3)) airct plot(airct) write(airct, file = "/home/paulo/Documentos/ctree.txt") Mas obtenho: "Error in cat(list(...), file, sep, fill, labels, append) : argumento 1 (tipo 'S4') não suportado por 'cat'" Tentei com writeLines, write.table, etc... inclusive adicionando print() ao objeto airct. Não funcionou. Alguém pode orientar? Estou errando algo simples creio. Grato, Paulo Nogueira Starzynski

3 3

Reorganizar base de dados
by Bernardo Rangel Tura 02 Ago '15

02 Ago '15

Bom dia a todos Estou precisando realizar uma operação de ordenação de um grupo de variáveis para o modelo que a equipe resolveu utilizar preciso transforma uma base de dados deste formato: idquestionario cod_cartao ord_cartao 1 1 33333 5 2 1 Morto 4 3 1 12123 2 4 1 22332 3 5 1 32332 5 6 1 11111 1 7 1 33132 5 8 1 33231 5 9 1 21333 3 10 10 23222 7 11 10 Morto 9 12 10 12233 4 13 10 33333 8 14 10 12111 2 15 10 22332 5 16 10 11333 3 17 10 32213 6 18 10 11111 1 Para este formato idquestionario 11111 33333 MORTO 12123 22332 32332 33132 33231 1 1 1 5 4 2 3 5 5 5 2 10 1 8 9 NA 5 NA NA NA 21333 23222 12233 12111 11333 32213 3 NA NA NA NA NA NA 7 4 2 3 6 Para deixar mais claro idquestionario = identificador, ord_cartao = posição que o cartão foi ordenado, cod_cartao = código do cartão (total 244 códigos) Cada pessoa responde 9 dos 244 códigos (3 fixos e 6 variáveis), cada linha inicia com o identificador cada código de cartão (cod_cartao) passa a ser uma nova variável e o seu valor naquela linha é o ranqueamento (ord_cartao) ou NA se o paciente não foi sorteado para este cartão. alguém pode me ajudar? Obrigado por antecipação []s Tura

2 2

Re: [R-br] Digest R-br, volume 56, assunto 3
by Diogo Braga Mendes 01 Ago '15

01 Ago '15

use a função *reshape* conforme código a seguir: reshape(dados, v.names = "ord_cartao", idvar = "idquestionario", timevar = "cod_cartao", direction = "*wide*") em que dados (codificando mortes por 99999) é: dados<-structure(list(idquestionario = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 10, 10, 10, 10, 10, 10, 10, 10, 10), cod_cartao = c(33333, 99999, 12123, 22332, 32332, 11111, 33132, 33231, 21333, 23222, 99999, 12233, 33333, 12111, 22332, 11333, 32213, 11111), ord_cartao = c(5, 4, 2, 3, 5, 1, 5, 5, 3, 7, 9, 4, 8, 2, 5, 3, 6, 1)), .Names = c("idquestionario", "cod_cartao", "ord_cartao"), row.names = c(NA, -18L), class = "data.frame") Diogo Braga Mendes Estatístico do IPAMB-Belém/PA. Mestre em Estatística-PPGME/UFPA. Bacharel em Estatística, CONRE 9558.

1 0