- R-br - DINF UFPR

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

Ajuste de regressão linear
by Andre Oliveira 02 Jul '15

02 Jul '15

Bom dia,este é um CMR hipotético. dados <- expand.grid(A=factor(paste("A", 1:3, sep="")),B=factor(paste("B", 1:4, sep="")), rep=1:4)dados$y <- rgamma(n=48,shape=2, scale = 2) summary(lm(y~rep, data = subset(dados,B=="B1")))plot(y ~resp) summary(lm(y~rep, data = subset(dados,B=="B2")))plot(y ~resp) summary(lm(y~rep, data = subset(dados,B=="B3")))plot(y ~resp) summary(lm(y~rep, data = subset(dados,B=="B4"))) Como posso automatizar estas regressões com um for? Alguma sugestão? Sem plot já está ótimo! Obrigado André Oliveira Souza

4 4

Gráfico de calor para uma matriz bidimensional
by Michelle Bau Graczyk 01 Jul '15

01 Jul '15

Boa noite, Eu preciso fazer um gráfico de calor de matrizes bidimensionais de correlacão entre volumes de 30 empresas. Eu tentei usar o heatmap.2 mas ele reorganiza as linhas e eu gostaria que ela ficasse na ordem que já estava. Vocês teria algum modelo no qual eu possa me espelhar ou alguma dica? Muuito obrigada, Michelle

2 1

Agregando ou reshape de dados?
by Pedro Emmanuel Alvarenga Americano do Brasil 30 Jun '15

30 Jun '15

Amigos de R, Há algumas horas estou empacado e ja tentei diversas coisas aqui, como o reshape e aggregate e não consigo. Bom a intensão é transformar os dados de tal forma que duas variáveis repetitivas sejam transformadas, de tal forma que as categorias dessas variáveis passem a ser as novas colunas e o seu conteúdo passe a ser True ou False para cada identificador, mais ou menos assim; Reparem que a linha 2 e 6 possuem o mesmo identificador Original: Unique_Identifier AdmissionReasonName_sec 1 {6841FBAD-EB26-4EB3-AD8C-3AE00E23A563} Trauma não-cirurgico 2 {D4BE0DF7-A143-4EE5-9B9F-659F198A817C} Neurológico 3 {A6240B7F-3571-473E-AEFF-944372B312A4} Infecção/Sepse 4 {946C597C-A79C-47A0-B0A6-16BEB9817571} Cardiovascular 5 {023FF83F-A65C-49FE-8E0C-D26D8BD82B40} Cardiovascular 6 {D4BE0DF7-A143-4EE5-9B9F-659F198A817C} Cardiovascular Formato desejado: Reparem que aqui o ID duplicado foi removido mas aparecem True nas variaveis Cardiovascular e Neurologico Unique_Identifier Trauma não-cirurgico Neurologico Infecção/Sepse Cardiovascular 1 {6841FBAD-EB26-4EB3-AD8C-3AE00E23A563} True False False False 2 {D4BE0DF7-A143-4EE5-9B9F-659F198A817C} False True False True 3 {A6240B7F-3571-473E-AEFF-944372B312A4} False False True False 4 {946C597C-A79C-47A0-B0A6-16BEB9817571} False False False True 5 {023FF83F-A65C-49FE-8E0C-D26D8BD82B40} False False False True Para montar o banco original s <- data.frame(Unique_Identifier=c("{6841FBAD-EB26-4EB3-AD8C-3AE00E23A563}","{D4BE0DF7-A143-4EE5-9B9F-659F198A817C}","{A6240B7F-3571-473E-AEFF-944372B312A4}","{946C597C-A79C-47A0-B0A6-16BEB9817571}","{023FF83F-A65C-49FE-8E0C-D26D8BD82B40}","{D4BE0DF7-A143-4EE5-9B9F-659F198A817C}"), AdmissionReasonName_sec=c("Trauma não-cirurgico","Neurológico","Infecção/Sepse","Cardiovascular","Cardiovascular","Cardiovascular")) s Abraço, Pedro Brasil

