Agosto 2014 - R-br

Matriz X de delineamento
by andrebvs＠bol.com.br 15 Ago '14

15 Ago '14

2 1

gráfico area
by luc.souza＠usp.br 14 Ago '14

14 Ago '14

Prezados, gostaria de fazer um gráfico no R, com os dados abaixo: d1 <- data.frame(X = 1:10, Y1 = c(-1.2,0,0,0,0,-30.3,-1.6,0,0,-2.4), Y2 = c(0,8.9,62.7,14,7.4,0,0,14.6,9.3,0)) O gráfico que gostaria de fazer é o mesmo que o gráfico de área do excel. Tentei fazer alguns gráficos, mas não ficou nem um pouco parecido. Primeiro tentei plotar duas linhas e depois preencher a área com uma cor. Não consegui. Alguém poderia me ajudar? Obrigado. Lucas Lucas F. de Souza Eng. Agrícola - CONAB M.Sc. Produção Vegetal - UENF Doutor em Agrometeorologia - ESALQ/USP tel: (19) 98160-7068

3 3

Integração numérica de valores tabulados
by Alexandre Castagna 14 Ago '14

14 Ago '14

Prezados, Preciso aplicar repetidas vezes uma integral em um distribuição de probabilidade e por isso preciso de um método muito rápido (são milhares de repetições). Porém não existe uma função matemática que descreva a distribuição, os valores estão tabulados. A função de integração que construí é muito lenta (0.2 segundos na minha máquina); queria reduzir essa velocidade em 10 ou 100 vezes. Apresento a função com os valores tabulados: angles <- log(c(0.001, 0.01, 0.1, 0.12589, 0.15849, 0.19953, 0.25119, 0.31623, 0.39811, 0.50119, 0.63096, 0.79433, 1, 1.2589, 1.5849, 1.9953, 2.5119, 3.1623, 3.9811, 5.0119, 6.3096, 7.9433, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 105, 110, 115, 120, 125, 130, 135, 140, 145, 150, 155, 160, 165, 170, 175, 180)*pi/180) mc.petold <- log(c(869105, 39180, 1.76661e3, 1.29564e3, 9.50172e2, 6.99092e2, 5.13687e2, 3.76373e2, 2.76318e2, 2.18839e2, 1.44369e2, 1.02241e2, 7.16082e1, 4.95803e1, 3.39511e1, 2.28129e1, 1.51622e1, 1.00154e1, 6.57957, 4.29530, 2.80690, 1.81927, 1.15257, 4.89344e-1, 2.44424e-1, 1.47151e-1, 8.60848e-2, 5.93075e-2, 4.20985e-2, 3.06722e-2, 2.27533e-2, 1.69904e-2, 1.31254e-2, 1.04625e-2, 8.48826e-3, 6.97601e-3, 5.84232e-3, 4.95306e-3, 4.29232e-3, 3.78161e-3, 3.40405e-3, 3.11591e-3, 2.91222e-3, 2.79696e-3, 2.68568e-3, 2.57142e-3, 2.47603e-3, 2.37667e-3, 2.32898e-3, 2.31308e-3, 2.36475e-3, 2.50584e-3, 2.66183e-3, 2.83472e-3, 3.03046e-3, 3.09206e-3, 3.15366e-3)) cum.beta.ph.p.petzold <- function(psi) { Integral <- integrate(function(x) { exp(approx(angles, mc.petold, xout = log(x))$y)*sin(x) }, lower = 1.745329e-05, upper = psi)$value*2*pi return(Integral) } cum.beta.ph.p.petzold(psi = pi) Onde psi deve ser entrado em radianos. A integral de "0" a pi é apenas aproximadamente 1, por conta de uma simplificação dos valores e pelo limite inferior de integração em 0.001º (1.745329e-05 rad), já que o valor é infinito em 0. A normalização para 1 é fácil de ser realizada. A função é aplicada milhares de vezes em um programa de simulação com números aleatórios, onde para tirar uma amostra aleatória dessa distribuição, eu uso: randon.n <- runif(1) optimize( function(x) abs(cum.beta.ph.p.petzold(x)-randon.n), interval = c(0, pi), tol = 1e-3)$minimum Agradeço se alguém puder oferecer uma sugestão. Obrigado, -- Alexandre Castagna Mourão e Lima Laboratório de Radioecologia e Mudanças Globais - UERJ Rua São Francisco Xavier, 524, Maracanã Pav. Haroldo Lisboa da Cunha, Sub-solo

1 1

Gráfico de Ishikawa
by Edson Lira 13 Ago '14

13 Ago '14

