Maio 2015 - R-br - DINF UFPR

erro na krigagem
by Hélio Gallo Rocha 16 Mai '15

16 Mai '15

Caros listeiros. Rodando krigagem de clima, com retirada de tendência polinomial de segunda ordem: KC.2nd <- krige.conv(gd,loc=gr,krige=krige.control(type="OK",obj=LIKFIT[[3]] , trend.d="2nd", trend.l=~gr0)) deu esse erro: Erro em solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : sistema é computacionalmente singular: condição recíproca número = 4.2826e-23 Além disso: Mensagens de aviso perdidas: 'overlayPointsWithPolygons' is deprecated. Use 'over' instead. See help("Deprecated") and help("sp-deprecated"). Agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

4 6

Dúvida sobre escrito em xlab
by Leonardo Monteiro 15 Mai '15

15 Mai '15

Senhores, boa noite. Como modifico o formato da legenda: t/ha para t ha-1? Sendo que o -1 deve ficar sobreescrito. obrigado, Leonardo

3 2

eficiência de mem?==?UTF-8?Q?ória no R
by Fatima do Nascimento Silva 15 Mai '15

15 Mai '15

prezados, boa noite. Venho por meio deste solicitar uma dica de como fazer a leitura(importação) de uma base de dados em formato txt, com caracteres separados por ; e com 29 colunas e 1.520.171 linhas no R, de uma maneira que não perca eficiência de memória. estou fazendo a leitura usando o comando: (dados <- read.table("SAJ_WIN1252.txt",header=T,dec=",",sep=";")) Como o arquivo tem 3Gb, o R demora uns 10 min para abri-lo, mas não consigo fazer mais nada porque o programa trava. Alguém teria dicas de como me ajudar? agradeço a atenção. Atenciosamente, Fátima Nascimento Estatística - UFRN Mestre em Ciência e Engª de Petróleo-PPGCEP/UFRN Doutoranda no Programa de Pós-Graduação em Ciência e Engª de Petróleo- PPGCEP/UFRN

3 6

Como fazer para dois graficos do ggplot ficarem lado a lado
by Pedro Emmanuel Alvarenga Americano do Brasil 12 Mai '15

12 Mai '15

Amigos de R, Estou mais uma vez batendo cabeça com o ggplot. Eu gostaria que doi graficos ficasse lado a lado na mesma janela, mas infelizmente ainda estou pensando com o plot.default() e não estou conseguindo. Alguma luz? d <- data.frame(ID=rep(1:10,100),outcome=c(rep(0,500),rep(1,500)),Sex=c(sample(c("F","M"),500,T,c(.8,.2)),sample(c("M","F"),500,T,c(.8,.2))),score=c(rnorm(500,5,10),rnorm(500,20,10)),grave=c("Sim","Não")) library("lme4") fit <- glmer(outcome ~ Sex + score + (1 | ID), data = d, family = binomial,subset=grave=="Sim") summary(fit) fit2 <- glmer(outcome ~ Sex + score + (1 | ID), data = d, family = binomial,subset=grave=="Não") summary(fit2) library("sjPlot") sjp.glmer(fit,type='fe') sjp.glmer(fit2,type='fe') Que a força esteja sempre com voces! Pedro Brasil

2 2

Primeiro contato
by Kamioka 12 Mai '15

12 Mai '15

Olá pessoal, boa noite. Já estudo data science/analysis há alguns meses usando R e python. Fiquei feliz em saber que temos uma lista brasileira. Eduardo

2 1

[OFF-TOPIC] Avaliação geral de modelos em SEM
by Vinícius Lionel Mateus 12 Mai '15

12 Mai '15