4 4

Re: [R-br] Método de Jacobi
by Tiago Souza Marçal 29 Jun '15

29 Jun '15

Da uma olhada neste site e verifique se os códigos te atendem. http://morotalab.org/Mrode2005/mme/mme.html#section0002 Att. Tiago. Wecsley Prates <woprates(a)gmail.com> escreveu: Olá a todos... Gostaria de saber se alguém tem o algoritmo do Método de Jacobi em linguagem do r para encontrar solução de equações lineares. Encontrei uma programação, mas está dando erro. Essa programação me fornece um erro do tipo Erro em a[i, j] <- a[i, j] - m * a[k, j] : substituto tem comprimento zero. Como resolvo??? Agradeço a atenção a = matriz dos coeficientes b <- vetor das constantes x <- chute inicial do vetor das variáveis for (k in 1:length(x)-1){ for (i in k+1:length(x)){ m <- a[i,k]/a[k,k] a[i,k] <- 0 for(j in k+1:length(x)){ a[i,j] <- a[i,j]-m*a[k,j] b[i] <- b[i]-m*b[k] } } } x[length(x)] <- b[length(x)]/a[length(x),length(x)] for (k in length(x)-1:1){ s <- 0 for (j in k+1:length(x)){ s <- s+a[k,j]*x[j] } x[k] <- (b[k]-s)/a[k,k] } -- * Wecsley O. Prates* *Doutorando em Estatística - Universidade Federal de Minas Gerais* _______________________________________________ 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.

3 3

Método de Jacobi
by Wecsley Prates 28 Jun '15

28 Jun '15

Olá a todos... Gostaria de saber se alguém tem o algoritmo do Método de Jacobi em linguagem do r para encontrar solução de equações lineares. Encontrei uma programação, mas está dando erro. Essa programação me fornece um erro do tipo Erro em a[i, j] <- a[i, j] - m * a[k, j] : substituto tem comprimento zero. Como resolvo??? Agradeço a atenção a = matriz dos coeficientes b <- vetor das constantes x <- chute inicial do vetor das variáveis for (k in 1:length(x)-1){ for (i in k+1:length(x)){ m <- a[i,k]/a[k,k] a[i,k] <- 0 for(j in k+1:length(x)){ a[i,j] <- a[i,j]-m*a[k,j] b[i] <- b[i]-m*b[k] } } } x[length(x)] <- b[length(x)]/a[length(x),length(x)] for (k in length(x)-1:1){ s <- 0 for (j in k+1:length(x)){ s <- s+a[k,j]*x[j] } x[k] <- (b[k]-s)/a[k,k] } -- * Wecsley O. Prates* *Doutorando em Estatística - Universidade Federal de Minas Gerais*

1 0

Sistema é computacionalmente singular - Matriz/Mahalanobis
by Andre Oliveira 26 Jun '15

26 Jun '15

Prezados boa noite,estou com problemas em uma matriz inversa! Alguma sugestão? Precisava muito desta matriz D2! require(biotools)setwd("/home/andre")dados <- read.csv("Multivariada.txt", sep="")Cov=cov(dados)D2=D2.dist(dados, Cov) # Não hpa inversa!Error in solve.default(cov, (x[i, ] - x[j, ])) : sistema é computacionalmente singular: condição recíproca número = 1.53112e-19 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

4 4

Dúvida glm (family)
by Karen Castillioni 25 Jun '15

25 Jun '15

Pessoal, estou tentando analisar dados de cobertura vegetal (5%, 10%, 15%, 20%..até 100%) e dados de biomassa (em gramas por m2) usando glm, pois nenhum dos dados apresenta distribuição normal, nem mesmo com transformação. Porém, estou em dúvida sobre qual família usar para cada variável. Alguém pode ajudar? Agradeço a ajuda.

3 3