Julho 2015 - R-br - DINF UFPR

Problemas com polygrid
by Arlindo Fabrício Corrêia 06 Jul '15

06 Jul '15

Boa tarde a todos, Pergunto se há alguma alteração na ferramenta "*polygrid* {geoR}" ou seria um problema que ocorre apenas em minha máquina. Já venho utilizando esta função há um tempo e atualmente estou com problemas para limitar um grid por meio de bordas utilizando o pacote. Pergunto também se é um problema do pacote ou a rotina deve ser alterada. Tentei gerar uma nova "borda" por meio dos comandos *SpatialPoints, SpatialPolygon* como apresentado no erro, porém, sem sucesso. Há uma outra solução aparente? Já fiz muitas tentativas e procurei exaustivamente e ainda não cheguei a uma solução. A seguir encaminho as linhas com o próprio exemplo da ferramenta. > #Examples> #> poly <- matrix(c(.2, .8, .7, .1, .2, .1, .2, .7, .7, .1), ncol=2)> plot(0:1, 0:1, type="n")> lines(poly)> poly.in <- polygrid(seq(0,1,l=11), seq(0,1,l=11), poly, vec=TRUE)Erro em (function (classes, fdef, mtable) : unable to find an inherited method for function ‘overlay’ for signature ‘"SpatialPoints"’> points(poly.in$xy)Erro em points(poly.in$xy) : objeto 'poly.in' não encontrado > Agradeço pela atenção. -- Arlindo Fabrício Corrêia Engenheiro Agrônomo Doutorando em Agronomia - UNIOESTE - Mal. Cdo. Rondon/PR

3 2

[off] z-score
by Rodolfo Timoteo da Silva 06 Jul '15

06 Jul '15

Pessoas, permitam-me fazer um flashback aproveitando o "expertise" em estatística da lista. Estou fazendo uma revisão em estatística, e só para confirmar em relação ao exercício da figura em anexo. x = 28 média = 24 desvio-padrão = 1.6 z = (28 - 24) / 1.6 z = 2.5 Assim, a proporção para ser maior que 28mpg seria, pela tabela da distribuição normal, (1 - 0.9938 = 0.0062 ou 0.62%). Correções são muito bem vindas. Antecipadamente P.S.: é possível calcular usando somente [1 - pnorm(28,24,1.6)] ou tem outro comando? -- --- Zhu

3 5

Função LOGIT
by Ricardo Gallo 05 Jul '15

05 Jul '15

Boa noite galera, Estou precisando de uma ajuda, estou com umas funções "logit" para distribuição binomial (0 e 1) usando GLM, e gostaria de aplicar a máxima verossimilhança no lugar. Segue o comando usado: mod1<-glm(XF ~ X1, data = dados, family = binomial(link="logit"))summary(mod1)exp(coef(mod1)) Acho que é algo simples, mas não estou conseguindo. Vlw - Ricardo Gallo - Engenheiro Florestal Doutorando em Ciência Florestal - UFV(31) 8213 5170 (claro)

2 2

Modelo Periódico Gama Autorregressivo
by Victor Eduardo 04 Jul '15

04 Jul '15

Pessoal, alguém aqui chegou a trabalhar com modelos gama autorregressivo no R? Tenho algumas dúvidas na minha programação. Alguém conhece algum fórum ou material que tenha discussão sobre esse assunto no R? Atenciosamente, Victor Eduardo

1 0

Transpor linhas em colunas
by Alexandre Loures 04 Jul '15

04 Jul '15

Bom dia! Tenho uma base em que há 1.241 ID's (FORMATO CHARACTER) sendo que cada ID repete 1.980 vezes. Por exemplo, ANO ID PREÇO 1989 0101 8 1990 0101 7 1991 0101 15 1989 0102 9 1990 0102 8 1991 0102 2 1989 0103 3 1990 0103 10 1991 0103 11 1989 0104 6 1990 0104 1 1991 0104 11 . . . Porém gostaria de transformar cada uma dessas ID's em colunas e que assumissem os valores da variável PREÇO. ANO ID_0101 ID_0102 ID_0103 ID_0104 ... 1989 8 9 3 6 1990 7 8 10 1 1991 15 2 11 11 Alguém poderia me ajudar? Desde já muito obrigado! -- *Alexandre Rodrigues Loures*

2 2

Ajuda com Kaplan-Meier
by "Ary Serpa Neto" 04 Jul '15

04 Jul '15

2 1

Beamer e RMarkdown
by David Feitosa 04 Jul '15

04 Jul '15