Caros liteiros, estou tentando usar a rotina abaixo: library(qcc) # Ativar o pacote qcc cause.and.effect(cause= list(Método.de.foto.cópia= c('Colocação do Original', 'Tempo de secagem'), Liquido= c('Qualidade', 'Qualidade quando novo', 'Sujeira'), Papel.de.cópia= c('Qualidade do papel', 'Foto sensibilidade'), Condições.ambientais= c('Mesa suja', 'Mãos sujas'), Papel.original= c('Dobras', 'Nitidez', 'Transparência'), Copiadora= c('Velocidade', 'Condições do enrolamento', 'Potência da lâmpada')), title=' Gestão da Qualidade \n Diagrama de Causas e Efeito', effect= 'Cópia de má qualidade', cex= c(0.7, 0.5, 1)) Como faço para excluir os pontos que separaram as palavras no gráfico, por exemplo como faço para aparecer Método de foto cópia e não Método.de.foto.cópia. A rotina exige o ponto entre as palavras das causas principais. Bom dia a todos. [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 2

ANOVA de fatorial com tratamento adicional
by helcio duarte pereira 12 Ago '14

12 Ago '14

Pessoal, vocês teriam outro pacote para indicar-me, além do ExpDes, que faz analise de variância de experimentos em esquema fatorial com tratamentos adicionais? Att. Helcio

2 2

odds ggplot2
by geovane barbosa 11 Ago '14

11 Ago '14

Ola pessoal boa tarde Estou tentando fazer o gráfico do odds de acordo com a tabela abaixo. var or lclucl A 1.5521.4011.719 B1.7891.4912.146 C 1.5611.2871.893 D1.4321.2371.658 E 1.2591.0251.548 F1.3631.1121.671 G1.5541.2581.920 H2.1571.7562.651 com o seguinte comando. dat <- read.table("tabela1.txt", header = T) ggplot(dat, aes(x = variaveis, y = or, ymin = lcl, ymax = ucl))+ ylab("Odds ratio & 95% CI") + geom_hline(aes(yintercept = 1)) + scale_color_discrete(name = "Lag") + xlab("") Erro em as.environment(where) : 'where' faltante Ele da o erro acima, o que eu faço? abraços

2 1

correção de erro package DEMEtics
by Marco Aurélio 11 Ago '14

11 Ago '14

Boa tarde, Alguém poderia me auxiliar a corrigir o erro abaixo. library(DEMEtics) D.Jost("genotipagemgerais", bias = "correct", object = FALSE, format.table = TRUE, pm = "pairwise", statistics = "CI", bt = 1000) Error in file(file, "rt") : cannot open the connection In addition: Warning message: In file(file, "rt") : cannot open file 'genotipagemgerais': No such file or directory Att., Marco Aurélio Desimoni Dias Zootecnista/UFLA Mestre em Produção e Nutrição de Não-Ruminantes/UFLA Doutorando em Produção e Nutrição de Não-Ruminantes/UFLA

3 3

Re: [R-br] Integral dupla no R
by Artur Lourival da Fonseca Machado 11 Ago '14

11 Ago '14

2 1

Integral dupla no R
by Jodavid Ferreira 11 Ago '14

11 Ago '14

estou precisando fazer uma integral de uma distribuição conjunta no R, e para isso é necessário fazer uma integral dupla. Alguém saberia me responder como posso calcular integral dupla no R? -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://rccpbjovem.com Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Bacharelando em Estatística - Universidade Federal da Paraíba - www.de.ufpb.br Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://blog.jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

1 0

Calcular distância entre cidades usando o nome (ASANTOS)
by João Augusto Batista 11 Ago '14

11 Ago '14

Olá Alexandre! estava imaginando > criar uma lista com dois objetos, na primeira com nomes de cidades que > corresponde a city from e um segundo objeto city to e a combinação > destes me dariam as distâncias. > Creio que existam 2 maneiras, uma por força bruta e outra mais inteligente. Veja nesta discussão <http://stackoverflow.com/questions/17361909/determing-the-distance-between-…> um código com ggmaps e busca por ZIP code nos EUA. Pode ser um bom ponto de partida. Do lado da força bruta, eu posso te dizer o meu passo a passo (não tenho mais o código, já tem tempo!). Quando tive esse problema, gerei um código que media a distância entre cidades a partir dos dados do IBGE e de uma fórmula simples, que você encontra na internet facilmente (existem duas, mas acho essa a melhor <http://staff.on.br/jlkm/geopath/>). Os arquivos shape do IBGE sobre os municípios possuem as coordenadas das cidades. Extrai essa informação em um data frame com 4 variáveis (coord x e y, nome da cidade e nome da UF). Depois, o script gerava um data frame com cada combinação possível de pares de cidades (sim, é gigantesco) e calculava a distância. Boa sorte e diga se deu certo! Abraços, *João Augusto P. B.* +55 21 992 453 599 FGV - Rio de Janeiro - Brazil LinkedIn Profile <http://www.linkedin.com/in/joaoaugustopb/en>

1 0