- R-br - DINF UFPR

Variável independente quantitativa em uma Análise de Discriminante Canonica
by Alisson Lucrecio 02 Mai '14

02 Mai '14

Caro Colegas, Boa noite. Um variável independente quantitativa (Tempo em dias) no modelo de análise de discriminante canonica será tratada como fator, fator ordenando ou outra forma? Obrigado. -- Alisson Lucrecio da Costa

1 0

Seleção de variáveis depedentes para Análise de Discriminate Canonica
by Alisson Lucrecio 02 Mai '14

02 Mai '14

Caro Colegas, Boa noite. Eu preciso fazer seleção de variáveis depedentes (repostas) para o modelo de Análise de Discrimiante Canonica. Qual é a forma mais fácil de fazer a seleção dessas variáveis? Obrigado -- Alisson Lucrecio da Costa

1 0

Kruskal-Wallis para grupos com apenas 1 elemento
by Letícia Pinheiro 02 Mai '14

02 Mai '14

Prezados, preciso fazer uma comparação entre valores de doses letais (DL50) para 12 populações diferentes, buscando saber se há diferenças entre as populações. Esses valores de DL50 foram encontrados a partir de uma análise Probit. Resumindo, tenho apenas 1 valor para cada população (1 valor para cada grupo) a ser comparada e gostaria de saber como usar o Kruskal-Wallis (é possível?) nesse caso. Obrigada. -- Letícia Cavalari Pinheiro

5 8

[OFF-TOPIC] Workshop Software Carpentry - FURG, Brasil
by Fernando Mayer 01 Mai '14

01 Mai '14

Amigos do R, entre os dias 12 e 13 de maio ocorrerá o primeiro workshop do Software Carpentry no Brasil, na Universidade Federal do Rio Grande (FURG). Quem tiver interesse em participar, entre em contato. Maiores informações em http://r-gaia-cs.github.io/2014-05-12-furg e na chamada abaixo: ##### Entre os dias 12 e 13 de maio de 2014, estaremos realizando um workshop do Software Carpentry na FURG, e gostariamos de convidar todos os interessados. O Software Carpentry é um projeto sem fins lucrativos e que faz parte do Mozilla Science Lab. O objetivo do projeto é ensinar conceitos básicos de computação para pesquisadores, para que eles possam fazer mais em menos tempo. No workshop serão abordados os seguintes temas - Introdução ao shell Unix - Introdução à programação usando o R - Controle de versões usando Git (e GitHub) Os workshops do Software Carpentry já foram realizados em diversos países como Austrália, Canadá, Estados Unidos e Reino Unido, e pela primeira vez será realizado em um país da América do Sul. Os instrutores são voluntários e aprenderam a usar estas ferramentas para otimizar seus trabalhos: Raniere Silva é estudante do curso de Matemática Aplicada e Computacional da UNICAMP, e Fernando Mayer é professor do Instituto de Matemática, Estatística e Física (IMEF) da FURG. O workshop do Software Carpentry é "hands-on", então os participantes devem trazer seus próprios laptops com os softwares necessários já instalados. Haverá um período no início do workshop para ajudar aqueles que tiveram alguma dificuldade na instalação dos programas, mas tente resolver seu problema antes procurando na internet ou enviando um email para os instrutores. Maiores informações sobre a instalação dos softwares, inscrição e detalhes sobre o workshop podem ser encontradas em http://r-gaia-cs.github.io/2014-05-12-furg. Para informações gerais sobre o Software Carpentry, acesse http://software-carpentry.org. ##### Atenciosamente, --- Fernando de Pol Mayer Laboratório de Estatística Ambiental - LEA Instituto de Matemática, Estatística e Física - IMEF Universidade Federal do Rio Grande - FURG URL: http://fernandomayer.github.io e-mail: fernando.mayer [@] {gmail.com, furg.br}

1 0

Matrizes e Classes
by David Feitosa 01 Mai '14

01 Mai '14

Bom dia! Eu gostaria de uma orientação sobre como criar uma matriz de um tipo de dado diferente. Preciso criar uma classe que armazene alguns atributos e, posteriormente, criar uma matriz desse tipo. Após a criação da matriz, atualizar os dados dos atributos de cada elemento da matriz. Já ouvi sugestões de "resolver o problema" usando uma matriz do tipo primitivo para cada atributo, mas creio que deva existir uma abordagem mais OO para o caso. Alguém poderia sugerir algo por favor? Não tenho ideia de como fazer. Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

2 5

Re: [R-br] [OFF-Topic] Distribuição Normal e teste de t
by Marcos Vital 01 Mai '14

01 Mai '14

Olá, Thales Bom, não estou respondendo a sua dúvida sobre Monte Carlo (não conheço bem o método), mas queria fazer uma observação sobre o pressuposto da normalidade. Como o seu teste é pareado, o cálculo do t é feito sobre as diferenças, e não sobre os valores brutos. Então acredito que você deva testar a normalidade das diferenças entre as variáveis x e y. Fiz o teste rapidamente com os dados que você incluiu no email, e o shapiro agora não rejeita a normalidade. Veja: dif<-x-y shapiro.test(dif) Shapiro-Wilk normality test data: dif W = 0.9716, p-value = 0.2683 Considere isso como possibilidade, então, pois é possível que você não precise sair do bom e velho teste t. Abraços Marcos Em 25 de abril de 2014 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Message: 2 > Date: Fri, 25 Apr 2014 08:55:04 +0200 > From: Tropidurus Torquatus <t.torquatus(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] [OFF-Topic] Distribuição Normal e teste de t > Message-ID: <E44FD9CC-9263-4D76-91B5-C084EDFB64C7(a)gmail.com> > Content-Type: text/plain; charset=windows-1252 > > Caros R-users, > > Coloquei como off-topic pois acredito que esta seja mais uma dúvida > conceitual do que de R ppd. > > Dado que tenho dois conjuntos de dados: > > x <- c ( 13.772 , 10.119 , 11.224 , 21.999 , 17.181 , 27.477 , 21.693 , > 26.706 , 14.576 , 19.689 , 12.494 , 19.511 , 16.785 , 16.551 , 20.780 , > 0.000 , 25.867 , 28.223 , 26.059 , 21.668 , 18.209 , 20.186 , 28.320 , > 26.946 , 16.218 , 16.814 , 16.859 , 0.000 , 0.000 , 3.986 , 1.625 , > 17.286 , 28.149 , 19.404 , 21.669 , 22.058 , 14.947 , 14.512 , 21.177 , > 29.997 , 20.958 , 15.123 , 28.932 , 23.807 , 30.589 , 10.208 , 20.076 , > 18.897 , 14.666 , 13.955) > > y <- c ( 13.902 , 10.102 , 11.301 , 22.081 , 17.234 , 27.540 , 21.741 , > 26.677 , 14.621 , 19.766 , 12.547 , 19.622 , 16.857 , 16.628 , 20.787 , > 0.000 , 25.943 , 28.313 , 26.166 , 21.724 , 18.255 , 20.307 , 28.387 , > 27.079 , 16.202 , 16.857 , 16.875 , 0.000 , 0.000 , 3.951 , 1.681 , > 17.353 , 28.273 , 19.422 , 21.724 , 22.164 , 14.902 , 14.458 , 21.174 , > 30.024 , 20.949 , 15.101 , 28.926 , 23.850 , 30.597 , 10.147 , 20.035 , > 18.873 , 14.691 , 13.892) > > As duas amostras estão pareadas. > As duas amostras falham para o teste de normalidade: > > > shapiro.test ( x ) > > Shapiro-Wilk normality test > > data: testetcompounds$tR...tM.UV.ABC > W = 0.9373, p-value = 0.01057 > > > shapiro.test ( y ) > > Shapiro-Wilk normality test > > data: testetcompounds$tR...tM.UV.ABC.Group10 > W = 0.9372, p-value = 0.01043 > > No livro Bayesian Computation with R > > Sugere que para distribuições não normais é necessário fazer o teste de > Monte Carlo > > Com a seguinte função > ======================================================= > tstatistic = function ( y , x ) > { > m = length ( y ) > n = length ( x ) > sp = sqrt( ( ( m - 1 ) * sd ( y ) ^ 2 + ( n - 1 ) * sd ( x ) ^ 2 ) / ( m > + n - 2 ) ) > t = ( mean ( y ) - mean ( x ) ) / ( sp * sqrt ( 1 / m + 1 / n) ) > return ( t ) > } > ======================================================== > Retornando > > t = -0.02242564 > > Em seguida o teste: > ======================================================== > alpha = 0.01 ; m = 50 ; n = 50 # sets alpha, m and n ( alpha is the > stated significance level, m and n are the samples sizes ) > N = 10000 # sets the number of > simulations > n.reject = 0 # counter of num. of > rejections > for ( i in 1 : N ) > { > x = rnorm ( m , mean = x , sd = x ) > > # simulates xs from population 1 > y = rnorm ( n , mean = y , sd = y ) > > # simulates ys from population 2 > t = tstatistic ( y , x ) # computes the t statistic > if ( abs ( t ) > qt ( 1 - alpha / 2 , n + m - 2) ) > n.reject = n.reject + 1 > > # reject if |t| exceeds critical > pt > } > true.sig.level = n.reject/N > > # est. is proportions of > rejections > ======================================================== > Retornando > > true.sig.level = 0.0049 > > > Problema Conceitual > > 1) Como interpretar esse resultado? > > > Porém ainda há uma diferença entre o t calculado por essa fórmula e o t > calculado com o comando t.stat (padrão do R) > > ================================================================ > t = t.test ( y , x , mu = 0 , paired = TRUE , var.equal = TRUE , > conf.level = 0.9999 ) > t > > Paired t-test > > data: y and x > t = 4.6437, df = 49, p-value = 2.594e-05 > alternative hypothesis: true difference in means is not equal to 0 > 99.99 percent confidence interval: > 0.003082235 0.067077765 > sample estimates: > mean of the differences > 0.03508 > > ================================================================ > > E finalmente? minha média da diferencia é conhecida e não é zero, assim > utilizei o seguinte comando: > > t = t.test ( y , x , mu = 0.87 , paired = TRUE , var.equal = TRUE , > conf.level = 0.9999, alternative = "less? ) > > Paired t-test > > data: y and x > t = -110.5231, df = 49, p-value < 2.2e-16 > alternative hypothesis: true difference in means is less than 0.87 > 99.99 percent confidence interval: > -Inf 0.06545182 > sample estimates: > mean of the differences > 0.03508 > > 2) Problema de Comando? como fazer a simulação de Monte Carlo na situação > acima? > > * Obs.: aceito sugestões de leitura sobre o assunto e comando alternativos > para comparar a médias > > > Muito obrigado e me desculpe pelo longo e-mail. > > > Thales > > > -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