Caros, boa tarde. Estou criando um documento e gostaria de usar RMarkdown no lugar do LaTex convencional. Tentei o abaixo seguindo as orientações do site e de alguns blogs. Deparei-me com dois "problemas": O 1º que consistia em gerar uma lista de segundo nível, consegui contornar adicionando mais dois espaços na indentação, pois a documentação dizia para usar apenas 2, mas só resolvi com 4. O 2º consiste em gerar a lista com os tópicos de cada slide. Quando escolho gerar apenas o documento como PDF, é gerado o sumário, mas quando mudo para o BEAMER, é gerado um slide em branco. Alguém teria ideia de como resolver? Segue o "código". Eu comento/descomento a linha do "pdf_document" e "beamer_presentation" para mudar o formato de documento final. Gostaria de usar a do beamer: --- title: "Título" author: "David Feitosa" output: # pdf_document: beamer_presentation: toc: true header-includes: - \usepackage[brazil]{babel} colortheme: "dolphin" --- # R Markdown 1 - A - B - C # R Markdown 2 - A - B - C # R Markdown 3 - A - B - C # R Markdown 4 - A - B - C # R Markdown 5 - A - B - C # R Markdown 6 - A - B - C Meu ambiente: R version 3.2.1 (2015-06-18) -- "World-Famous Astronaut" Platform: x86_64-w64-mingw32/x64 (64-bit) RStudio Version 0.99.442 knitr 1.10.5 pdfTeX, Version 3.14159265-2.6-1.40.15 (MiKTeX 2.9 64-bit) pandoc 1.14.0.1 Compiled with texmath 0.8.2, highlighting-kate 0.6. Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

1 0

Pacote que trabalhe a distribuição gama com parâmetro location
by Victor Eduardo 04 Jul '15

04 Jul '15

Pessoal, vocês indicam algum pacote que trabalhe com a distribuição gama que tenha o parâmetro location? Atenciosamente, Victor Eduardo

5 5

Identificação de pontos críticos associados a uma legenda em ggplot
by Iábita Fabiana 04 Jul '15

04 Jul '15

Boa Tarde! Estou fazendo alguns gráficos utilizando o ggplot e gostaria de identificar alguns pontos no gráfico, bem como adicionar uma legenda desses pontos, alguém poderia me ajudar sobre como poderia associar os pontos a uma legenda? E se há uma maneira de chamar esses pontos no gráfico de um modo mais prático? Lines <- " x y 0 34 1 40 2 58 3 78 4 108 5 126 6 130 7 147 8 157 9 162 10 168 11 175 12 165 13 188 14 232 15 250 16 278 17 284 18 280 19 310 20 312 21 320 22 352 23 380 24 386 25 426 26 422 27 440 28 455 29 468 30 480 31 486 32 485 33 490 34 492 35 500 36 490 37 505 38 504 39 505 40 500 41 503 42 492 43 486 44 490 45 496 46 500 47 507 48 500 49 492 " " dados <- read.table(textConnection(Lines), header=TRUE); closeAllConnections() dados require(ggplot) qplot(x,y,data=dados)+ geom_smooth( method = "nls", formula='y~b1/(1+exp(-b2-b3*x))' , data=dados,correlation=corAR1(form=~x), start=c(b1=514,b2=-3,b3=0.1716),se = FALSE, linetype = 1.5,lwd=1)+ annotate("point", x = 13, y = 220, shape=8, size=6,label ="some text") + annotate("text", x = 14, y = 210, shape=15, size=4,label = "P4") Obrigada. *Iábita * Doutoranda em Estatística e Experimentação agronômica na ESALQ/USP Mestre em Estatística e Experimentação Agropecuária na UFLA Licenciada em Matemática - UFPI

2 2

Re: [R-br] Função LOGIT
by Ricardo Gallo 04 Jul '15

04 Jul '15

Walmes, Então não sei se estou fazendo besteira, mas pelo Stata nos fornece modelos usando a função de máxima verossimilhança. Não entendo muito desta análise, vou usa-la para análise de sobrevivência e como o trabalho todo fiz o uso do R gostaria de fazer estas análises por ele também. Tem alguma sugestão? Posso usar o exemplo passado mesmo? Obrigado. Em 04/07/2015 11:07, Walmes Zeviani <walmeszeviani(a)gmail.com> escreveu: Ricardo, 1) seu exemplo não é reproduzível. 2) seu objetivo não foi claro. Você gostaria de fazer a estimação de um modelo glm binomial com função de ligação logit por máxima verossimilhança? À disposição. Walmes. _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0