Olá a todos, Alguém tem uma boa sugestão de bibliografia para avaliação de "model fit" em SEM? Tenho uma amostra relativamente grande (454 x 25) e que não segue a multinormal. Minhas referências básicas têm sido Kline (2014) e Ralph (1996). Gostaria de investigar a magnitude e fontes de /misfit/ entre minhas matrizes de covariância, oriundas da amostra e impostas pelo modelo. Desde já, muito obrigado. Abs, Vinícius -- Best regards, Vinícius Lionel Mateus, MSc (http://lattes.cnpq.br/6501001637020665) Chemistry PhD Student Atmospheric Chemistry Laboratory - Dep. Chemistry PUC-Rio - Pontifical Catholic University of Rio de Janeiro Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brazil CEP: 22451-900 Telefone: (+45) 26 28 28 51 (+55) (21) 3527-1327 (+55) (21) 993 - 588 - 051 Skype: vinicius.lionel http://www.qui.puc-rio.br/index.html

1 0

Knitr iniciante
by Carlos Andrade 12 Mai '15

12 Mai '15

Prezados amigos da lista. Baixei e instalei o pacote knitr. Estou usando o RStudio versão 99.9.9 no Slackware 14.1, o R é da versão 3.1.2. Em uma pasta que fiz alguns exercicios e salvei um gráfico no formato png iniciei o RStudio e abri um arquivo testeKnitr.Rmd. Executei knitHTML e funcionou. Substitui o plot(cars) do arquivo demo que funcionou por plot(luglikp1) com um arquivo que criei. Executei knitHTML e recebi a mensagem de erro: ## Error in plot(luglikp1): objeto 'luglikp1' não encontrado Também deu erro no Windows 7, R version 3.2.0 Que ajustes estão faltando para funcionar corretamente? -- Atenciosamente, Prof. Carlos A. S. de Andrade LAPEA - Laboratório de Pesquisa em Economia Aplicada e Engenharia de Produção Universidade Federal de Campina Grande. Centro de Humanidades Unidade Acadêmica de Economia

3 2

javareconf + windows + R 3.2
by Leonard de Assis 11 Mai '15

11 Mai '15

Olá, O manual do R (http://cran.r-project.org/doc/manuals/r-release/R-admin.html#Java-support) diz o seguinte: <PRINT> Note that this unfortunately depends on the exact version of the JRE/JDK installed, and so may need updating if the Java installation is updated. This can be done by running R CMD javareconf which updates settings in both etc/Makeconf and R_HOME/etc/ldpaths. See R CMD javareconf --help for details. </PRINT> Contudo, ao executar tal comando no Windows, recebo a seguinte msg: <PRINT> C:\Program Files\R\R-3.2.0\bin>R CMD javareconf 'javareconf' não é reconhecido como um comando interno ou externo, um programa operável ou um arquivo em lotes. </PRINT> Ao executar o Help do comando, deparo com o seguinte: <PRINT> C:\Program Files\R\R-3.2.0\bin>R --arch x64 --help Usage: Rterm [options] [< infile] [> outfile] [EnvVars] Start R, a system for statistical computation and graphics, with the specified options EnvVars: Environmental variables can be set by NAME=value strings Options: -h, --help Print usage message and exit --version Print version info and exit --encoding=enc Specify encoding to be used for stdin --encoding enc ditto --save Do save workspace at the end of the session --no-save Don't save it --no-environ Don't read the site and user environment files --no-site-file Don't read the site-wide Rprofile --no-init-file Don't read the .Rprofile or ~/.Rprofile files --restore Do restore previously saved objects at startup --no-restore-data Don't restore previously saved objects --no-restore-history Don't restore the R history file --no-restore Don't restore anything --vanilla Combine --no-save, --no-restore, --no-site-file, --no-init-file and --no-environ --max-mem-size=N Set limit for memory to be used by R --max-ppsize=N Set max size of protect stack to N -q, --quiet Don't print startup message --silent Same as --quiet --slave Make R run as quietly as possible --verbose Print more information about progress --internet2 Use Internet Explorer settings for proxies etc. --args Skip the rest of the command line --ess Don't use getline for command-line editing and assert interactive use -f file Take input from 'file' --file=file ditto -e expression Use 'expression' as input One or more -e options can be used, but not together with -f or --file An argument ending in .RData (in any case) is taken as the path to the workspace to be restored (and implies --restore) Or: R CMD command args where 'command' is one of: INSTALL Install add-on packages REMOVE Remove add-on packages SHLIB Make a DLL for use with dynload BATCH Run R in batch mode build Build add-on packages check Check add-on packages Rprof Post process R profiling files Rdconv Convert Rd format to various other formats Rdiff difference R output files Rd2pdf Convert Rd format to PDF Rd2txt Convert Rd format to pretty text Stangle Extract S/R code from vignette Sweave Process vignette documentation config Obtain configuration information about R open Open a file via Windows file associations texify Process a latex file Use R CMD command --help for usage information for each command. </PRINT> A questão é a seguinte: Foi abandonado o suporte R ao java na versão 3.2, visto que tal comando agora não existe? Leonard

3 6

Java 8 e bibliotecas
by Andre Oliveira 11 Mai '15

11 Mai '15

Pessoal,configurei tudo aqui a xlsx, gdata, XLConnect entre outras que dependem do java, deu tudo certo, mas tive que colocar java 8 por causa do site da caixa e ai voltou tudo a zero e não consigo instalar de novo. Alguma dica? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

2 2

Re: [R-br] Estimativas de parâmetros de uma Weibull utilizando duas metodologias: Newton-Raphson "vs" Nelder Mead (estimativas muito diferentes) o que pode estar acontecendo??
by Romero Luiz M. Sales Filho 11 Mai '15

11 Mai '15

Olá Wagner, segui o seu conselho e coloquei para o caso Nelder-Mead as estimativas iniciais como sendo as estimativas encontradas no caso de Newton-Raphson, o valor foi realmente quase em cima! Entaõ o problema realmente era esse, a seleção dos pontos iniciais. Fica então uma dúvida: Quando tenho dist. de probabilidades que não conseguimos estimar via Newton Raphson, por exemplo, como podemos ter uma idéias dos melhores pontos iniciais a serem utilizados?? Tem como encontrar alguma forma de se encontrar algo interessante ou é somente por tentativa e erro? Elias, muito obrigado também por compartilhar informação! Acho que vc deu uma dica de como encontrar os melhores pontos iniciais, né? Infelizmente não entendi... vc poderia me explicar melhor... Ah, outra coisa, com relação ao ks.test tbm tinha dúvidas se ele realmente era adequado, uma vez que não tenho os parâmetros e sim estimativas desses, aí gero uma amostra utilizando essas estimativas e isso talvez não fosse o mais indicado pq se o algoritomo tivesse convergido para valores errados de parâmetros e mesmo assim a distribuição osse Weibull o ajuste seria excelente quando na verdade os parâmetros estariam errados, né isso (por favor, me corrija se eu estiver errado)??... Então vc colocou uma abordagem utilizando bootstrap, essa seria uma forma alternativa de utilização do ks.test sendo mais apropriada para esse tipo de problema em que não conheço os parâmetros?? Muito obrigado à todos!! Abraço, Romero. -- *Profº Romero Luiz Mendonça Sales Filho* *Estatístico / Mestre em Engª de ProduçãoUFRPE - UAGCel.: (81) 9667-4477* *romero.sfilho(a)gmail.com* <romero.sfilho(a)gmail.com> *romero_sfilho(a)yahoo.com.br* <romero_sfilho(a)yahoo.com.br>

3 3