2 1

Seleção de modelos não-aninhados com GAMM (mgcv)
by Pavel Dodonov 01 Mai '14

01 Mai '14

Saudações... Eu estou querendo comparar três modelos aditivos mistos (generalized additive mixed models), com distribuição de Poisson e um fator aleatório. Os modelos são 1)modelo nulo (intercepto), 2) modelo com uma variávei explanatória categórica, 3) modelo aditivo com uma variável explanatória contínua. Os modelos 2) e 3) portanto não são aninhados... O único jeito que conheço para comparar modelos não-aninhados é com AIC. Mas não sei como (ou se é possível) calcular o AIC para um modelo GAMM... Alguma sugestão? Estou usando os pacotes mgcv para fazer a análise e tentei usar o bbmle para calcular os AIC. Abs... - Pavel -- ------- Pavel Dodonov Biologist, PhD candidate in Ecology and Natural Resources, São Carlos Federal University (UFSCar), SP, Brazil Publications: http://ufscar-br.academia.edu/pdodonov Lattes CV: http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4234326T5 *"The highest function of ecology is understanding consequences." - Frank Herbert, Dune, 1965*

1 0

Legenda externa ao gráfico
by Helder Gramacho 30 Abr '14

30 Abr '14

Bom dia pessoal, Gostaria de uma orientação sobre como elaborar uma legenda como mostrado na figura abaixo, uma barra externa ao gráfico, pois com o código que estou executando, a legenda está dividida em classes e posicionada internamente. #Download do arquivo no dropbox (copiado do Alexandre Santos) links <- c("https://www.dropbox.com/s/o914ov34mv6r0sc/cadmio.csv") tokens <- gsub("^.*/s/","",dirname(links)) fileNames <- basename(links) newLinks <- file.path("http://dl.dropbox.com/s", tokens, fileNames); newLinks for (a in newLinks) { tryCatch(download.file(a, dest=basename(a), mode='wb'), error=function(...) print("Falha no download!"))} cd<-read.table(file="cadmio.csv", sep=",", header=T,dec=".") attach(cd) cd<-as.geodata(cbind(cd$Coord1,cd$Coord2,cd$data)) points(cd, pt.div= "quint", x.leg=c(0,2.1), y.leg=c(3.8,5.8)) [image: Imagem inline 1] Desde já agradeço, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

3 6

Fwd: Mover coordenadas em shapefile
by ASANTOS 30 Abr '14

30 Abr '14

-------- Mensagem original -------- Assunto: Mover coordenadas em shapefile Data: Tue, 29 Apr 2014 15:11:34 -0400 De: ASANTOS <alexandresantosbr(a)yahoo.com.br> Para: r-br(a)listas.c3sl.ufpr.br Bom tarde Pessoal, Estou tentando mover as coordenadas de um arquivo shapefile sem sucesso, preciso mover as coordenadas do meu polígono (que esta em UTM), 1032m à esquerda e 280m para cima, tentei fazer operações trigonométricas sobre os pontos da extensão *bdf, mas não consegui a visualização em formato *shp, alguém conhece alguma função em algum pacote para fazer isto ou teria alguma sugestão, segue CRM dos dados: #Download do controno da área em shapefile links <- c( "https://www.dropbox.com/s/1oz5lhnk96jih3w/Tamandua288.dbf", "https://www.dropbox.com/s/6b3uygyrjhq1bby/Tamandua288.shp", "https://www.dropbox.com/s/wntluhn8bi7us2f/Tamandua288.shx") tokens <- gsub("^.*/s/","",dirname(links)) fileNames <- basename(links) newLinks <- file.path("http://dl.dropbox.com/s", tokens, fileNames); newLinks for (a in newLinks) { tryCatch(download.file(a, dest=basename(a), mode='wb'), error=function(...) print("Falha no download!"))} ## Abrir e visualizar contorno require(rgdal) plotRGB(landa) contorno_line <- readOGR(dsn="/Documents and Settings/ASANTOS/Desktop/Percevejo GIS", layer="Tamandua288") CRS.new <- CRS("+proj=utm +zone=23 +south +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0") contorno_line.WGS84 <- spTransform(contorno_line, CRS.new) projection(landa) projection(contorno_line.WGS84) lines(contorno_line.WGS84,col="red") Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 4

Problema com abordagem inflacionada de zeros
by ASANTOS 29 Abr '14

29 Abr '14

Boa dia Pessoal, Tenho alguns dados que estou com dificuldade de encontrar uma abordagem de análise, tenho uma variável n.insetos (são dados de contagem de insetos, mas tenho várias parcelas com 0 insetos) e a variável resposta n.ataques (são dados de contagem do número de plantas atacadas, também com zeros), o meu problema é que: quando não tenho insetos, não tenho ataque, mas também tenho as vezes insetos, mas não tenho ataque. Tentei pela abordagem de Poisson inflacionado por zeros, sem sucesso (minha variável explicativa tem inúmeros zeros e é uma informação com implicação biológica que não pode ser descartada). Alguém saberia me dizer que abordagem eu poderia tentar para modelar essa relação? Segue CRM com dados parciais e a tentativa de ajuste inflacionada por zeros: #------------------------------------------------------------------------------- # Definições da sessão. rm(list=ls()) require(pscl) require(VGAM) #------------------------------------------------------------------------------- #Dados n.insetos<-c(4,3,1,2,2,2,1,4,1,1,1,0,2,2,3,3,2,6,2,3,3,4,3,0,2,5,6,3,3,3,5,2,1, 2,1,1,0,2,2,1,3,6,3,1,2,0,3,1,0,4,2,2,3,2,0,2,2,2,1,2,4,4,3,5,2,1,5,2,3,2,5,3,3, 3,4,3,2,5,3,4,1,2,2,2,1,2,3,1,4,2,4,4,1,2,2,9,3,1,2,2,2,1,3,3,3,1,1,2,5,2,1,1,0, 5,0,4,0,3,2,2,3,4,1,1,2,1,1,2,3,3,2,5,2,2,2,3,1,2,6,1,2,3,1,2) n.ataques<-c(0,1,3,0,1,2,1,1,2,2,0,0,1,0,1,1,0,0,2,0,1,0,2,0,0,1,3,0,1,0,0,1,1, 1,0,1,0,0,1,1,3,1,1,1,2,0,1,1,0,0,0,1,0,0,0,1,1,0,0,1,2,0,1,2,0,1,1,2,0,1,0,0,1, 2,1,1,3,1,3,2,0,1,0,0,1,0,2,2,2,2,2,0,1,0,2,3,2,2,3,0,1,2,1,1,0,1,2,1,2,0,1,1,0, 1,0,0,0,1,0,0,1,0,1,0,0,2,2,3,1,0,2,1,4,1,1,0,0,3,0,2,0,2,1,0) plot(n.insetos, n.ataques) #------------------------------------------------------------------------------- #------------------------------------------------------------------------------- # Modelo completo. compl.mod <- zeroinfl(n.ataques~n.insetos|n.insetos) coef(compl.mod) summary(compl.mod) #------------------------------------------------------------------------------- Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 1

for e If
by Evelyn Alecrim 27 Abr '14

27 Abr '14

Boa tarde, Tenho a seguinte situação, tenho três variáveis: Sit, d1 e d2. Gostaria que um valor (k) fosse calculado levando em consideração o valor de Sit, ou seja, para cada Sit uma fórmula diferente seria usada, como abaixo. Sit=c(3, 2, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1) d1=c(1.27, 5.80, 2.42, 3.11, 5.25, 4.74, 5.42, 5.14, 5.27, 5.61, 6.09, 6.38, 6.34, 8.60, 7.36, 9.96, 10.86, 11.35, 8.94, 9.19, 8.66, 8.98, 9.86, 9.94, 11.59, 15.61, 6.28, 7.83, 7.94, 7.67, 8.62, 9.43, 9.00) d2=c(9.12, 10.81, 10.21, 7.11, 7.23, 7.62, 7.69, 7.59, 5.73, 5.39, 5.35, 5.42, 5.32, 3.84, 2.89, 0.70, 1.65, 2.32, 4.01, 4.16, 4.32, 4.74, 5.12, 3.01, 1.91, 7.33, 3.81, 2.27, 2.22, 4.44, 4.76, 7.35, 7.35) for(i in 1:seq_along(Sit)){ if(Sit==1){ k<-(100 + d1 ^ 2 - d2 ^ 2) / 20 } else if(Sit==2){ k<-(10 + (100 + d1 ^ 2 - d2 ^ 2) / 20) } else { k<-(400 + d1 ^ 2 - d2 ^ 2) / 40 } print(k) } No entanto, o R me retorna a seguinte mensagem: [1] 7.960963 7.919597 7.540307 8.978000 9.382240 9.110080 9.256007 9.220288 9.873500 10.060500 10.211640 10.283200 10.297330 11.480360 11.145437 [16] 12.467790 12.880427 13.086003 11.596087 11.678762 11.408330 11.454320 11.775130 12.243588 13.267000 14.748580 10.623057 11.403900 11.452880 10.977882 [31] 11.291170 10.872560 10.674437 Mensagens de aviso perdidas: 1: In 1:seq_along(Sit) : expressão numérica tem 33 elementos: somente o primeiro usado 2: In if (Sit == 1) { : a condição tem comprimento > 1 e somente o primeiro elemento será usado 3: In if (Sit == 2) { : a condição tem comprimento > 1 e somente o primeiro elemento será usado Parece que ele está calculando usando apenas a terceira fórmula. Alguém pode dar uma luz? Obrigada, Emily

7 8

Estimador de Máxima verossimilhança para parâmetros da dist. q-exponencial
by Romero Luiz M. Sales Filho 27 Abr '14

27 Abr '14

Caros amigos, estou precisando encontrar os estimadores para os parâmetros da dist. q-ponencial, cuja pdf é dada da seguinte forma: [image: { (2-q) \lambda e_q^{-\lambda x}}] onde, [image: e_q(x) = [1+(1-q)x]^{1 \over 1-q}]. Estou tendo problemas pois utilizando o pacote Optim do R, não ocorre convergência em alguns métodos, na verdade só há convergência usando o método SANN, e mesmo assim as estimativas estão muito diferentes dos parâmetros que foram utilizados para gerar uma amostra aleatória da dist. q-exp. O código utilizado no R é o seguinte: ## Gerando Valores de uma q-exponencial n<-10000 u<-runif(n) z<-1.8 beta<- 1.2 dados1 = (1/(beta*(z-2)))*((u^((1-z)/(2-z))-1)/((1-z)/(2-z))) ## Definindo a função de verossimilhança vero <- function(par,x){ q = par[1] lambda = par[2] ((sum(log(1 - ((1 - q)*x*lambda))))/(1 - q)) + ((n)*log(2 - q)) + (n*log(lambda)) } ## Otimizando os parâmetros - o método SANN sempre converge. optim(par=c(1.8,4),fn=vero, method="SANN",x=dados1 ) E a saída é a seguinte: $par [1] 1.000014 3.807168 $value [1] -284366598 $counts function gradient 10000 NA $convergence [1] 0 $message NULL Houve 50 ou mais avisos (use warnings() para ver os primeiros 50) Alguém pode me ajudar nessa questão, uma vez que as estimativas dos parâmetros estão muito diferentes dos parâmetros originais? Obrigado! Romero.

5 4

Erro: estouro da pilha de contexto (O que fazer?)
by Alexandre Loures 26 Abr '14

26 Abr '14

Boa noite! Estou usando a rotina abaixo para construir uma dummy para as fronteiras dos estados brasileiros. Porém, quando tendo rodar o R reporta o seguinte erro: Erro: estouro da pilha de contexto na linha 18 O que devo fazer? mydata$contiguity<-ifelse((origin=="ac") & (destination=="am"),1, ifelse((origin=="am") & (destination=="ro"),1,ifelse((origin=="am") & (destination=="rr"),1,ifelse((origin=="am") & (destination=="pa"),1,ifelse((origin=="am") & (destination=="mt"),1, ifelse((origin=="ro") & (destination=="mt"),1,ifelse((origin=="ro") & (destination=="am"),1, ifelse((origin=="rr") & (destination=="am"),1,ifelse((origin=="rr") & (destination=="pa"),1, ifelse((origin=="pa") & (destination=="mt"),1,ifelse((origin=="pa") & (destination=="am"),1,ifelse((origin=="pa") & (destination=="am"),1,ifelse((origin=="pa") & (destination=="rr"),1, ifelse((origin=="pa") & (destination=="ap"),1,ifelse((origin=="pa") & (destination=="ma"),1,ifelse((origin=="pa") & (destination=="to"),1,ifelse((origin=="pa") & (destination=="mt"),1, ifelse((origin=="ap") & (destination=="pa"),1, ifelse((origin=="ma") & (destination=="pa"),1,ifelse((origin=="ma") & (destination=="to"),1,ifelse((origin=="ma") & (destination=="pi"),1, ifelse((origin=="pi") & (destination=="ma"),1,ifelse((origin=="pi") & (destination=="to"),1,ifelse((origin=="pi") & (destination=="pe"),1,ifelse((origin=="pi") & (destination=="ce"),1, ifelse((origin=="ce") & (destination=="pi"),1,ifelse((origin=="ce") & (destination=="pe"),1,ifelse((origin=="ce") & (destination=="pb"),1,ifelse((origin=="ce") & (destination=="rn"),1, ifelse((origin=="rn") & (destination=="ce"),1,ifelse((origin=="rn") & (destination=="pb"),1, ifelse((origin=="pb") & (destination=="rn"),1,ifelse((origin=="pb") & (destination=="ce"),1,ifelse((origin=="pb") & (destination=="pe"),1, ifelse((origin=="pe") & (destination=="pb"),1,ifelse((origin=="pe") & (destination=="ce"),1,ifelse((origin=="pe") & (destination=="pi"),1,ifelse((origin=="pe") & (destination=="ba"),1, ifelse((origin=="pe") & (destination=="al"),1, ifelse((origin=="al") & (destination=="pe"),1,ifelse((origin=="al") & (destination=="ba"),1,ifelse((origin=="al") & (destination=="se"),1, ifelse((origin=="se") & (destination=="al"),1,ifelse((origin=="se") & (destination=="ba"),1, ifelse((origin=="ba") & (destination=="se"),1,ifelse((origin=="ba") & (destination=="al"),1,ifelse((origin=="ba") & (destination=="pe"),1,ifelse((origin=="ba") & (destination=="pi"),1, ifelse((origin=="ba") & (destination=="to"),1,ifelse((origin=="ba") & (destination=="go"),1,ifelse((origin=="ba") & (destination=="mg"),1,ifelse((origin=="ba") & (destination=="es"),1, ifelse((origin=="to") & (destination=="pa"),1,ifelse((origin=="to") & (destination=="mt"),1,ifelse((origin=="to") & (destination=="go"),1,ifelse((origin=="to") & (destination=="ba"),1, ifelse((origin=="to") & (destination=="pi"),1,ifelse((origin=="to") & (destination=="ma"),1, ifelse((origin=="mt") & (destination=="ro"),1,ifelse((origin=="mt") & (destination=="am"),1,ifelse((origin=="mt") & (destination=="pa"),1,ifelse((origin=="mt") & (destination=="to"),1, ifelse((origin=="mt") & (destination=="go"),1,ifelse((origin=="mt") & (destination=="ms"),1, ifelse((origin=="ms") & (destination=="mt"),1,ifelse((origin=="ms") & (destination=="go"),1,ifelse((origin=="ms") & (destination=="mg"),1,ifelse((origin=="ms") & (destination=="sp"),1, ifelse((origin=="ms") & (destination=="pr"),1, ifelse((origin=="go") & (destination=="mg"),1,ifelse((origin=="go") & (destination=="ba"),1,ifelse((origin=="go") & (destination=="to"),1,ifelse((origin=="go") & (destination=="mt"),1, ifelse((origin=="go") & (destination=="ms"),1,ifelse((origin=="go") & (destination=="df"),1, ifelse((origin=="mg") & (destination=="sp"),1,ifelse((origin=="mg") & (destination=="rj"),1,ifelse((origin=="mg") & (destination=="es"),1,ifelse((origin=="mg") & (destination=="ba"),1, ifelse((origin=="mg") & (destination=="go"),1,ifelse((origin=="mg") & (destination=="ms"),1, ifelse((origin=="sp") & (destination=="rj"),1,ifelse((origin=="sp") & (destination=="mg"),1,ifelse((origin=="sp") & (destination=="ms"),1,ifelse((origin=="sp") & (destination=="pr"),1, ifelse((origin=="rj") & (destination=="mg"),1,ifelse((origin=="rj") & (destination=="mg"),1,ifelse((origin=="rj") & (destination=="es"),1, ifelse((origin=="es") & (destination=="ba"),1,ifelse((origin=="es") & (destination=="mg"),1,ifelse((origin=="es") & (destination=="rj"),1, ifelse((origin=="pr") & (destination=="sp"),1,ifelse((origin=="pr") & (destintaion=="sc"),1,ifelse((origin=="pr") & (destination=="ms"),1, ifelse((origin=="sc") & (destination=="pr"),1,ifelse((origin=="sc") & (destination=="rs"),1, ifelse((origin=="rs") & (destination=="sc"),1, ifelse((origin=="df") & (destination=="go"),1,0)))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))))) -- *Alexandre Rodrigues Loures*

2 2

Criar Dummy
by Alexandre Loures 25 Abr '14

25 Abr '14

Olá pessoal! Estou trabalhando numa base de dados e preciso criar algumas dummies, sendo que comando que geralmente utilizo é: mydata$dac<-ifelse(mydata$origin=="ac",1,0) Mas como preciso criar várias dummies, gostaria de saber se há como fazer isso de uma só vez? Isto é, utilizando um único comando. As dummies são para origin, destination e ano. Ou seja: Dummies para as regiões de origem (origin) Dummies para as regiões de destino (destination) Dummies para as regiões de origem com relação ao ano (Por exemplo: origin ac e ano 1998 = 1, ...) Dummies para as regiões de destino com relação ao ano (Por exemplo: destination ac e ano 1998 = 1, ...) Dummy para ano > str(mydata) 'data.frame': 4158 obs. of 9 variables: $ origin : Factor w/ 27 levels "ac","al","am",..: 1 1 1 1 1 1 1 1 1 1 ... $ destination: Factor w/ 78 levels "ac","africa",..: 3 3 5 5 7 7 13 13 17 17 ... $ ano : int 1998 1999 1999 1998 1999 1998 1999 1998 1998 1999 ... $ pib_i : num 1694329 1866542 1866542 1694329 1866542 ... $ pop_i : int 514050 527937 527937 514050 527937 514050 527937 514050 514050 527937 ... $ value : num 785 115 32290 0 0 ... $ pib_j : num 6676333 6935080 13933920 13198510 2128307 ... $ pop_j : num 2688117 2713203 2580860 2520684 439781 ... $ dist : num 5039 5039 1445 1445 NA ... Desde já muito obrigado! -- *Alexandre Rodrigues Loures*

3 5

Instalar pacotes no software R
by Jodavid Ferreira 25 Abr '14

25 Abr '14

Como posso instalar todos os pacotes de uma vez no software R no linux? -- 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 Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

5 6

[OFF-Topic] Distribuição Normal e teste de t
by Tropidurus Torquatus 25 Abr '14

25 Abr '14

Caros R-users, Coloquei como off-topic pois acredito que esta seja mais uma dúvida conceitual do que de R ppd. Dado que tenho dois conjuntos de dados: x <- c ( 13.772 , 10.119 , 11.224 , 21.999 , 17.181 , 27.477 , 21.693 , 26.706 , 14.576 , 19.689 , 12.494 , 19.511 , 16.785 , 16.551 , 20.780 , 0.000 , 25.867 , 28.223 , 26.059 , 21.668 , 18.209 , 20.186 , 28.320 , 26.946 , 16.218 , 16.814 , 16.859 , 0.000 , 0.000 , 3.986 , 1.625 , 17.286 , 28.149 , 19.404 , 21.669 , 22.058 , 14.947 , 14.512 , 21.177 , 29.997 , 20.958 , 15.123 , 28.932 , 23.807 , 30.589 , 10.208 , 20.076 , 18.897 , 14.666 , 13.955) y <- c ( 13.902 , 10.102 , 11.301 , 22.081 , 17.234 , 27.540 , 21.741 , 26.677 , 14.621 , 19.766 , 12.547 , 19.622 , 16.857 , 16.628 , 20.787 , 0.000 , 25.943 , 28.313 , 26.166 , 21.724 , 18.255 , 20.307 , 28.387 , 27.079 , 16.202 , 16.857 , 16.875 , 0.000 , 0.000 , 3.951 , 1.681 , 17.353 , 28.273 , 19.422 , 21.724 , 22.164 , 14.902 , 14.458 , 21.174 , 30.024 , 20.949 , 15.101 , 28.926 , 23.850 , 30.597 , 10.147 , 20.035 , 18.873 , 14.691 , 13.892) As duas amostras estão pareadas. As duas amostras falham para o teste de normalidade: > shapiro.test ( x ) Shapiro-Wilk normality test data: testetcompounds$tR...tM.UV.ABC W = 0.9373, p-value = 0.01057 > shapiro.test ( y ) Shapiro-Wilk normality test data: testetcompounds$tR...tM.UV.ABC.Group10 W = 0.9372, p-value = 0.01043 No livro Bayesian Computation with R Sugere que para distribuições não normais é necessário fazer o teste de Monte Carlo Com a seguinte função ======================================================= tstatistic = function ( y , x ) { m = length ( y ) n = length ( x ) sp = sqrt( ( ( m - 1 ) * sd ( y ) ^ 2 + ( n - 1 ) * sd ( x ) ^ 2 ) / ( m + n - 2 ) ) t = ( mean ( y ) - mean ( x ) ) / ( sp * sqrt ( 1 / m + 1 / n) ) return ( t ) } ======================================================== Retornando t = -0.02242564 Em seguida o teste: ======================================================== alpha = 0.01 ; m = 50 ; n = 50 # sets alpha, m and n ( alpha is the stated significance level, m and n are the samples sizes ) N = 10000 # sets the number of simulations n.reject = 0 # counter of num. of rejections for ( i in 1 : N ) { x = rnorm ( m , mean = x , sd = x ) # simulates xs from population 1 y = rnorm ( n , mean = y , sd = y ) # simulates ys from population 2 t = tstatistic ( y , x ) # computes the t statistic if ( abs ( t ) > qt ( 1 - alpha / 2 , n + m - 2) ) n.reject = n.reject + 1 # reject if |t| exceeds critical pt } true.sig.level = n.reject/N # est. is proportions of rejections ======================================================== Retornando true.sig.level = 0.0049 Problema Conceitual 1) Como interpretar esse resultado? Porém ainda há uma diferença entre o t calculado por essa fórmula e o t calculado com o comando t.stat (padrão do R) ================================================================ t = t.test ( y , x , mu = 0 , paired = TRUE , var.equal = TRUE , conf.level = 0.9999 ) t Paired t-test data: y and x t = 4.6437, df = 49, p-value = 2.594e-05 alternative hypothesis: true difference in means is not equal to 0 99.99 percent confidence interval: 0.003082235 0.067077765 sample estimates: mean of the differences 0.03508 ================================================================ E finalmente… minha média da diferencia é conhecida e não é zero, assim utilizei o seguinte comando: t = t.test ( y , x , mu = 0.87 , paired = TRUE , var.equal = TRUE , conf.level = 0.9999, alternative = "less” ) Paired t-test data: y and x t = -110.5231, df = 49, p-value < 2.2e-16 alternative hypothesis: true difference in means is less than 0.87 99.99 percent confidence interval: -Inf 0.06545182 sample estimates: mean of the differences 0.03508 2) Problema de Comando… como fazer a simulação de Monte Carlo na situação acima? * Obs.: aceito sugestões de leitura sobre o assunto e comando alternativos para comparar a médias Muito obrigado e me desculpe pelo longo e-mail. Thales

1 0

[Off] Pacotes/Artigos na Journal of Statistical Software
by FHRB Toledo 24 Abr '14

24 Abr '14

Recomendações do editor... interessante! http://www.r-bloggers.com/publishing-an-r-package-in-the-journal-of-statist… att, FH

1 0

gravando png
by Hélio Gallo Rocha 23 Abr '14

23 Abr '14

Caros, bom dia. Estou tentando aumentar a resolução de um plot utilizando o comando 'png'. Ele gera um arquivo, mas com 0kb e trava o visualizador de plot, ou seja, se tentar plotar outra imagem, ele não gera, tenho de fechar o trabalho Estou usando o comando png(filename = "k:/2014/imas/image.png", width = 1890, height = 2550, bg = "white", res = 300, restoreConsole = TRUE) Se alguém puder auxiliar, agradeço desde já -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Help
by Bruce Kelly 22 Abr '14

22 Abr '14

Olá a todos, estou com um pequeno problema, estou usando o pacote 'rela' usando a função paf para aplica a técnica de analise fatorial, pois esta roda com NA, diferentemente da factanal. Bem quando tenho meu banco de dados .csv, tendo rodar mas ele da a seguinte mensagem de erro: Your dataset is not a numeric object, gostaria de tornar esse banco numérico, como posso fazer isso. Obrigado pela atenção de todos. Bruce Kelly

3 4

R activity around the world
by Fernando Souza 22 Abr '14

22 Abr '14

Caros colegas, foi publicado hoje no R-bloggers, que mostra o nível de atividade do R ao redor do mundo, mostrando em um mapa interativo onde há mais atividades e tipo de atividades. Muito legal! http://blog.rapporter.net/2014/04/r-activity-around-world.html?utm_source=f…

1 0

Retirar strings
by Edson Lira 22 Abr '14

22 Abr '14

Tenho a variável idade em um banco de dados da forma abaixo: 1 ano(s) 10 ano(s) 11 ano(s) 12 ano(s) 13 ano(s) 14 ano(s) Quero retirar a palavra ano(s), transformando-a de fator para numérica Olhei o help da ??stringsAsFactors, não consegui encontrar resposta. [ ]'s. Edson Lira Estatístico Manaus-Amazonas

5 7

Importar dados
by Alexandre Loures 21 Abr '14

21 Abr '14

Boa tarde pessoal! Tenho uma base de dados em .csv em que algumas variáveis são numéricas, todavia quando importo para o R duas dessas variáveis estão sendo importados como factor (pib_j e dist). Alguém saberia me dizer como faço para importar mantendo as classes das variáveis? Pois quando tendo rodar um ols com o log das variáveis dá erro. A solução encontrada por mim foi alterar variável por variável, mas será que há alguma forma de importar já dizendo para manter a classe das variáveis? Desde já muito obrigado! > dados<-read.csv("states.csv",header=TRUE,sep=";",dec=",") > > str(dados) 'data.frame': 1404 obs. of 9 variables: $ origin : Factor w/ 27 levels "ac","al","am",..: 1 1 1 1 1 1 1 1 1 1 ... $ destination: Factor w/ 27 levels "ac","al","am",..: 2 2 3 3 4 4 5 5 6 6 ... $ ano : int 1998 1999 1999 1998 1999 1998 1999 1998 1998 1999 ... $ pib_i : num 1694329 1866542 1866542 1694329 1866542 ... $ pop_i : int 514050 527937 527937 514050 527937 514050 527937 514050 514050 527937 ... $ value : int 785 115 32290 0 0 0 4259 24289 5318 6059 ... $ pib_j : Factor w/ 54 levels "10172207.9512881",..: 43 45 9 7 21 18 32 30 16 19 ... $ pop_j : int 2688117 2713203 2580860 2520684 439781 420834 12993011 12851268 7013376 7106605 ... $ dist : Factor w/ 343 levels ".","1004","1014",..: 281 281 38 38 1 1 259 259 295 295 ... > reg<-lm(value~log(pib_i)+log(pib_j)+log(pop_i)+log(pop_j)+log(dist),data=dados) Erro em Math.factor(c(43L, 45L, 9L, 7L, 21L, 18L, 32L, 30L, 16L, 19L, : log not meaningful for factors -- *Alexandre Rodrigues Loures*

3 3

Soundex
by Fátima Lima Paula 21 Abr '14

21 Abr '14

Olá Daniel Marcelino. Estou terminando minha tese e como usei o soundex de sua autoria para fazer o linkage, gostaria de saber se no anexo 1 que terá os comandos do R mais importantes que utilizei, posso colocar o soundex (naturalmente citando o seu nome). Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

1 0

função para cálculo da área abaixo da curva de progresso da doença
by helcio duarte pereira 21 Abr '14

21 Abr '14

Vocês sabem de outra função sem ser a audpc() do pacote agricolae para cálculo de área abaixo da curva de progresso da doença, em estudos de severidade de doenças foliares?

2 1

Re: [R-br] Link do R para android
by Gilenio Borges Fernandes 21 Abr '14

21 Abr '14

Prezados(as) Circulou na lista um link do R para android. Alguém poderia reenviar? Grato Gilenio Fernandes Em 21 de abril de 2014 12:00, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. função para cálculo da área abaixo da curva de progresso da > doença (helcio duarte pereira) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 20 Apr 2014 16:11:28 -0700 (PDT) > From: helcio duarte pereira <hhelciopassos(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] função para cálculo da área abaixo da curva de > progresso da doença > Message-ID: > <1398035488.68828.YahooMailNeo(a)web122102.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > > Vocês sabem de outra função sem ser a audpc() do pacote agricolae para > cálculo de área abaixo da curva de progresso da doença, em estudos de > severidade de doenças foliares? > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140420/d6a1de8c/atta… > > > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 40, assunto 20 > ***************************************** > -- Gilenio Borges Fernandes Professor Associado IV Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6340/6341/6337 Fax: (071)3283-6336 URL: http://lattes.cnpq.br/6764860618464860

2 1

Integrando o R com o QGIS.
by Paulo Pimenta 19 Abr '14

19 Abr '14

Sou um usuário do R há algum tempo e um novo usuário no QGIS, por isso ainda estou tendo dificuldades para integrar ambas as ferramentas. Como estou a usar a versão do QGIS 2.x "Valmiera" não há necessidade de usar o pacote SEXTANTE, pois o R já vem integrado por meio do pacote PROCESSING. Meu problema esta na hora de abrir e rodar os scripts. sempre que faço isso uma mensagem é mostrada dizendo que ele não esta configurado corretamente ou que não esta instalado no meu computador. Gostaria de saber se o problema esta nos PATH indicando o caminho para o QGIS do R ou em alguma configuração do QGIS/R. Agradeço desde já a todos!!

2 1

Ajuda quasi likelihood
by Wagner Bonat 18 Abr '14

18 Abr '14

Caros, Alguem tem ideia do que acontece nesta situacao quasi(link = "log", variance = "mu") Family: quasi Link function: log Resultado perfeito, agora preciso deixar generico para usar a funcao quasi dentro de outra funcao para obter valores iniciais, entao pensei algo assim link = "log" variance = "mu" quasi(link = link, variance = variance) Error in quasi(link = link, variance = variance) (from quasi.R@3553CaW#69) : 'variance' "NA" is invalid: possible values are "mu(1-mu)", "mu", "mu^2", "mu^3" and "constant" Alguma ideia ?? Obrigado -- Wagner Hugo Bonat LEG - Laboratório de Estatística e Geoinformação UFPR - Universidade Federal do Paraná

4 6

Desenho experimental- fatorial com medidas repetidas no tempo
by Fernando Souza 18 Abr '14

18 Abr '14

Caros amigos, estou com a seguinte dúvida Estou planejando um projeto de pesquisa onde tenho um delineamento inteiramente casualizado com fatorial. Os tratamentos são duas fontes de leguminosas (A e B) e 4 tratamentos (1,2,3,4) com 5 repetições. Eu pretendo fazer uma avaliação comportamental e para isso será realizado observações a cada 10 minutos durante um período de 24 horas (ou seja 144 observações por animais) a minha dúvida é como ficará a tabela de análise de variância desse experimento? Dic com fatorial e medidas repetidas no tempo? Qual seria a melhor forma de analisá-lo. Não há como blocar os animais pois cada grupo será abatido no final do período experimental.

1 0

Falha ao instalar pacote gmp
by Fernando Souza 18 Abr '14

18 Abr '14

Caros amigos, Estou com o seguinte problema, estou tentando instalar um pacote e ele depende de outro chamado "gmp" http://cran.r-project.org/web/packages/gmp/index.html. Acontece que não estou conseguindo fazer a instalação este pacote de forma automática ou manual. O R retorna a seguinte mensage Installing package into ‘/home/fernando/R/i686-pc-linux-gnu-library/3.1’ (as ‘lib’ is unspecified) * installing *source* package ‘gmp’ ... ** package ‘gmp’ successfully unpacked and MD5 sums checked creating cache ./config.cache checking for __gmpz_ui_sub in -lgmp... no configure: error: GNU MP not found, or not 4.1.4 or up, see http://gmplib.org ERROR: configuration failed for package ‘gmp’ * removing ‘/home/fernando/R/i686-pc-linux-gnu-library/3.1/gmp’ Warning in install.packages : installation of package ‘/home/fernando/Downloads/gmp_0.5-3.tar.gz’ had non-zero exit status Alguem sabe como resolver isto? Abraços

2 2

Análises de regressão
by marcelo claro de souza 16 Abr '14

16 Abr '14

Bom dia Por favor, alguém poderia me informar os comandos para rodas as análises de regressão polinomial de 2 grau, regressão logarítmica e regressão exponencial? Muito obrigado pela ajuda. Abraço Marcelo

3 2

mistura sen cos para arimax
by Vinicius Brito Rocha 15 Abr '14

15 Abr '14

Senhores, estou trabalhando com uma consultoria onde eles realizam previsões usando o R via sarima. Conversando com um dos consultores ele me disse que estão colocando mistura de sen e cos como input do modelo para reforçar a sazonalidade. honestamente não consegui pensar como eles estão fazendo. Alguém sabe como fazer isso? -- *Vinicius Brito Rocha, PhD. Student* *Estatístico e Atuário * * M.Sc. Engenharia de Produção/PO*

1 0

Re: [R-br] Ler .MAP no R - Shapefile
by Éder Comunello 15 Abr '14

15 Abr '14

Mauro, bom dia! Use o segundo bloco de código informado no email pra criar um arquivo readMap.R. Melhor criar o arquivo no diretório de trabalho (que pode ser definido por setwd()). Éder Comunello ---------- Mensagem encaminhada ---------- De: Mauro Sznelwar <sznelwar(a)uol.com.br> Data: 13 de abril de 2014 19:42 Assunto: Re: [R-br] Ler .MAP no R - Shapefile Para: ecomunel(a)gmail.com Como eu carrego esta rotina? source("readMap.R")

3 6

Como faço para descobrir o poder de um teste de hipótese
by Jodavid Ferreira 15 Abr '14

15 Abr '14

Olá, Estou gerando duas amostras aleatórias de poisson: amostra1<-rpois(500,lambda=2.5) amostra2<-rpois(500,lambda=3.5) e aplicando o teste de wilcoxon e o teste t para dados pareados. wilcox.test(amostra1,amostra2, paired=T) t.test(amostra1,amostra2, paired=T) e gostaria de obter o poder de cada um dos testes acima. -- 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 Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

3 2

lme
by ana paula coelho madeira 15 Abr '14

15 Abr '14

Boa tarde. Tenho um conjunto de dados avaliando o número de ácaros em função de tempo, tratamento e planta. Considerando o tempo e a planta como efeitos aleatórios e tratamento(ácaros) como efeito fixo, escrevi o modelo: model <- lme(n_individ~acaro,random=~time|plant) Erro em lme.formula(n_individ ~ acaro, random = ~time | plant) : nlminb problem, convergence error code = 1 message = iteration limit reached without convergence (10) No entanto apareceu a mensagem de erro. Alguém sabe me informar o motivo desse erro? Desde já agradeço, Att., Ana Paula

3 2

Criar uma nova variável a partir de duas outras já existentes
by Alexandre Loures 14 Abr '14

14 Abr '14

Boa noite pessoal! Mais uma vez recorro à vocês. Desta vez minha dúvida é: tenho uma base com as variáveis ano, origin, destination e quero criar uma quarto variável pib. Sendo assim, gostaria de uma rotina da seguinte forma: quando o ano for igual a 1997 e destination igual a africa a variável pib assumirá valor igual à 1.289.989,00 e assim sucessivamente para os demais anos (1997 à 2013) e destination (51 países) - sendo que os valores para o pib mudarão de ano para ano e de país para país. Alguém poderia me ajudar? Desde já muito obrigado! -- *Alexandre Rodrigues Loures*

1 0

Off-topic: R App para Android
by Éder Comunello 14 Abr '14

14 Abr '14

Senhores, bom dia! Segue o link, caso tenham curiosidade em ver e testar uma versão do R portada para dispositivos Android (tablets e celulares). Ainda é bastante incipiente, mas tem alguma funcionalidade básica. https://play.google.com/store/apps/details?id=com.appsopensource.R Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

2 2

figura
by Hélio Gallo Rocha 14 Abr '14

14 Abr '14

Caros da Lista, bom dia. Teria como inserir uma figura em num "image"? Grato, -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

3 4

Analisando os dados da amostra do IBGE
by Rodrigo Coster 14 Abr '14

14 Abr '14

Caros, Estou com uma dúvida mais metodológica do que computacional: Qual é o delineamento amostral da amostra do IBGE que ocorre paralelo ao censo? Dentre as variáveis, temos: V0011 - Area de Ponderação V0300 - Controle (Pelo que eu consegui entender, é dado um código de controle para cada residência, mas não achei isso documentado em algum lugar) V0010 - Peso Pelo que vi, a amostra é estratificada (dado que a cidade é dividida em AEDs) e, dentro de cada estrato, é por conglomerado (visto que são entrevistados todos os moradores das casas). Estou certo? Se sim, como declaro esse delineamento no pacote survey? svydesign( ids = ~ V0300, strata = ~ V0011, weights = ~ V0010, data = dados) ?

2 3

Programa fecha sozinho
by Fátima Lima Paula 12 Abr '14

12 Abr '14

Prezados, de uma hora para outra o windows diz que surgiu um problema e que o programa não está respondendo e fecha tudo. Alguém sabe o que pode ser? Minha versão é a 3.0.3. Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

3 4

Ler .MAP no R - Shapefile
by Sérgio Henrique almeida da silva ju 12 Abr '14

12 Abr '14

Olá Alguém sabe como faço para ler arquivos .MAP no R? Esse formato é usado pelo tabwin para ler shapefiles. Sei que existe o script abaixo http://dev-pg2.datasus.gov.br/DATASUS/tabwin/rx/Importa_mapa.htm Porém quando faço > mapa = read.map("BR.MAP") Erro em readChar(zz, 25) (from #33) : entrada UTF-8 inválida em readChar() Não roda. Abraços -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 968463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

2 5

Erro na leitura de base em SQL
by Edson Lira 11 Abr '14

11 Abr '14

Caros amigos, estou lendo uma base de dados em SQL com a rotina abaixo #DOAÇÃO GERAL channel<-odbcConnect("saddw",uid="saddw",case="tolower") doacao.dados<- sqlQuery(channel, "select doacao_dw.numerobolsa, doacao_dw.datahoradoacao,doacao_dw.tipodoacao, doacao_dw.tipodoacaodirigida,doacao_dw.localcoleta,doacao_dw.localespecificocoleta, doacao_dw.apto, doacao_dw.codigoinaptidao,doador_dw.tipagemsanguinea,doador_dw.rh, doador_dw.idade, doador_dw.sexo, doador_dw.data_primeira_doacao, doador_dw.data_ultima_doacao,bairroresidencia, bairro_nao_tabelado,doador_dw.qtde_doacoes,doador_dw.municipioorigem, doador_dw.naturalidade, doador_dw.numerogente from saddw.doador_dw, saddw.doacao_dw where doador_dw.numerogente = doacao_dw.numerogente and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') >= to_date('01/01/1992','DD/MM/YYYY') and to_char(doacao_dw.datahoradoacao,'DD/MM/YYYY') <= to_date('31/12/1999','DD/MM/YYYY')" ) Que está me gerando a mensagem abaixo: Erro em as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format Os anos com problema são: 1992, 1995 e 1996 - rodei cada ano separadamente. Como consequência, não gera o data frame doacao.dados Alguém teria alguma sugestão para corrigir esse erro? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

1 0

R+SIG
by Emerson Cotta Bodevan 11 Abr '14

11 Abr '14

Pessoal, alguém poderia me indicar um bom material sobre uso do R com SIG. Como ler e trabalhar com arquivos shape. Como converter de um sistema de coordenadas geográfica para outro. (p. ex. mudar do South America 69 para o SIRGAS 2000, ou ainda do South America para o SIRGAS UTM 23, etc). Agradeço qualquer ajuda. Att., *Emerson*

3 3

Problema "Desconhecido no R"
by Rodrigo Muniz 10 Abr '14

10 Abr '14

Caros membros da lista, gostaria mais uma vez de solicitar a ajuda de vocês. Estou tendo um problema, que não sei identificar o que é. O problema é o seguinte: Tenho um data frame 20 colunas chamado "Carbono", onde um das colunas é um fator identificando os tratamentos do experimento, e tenho um outra coluna com as classes texturais do solo. Eu fiz uma análise por classe textural, então extraí um outro data frame chamado "Arenosa", utilizando a função "subset" somente para a classe arenosa, desta forma neste segundo data frame não há todos os tratamentos que estão presentes no data frame "Carbono". Porém qualquer gráfico que eu faça como por exemplo, plot(tapply(arenosa, arenosa$Tratamento, mean), aparece no gráfico todos os tratamentos do data frame "Carbono". Esse é um exemplo, mas tenho tido problemas recorrentes referente a isto. Alguém tem alguma ideia do que possa ser? A título de informação, meu sistema operacional é Windows 7 e eu utilizo o Rstudio. Desde já agradeço... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

5 7

Repetição dentro da parcela
by ROSIANA RODRIGUES ALVES 10 Abr '14

10 Abr '14

Olá Pessoal, Estou com um problema em que minha parcela é um aquario com 15 peixes, com 3 tratamentos e 3 aquários por tratamento. Como faço para identificar no meu modelo essa repetição dentro da parcela? Att, Rosiana Rodrig ues Alves Estat ística e Experi mentação Agropecuária Em brapa Pesca e Aquicu ltura Empresa Brasileira de Pesquisa Agropecuária (Embrapa) Palmas/TO +55 (63) 3229-7800/7850 | Skype: cnpasa.rosiana http://cnpasa.sede.embrapa.br | twitter.com/embrapa Confira também: www.facebook.com/agrosustentavel ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure.

2 1

Somar intervalo de valores nas colunas
by Alexandre Loures 10 Abr '14

10 Abr '14

Olá pessoal! Tenho uma base de dados com 216 observações de 10 variáveis. Sendo que essas observações são mensais, mas tenho que trabalhar com o total anual. Ou seja, eu tenho dados para os meses de 1997, 1998, 1999, 2000, ..., 2014. Mas quero os totais desses anos em cada coluna. Alguém saberia como fazer isso? Desde já muito obrigado! -- *Alexandre Rodrigues Loures*

6 11

Caracteres especiais como string
by Rodrigo Muniz 10 Abr '14

10 Abr '14

Caros membros da lista, gostaria mais uma vez de solicitar a ajuda de vocês. Estou fazendo um gráfico boxplot onde eu gostaria de nomear os rótulos com o argumento ylab="δ13C", porém o resultado no plot é "d13C". Outra coisa é que o número "13" precisa estar em sobrescrito. Como eu poderia fazer esses ajustes? Desde já agradeço... -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

4 4

Re: [R-br] Somar intervalo de valores nas colunas
by Robert Iquiapaza 08 Abr '14

08 Abr '14

Alexandre, O objeto que enviou não parece ter a estrutura que vc falou. Mas se seus dados estiverem ordenados por meses, > sp[c(1,5,89,200,216),] Sp1 Sp2 Sp3 Sp4 1 101.20315 96.83025 95.74570 104.38100 5 99.13003 98.05711 98.44344 97.72637 89 98.88583 101.97051 99.31058 102.20602 200 95.72225 100.38504 101.06742 102.75147 216 100.06324 100.50797 100.69345 101.10212 O código a seguir pode funcionar> mes=rep(1:12,(2014-1996)) > ano=rep(1997:2014,each = 12) > (Tot <-aggregate(sp, by=list(ano), FUN=sum, na.rm=TRUE)) Group.1 Sp1 Sp2 Sp3 Sp4 1 1997 1203.948 1197.491 1203.537 1207.432 2 1998 1201.671 1225.573 1191.417 1192.366 3 1999 1206.173 1206.108 1191.918 1218.957 4 2000 1199.445 1214.414 1197.911 1190.456 Att From: Alexandre Loures Sent: Tuesday, April 08, 2014 8:37 PM To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Somar intervalo de valores nas colunas Olá Rodrigo Coster! É isto mesmo que você solicitou? > dput(head(sp)) structure(c(38795229L, 63924605L, 0L, 0L, 0L, 0L, 33680798L, 51477786L, 0L, 0L, 0L, 0L, 32400472L, 39645654L, 0L, 0L, 0L, 0L, 68022854L, 56066200L, 0L, 0L, 0L, 0L, 31881633L, 67045744L, 0L, 0L, 0L, 0L, 452118806L, 500158705L, 0L, 0L, 0L, 0L, 21795038L, 13562364L, 0L, 0L, 0L, 0L, 2730838L, 4259594L, 0L, 0L, 0L, 0L, 103132447L, 85622631L, 0L, 0L, 0L, 0L, 60353896L, 61121454L, 0L, 0L, 0L, 0L, 32403301L, 40519152L, 0L, 0L, 0L, 0L, 116570336L, 110079499L, 0L, 0L, 0L, 0L, 166832925L, 149309546L, 0L, 0L, 0L, 0L, 78060254L, 101650482L, 0L, 0L, 0L, 0L, 25985390L, 36699369L, 0L, 0L, 0L, 0L, 69222831L, 80306324L, 0L, 0L, 0L, 0L, 4392027L, 4567085L, 0L, 0L, 0L, 0L, 33752865L, 41537988L, 0L, 0L, 0L, 0L, 80216802L, 88982746L, 0L, 0L, 0L, 0L, 1489612L, 1465175L, 0L, 0L, 0L, 0L, 60973115L, 54078115L, 0L, 0L, 0L, 0L, 28441802L, 50033001L, 0L, 0L, 0L, 0L, 60420736L, 17987780L, 0L, 0L, 0L, 0L, 28824716L, 48507914L, 0L, 0L, 0L, 0L, 444712L, 328451L, 0L, 0L, 0L, 0L, 786654L, 838397L, 0L, 0L, 0L, 0L, 12933000L, 10051098L, 0L, 0L, 0L, 0L, 39827569L, 45256447L, 0L, 0L, 0L, 0L, 58803226L, 76290276L, 0L, 0L, 0L, 0L, 46468076L, 38537192L, 0L, 0L, 0L, 0L, 53314093L, 33200142L, 0L, 0L, 0L, 0L, 34642085L, 23556561L, 0L, 0L, 0L, 0L, 119213241L, 121005345L, 0L, 0L, 0L, 0L, 274982262L, 244713051L, 0L, 0L, 0L, 0L, 74376824L, 35888744L, 0L, 0L, 0L, 0L, 2144991L, 2440378L, 0L, 0L, 0L, 0L, 114384520L, 74615518L, 0L, 0L, 0L, 0L, 54130383L, 62022477L, 0L, 0L, 0L, 0L, 3340357L, 4785296L, 0L, 0L, 0L, 0L, 3396350L, 3538952L, 0L, 0L, 0L, 0L, 7356480L, 6493192L, 0L, 0L, 0L, 0L, 30067652L, 34811828L, 0L, 0L, 0L, 0L, 10236714L, 14744233L, 0L, 0L, 0L, 0L, 21064904L, 28239710L, 0L, 0L, 0L, 0L, 5987416L, 9771117L, 0L, 0L, 0L, 0L, 16779486L, 11925052L, 0L, 0L, 0L, 0L, 9060866L, 13985017L, 0L, 0L, 0L, 0L, 15448503L, 24449350L, 0L, 0L, 0L, 0L, 8407209L, 20023742L, 0L, 0L, 0L, 0L, 38288341L, 43777071L, 0L, 0L, 0L, 0L, 610282312L, 508090709L, 0L, 0L, 0L, 0L), .Dim = c(6L, 51L), .Dimnames = list( NULL, c("sp1", "sp2", "sp3", "sp4", "sp5", "sp6", "sp7", "sp8", "sp9", "sp10", "sp11", "sp12", "sp13", "sp14", "sp15", "sp16", "sp17", "sp18", "sp19", "sp20", "sp21", "sp22", "sp23", "sp24", "sp25", "sp26", "sp27", "sp28", "sp29", "sp30", "sp31", "sp32", "sp33", "sp34", "sp35", "sp36", "sp37", "sp38", "sp39", "sp40", "sp41", "sp42", "sp43", "sp44", "sp45", "sp46", "sp47", "sp48", "sp49", "sp50", "sp51"))) Em 08-04-2014 18:37, Rodrigo Coster escreveu: Sem saber como são os teus dados fica complicado a gente de dar um comando que certamente funcionará... da um dput(head(dados)) e nos envia o que aparecer na tela 2014-04-08 18:34 GMT-03:00 Daniel C Bezerra <danielcbezerra(a)gmail.com>: apply(nomedobanco,2,sum) 2014-04-08 18:07 GMT-03:00 Alexandre Loures <alexandre.loures(a)ymail.com>: Olá pessoal! Tenho uma base de dados com 216 observações de 10 variáveis. Sendo que essas observações são mensais, mas tenho que trabalhar com o total anual. Ou seja, eu tenho dados para os meses de 1997, 1998, 1999, 2000, ..., 2014. Mas quero os totais desses anos em cada coluna. Alguém saberia como fazer isso? Desde já muito obrigado! -- Alexandre Rodrigues Loures _______________________________________________ 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. _______________________________________________ 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. _______________________________________________ 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. -- Alexandre Rodrigues Loures -------------------------------------------------------------------------------- _______________________________________________ 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

Trocar caracteres específicos repetidos
by Waldelene Moura 07 Abr '14

07 Abr '14

Olá Sou iniciante no R, preciso retirar de uma string um caractere específico que esteja repetido, trocando estes por apenas um, Ex: a = "23345556333388" gostaria de retirar o 3 especificamente que estão repetidos, logo a string ficaria a = "2345556388" - só altero onde o 3 aparece, os outros continuam como estão. outro exemplo texto = "BOLLA CARRO LLATA CASSINO " aplicando somente para a letra L texto = "BOLA CARRO LATA CASSINO" Alguém poderia me dar uma dica de qual função usar? Waldelene Moura Analista de Sistemas Universidade Federal do Maranhão - PROEN

2 2

table with >= 2^31 elements
by Roney Fraga Souza 07 Abr '14

07 Abr '14

Pessoal, Utilizando o pacote survey resolvi estratificar a amostra de uma base de dados de com as seguintes dimensões: > dim(mydata) [1] 20274412 4 com o seguinte comando: monom00.svy <- svydesign(id=~V0300, strata=~AREAP, weights=~P001p, data=mydata) obtive a seguinte mensagem de erro: Error in table(ids[, 1], strata[, 1]) : attempt to make a table with >= 2^31 elements Com o comando traceback() tenho: 6: stop("attempt to make a table with >= 2^31 elements") 5: table(ids[, 1], strata[, 1]) 4: is.data.frame(x) 3: rowSums(table(ids[, 1], strata[, 1]) > 0) 2: svydesign.default(id = ~V0300, strata = ~AREAP, weights = ~P001p, data = mydata) 1: svydesign(id = ~V0300, strata = ~AREAP, weights = ~P001p, data = mydata) Alguém tem alguma sugestão de como lidar com isso? Informações úteis: > .Machine$integer.max [1] 2147483647 > sessionInfo() R version 3.1.0 beta (2014-03-28 r65330) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=pt_BR.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=en_US.UTF-8 LC_PAPER=pt_BR.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] splines grid tools stats graphics grDevices utils datasets methods base other attached packages: [1] plotrix_3.5-5 ggplot2_0.9.3.1 doBy_4.5-10 MASS_7.3-31 pastecs_1.3-18 [6] boot_1.3-11 psych_1.4.3 weights_0.80 gdata_2.13.3 Hmisc_3.14-3 [11] Formula_1.1-1 survival_2.37-7 lattice_0.20-27 survey_3.29-5 reshape_0.8.4 [16] plyr_1.8.1 quantreg_5.05 SparseM_1.03 ETLUtils_1.2 ffbase_0.11.3 [21] ff_2.2-12 bit_1.1-11 sqldf_0.4-7.1 RSQLite.extfuns_0.0.1 RSQLite_0.11.4 [26] DBI_0.2-7 gsubfn_0.6-5 proto_0.3-10 vimcom.plus_0.9-93 setwidth_1.0-3 [31] colorout_1.0-2 loaded via a namespace (and not attached): [1] chron_2.3-45 cluster_1.15.2 colorspace_1.2-4 dichromat_2.0-0 digest_0.6.4 [6] fastmatch_1.0-4 gtable_0.1.2 gtools_3.3.1 labeling_0.2 latticeExtra_0.6-26 [11] lme4_1.1-5 Matrix_1.1-3 minqa_1.2.3 munsell_0.4.2 nlme_3.1-117 [16] RColorBrewer_1.0-5 Rcpp_0.11.1 RcppEigen_0.3.2.1.1 reshape2_1.2.2 scales_0.2.3 [21] stringr_0.6.2 tcltk_3.1.0 Atenciosamente Roney

1 0

Carta de Controle
by Edson Lira 07 Abr '14

07 Abr '14

Caros amigos, recorro novamente para perguntar-lhes o seguinte: Com a rotina abaixo: data(pistonrings) attach(pistonrings) # 40 sample of 5 obs each qcc.groups(diameter, sample) # some obs are removed, the result is still a 40x5 matrix but with NAs added qcc.groups(diameter[-c(1,2,50,52, 199)], sample[-c(1,2,50,52, 199)]) old <- qcc.options() # save defaults qcc.options("cex.stats") # get a single parameter qcc.options("cex.stats"=1.2) # change parameters qcc.options(bg.margin="azure2") qcc.options("violating.runs" = list(pch = 15, col = "purple")) qcc.options("beyond.limits" = list(pch = 15, col = "orangered")) par(center=1,sd=0.5) qcc(rnorm(100),type = "xbar", center=1,std.dev = 0.5) # see the results qcc.options(old) # restore old default Gero o gráfico de xbarra para os dados do banco pistonrings, já consegui alterar os tamanhos das fontes dos labels de x e de y, não consegui alterar o tamanho da fonte dos limites "LCL" e "UCL", alguém tem alguma sugestão? Muito obrigado. [ ]'s. Edson Lira Estatístico Manaus-Amazonas

1 0