- R-br - DINF UFPR

Classificação de distâncias a partir de um raio em uma matriz de distância
by ASANTOS 07 Jan '15

07 Jan '15

Caros listeiros, Tenho um objeto com coordenadas e um identificador (da) que representa um conjunto de pontos. Calculei a matriz de distância entre todos os pontos e gostaria de a partir de um raio arbitrário, criar uma nova coluna de classificação dos pontos (em da) que estão contidos nesta distância, ou seja, quais pontos estão próximos uns dos outros em uma distância escolhida, segue CRM: #Pacotes require(spatstat) require(sp) ## Gera um processo de pontos xypois=rpoispp(25,win=owin(c(0,1),c(0,1)))## ~25 pontos aleatórios xd=xypois$x yd=xypois$y cod<-seq(1,length(yd)) da<-cbind(xd, yd, cod) ## Crio um vetor com coordenada x e y e um identificador ##Transformando em um objeto da classe espacial pontos<- SpatialPoints(cbind(da[,1],da[,2])) pontos<- SpatialPoints(list(da[,1],da[,2])) pontos<- SpatialPoints(data.frame(da[,1],da[,2])) ##Cria a matriz de distâncias res<-spDists(pontos) res Aqui eu calculei a matriz de distância, mas agora preciso de alguma função que à partir de um raio, como R=0.1, faça a classificação que eu preciso, alguém teria alguma sugestão ou conhecimento de alguma função que faça isto? 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 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

2 2

Homogeneidade de variancias
by Andre Oliveira 06 Jan '15

06 Jan '15

Pessoal boa tarde, dado meu modelo em fatorial. modelo=aov(respsota~factor1*factor2) Como avaliar homogeneidade de variâncias? O que fiz aqui foi ..par(mfrow=c(2,2)) boxplot(resp~fat1) boxplot(resp~fat2) boxplot(resp~fat1*fat2) interaction.plot(fat1,fat2,resp,ylab="Médias",main="Interação Fat1 vs Fat2") Mas está muito subjetivo. 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

6 7

Suavização pela mediana (median polish)
by Felinto COSTA 06 Jan '15

06 Jan '15

Prezados. Por razões didáticas gostaria de decompor a variável regionalizada "v" no exemplo abaixo em duas partes: (i) uma parte atribuída a uma tendência associada às suas coordenadas (x,y) (ii) outra que seriam seus resíduos suavizados pela mediana (median polish) library(geoR) dados.mat=matrix(nrow=8, ncol=3,byrow=F, c(486259,483200,487021,478100,482646,481735,480789,484181, 7428046,7428217,7424717,7420978,7426862,7427070,7427206,7427091, 159,350,220,265,450,110,224,530)) colnames(dados.mat)=c("x","y","v") dados.mp=medpolish(dados.mat) dados.res=as.vector(dados.mp$residuals[,3]) dados.mat.mp[,3]=dados.res dados.mat.mp geo.dados.mat.mp=as.geodata(dados.mat.mp,coords.col = 1:2,data.col=3) Minhas dúvidas são: (i) a "extração" dos resíduos está correta? (ii) numa análise geoestatística, ao se modelar os resíduos (semivariograma experimental > ajuste dos parâmetros de uma função teórica) e usar esta função para estimar os valores numa grade de pontos (krigagem ordinária) ainda deve ser "adicionada" a parte atribuída à tendência removida dos dados experimentais no começo? (ou teoricamente não há razão para isso?) Antecipadamente agradeço a colaboração e desejo a todos um ótimo 2015. Felinto

1 0

Correlação
by Gustavo Faria 06 Jan '15

06 Jan '15

Boa tarde, estou com dificuldade (iniciantes ¬¬ hehe) em fazer uma correlação com dados qualitativos, meus dados são 3 metodologias de identificação de organismos, e eu gostaria de saber o quão parecido e/ou divergente eles foram nas identificações realizadas, como poderia montar essa matriz? ja que tenho somente as identificações realizadas pelas 3 metodologias.... -- --- Gustavo de Faria Silva Biólogo herpetólogo Crbio 87421/04-P Pós-Graduado em Consultoria e Licenciamento Ambiental Mestrando em Zoologia de Vertebrados Telefone: (31) 9579-8473 (Vivo)

3 3

Tratamento de colunas
by Angélica Jaconi 06 Jan '15

06 Jan '15

Prezados, Eu tenho uma planilha ocm numeros. Essa tem o cabecalho (head) indicando as repeticoes, por exemplo: A1, A2, A3...A10, média. Gostaria de uma ideia de como tratar essas colunas iniciadas em "A" como uma única coluna para aplicar posteriores tratamento de dados. Desde ja obrigada pela atencao AJ

2 1

Letras gregas dentro do caption - xtable
by Andre Oliveira 04 Jan '15

04 Jan '15

Pessoas bom dia,estou tentando colocar, em grego, o lambda dentro do caption da tabela gerada pelo xtable, mas sem sucesso. Já olhei no ?plotmath tentei $\lambda$ ~lambda e não consegui. Alguma sugestão? <<keep.source=TRUE,results=tex>>=x=seq(0:12)p=dpois(0:12, 3)dados=data.frame(x,p)require(xtable)x=xtable(dados,caption="Distribuição Poisson ~lambda = 3",align="lll")print(x,caption.placement='top',table.placement="H",size="scriptsize",include.rownames=T, floating=TRUE)@ 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

ANOVA medidas repetidas no R
by Roberto Guevara 03 Jan '15

03 Jan '15

Caros usuários de R, minha dúvida na verdade é mais sobre estatística do que o R em si. Estou utilizando uma ANOVA de medidas repetidas com três fatores sendo que há repetição somente para dois destes fatores. O código e o resultado que estou obtendo segue o modelo abaixo: > summary(aov(dependente ~ fator1 * fator2 * fator3 + + Error(individuo/(fator2 * fator3)))) Error: individuo Df Sum Sq Mean Sq F value Pr(>F) fator1 1 80616 80616 0.260 0.645 fator2 1 809787 809787 2.611 0.205 Residuals 3 930281 310094 Error: individuo:fator2 Df Sum Sq Mean Sq F value Pr(>F) fator2 2 181649 90824 1.491 0.282 fator1:fator2 2 186811 93405 1.533 0.273 Residuals 8 487309 60914 Error: individuo:fator3 Df Sum Sq Mean Sq F value Pr(>F) fator3 2 474691 237345 14.526 0.00502 ** fator1:fator3 2 138023 69012 4.224 0.07163 . fator2:fator3 2 184239 92120 5.638 0.04189 * Residuals 6 98035 16339 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Error: individuo:fator2:fator3 Df Sum Sq Mean Sq F value Pr(>F) fator2:fator3 4 58378 14594 0.374 0.824 fator1:fator2:fator3 4 198072 49518 1.270 0.323 Residuals 16 624020 39001 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Residuals 51 2676429 52479 Todas as cinco tabelas são imprescindíveis para fins de avaliar o efeito dos fatores sobre a variável dependente? Como diferencio o que cada uma expõe? Por exemplo, o resultado da primeira tabela indica o efeito dos fatores 1 e 2 isoladamente ao retirar o efeito dos indivíduos? Qual a diferença entre os resultados expostos para o fator2 na primeira e na segunda tabela? Desde já, muito obrigado pela atenção. -- Roberto Guevara Ferreira Lima Biólogo - CRBio 6ª Região: 73146/06-D Mestre em Zoologia Laboratório de Ecologia e Zoologia de Vertebrados Instituto de Ciências Biológicas Universidade Federal do Pará Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110 Currículo Lattes: http://lattes.cnpq.br/7363382159007600

1 1

Problema com Dates no Macbook pro
by Edson Lira 03 Jan '15

03 Jan '15

Estou com problemas na conversão de datas. Quero separar dia, mês e ano. Abaixo o link de uma amostra dos dados com as variáveis de datas Download TESTE DATA.csv | | | | | | | | | | | Download TESTE DATA.csvDownload TESTE DATA.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Link do Script Download CB.R | | | | | | | | | | | Download CB.RDownload CB.R free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Feliz 2015 para os amigos!!! Edson Lira Estatístico Manaus-Amazonas

2 2

Espaço entre texto e imagens - Sweave
by Andre Oliveira 30 Dez '14

30 Dez '14

Pessoas, boa noite. Depois de muita pesquisas eu recorro a vocês. O que posso colocar no preâmbulo de meu .rnw para reduzir os espaços entre imagens e texto? obrigado 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 7

Selecionando níveis dentro de fator
by ASANTOS 29 Dez '14

29 Dez '14

Caros listeiros, Tenho um tratamento constituído de 26 letras e gostaria de criar um vetor arbitrário para separar alguns níveis dentro deste tratamento, porém não quero utilizar a função & para cada nível (pois no meu banco de dados original ficaria muito longo), mas especificar diretamente um lista com eles sendo: samp.niv<-c("a","k","o","p","r","s","t")##Níveis escolhidos dentro de trat Porém, segundo o meu CRM se faço uma relação de igualdade dentro do banco de dados, no meu caso da, não funciona, alguém poderia me ajudar, segue CRM: ##Criando um banco artificial rm(list=ls(all=TRUE)) da <- expand.grid(trat=rep(paste(letters[1:26]),100), tempo=1:10) da$y1 <- rnorm(da$trat) samp.niv<-c("a","k","o","p","r","s","t")##Níveis escolhidos dentro de trat da.1<-da[da[,1]==samp.niv,]## Selecionando os níveis de interesse 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 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

2 2

Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
by Luciane Maria Pilotto 22 Dez '14

22 Dez '14

Não consegui usar o comando svyolr, está dando os erros abaixo. # Utilizando as variáveis criadas q11 e q13 que foram transformadas em fator: > id$q11 <- as.factor(id$q11) > id$q13 <- as.factor(id$q13) > dstrat<-svydesign(ids=~1, probs=NULL, strata = NULL, variable = NULL, fpc=NULL, + data = NULL, nest = FALSE, weights=~peso,pps=FALSE) > dstrat Independent Sampling design (with replacement) svydesign(ids = ~1, probs = NULL, strata = NULL, variable = NULL, fpc = NULL, data = NULL, nest = FALSE, weights = ~peso, pps = FALSE) > svyolr(q13~q11,design=dstrat, method="logistic") Erro em svyolr.survey.design2(q13 ~ q11, design = dstrat, method = "logistic") : response must be a factor # Utilizando as variáveis quest_11 e quest_13 que são fatores no banco: > svyolr(quest_13~quest_11,design=dstrat, method="logistic") Erro em tapply(1:NROW(x), list(factor(strata)), function(index) { : arguments must have same length ========================================================== Eu tb poderia usar os pesos amostrais nas tabelas de contingencia, no entanto, tb não consigo usar os pesos nestas. Utilizei o comando crosstab (crosstab(q11, q13, peso, expected = TRUE,chisq = TRUE, plot = FALSE)). Dicas com outros comandos que utilizam os pesos amostrais são bem vindas!!! __________________________________________________ Luciane Maria Pilotto Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS -------------------------------------------- Em ter, 25/11/14, Leonardo Ferreira Fontenelle <leonardof(a)leonardof.med.br> escreveu: Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais Para: r-br(a)listas.c3sl.ufpr.br Data: Terça-feira, 25 de Novembro de 2014, 17:09 Não serve a http://r-survey.r-forge.r-project.org/survey/html/svyolr.html? Leonardo Ferreira Fontenelle Em Ter 25 nov. 2014, às 15:36, lutipilotto(a)yahoo.com.br escreveu: Olá, estou tentando rodar regressão logística ordinal e não estou conseguindo usar os pesos amostrais, está dando o seguinte erro: > m11 <- polr(q13 ~ q11 , method="logistic", weights = peso, data=id, Hess=TRUE) Aviso: glm.fit: fitted probabilities numerically 0 or 1 occurred Erro em optim(s0, fmin, gmin, method = "BFGS", ...) : valor inicial em vmmin não é finito *A variável q13 possui 5 ctegorias (1,2,3,4,5) e a variável q11 possui 4 categorias (1,2,3,4). Tem como solicitar o qui-quadrado, o número de observações e o pseudo R2 no mesmo comando utilizado para reg. log. ord.? Como utilizar o peso amostral nas tabelas de contingência? Estou usando os seguintes comandos: table(q11,q13) rowPercents(table(q11,q13)) ou CrossTable(q11, q13, digits=2, expected=TRUE, prop.r=TRUE, prop.c=TRUE, prop.t=FALSE,prop.chisq=FALSE, chisq = FALSE, format=c("SPSS"), dnn = c("local atend","satisfação")) Atenciosamente, Luciane M. Pilotto Doutoranda em Saúde Bucal Coletiva/UFRGS Enviado do Email do Windows _______________________________________________ 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. -----Anexo incorporado----- _______________________________________________ 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.

4 15

Problema em biblioteca do R para uso com Rstudio no Ubuntu
by ASANTOS 20 Dez '14

20 Dez '14

Caros listeiros, Estou usando o R version 2.15.0 (2012-03-30) no Ubuntu em função de pacotes pessoais e outros não atualizados para as versões 3.x, no entanto, ao tentar instalar o Rstudio para utilizar o R, esse não funciona com o erro R shared library (/usr/local/lib/R/lib/libR.so) not found. Minha biblioteca do R esta em /usr/local/bin/R e gostaria de saber o que devo fazer para o Rstudio funcionar, 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 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

3 3

simetria e curtose, limites - distribuição normal
by Luiz Roberto Martins Pinto 17 Dez '14

17 Dez '14

Caros, Considerem x=rnorm(1000). As funções skewness(x) # retorna valor simetria. kurtosis(x) # retorna valor curtose. A distribuição normal teórica é simetrica, com skewness(x)=0, e mesocúrtica, com kurtosis(x)=0. Preciso saber os intervalos para skewness(x) e kurtosis(x), de modo que o vetor x seja considerado 'normalmente distribuído'. Luiz Roberto. Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

2 3

parallel no R
by Fernanda De Bastiani 17 Dez '14

17 Dez '14

Caros, alguém pode me indicar algum material (por ex. tutorial) que tenha exemplos de como fazer paralelização de função no R, além do manual do pacote parallel? Desde já muito obrigada,Fernanda

4 3

Mestrado em Modelagem Matemática e Computacional
by Wagner Tassinari 15 Dez '14

15 Dez '14

Gostaria de anunciar, que o Programa em Modelagem Matemática e Computacional (PPGMMC) está com as inscrições abertas para o Processo Seletivo de ingresso para o ano de 2015 . As inscrições são gratuitas e estão abertas do dia 19 de novembro de 2014 a 13 de janeiro de 2015. Conforme consta no edital, o endereço de inscrição é: Secretaria do PROPPG - UFRRJ – Campus Seropédica, sala 115, Prédio P1 BR 465 Km 7, Seropédica, RJ CEP.: 23897-000 Tel./Fax: (21) 26814707 As inscrições podem ser efetuadas pessoalmente (das 9:00 às 11:30 horas e das 13:30 às 16:00 horas) ou através de postagem. Quaisquer dúvidas ou informações os candidatos podem entrar em contato através dos emails: ppgmmc.ufrrj(a)gmail.com ppgmmc(a)ufrrj.br Site: Home Telefones: (21) 2681-4707 (21) 2681-4600 | | | | | | | | | HomeNotícias do MMC-UFRRJ | | | | Visualizar em r1.ufrrj.br | Visualizado por Yahoo | | | | |

1 0

Cálculo do deviance das variáveis explicativas em modelo de Poisson inflacionado de zeros
by ASANTOS 15 Dez '14

15 Dez '14

Bom dia Pessoal, Estou querendo saber como realizar o cálculo do deviance das variáveis explicativas em um modelo de Poisson inflacionado de zeros com a função zeroinfl() do pacote pscl, pois para um modelo do tipo glm eu faria uma anova(modelo, test="Chi"), mas um objeto criado pela função zeroinfl() não permite anova() e o teste de Vuong serve apenas para comparação entre modelos, existe alguma opção? Segue CRM: #------------------------------------------------------------------ # Definições da sessão. rm(list=ls()) require(pscl) #------------------------------------------------------------------ # Dados artificiais. da <- expand.grid(trat=gl(4,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) betas <- c(0.1,0.9,0.6,0.3,0.7) eta <- X%*%betas y1 <- rpois(da$trat, lambda=exp(eta)) y2 <- rbinom(y1, size=1, prob=0.7) da$y <- y1*y2 str(da) #------------------------------------------------------------------ # Ajuste do modelo. m0 <- zeroinfl(y~trat+tempo|trat, data=da) summary(m0) #------------------------------------------------------------------ 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 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

2 2

Numeração de figuras Sweave na classe artigo
by Andre Oliveira 14 Dez '14

14 Dez '14

Pessoal,bom dia. Minhas figuras estão saindo assim. Figura 1Figura 2 Gostaria que fosse impressa assim. Figura 1.1Figura 1.2 Estou trabalhando assim. \SweaveOpts{width=15,height=12}\setkeys{Gin}{width=0.8\textwidth}\begin{figure}[H]\centering<<fig=T,keep.source=TRUE>>= par(mfrow=c(1,2))curve(dnorm(x,mean=1,sd=sqrt(2)),lwd=2,from=-6,to=17,ylab="Densidade")curve(dnorm(x,mean=10,sd=sqrt(2)),col=2,lwd=2,add=T)legend('topright',legend=c(expression(mu==1),expression(mu==10)),text.col=c(1,2),cex=1.5,bty="n")title("Comparando a distribuição normal com médias diferentes")curve(dnorm(x,mean=4,sd=sqrt(3)),lwd=2,from=-15,to=25,ylab="Densidade")curve(dnorm(x,mean=4,sd=sqrt(20)),col=2,lwd=2,add=T)legend('topright',legend=c(expression(sigma^2==3),expression(sigma^2==20)),text.col=c(1,2),cex=1.5,bty="n")title("Comparando a distribuição normal com variâncias diferentes")@\caption{Curva normal}\label{curvas}\end{figure} Tem como acertar isto? obrigado 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

1 0

Contastes não lineares
by Leonardo Ferreira Fontenelle 12 Dez '14

12 Dez '14

Boa tarde a todos! Procurei na web, mas ainda não achei a resposta para esta dúvida. Tenho um objeto *d* criado pela função *svydesign()*, do pacote *survey*. Duas variáveis desse objeto são *x* (uma variável contínua de 0 a 1) e *y* (uma variável dicotômica, 0 ou 1). Depois de ajustar o seguinte modelo: m <- svyglm(formula = y ~ x, design = subset(d, complete.cases(x, y)), family = quasibinomial) Eu gostaria de estimar a proporção (não as chances) de *y* quando *x* é igual a 0; o mesmo para quando *x* é igual a 1; e então calcular a diferença ou a razão entre essas proporções. Pelo que li, *svycontrast()* é a função mais flexível para estimar contrastes, inclusive não lineares. No entanto, ainda não consegui descobrir qual é a sintaxe para eu convencer o *svycontrast()* a estimar o que estou procurando. Alguém poderia me ajudar com isso? Obrigado! Leonardo Ferreira Fontenelle[1] Links: 1. http://lattes.cnpq.br/9234772336296638

1 1

Re: [R-br] Barplot múltiplo com dados de proporção de variáveis categóricas
by Ivan Bezerra Allaman 11 Dez '14

11 Dez '14

Há relação entre as variáveis? Se houver você pode fazer: tab <- ftable(d) mosaicplot(tab) O resto da confecção você tem capacidade de fazer! (s,f,p) Allaman

1 0

Barplot múltiplo com dados de proporção de variáveis categóricas
by ASANTOS 10 Dez '14

10 Dez '14

Bom dia Pessoal, Tenho 3 variáveis categóricas a, b e c, sendo a com três fatores, b com dois e c com quatro. Eu consigo representar a proporção delas em um barplot individualmente, porém gostaria de colocar as três variáveis em um mesmo barplot, isso é possível? Segue CRM: d <- data.frame(a=factor(sample(c('doente','sadio','indeterminado'), 100, rep=TRUE)), b=factor(sample(c('sim','nao'), 100, rep=TRUE)), c=factor(sample(c('ruim','regular','bom','otimo'), 100, rep=TRUE))) a.tb <- prop.table(table(d$a)) barplot(a.tb) b.tb <- prop.table(table(d$b)) barplot(b.tb) c.tb <- prop.table(table(d$c)) barplot(c.tb) 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 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

1 0

Teste
by Fabiana Goncalves Barbosa 07 Dez '14

07 Dez '14

Tenho os seguintes dados indep inter 288 1887 187 2195 377 111 191 154 165 379 350 156 189 153 100 250 173 186 176 122 431 108 249 1841 213 426 314 252 122 184 635 138 168 169 151 166 154 117 109 141 100 313 1030 125 350 116 329 726 253 420 218 245 123 111 108 135 313 124 275 118 192 100 197 246 222 247 103 567 557 558 299 211 102 184 116 180 119 609 259 155 151 634 343 472 299 454 242 381 167 249 162 248 145 124 126 169 125 151 157 156 208 122 161 106 139 104 101 434 NA 190 NA 171 NA 122 NA 117 NA 2016 NA 349 NA 309 NA 271 NA 255 NA 259 NA 239 NA 158 NA 199 NA 200 NA 168 NA 1368 NA 622 NA 530 NA 373 NA 271 NA 121 NA 188 NA 123 NA 194 NA 125 NA 1811 NA 403 NA 322 NA 229 NA 260 NA 215 NA 151 NA 129 NA 135 NA 104 NA 543 NA 460 NA 286 NA 265 NA 248 NA 164 NA 130 NA 125 NA 110 NA 192 NA 122 NA 103 NA 361 NA 314 NA 262 NA 193 NA 155 NA 513 NA 301 NA 243 NA 177 NA 155 NA 188 NA 193 NA 122 NA 108 NA 109 NA 191 NA 324 NA 163 Que teste uso para verificar se existe diferença entre as duas amostras?

3 4

extracao da Longitude correta para dados de SST
by Humberto 06 Dez '14

06 Dez '14

Boa noite pessoal, Estou extraindo de um arquivo ncdf os valores de latitude (lat) longitude( lon) e temperatura da superfície da agua do mar( SST), entretanto a longitude (lon) ela está na projeção 0 a 360 graus. Minha pergunta é como faço para converter ela para -180 a 180 graus? Segue o script que estou utilizando juntamente com o link dos dados. Desde já agradeço a ajuda Humberto lon lat data sst Lon como deveria ser 299.5 10.5 11/30/1981 27.85 -60.000 300.5 10.5 11/30/1981 27.96 -59.500 301.5 10.5 11/30/1981 28.02 -58.500 302.5 10.5 11/30/1981 28.02 -57.500 303.5 10.5 11/30/1981 27.97 -56.500 304.5 10.5 11/30/1981 27.84 -55.500 305.5 10.5 11/30/1981 27.81 -54.500 306.5 10.5 11/30/1981 27.78 -53.500 307.5 10.5 11/30/1981 27.75 -52.500 Link para acesso dos dados https://www.dropbox.com/s/xnab2g365mduykx/SST_81_14.nc?dl=0 Script: require(ncdf) #'ncdf' versao 1.6.6 require(reshape) require(chron) ###sst nc <- open.ncdf("SST_81_14.nc") nc lat <- get.var.ncdf(nc, "lat") lon <- get.var.ncdf(nc, "lon") tunits <- att.get.ncdf(nc,"time","units") # print tunits tunits$value # criar tempo "real" Times time<- get.var.ncdf(nc,"time") Time<-chron(time,origin=c(month=1, day=1, year=1900)) ## array da variavel sst wm <- get.var.ncdf(nc, "sst") str(wm) ## atribuindo dimnames conforme lon, lat e tempo dimnames(wm)[[1]] <- as.character(lon) dimnames(wm)[[2]] <- as.character(lat) dimnames(wm)[[3]] <- as.character(Time) str(wm) ## melt para reestruturação dos dados mwm <- melt(wm) head(mwm)

6 6

RCurl: passando parâmetros e executando função javascript
by Éder Comunello 04 Dez '14

04 Dez '14

Colegas, boa tarde! Estou precisando de ajuda com relação ao pacote *RCurl*. Tomando por exemplo o endereço < http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001>, ao acessá-lo via *browser *pode-se notar dois botões que executam trechos de código em *javascript *e geram links para o *download *de arquivos de dados: * Botão [Arquivo Access] remete ao comando "javascript:criarArq(2254001,2)" * Botão [Arquivo Texto] remete ao comando "javascript:criarArq(2254001,1)" Estou tentando usar o *RCurl *para passar diretamente os parâmetros para criação do *link *para o "Arquivo Texto". Isto feito, posteriormente pretendo usar o *link *gerado para *download* dos dados. Contudo, não estou tendo sucesso e ficarei muito grato por qualquer ajuda. Abaixo segue o código que tentei aplicar: ### <BEGIN> ### library(RCurl) myURL <- 'http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001' ### URL resultante da operação via browser: ### myURL2 <- ' http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001&CriaArq=true&TipoArq=1 ' myCurl <- getCurlHandle(cookiefile = "myCookies.txt", useragent="Mozilla/5.0", verbose = T, autoreferer = T, followlocation = T) myParams <- list(Codigo=2254001, CriaArq='true', TipoArq=1) ### Codigo: ID da estação ### CriaArq: true (executa função?), ### TipoArq: 1 (tipo texto) e 2 (Access) postForm(myURL, .params=myParams, curl=myCurl) getURL(myURL, curl=myCurl) ### Operação falha, retornando página original! ### Em caso de sucesso embutiria no código um link no formato: ### <a href="ARQ/A20131008-164720-14/CHUVAS.ZIP"> ### <END> ### -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

4 3

Anisotropia, verificação gráfica
by Felinto COSTA 04 Dez '14

04 Dez '14

Prezados. Para o traçado da elipse ("rose diagram") tal como sugerido por Isaaks e Srivastava devo usar os alcances com cor=0.05 (practical ranges) para as diversas direções estudadas ou os alcances (phi) dos patamares (sigmasq) calculados por variofit ou likft (GeoR)? Antecipadamente agradeço a atenção dos colegas. FCosta

2 2

Formatação de várias variáveis
by Alessandro Corrêa 03 Dez '14

03 Dez '14

Prezados colegas, Estou com um banco de dados (x) no seguinte formato sexo idade tempo Q1 Q2 Q3 Q4 Q5 Q6 ...Q63 1 44 24.0 4 3 3 3 2 4 0 62 28.0 1 5 2 2 3 1 1 52 6.0 4 5 4 2 5 4 0 51 27.0 4 4 4 3 3 3 1 32 3.5 1 3 2 1 4 3 0 33 5.0 1 5 1 1 1 1 As variáveis de Q1 a Q6 são variáveis na escala de Likert1 = Discordo totalmente2 = Discordo 3 = Neutro 4 = Concordo 5 = Concordo totalmente Gostaria de saber se é possível criar uma rotina ou se há uma função capaz deformatar as variáveis de Q1 a Q63 em fatores com as etiquetas referentes a cada nível. Se eu fosse fazer uma a uma, repetiria 63 vezes o seguinte comando para cada variável: x$Q1<-factor(x$Q1, label=c("Discordo totalmente", "Discordo", "Neutro", "Concordo", "Concordo totalmente"), levels=1:5) Acredito que haja uma maneira mais prática. Obrigado desde já. Alessandro

4 7

Numero de dias acima de um valor
by Jefferson Ferreira-Ferreira 02 Dez '14

02 Dez '14

Olá amigos!! Estou com uma questão que não estou sabendo como resolver. Eu tenho uma série de dados de nivel d'água (3 vezes por dia). O que eu quero saber é em qual data começa a enchente. Como o nível sofre algumas "pequenas enchentes" que de fato não representam o real início da enchente, o que eu preciso é saber em qual data começa o seguinte evento: nível > 2cm por 45 dias consecutivos. Isso porque se o nível permanecer por menos de 45 dias consecutivos acima de 2cm eu não estou considerando isso como o início da enchente, mas como uma subida intermitente da água. Outro complicador é que preciso saber essa informação para cada um dos 10 equipamento (datalogger), discriminado na culona que chamei de logger. Ou seja, a data (em dia juliano) do início da enchente para cada logger. Em anexo envio os dados de apenas um dos meus aparelhos (logger == 1) Alguma ideia?

4 8

salvando codigos carregados
by Hélio Gallo Rocha 02 Dez '14

02 Dez '14

Bom dia a todos da lista. Vendo alguns posts vi um sobre semente, onde o listeiro Alejandro onde salva uma simulação e reinicia a partir de um ponto com .Random.seed Haveria algo parecido em se salvar relatórios e dados gerados a partir de códigos rodados numa seção e abrir a partir daquele ponto em outra seção? Agradeço desde já. -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Re: [R-br] Digest R-br, volume 48, assunto 3
by Fabiana Goncalves Barbosa 02 Dez '14

02 Dez '14

Olá Cesar A variável "total" é uma variável de contagem. Em 2 de dezembro 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. Re: Teste (David Feitosa) > 2. Re: Teste (Cesar Rabak) > 3. Re: Teste (Cesar Rabak) > 4. Re: Teste (Cesar Rabak) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Mon, 1 Dec 2014 20:18:17 -0300 > From: David Feitosa <contato(a)davidfeitosa.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Teste > Message-ID: > < > CAL1rLto7Dn-1HFEmxx70vh56D6n96YBdR-+1oTAyhfxf8Pop9Q(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Diferença em que sentido? > > > > Atenciosamente, > > David Feitosa > > (\_(\ > (=°;°) > (("")("") > > Em 1 de dezembro de 2014 20:12, Fabiana Goncalves Barbosa < > fabibarbos(a)gmail.com> escreveu: > > > Tenho os seguintes dados > > > > indep inter > > 288 1887 > > 187 2195 > > 377 111 > > 191 154 > > 165 379 > > 350 156 > > 189 153 > > 100 250 > > 173 186 > > 176 122 > > 431 108 > > 249 1841 > > 213 426 > > 314 252 > > 122 184 > > 635 138 > > 168 169 > > 151 166 > > 154 117 > > 109 141 > > 100 313 > > 1030 125 > > 350 116 > > 329 726 > > 253 420 > > 218 245 > > 123 111 > > 108 135 > > 313 124 > > 275 118 > > 192 100 > > 197 246 > > 222 247 > > 103 567 > > 557 558 > > 299 211 > > 102 184 > > 116 180 > > 119 609 > > 259 155 > > 151 634 > > 343 472 > > 299 454 > > 242 381 > > 167 249 > > 162 248 > > 145 124 > > 126 169 > > 125 151 > > 157 156 > > 208 122 > > 161 106 > > 139 104 > > 101 434 > > NA 190 > > NA 171 > > NA 122 > > NA 117 > > NA 2016 > > NA 349 > > NA 309 > > NA 271 > > NA 255 > > NA 259 > > NA 239 > > NA 158 > > NA 199 > > NA 200 > > NA 168 > > NA 1368 > > NA 622 > > NA 530 > > NA 373 > > NA 271 > > NA 121 > > NA 188 > > NA 123 > > NA 194 > > NA 125 > > NA 1811 > > NA 403 > > NA 322 > > NA 229 > > NA 260 > > NA 215 > > NA 151 > > NA 129 > > NA 135 > > NA 104 > > NA 543 > > NA 460 > > NA 286 > > NA 265 > > NA 248 > > NA 164 > > NA 130 > > NA 125 > > NA 110 > > NA 192 > > NA 122 > > NA 103 > > NA 361 > > NA 314 > > NA 262 > > NA 193 > > NA 155 > > NA 513 > > NA 301 > > NA 243 > > NA 177 > > NA 155 > > NA 188 > > NA 193 > > NA 122 > > NA 108 > > NA 109 > > NA 191 > > NA 324 > > NA 163 > > > > Que teste uso para verificar se existe diferença entre as duas amostras? > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141201/98b0d5a3/atta… > > > > ------------------------------ > > Message: 2 > Date: Mon, 1 Dec 2014 21:20:40 -0200 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Teste > Message-ID: > <CAKrF98=_+f4DZ1RwUKn8GxHhHgu6nzdfst6cFUwB-= > L8h0osdg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Fabiana, > > O teste estatístico mais apropriado para seus dados depende de uma série de > premissas que se põe a respeito de como eles são produzidos. > > A variável "total" é uma variável de contagem ou é uma variável de outra > natureza? > > Comentando as respostas que você já obteve, afirmo que com os dados [que > você mostrou] sem nenhuma transformação um teste t (ou mesmo o F como > propõe o Fernando) seriam uma "loteria" pois esses testes são para detetar > diferenças nas médias e essa estatística é muito afetada por valores > extremos (cuja distribuição dos seus dados mostra que ela é assimétrica). > > O inst que tenha a amostra com o maior valor "levaria" vantagem pois o > máximo é 7,5 vezes a média enquanto o menor valor é apenas 0,34 vezes ela. > > Recomendo que você faça uma análise exploratória dos seus dados e > identifique face ao processo que os geram qual seria a melhor forma de > comprovar a hipótese experimental que você que verificar com o teste > estatístico. > > meus 0,019999. > > -- > Cesar Rabak > > > 2014-12-01 16:03 GMT-02:00 Fernando Antonio de souza < > nandodesouza(a)gmail.com > >: > > > Teste t ou F. > > Em 01/12/2014 13:14, "Luciano F. Sgarbi" <luciano.f.sgarbi(a)gmail.com> > > escreveu: > > > > Fabiana. > >> Não seria o caso de um teste t? > >> Abraço > >> Em 01/12/2014 11:27, "Fabiana Goncalves Barbosa" <fabibarbos(a)gmail.com> > >> escreveu: > >> > >>> Olá Pessoal > >>> > >>> Tenho os seguintes dados: inst = variavel independente e total = > >>> variavel dependente. Sendo que, 54 valores são para "indep" e 119 para > >>> "inter". > >>> > >>> Quero verificar se há diferença significativa no "total" nos dois > grupos > >>> "indep" e "inter". Que teste usar? > >>> > >>> Segue os dados: > >>> > >>> inst total indep 288 indep 187 indep 377 indep 191 indep 165 > >>> indep 350 indep 189 indep 100 indep 173 indep 176 indep 431 indep > >>> 249 indep 213 indep 314 indep 122 indep 635 indep 168 indep 151 > >>> indep 154 indep 109 indep 100 indep 1030 indep 350 indep 329 > indep > >>> 253 indep 218 indep 123 indep 108 indep 313 indep 275 indep 192 > >>> indep 197 indep 222 indep 103 indep 557 indep 299 indep 102 indep > >>> 116 indep 119 indep 259 indep 151 indep 343 indep 299 indep 242 > >>> indep 167 indep 162 indep 145 indep 126 indep 125 indep 157 indep > >>> 208 indep 161 indep 139 indep 101 inter 1887 inter 2195 inter 111 > >>> inter 154 inter 379 inter 156 inter 153 inter 250 inter 186 inter > >>> 122 inter 108 inter 1841 inter 426 inter 252 inter 184 inter 138 > >>> inter 169 inter 166 inter 117 inter 141 inter 313 inter 125 inter > >>> 116 inter 726 inter 420 inter 245 inter 111 inter 135 inter 124 > >>> inter 118 inter 100 inter 246 inter 247 inter 567 inter 558 inter > >>> 211 inter 184 inter 180 inter 609 inter 155 inter 634 inter 472 > >>> inter 454 inter 381 inter 249 inter 248 inter 124 inter 169 inter > >>> 151 inter 156 inter 122 inter 106 inter 104 inter 434 inter 190 > >>> inter 171 inter 122 inter 117 inter 2016 inter 349 inter 309 > inter > >>> 271 inter 255 inter 259 inter 239 inter 158 inter 199 inter 200 > >>> inter 168 inter 1368 inter 622 inter 530 inter 373 inter 271 > inter > >>> 121 inter 188 inter 123 inter 194 inter 125 inter 1811 inter 403 > >>> inter 322 inter 229 inter 260 inter 215 inter 151 inter 129 inter > >>> 135 inter 104 inter 543 inter 460 inter 286 inter 265 inter 248 > >>> inter 164 inter 130 inter 125 inter 110 inter 192 inter 122 inter > >>> 103 inter 361 inter 314 inter 262 inter 193 inter 155 inter 513 > >>> inter 301 inter 243 inter 177 inter 155 inter 188 inter 193 inter > >>> 122 inter 108 inter 109 inter 191 inter 324 inter 163 > >>> Grata > >>> > >>> > >>> > >>> > >>> _______________________________________________ > >>> 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141201/f257fc3c/atta… > > > > ------------------------------ > > Message: 3 > Date: Mon, 1 Dec 2014 21:28:48 -0200 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Teste > Message-ID: > <CAKrF98nABg0== > h+mHEX14i9xTWecm8kTAdhj+qC+-QaNResJtw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > David, > > Via de regra considera-se mais "seguro" usar testes bicaudais para > verificações entre tratamentos para a resposta a uma variável numérica > (como parece ser o caso aqui) por isso o "sentido" da diferença não deveria > ser uma preocupação do pesquisador (a menos que haja uma razão dada pela > natureza do problema muito grande para descartar um dos lados da cauda). > > []s > > > 2014-12-01 21:18 GMT-02:00 David Feitosa <contato(a)davidfeitosa.com>: > > > Diferença em que sentido? > > > > > > > > Atenciosamente, > > > > David Feitosa > > > > (\_(\ > > (=°;°) > > (("")("") > > > > Em 1 de dezembro de 2014 20:12, Fabiana Goncalves Barbosa < > > fabibarbos(a)gmail.com> escreveu: > > > >> Tenho os seguintes dados > >> > >> indep inter > >> 288 1887 > >> 187 2195 > >> 377 111 > >> 191 154 > >> 165 379 > >> 350 156 > >> 189 153 > >> 100 250 > >> 173 186 > >> 176 122 > >> 431 108 > >> 249 1841 > >> 213 426 > >> 314 252 > >> 122 184 > >> 635 138 > >> 168 169 > >> 151 166 > >> 154 117 > >> 109 141 > >> 100 313 > >> 1030 125 > >> 350 116 > >> 329 726 > >> 253 420 > >> 218 245 > >> 123 111 > >> 108 135 > >> 313 124 > >> 275 118 > >> 192 100 > >> 197 246 > >> 222 247 > >> 103 567 > >> 557 558 > >> 299 211 > >> 102 184 > >> 116 180 > >> 119 609 > >> 259 155 > >> 151 634 > >> 343 472 > >> 299 454 > >> 242 381 > >> 167 249 > >> 162 248 > >> 145 124 > >> 126 169 > >> 125 151 > >> 157 156 > >> 208 122 > >> 161 106 > >> 139 104 > >> 101 434 > >> NA 190 > >> NA 171 > >> NA 122 > >> NA 117 > >> NA 2016 > >> NA 349 > >> NA 309 > >> NA 271 > >> NA 255 > >> NA 259 > >> NA 239 > >> NA 158 > >> NA 199 > >> NA 200 > >> NA 168 > >> NA 1368 > >> NA 622 > >> NA 530 > >> NA 373 > >> NA 271 > >> NA 121 > >> NA 188 > >> NA 123 > >> NA 194 > >> NA 125 > >> NA 1811 > >> NA 403 > >> NA 322 > >> NA 229 > >> NA 260 > >> NA 215 > >> NA 151 > >> NA 129 > >> NA 135 > >> NA 104 > >> NA 543 > >> NA 460 > >> NA 286 > >> NA 265 > >> NA 248 > >> NA 164 > >> NA 130 > >> NA 125 > >> NA 110 > >> NA 192 > >> NA 122 > >> NA 103 > >> NA 361 > >> NA 314 > >> NA 262 > >> NA 193 > >> NA 155 > >> NA 513 > >> NA 301 > >> NA 243 > >> NA 177 > >> NA 155 > >> NA 188 > >> NA 193 > >> NA 122 > >> NA 108 > >> NA 109 > >> NA 191 > >> NA 324 > >> NA 163 > >> > >> Que teste uso para verificar se existe diferença entre as duas amostras? > >> > >> _______________________________________________ > >> 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141201/d73c49f0/atta… > > > > ------------------------------ > > Message: 4 > Date: Mon, 1 Dec 2014 21:29:15 -0200 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Teste > Message-ID: > < > CAKrF98nTx_GZraqEdnUuv+ZS0RL2gA0H9Hmz5yZ8B3FdE18iyg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Fabiana, > > Os dados desse post são os mesmos do post original, apenas a forma de > organizá-los está diferente. A teoria para testá-los é a mesma. > > []s > > > 2014-12-01 21:12 GMT-02:00 Fabiana Goncalves Barbosa <fabibarbos(a)gmail.com > >: > > > Tenho os seguintes dados > > > > indep inter > > 288 1887 > > 187 2195 > > 377 111 > > 191 154 > > 165 379 > > 350 156 > > 189 153 > > 100 250 > > 173 186 > > 176 122 > > 431 108 > > 249 1841 > > 213 426 > > 314 252 > > 122 184 > > 635 138 > > 168 169 > > 151 166 > > 154 117 > > 109 141 > > 100 313 > > 1030 125 > > 350 116 > > 329 726 > > 253 420 > > 218 245 > > 123 111 > > 108 135 > > 313 124 > > 275 118 > > 192 100 > > 197 246 > > 222 247 > > 103 567 > > 557 558 > > 299 211 > > 102 184 > > 116 180 > > 119 609 > > 259 155 > > 151 634 > > 343 472 > > 299 454 > > 242 381 > > 167 249 > > 162 248 > > 145 124 > > 126 169 > > 125 151 > > 157 156 > > 208 122 > > 161 106 > > 139 104 > > 101 434 > > NA 190 > > NA 171 > > NA 122 > > NA 117 > > NA 2016 > > NA 349 > > NA 309 > > NA 271 > > NA 255 > > NA 259 > > NA 239 > > NA 158 > > NA 199 > > NA 200 > > NA 168 > > NA 1368 > > NA 622 > > NA 530 > > NA 373 > > NA 271 > > NA 121 > > NA 188 > > NA 123 > > NA 194 > > NA 125 > > NA 1811 > > NA 403 > > NA 322 > > NA 229 > > NA 260 > > NA 215 > > NA 151 > > NA 129 > > NA 135 > > NA 104 > > NA 543 > > NA 460 > > NA 286 > > NA 265 > > NA 248 > > NA 164 > > NA 130 > > NA 125 > > NA 110 > > NA 192 > > NA 122 > > NA 103 > > NA 361 > > NA 314 > > NA 262 > > NA 193 > > NA 155 > > NA 513 > > NA 301 > > NA 243 > > NA 177 > > NA 155 > > NA 188 > > NA 193 > > NA 122 > > NA 108 > > NA 109 > > NA 191 > > NA 324 > > NA 163 > > > > Que teste uso para verificar se existe diferença entre as duas amostras? > > > > _______________________________________________ > > 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. > > > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141201/8116ae2c/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 48, assunto 3 > **************************************** >

1 0

Teste
by Fabiana Goncalves Barbosa 01 Dez '14

01 Dez '14

Olá Pessoal Tenho os seguintes dados: inst = variavel independente e total = variavel dependente. Sendo que, 54 valores são para "indep" e 119 para "inter". Quero verificar se há diferença significativa no "total" nos dois grupos "indep" e "inter". Que teste usar? Segue os dados: inst total indep 288 indep 187 indep 377 indep 191 indep 165 indep 350 indep 189 indep 100 indep 173 indep 176 indep 431 indep 249 indep 213 indep 314 indep 122 indep 635 indep 168 indep 151 indep 154 indep 109 indep 100 indep 1030 indep 350 indep 329 indep 253 indep 218 indep 123 indep 108 indep 313 indep 275 indep 192 indep 197 indep 222 indep 103 indep 557 indep 299 indep 102 indep 116 indep 119 indep 259 indep 151 indep 343 indep 299 indep 242 indep 167 indep 162 indep 145 indep 126 indep 125 indep 157 indep 208 indep 161 indep 139 indep 101 inter 1887 inter 2195 inter 111 inter 154 inter 379 inter 156 inter 153 inter 250 inter 186 inter 122 inter 108 inter 1841 inter 426 inter 252 inter 184 inter 138 inter 169 inter 166 inter 117 inter 141 inter 313 inter 125 inter 116 inter 726 inter 420 inter 245 inter 111 inter 135 inter 124 inter 118 inter 100 inter 246 inter 247 inter 567 inter 558 inter 211 inter 184 inter 180 inter 609 inter 155 inter 634 inter 472 inter 454 inter 381 inter 249 inter 248 inter 124 inter 169 inter 151 inter 156 inter 122 inter 106 inter 104 inter 434 inter 190 inter 171 inter 122 inter 117 inter 2016 inter 349 inter 309 inter 271 inter 255 inter 259 inter 239 inter 158 inter 199 inter 200 inter 168 inter 1368 inter 622 inter 530 inter 373 inter 271 inter 121 inter 188 inter 123 inter 194 inter 125 inter 1811 inter 403 inter 322 inter 229 inter 260 inter 215 inter 151 inter 129 inter 135 inter 104 inter 543 inter 460 inter 286 inter 265 inter 248 inter 164 inter 130 inter 125 inter 110 inter 192 inter 122 inter 103 inter 361 inter 314 inter 262 inter 193 inter 155 inter 513 inter 301 inter 243 inter 177 inter 155 inter 188 inter 193 inter 122 inter 108 inter 109 inter 191 inter 324 inter 163 Grata

4 3

Raiz Quadrada do Erro Médio Quadrático e Teste de Médias
by Helder Gramacho 30 Nov '14

30 Nov '14

Olá Pessoal, boa noite! A seguir tem as explicações do problema e o código para reproduzir os resultados, procurei detalhar bem para facilitar a vida de quem se dispõe a ajudar. No trabalho disponível no link abaixo: http://www.lsie.unb.br/rbc/index.php/rbc/article/download/433/428 os autores validam o resultado de uma interpolação utilizando pontos de controle e calculando a Raiz Quadrada do Erro Médio Quadrático (REMQ) dado pela seguinte equação: [image: Imagem inline 2] Onde VMLTC é o Valor do Metro Linear de Testada Corrigida. Repeti o procedimento para 4 diferentes superfícies de interpolação e calculei o REMQ para cada uma como mostrado abaixo: #Download do arquivo no dropbox links <- c("https://www.dropbox.com/s/1grlssvxpihtcmt/valid.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!"))} valid<-read.table(file="valid.csv",sep=",",header=T,dec=".") # Função que calcula o REMQ (ve=valor estimado vr=valor real) remq<-function(ve,vr) { r<-ve-vr r2<-r^2 soma<-sum(r2) n<-length(r2) remq<-sqrt((1/(n-1))*soma) result<-c("Raiz do Erro Médio Quadrático"=remq) return(result) } # Calculando o REMQ remq.orig<-remq(valid[,2],valid[,1]) remq.mod.1.1<-remq(valid[,3],valid[,1]) remq.mod.2.1<-remq(valid[,4],valid[,1]) remq.mod.2.1.sar<-remq(valid[,5],valid[,1]) resumo.remq<-rbind(orig=remq.orig,mod.1.1=remq.mod.1.1,mod.2.1=remq.mod.2.1,mod.sar=remq.mod.2.1.sar);resumo.remq Agora gostaria de saber se as diferenças entre os REMQ são significativas estatisticamente e é aqui que está a minha dúvida, gostaria que avaliassem se o procedimento abaixo está correto, utilizo as diferenças entre os pontos de controle e os pontos obtidos nas superfícies para realizar a análise de variância e em seguida o Teste de Tukey. # Preparando para Teste de Médias m1<-valid[,2]-valid[,1] m2<-valid[,3]-valid[,1] m3<-valid[,4]-valid[,1] m4<-valid[,5]-valid[,1] dados.mod<-cbind(m1,m2,m3,m4);dados.mod dados.org <-data.frame(dif.mod=c(dados.mod[,1],dados.mod[,2],dados.mod[,3],dados.mod[,4]),Tipo=factor(c(rep("m1",20),rep("m2",20),rep("m3",20),rep("m4",20))),observ=factor(rep(1:20,4))) # Realizando a análise de variância tapply(dados.org$dif.mod,dados.org$Tipo,mean) ajuste<-aov(dados.org$dif.mod~dados.org$Tipo+dados.org$observ) summary(ajuste) posthoc<-TukeyHSD(x=ajuste,'dados.org$Tipo',conf.level=0.95);posthoc E então? Está correto utilizar as diferenças entre os pontos de controle e os pontos obtidos na superfície para realizar o teste de médias? Observe que os REMQ são bem diferentes, mas no teste médias as diferenças não são significativas, é aí que está a dúvida. Não sei se da forma como realizei o teste de médias está realmente avaliando as diferenças entre os REMQ. Desde já agradeço toda ajuda, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 2

Krigagem lognormal
by Helder Gramacho 27 Nov '14

27 Nov '14

Bom dia, Quando numa variável foi aplicada a transformação logarítmica como é possível realizar a krigagem lognormal no geo-r? e após a obtenção da estimativa como deve ser feita a transformação para a escala original? Alguém possui um exemplo de como isso pode ser realizado? Desde já muito obrigado, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

3 4

Arrumando um banco de dados
by Simone D. Sartorio 26 Nov '14

26 Nov '14

Pessoal, boa noite! Gostaria de uma ajuda, talvez seja bem fácil para aqueles que programam... 1) Tenho um banco de dados ("dados_prod") grande (110mil linhas) e preciso arrumar uma variavel que irei utilizar...a variavel é "dia_ciclo". O problema é que preciso q esta variavel (dia_ciclo), dentro de cada "n_injeção" esteja entre 1 a 14 (Ou seja, da injeção de n. 0 a 17, cada injeção tem 14 dias). O problema seria fácil se todas as minhas 240 vacas tivessem o mesmo número de tempos observados no inicio do experimento...mas não temos...alguma sugestão para resolver isso? 2) Deste mesmo banco de dados, um outro problema, seria retirar todas as injeções "n_injecao= -1". Estive pensando aqui, resolvendo o problema 2) primeiro o 1) fica bem facil... Alguém pode me ajudar nisso?Desde já, muito obrigadaSimone Exemplo: Dados de 1 vaca!> (vaca1731<- subset(dados_prod, dados_prod$n_vaca==1731)); dim(vaca1731) n_vaca n_lactacao trat pri_multiparas n_injecao dia_CICLO dia_exp producao saude1 1731 3 1 2 -1 -7 -21 46.0 healthy2 1731 3 1 2 -1 -6 -20 49.7 healthy3 1731 3 1 2 -1 -5 -19 50.2 healthy4 1731 3 1 2 -1 -4 -18 54.3 healthy5 1731 3 1 2 -1 -3 -17 47.4 heat6 1731 3 1 2 -1 -2 -16 49.3 healthy7 1731 3 1 2 -1 -1 -15 46.7 healthy8 1731 3 1 2 -1 0 -14 45.1 healthy9 1731 3 1 2 0 -13 -13 47.8 healthy10 1731 3 1 2 0 -12 -12 45.5 healthy11 1731 3 1 2 0 -11 -11 45.2 healthy12 1731 3 1 2 0 -10 -10 47.4 enrollment13 1731 3 1 2 0 -9 -9 48.3 enrollment14 1731 3 1 2 0 -8 -8 46.4 enrollment15 1731 3 1 2 0 -7 -7 47.9 enrollment16 1731 3 1 2 0 -6 -6 50.3 enrollment17 1731 3 1 2 0 -5 -5 50.8 enrollment18 1731 3 1 2 0 -4 -4 47.9 enrollment19 1731 3 1 2 0 -3 -3 46.2 heat20 1731 3 1 2 0 -2 -2 50.8 healthy21 1731 3 1 2 0 -1 -1 45.8 healthy22 1731 3 1 2 0 0 0 50.5 healthy23 1731 3 1 2 1 1 1 49.0 healthy24 1731 3 1 2 1 2 2 48.8 healthy25 1731 3 1 2 1 3 3 53.9 healthy26 1731 3 1 2 1 4 4 45.4 healthy27 1731 3 1 2 1 5 5 38.6 healthy28 1731 3 1 2 1 6 6 41.3 healthy29 1731 3 1 2 1 7 7 50.2 healthy30 1731 3 1 2 1 8 8 50.2 healthy31 1731 3 1 2 1 9 9 50.4 healthy32 1731 3 1 2 1 10 10 49.9 healthy33 1731 3 1 2 1 11 11 50.1 healthy34 1731 3 1 2 1 12 12 52.0 healthy35 1731 3 1 2 1 13 13 47.8 healthy36 1731 3 1 2 1 0 14 50.0 healthy37 1731 3 1 2 2 1 15 49.7 healthy38 1731 3 1 2 2 2 16 50.4 healthy39 1731 3 1 2 2 3 17 48.9 healthy40 1731 3 1 2 2 4 18 49.5 healthy41 1731 3 1 2 2 5 19 48.8 healthy42 1731 3 1 2 2 6 20 56.3 healthy43 1731 3 1 2 2 7 21 52.1 healthy44 1731 3 1 2 2 8 22 51.6 healthy45 1731 3 1 2 2 9 23 50.5 healthy46 1731 3 1 2 2 10 24 46.0 healthy47 1731 3 1 2 2 11 25 53.7 healthy48 1731 3 1 2 2 12 26 51.7 healthy49 1731 3 1 2 2 13 27 49.1 healthy50 1731 3 1 2 2 0 28 51.9 healthy51 1731 3 1 2 3 1 29 51.4 healthy52 1731 3 1 2 3 2 30 43.6 healthy...259 1731 3 1 2 17 13 237 19.9 lameness260 1731 3 1 2 17 0 238 18.8 lameness[1] 259 9 ***---***---***---***---***---***---***---***---***---***---***---***---**** Simone Daniela Sartorio* Professora Adjunta I no CCA/UFSCar, Araras/SP.* Doutora e Mestre em Estatística e Experimentação Agronômica - ESALQ/USP.* Licenciada em Matemática - UNESP/Rio Claro.* Contatos: * tel: (19)3543-2614, ramal 2664* cel: (19)98182-0586 - Tim***---***---***---***---***---***---***---***---***---***---***---***---*** Tenha um bom dia! ;)

2 1

(sem assunto)
by Thaís Brenda Martins 26 Nov '14

26 Nov '14

Boa noite prezad@s, Desenvolvi uma interface gráfica no R com a biblioteca Tcl/Tk. Gostaria que o usuário inserisse os elementos de uma amostra populacional na caixa de entrada (tkentry) ou importasse esses dados de um txt por exemplo. A primeira parte já fiz, funciona corretamente. Mas a segunda não consegui. O que tentei fazer foi: - criei um botão pra que o usuário escolha o arquivo. - Na escopo da função desse botão precisaria jogar os valores do arquivo em forma de uma string na caixa de entrada (onde o usuário digita os elementos quando não quer importar). Para quando o usuário clicar no botão de realizar teste não seja impedido de prosseguir porque o campo esta vazio? Alguém que trabalhe com está biblioteca poderia me dar uma ajuda? Att., Thaís Brenda Martins

3 3

Resolução de grafico word - pdf
by Leonardo Monteiro 26 Nov '14

26 Nov '14

Pessoal, não é uma pergunta relacionada ao R, mas eu estou em fase de termino da tese, e quando passo o work para pdf, os graficos e mapas perdem mto a qualidade (Fica com aparencia de borrado). Eu estou gerando ambos em .tiff com 300 dpi. Por favor, se alguem tiver uma sugestão, eu agradeço. Leonardo

6 5

Facebook
by Vinicius Brito Rocha 25 Nov '14

25 Nov '14

Pessoal, ninguém nunca extraiu dados do facebook com o R? Perguntei mas ninguém respondeu abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário * *M.Sc. Engenharia de Produção/PO* *Ph.D. Estudante Sistemas Computacionais*

6 7

Modelo misto - Coeficiente parcial de determinação
by Fernando Souza 25 Nov '14

25 Nov '14

caros amigos, Estou tentando calcular o coeficiente parcial de determinação, para um efeito de modelo misto, conforme recomendação do livro do Kutner et al 2004 (Applied linear statistical models, 5º edição; para edições anteriores a 4º o autor principal era o Neter). Nesse livro ele calcula ess R2parcial com: R2parcial = (SSE(X1)-SSE(X1,X2)) / SSE(X1) ; onde SSE(X1) é a soma de quadrado do erro do modelo reduzido (modelo sem a variável de interesse) e SSE(X1,X2) a soma de quadrado mo modelo completo (contendo todas as variáveis) Esse R2 mede a porcentagem de redução da variação da variável resposta quando a variável de interesse é adicionada ao modelo. Ele propõe para modelo linear generalizado. Desconfio que para modelo misto não seja realizado da mesma forma devido ao modelo aleatório. Na função abaixo utilizei a função resid() para obtenção dos resíduos e realizar as soma de quadrado, no entanto esse resíduo é resultado do modelo de efeito fixo e aleatório. Há uma forma de remover o resíduo do modelo fixo dos modelos por meio do comando : modelo.lme$residuals[,1]. minha dúvida é como eu devo me proceder para obter esta estatística. R2parcial <- function(modelocompleto,modeloreduzido){ SSR.compl<-sum(resid(modelocompleto)^2) SSR.reduz<- sum(resid(modeloreduzido)^2) R2P <- (SSR.reduz-SSR.compl)/SSR.reduz return(paste("R2 parcial:",round(R2P,3))) } #um pequeno exemplo. somente para demonstração, não considerei os ajustes necessários a esses dados. library(nlme) #o banco de dados Gasoline, está presente no pacote nlme modelocompleto <- lme(yield~endpoint+API+vapor+ASTM,random=~1|Sample,data=Gasoline) modeloreduzido <-lme(yield~API+vapor+ASTM,random=~1|Sample,data=Gasoline) # reduçao gerada pela variavel endpoint R2parcial(modelocompleto,modeloreduzido)

1 0

Re: [R-br] Diferença entre métodos
by Ari Clecius 24 Nov '14

24 Nov '14

Boa noite, alguém poderia me ajudar o porque da diferença? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 Em Sábado, 22 de Novembro de 2014 10:06, Ari Clecius <ari072000(a)yahoo.com.br> escreveu: Bom dia, estou rodando um experimento do tipo Box Behnken e usei a função reg e comparei com o obtido pela função rsm, o problema é que o que os dois consideram significativo difere, aonde posso estar errando? Testei também a função AIC para o modelo que ele fornece e dá diferente pela função AIC(). Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #CMR library(DoE.wrapper)Design.1_2<- bbd.design(nfactors= 3 , ncenter= 4 ,randomize= FALSE,seed=7480, factor.names=list( A=c(0,400),B=c(600,2200),C=c(1,10) ) ) res=c(73.450,75.725,63.980,67.670,76.070,77.345,66.005,67.335,84.740,81.595,74.830,70.325,72.460,73.170,73.925,73.800)length(res)#Medidas resumoDesign.1_2<- add.response(Design.1_2,res, replace=FALSE)plancoded<- code.design(Design.1_2)plancoded LinearModel.1 <- aov(res ~(x1+x2+x3)^2+ I(x1^2) + I(x2^2) + I(x3^2), data=plancoded)sreg = step(LinearModel.1,direction="both") #Melhor modelo#Step: AIC=14.07#res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2) # Df Sum of Sq RSS AIC#<none> 16.074 14.074#+ x1:x2 1 0.501 15.574 15.568#+ x2:x3 1 0.462 15.612 15.607#+ x1:x3 1 0.001 16.073 16.073#- x1 1 9.181 25.255 19.303#- I(x2^2) 1 9.310 25.384 19.385#- I(x3^2) 1 36.195 52.269 30.941#- x2 1 79.223 95.297 40.551#- I(x1^2) 1 86.793 102.867 41.774#- x3 1 212.747 228.821 54.566 summary(sreg) # Df Sum Sq Mean Sq F value Pr(>F) #x1 1 9.18 9.18 5.140 0.04959 * #x2 1 79.22 79.22 44.357 9.27e-05 ***#x3 1 212.75 212.75 119.118 1.72e-06 ***#I(x1^2) 1 86.79 86.79 48.596 6.53e-05 ***#I(x2^2) 1 9.31 9.31 5.213 0.04831 * #I(x3^2) 1 36.20 36.20 20.266 0.00149 ** #Residuals 9 16.07 1.79 #---#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Design.1_2.coded<-plancodedrsmModel.2 <- rsm(res ~ FO(x1, x2, x3) + TWI(x1, x2, x3) + PQ(x1, x2, x3), data=Design.1_2.coded)summary(rsmModel.2) # Estimate Std. Error t value Pr(>|t|) #(Intercept) 73.33875 0.79347 92.4275 1.081e-10 ***#x1 1.07125 0.56107 1.9093 0.104810 #x2 -3.14688 0.56107 -5.6087 0.001370 ** #x3 -5.15688 0.56107 -9.1911 9.350e-05 ***#x1:x2 0.35375 0.79347 0.4458 0.671360 #x1:x3 0.01375 0.79347 0.0173 0.986736 #x2:x3 -0.34000 0.79347 -0.4285 0.683246 #x1^2 -4.65813 0.79347 -5.8705 0.001081 ** #x2^2 1.52563 0.79347 1.9227 0.102878 #x3^2 3.00813 0.79347 3.7911 0.009062 ** #---#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Multiple R-squared: 0.9664, Adjusted R-squared: 0.916 #F-statistic: 19.17 on 9 and 6 DF, p-value: 0.0009315 LinearModel.2 <- aov(res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2), data=plancoded)summary(LinearModel.2 )AIC(LinearModel.2)#[1] 61.47997 #valor obtido anteriormente Step: AIC=14.07

1 0

Ler arquivo .xls - RDOBC e gdata
by Andre Oliveira 24 Nov '14

24 Nov '14

Bom dia,sempre usei o gdata no linux e o RDOBC e gdata no win para ler arquivos .xls, no entanto venho tendo problemas para colocar o gdata e o xlsx mesmo no linux, e quando passo algum trabalho para alunos sempre dá problemas com rjava etc etc. Diante tudo isto recorro a lista e pergunto existe alguma forma mais simples de fazer o serviço de ler um .xls em linux e windows? 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

grafico de efeitos espaciais - est. espacial
by Edson Ferreira 22 Nov '14

22 Nov '14

Paz e Bem !! estou com o seguinte codigo, para modelos CAR com dados de contagem (Poisson) e gostaria de plotar no grid criado um grafico em tom de cinzacom os efeitos espaciais e um outro com a média de cada regiãonão estou conseguindo.neste modelo qual a saída q mostra o valor estimado pro rho.?desde ja obrigado ############################ CARBayes Leroux vs CAR ####################################### ##################################################################### rm(list=ls()) require(spatstat) require(lattice) require(MASS) require(coda) require (MCMCpack) require(plotrix) require(CARBayes) ########## simulação x.easting= 1:10 x.northing <- 1:10 Grid <- expand.grid(x.easting, x.northing) n <- nrow(Grid) distance <-array(0, c(n,n)) W <-array(0, c(n,n)) for(i in 1:n) { for(j in 1:n) { temp <- (Grid[i,1] - Grid[j,1])^2 + (Grid[i,2] - Grid[j,2])^2 distance[i,j] <- sqrt(temp) if(temp==1) W[i,j] <- 1 } } x1= rnorm(n) # variavel 1 x2= rnorm(n) # variavel 2 E= rep(50,n) # valor esperado para cada regiao theta= rnorm(n, sd=0.05) # phi <- mvrnorm(n=1, mu=rep(0,n), Sigma=0.4 * exp(-0.1 * distance)) risk <- exp(-0.2 + 0.1 * x1 + 0.1*x2 + theta + phi) fitted <- E * risk Y <- rpois(n=n, lambda=fitted) formula <- Y ~ x1 + x2 + offset(log(E)) model <- poisson.lerouxCAR(formula=formula, W=W, burnin=50, n.sample=1000) model$formula model$samples model$fitted.values model$random.effects model$residuals model$W.summary model$DIC model$p.d model$summary.results model$model model$accept AttEdson

1 0

Diferença entre métodos
by Ari Clecius 22 Nov '14

22 Nov '14

Bom dia, estou rodando um experimento do tipo Box Behnken e usei a função reg e comparei com o obtido pela função rsm, o problema é que o que os dois consideram significativo difere, aonde posso estar errando? Testei também a função AIC para o modelo que ele fornece e dá diferente pela função AIC(). Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)88412345 (085)33669042 #CMR library(DoE.wrapper)Design.1_2<- bbd.design(nfactors= 3 , ncenter= 4 ,randomize= FALSE,seed=7480, factor.names=list( A=c(0,400),B=c(600,2200),C=c(1,10) ) ) res=c(73.450,75.725,63.980,67.670,76.070,77.345,66.005,67.335,84.740,81.595,74.830,70.325,72.460,73.170,73.925,73.800)length(res)#Medidas resumoDesign.1_2<- add.response(Design.1_2,res, replace=FALSE)plancoded<- code.design(Design.1_2)plancoded LinearModel.1 <- aov(res ~(x1+x2+x3)^2+ I(x1^2) + I(x2^2) + I(x3^2), data=plancoded)sreg = step(LinearModel.1,direction="both") #Melhor modelo#Step: AIC=14.07#res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2) # Df Sum of Sq RSS AIC#<none> 16.074 14.074#+ x1:x2 1 0.501 15.574 15.568#+ x2:x3 1 0.462 15.612 15.607#+ x1:x3 1 0.001 16.073 16.073#- x1 1 9.181 25.255 19.303#- I(x2^2) 1 9.310 25.384 19.385#- I(x3^2) 1 36.195 52.269 30.941#- x2 1 79.223 95.297 40.551#- I(x1^2) 1 86.793 102.867 41.774#- x3 1 212.747 228.821 54.566 summary(sreg) # Df Sum Sq Mean Sq F value Pr(>F) #x1 1 9.18 9.18 5.140 0.04959 * #x2 1 79.22 79.22 44.357 9.27e-05 ***#x3 1 212.75 212.75 119.118 1.72e-06 ***#I(x1^2) 1 86.79 86.79 48.596 6.53e-05 ***#I(x2^2) 1 9.31 9.31 5.213 0.04831 * #I(x3^2) 1 36.20 36.20 20.266 0.00149 ** #Residuals 9 16.07 1.79 #---#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Design.1_2.coded<-plancodedrsmModel.2 <- rsm(res ~ FO(x1, x2, x3) + TWI(x1, x2, x3) + PQ(x1, x2, x3), data=Design.1_2.coded)summary(rsmModel.2) # Estimate Std. Error t value Pr(>|t|) #(Intercept) 73.33875 0.79347 92.4275 1.081e-10 ***#x1 1.07125 0.56107 1.9093 0.104810 #x2 -3.14688 0.56107 -5.6087 0.001370 ** #x3 -5.15688 0.56107 -9.1911 9.350e-05 ***#x1:x2 0.35375 0.79347 0.4458 0.671360 #x1:x3 0.01375 0.79347 0.0173 0.986736 #x2:x3 -0.34000 0.79347 -0.4285 0.683246 #x1^2 -4.65813 0.79347 -5.8705 0.001081 ** #x2^2 1.52563 0.79347 1.9227 0.102878 #x3^2 3.00813 0.79347 3.7911 0.009062 ** #---#Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 #Multiple R-squared: 0.9664, Adjusted R-squared: 0.916 #F-statistic: 19.17 on 9 and 6 DF, p-value: 0.0009315 LinearModel.2 <- aov(res ~ x1 + x2 + x3 + I(x1^2) + I(x2^2) + I(x3^2), data=plancoded)summary(LinearModel.2 )AIC(LinearModel.2)#[1] 61.47997 #valor obtido anteriormente Step: AIC=14.07

1 0

Erro na instalação do R-2.15.0 no Ubuntu 14.04
by ASANTOS 21 Nov '14

21 Nov '14

Boa tarde Pessoal, Por motivos de pacotes que os autores não passaram para R3.0, no ubuntu 14.04 utilizo o R-2.15.0, no entanto, precisei formatar o computador e tentando instalar novamente o R-2.15.0, abaixado no com o arquivo R-2.15.0.tar.gz, não funciona, sendo que após a descompactação do arquivo, na pasta criada que fiz os passos: santosa@santosa-K43E:~$ cd Documentos santosa@santosa-K43E:~/Documentos$ cd R-2.15.0 santosa@santosa-K43E:~/Documentos/R-2.15.0$ ./configure santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install E recebi o erro no checking: configure: error: No F77 compiler found e quando tentei fazer o make install :santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install make: *** Sem regra para processar o alvo `install'. Pare. Alguém sabe se existe alguma incompatibilidade entre o R-2.15.0 e o ubuntu 14.04 ou estou fazendo algo errado? Obrigado, segue uma lista do resultado completo com vários no que recebi no processo de checking usando ./configure: santosa@santosa-K43E:~$ cd Documentos santosa@santosa-K43E:~/Documentos$ cd R-2.15.0 santosa@santosa-K43E:~/Documentos/R-2.15.0$ ./configure checking build system type... i686-pc-linux-gnu checking host system type... i686-pc-linux-gnu loading site script './config.site' loading build-specific script './config.site' checking for pwd... /bin/pwd checking whether builddir is srcdir... yes checking for working aclocal... missing checking for working autoconf... missing checking for working automake... missing checking for working autoheader... missing checking for gawk... no checking for mawk... mawk checking whether ln -s works... yes checking for bison... no checking for byacc... no checking for ar... ar checking for a BSD-compatible install... /usr/bin/install -c checking for sed... /bin/sed checking for which... /usr/bin/which checking for less... /usr/bin/less checking for gtar... no checking for gnutar... no checking for tar... /bin/tar checking for tex... no checking for pdftex... no configure: WARNING: you cannot build PDF versions of the R manuals checking for pdflatex... no configure: WARNING: you cannot build PDF versions of vignettes and help pages checking for makeindex... no checking for makeinfo... no configure: WARNING: you cannot build info or HTML versions of the R manuals checking for texi2dvi... no checking for kpsewhich... no checking for unzip... /usr/bin/unzip checking for zip... /usr/bin/zip checking for gzip... /bin/gzip checking for bzip2... /bin/bzip2 checking for firefox... /usr/bin/firefox using default browser ... /usr/bin/firefox checking for acroread... no checking for acroread4... no checking for xdg-open... /usr/bin/xdg-open checking for notangle... false checking for pkg-config... /usr/bin/pkg-config checking for gcc... gcc checking whether the C compiler works... yes checking for C compiler default output file name... a.out checking for suffix of executables... checking whether we are cross compiling... no checking for suffix of object files... o checking whether we are using the GNU C compiler... yes checking whether gcc accepts -g... yes checking for gcc option to accept ISO C89... none needed checking how to run the C preprocessor... gcc -E checking for grep that handles long lines and -e... /bin/grep checking for egrep... /bin/grep -E checking whether gcc needs -traditional... no checking how to run the C preprocessor... gcc -E checking for gfortran... no checking for f95... no checking for fort... no checking for xlf95... no checking for ifort... no checking for ifc... no checking for efc... no checking for pgf95... no checking for lf95... no checking for gfortran... no checking for ftn... no checking for g95... no checking for f90... no checking for xlf90... no checking for pgf90... no checking for pghpf... no checking for epcf90... no checking for g77... no checking for f77... no checking for xlf... no checking for frt... no checking for pgf77... no checking for cf77... no checking for fort77... no checking for fl32... no checking for af77... no checking for fc... no configure: error: No F77 compiler found santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install make: *** Sem regra para processar o alvo `install'. Pare. -- ====================================================================== 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 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

4 6

Re: [R-br] Erro na instalação do R-2.15.0 no Ubuntu 14.04
by Alisson Lucrécio 20 Nov '14

20 Nov '14

Voce baixou o arquivo .deb ou .taz? Segue umas sugestões que acho pertinente. É por isso, que a partir de agora salvo tmb minha library, assim mesmo que tenha o R-3.1.2, se tiver minha biblioteca salva posso carregar os pacotes dela. Vc já tenteu instalar com biocLite? source("http://bioconductor.org/biocLite.R") biocLite("biology") Att. On 20 November 2014 15:22, Alisson Lucrécio < alisson.lucrecio(a)ifgoiano.edu.br> wrote: > É por isso, que a partir de agora salvo tmb minha library, assim mesmo que > tenha o R-3.1.2, se tiver minha biblioteca salva posso carregar os pacotes > dela. > > Vc já tenteu instalar com biocLite? > > > On 20 November 2014 15:05, ASANTOS [via R-br] < > ml-node+s2285057n4663561h27(a)n4.nabble.com> wrote: > >> Boa tarde Pessoal, >> >> Por motivos de pacotes que os autores não passaram para >> R3.0, no ubuntu 14.04 utilizo o R-2.15.0, no entanto, precisei formatar >> o computador e tentando instalar novamente o R-2.15.0, abaixado no com o >> arquivo R-2.15.0.tar.gz, não funciona, sendo que após a descompactação >> do arquivo, na pasta criada que fiz os passos: >> >> santosa@santosa-K43E:~$ cd Documentos >> santosa@santosa-K43E:~/Documentos$ cd R-2.15.0 >> santosa@santosa-K43E:~/Documentos/R-2.15.0$ ./configure >> santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install >> >> >> E recebi o erro no checking: configure: error: No F77 >> compiler found e quando tentei fazer o make install >> :santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install >> make: *** Sem regra para processar o alvo `install'. Pare. >> >> Alguém sabe se existe alguma incompatibilidade entre >> o R-2.15.0 e o ubuntu 14.04 ou estou fazendo algo errado? >> Obrigado, segue uma lista do resultado completo com vários no que recebi >> no processo de checking usando ./configure: >> >> santosa@santosa-K43E:~$ cd Documentos >> santosa@santosa-K43E:~/Documentos$ cd R-2.15.0 >> santosa@santosa-K43E:~/Documentos/R-2.15.0$ ./configure >> checking build system type... i686-pc-linux-gnu >> checking host system type... i686-pc-linux-gnu >> loading site script './config.site' >> loading build-specific script './config.site' >> checking for pwd... /bin/pwd >> checking whether builddir is srcdir... yes >> checking for working aclocal... missing >> checking for working autoconf... missing >> checking for working automake... missing >> checking for working autoheader... missing >> checking for gawk... no >> checking for mawk... mawk >> checking whether ln -s works... yes >> checking for bison... no >> checking for byacc... no >> checking for ar... ar >> checking for a BSD-compatible install... /usr/bin/install -c >> checking for sed... /bin/sed >> checking for which... /usr/bin/which >> checking for less... /usr/bin/less >> checking for gtar... no >> checking for gnutar... no >> checking for tar... /bin/tar >> checking for tex... no >> checking for pdftex... no >> configure: WARNING: you cannot build PDF versions of the R manuals >> checking for pdflatex... no >> configure: WARNING: you cannot build PDF versions of vignettes and help >> pages >> checking for makeindex... no >> checking for makeinfo... no >> configure: WARNING: you cannot build info or HTML versions of the R >> manuals >> checking for texi2dvi... no >> checking for kpsewhich... no >> checking for unzip... /usr/bin/unzip >> checking for zip... /usr/bin/zip >> checking for gzip... /bin/gzip >> checking for bzip2... /bin/bzip2 >> checking for firefox... /usr/bin/firefox >> using default browser ... /usr/bin/firefox >> checking for acroread... no >> checking for acroread4... no >> checking for xdg-open... /usr/bin/xdg-open >> checking for notangle... false >> checking for pkg-config... /usr/bin/pkg-config >> checking for gcc... gcc >> checking whether the C compiler works... yes >> checking for C compiler default output file name... a.out >> checking for suffix of executables... >> checking whether we are cross compiling... no >> checking for suffix of object files... o >> checking whether we are using the GNU C compiler... yes >> checking whether gcc accepts -g... yes >> checking for gcc option to accept ISO C89... none needed >> checking how to run the C preprocessor... gcc -E >> checking for grep that handles long lines and -e... /bin/grep >> checking for egrep... /bin/grep -E >> checking whether gcc needs -traditional... no >> checking how to run the C preprocessor... gcc -E >> checking for gfortran... no >> checking for f95... no >> checking for fort... no >> checking for xlf95... no >> checking for ifort... no >> checking for ifc... no >> checking for efc... no >> checking for pgf95... no >> checking for lf95... no >> checking for gfortran... no >> checking for ftn... no >> checking for g95... no >> checking for f90... no >> checking for xlf90... no >> checking for pgf90... no >> checking for pghpf... no >> checking for epcf90... no >> checking for g77... no >> checking for f77... no >> checking for xlf... no >> checking for frt... no >> checking for pgf77... no >> checking for cf77... no >> checking for fort77... no >> checking for fl32... no >> checking for af77... no >> checking for fc... no >> configure: error: No F77 compiler found >> santosa@santosa-K43E:~/Documentos/R-2.15.0$ make install >> make: *** Sem regra para processar o alvo `install'. Pare. >> >> -- >> ====================================================================== >> 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:[hidden email] >> <http://user/SendEmail.jtp?type=node&node=4663561&i=0> >> [hidden email] >> <http://user/SendEmail.jtp?type=node&node=4663561&i=1> >> Lattes: http://lattes.cnpq.br/1360403201088680 >> ====================================================================== >> >> >> --- >> Este email foi escaneado pelo Avast antivírus. >> http://www.avast.com >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4663561&i=2> >> 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. >> >> >> ------------------------------ >> If you reply to this email, your message will be added to the >> discussion below: >> >> http://r-br.2285057.n4.nabble.com/R-br-Erro-na-instalacao-do-R-2-15-0-no-Ub… >> To unsubscribe from R-br, click here >> <http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >> . >> NAML >> <http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> > > > > -- > Alisson Lucrecio da Costa > -- Alisson Lucrecio da Costa

1 1

Dificuldade na utilização de um loop duplo para a criação de várias amostras
by Thiago Touya 20 Nov '14

20 Nov '14

Boa tarde! Senhores, peço sua ajuda para criar um código que gere um número 'n' de amostras de tamanho 'n'. Atualmente, consigo fazer apenas uma amostra de tamanho 'n', mas gostaria de fazer de um modo prático a geração de outras 'n-1' amostras. No final gostaria de armazenar tudo numa matriz 'nxn'. Segue o código que estou usando para a geração de uma amostra: ## Tamanho da Amostra n=100 ## Parâmetros 1 a <- c(4.5,1.5) b <- c(3.0,9.0) p <- c(0.5,0.5) ## Função Acumulada Inversa fmenos1_kuma = function(x,a,b){ (1-(1-x)^(1/b))^(1/a) } ## Mistura mkuma = function(x){ (a[1]*b[1]*(x^(a[1]-1))*(1-x^a[1])^(b[1]-1))*p[1] + (a[2]*b[2]*(x^(a[2]-1))*(1-x^a[2])^(b[2]-1))*p[2] } ## Simulação da amostra u1 <- runif(n) u2 <- runif(n) x <- numeric(n) for (i in 1:n){ if(u1[i] < p[1]) x[i]= fmenos1_kuma(u2[i],a[1],b[1]) if(u1[i] >= p[1]) x[i]= fmenos1_kuma(u2[i],a[2],b[2]) } O que eu pensei em fazer foi adicionar um outro 'For' no passo da simulação da amostra, mas não consegui alimentar um objeto 'y' com os valores gerados para cada amostra. Obrigado. -- Att, Thiago Morais de Carvalho (61 - 8169-5094)

4 5

Trocar valores de colunas condicionalmente
by Thiago V. dos Santos 20 Nov '14

20 Nov '14

Olá pessoal, Considerem o seguinte data.frame: df <- structure(list(date = structure(c(1251350100.288, 1251351900, 1251353699.712, 1251355500.288, 1251357300, 1251359099.712), class = c("POSIXct", "POSIXt")), mix.ratio.csi = c(442.78316237477, 436.757082063885, 425.742872761246, 395.770804307671, 386.758335309866, 392.115887652156), mix.ratio.licor = c(447.141491945547, 441.319548211994, 430.854166343173, 402.232640566763, 393.683007533694, 398.388336602215), ToKeep = c(FALSE, FALSE, TRUE, TRUE, TRUE, TRUE)), .Names = c("date", "value1", "value2", "ToKeep"), index = structure(integer(0), ToKeep = c(1L, 2L, 8L, 52L, 53L, 54L, 55L, 85L, 86L, 87L, 88L, 89L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 102L, 103L, 105L, 106L, 192L, 193L, 220L, 223L, 225L, 228L, 229L, 260L, 263L, 264L, 265L, 266L, 267L, 305L, 306L, 307L, 308L, 309L, 310L, 311L, 312L, 313L, 314L, 315L, 352L, 353L, 354L, 375L, 376L, 378L, 379L, 380L, 383L, 411L, 412L, 413L, 414L, 415L, 416L, 418L, 419L, 445L, 453L, 463L, 464L, 465L, 466L, 467L, 468L, 497L, 504L, 547L, 548L, 549L, 586L, 589L, 630L, 631L, 632L, 633L, 634L, 635L, 636L, 644L, 645L, 646L, 647L, 648L, 649L, 650L, 651L, 674L, 675L, 676L, 677L, 678L, 682L, 687L, 690L, 691L, 724L, 725L, 726L, 727L, 728L, 729L, 730L, 731L, 732L, 733L, 734L, 735L, 736L, 739L, 740L, 741L, 742L, 768L, 771L, 772L, 773L, 774L, 775L, 776L, 777L, 778L, 779L, 3L, 4L, 5L, 6L, 7L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L, 17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L, 44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 56L, 57L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L, 90L, 91L, 101L, 104L, 107L, 108L, 109L, 110L, 111L, 112L, 113L, 114L, 115L, 116L, 117L, 118L, 119L, 120L, 121L, 122L, 123L, 124L, 125L, 126L, 127L, 128L, 129L, 130L, 131L, 132L, 133L, 134L, 135L, 136L, 137L, 138L, 139L, 140L, 141L, 142L, 143L, 144L, 145L, 146L, 147L, 148L, 149L, 150L, 151L, 152L, 153L, 154L, 155L, 156L, 157L, 158L, 159L, 160L, 161L, 162L, 163L, 164L, 165L, 166L, 167L, 168L, 169L, 170L, 171L, 172L, 173L, 174L, 175L, 176L, 177L, 178L, 179L, 180L, 181L, 182L, 183L, 184L, 185L, 186L, 187L, 188L, 189L, 190L, 191L, 194L, 195L, 196L, 197L, 198L, 199L, 200L, 201L, 202L, 203L, 204L, 205L, 206L, 207L, 208L, 209L, 210L, 211L, 212L, 213L, 214L, 215L, 216L, 217L, 218L, 219L, 221L, 222L, 224L, 226L, 227L, 230L, 231L, 232L, 233L, 234L, 235L, 236L, 237L, 238L, 239L, 240L, 241L, 242L, 243L, 244L, 245L, 246L, 247L, 248L, 249L, 250L, 251L, 252L, 253L, 254L, 255L, 256L, 257L, 258L, 259L, 261L, 262L, 268L, 269L, 270L, 271L, 272L, 273L, 274L, 275L, 276L, 277L, 278L, 279L, 280L, 281L, 282L, 283L, 284L, 285L, 286L, 287L, 288L, 289L, 290L, 291L, 292L, 293L, 294L, 295L, 296L, 297L, 298L, 299L, 300L, 301L, 302L, 303L, 304L, 316L, 317L, 318L, 319L, 320L, 321L, 322L, 323L, 324L, 325L, 326L, 327L, 328L, 329L, 330L, 331L, 332L, 333L, 334L, 335L, 336L, 337L, 338L, 339L, 340L, 341L, 342L, 343L, 344L, 345L, 346L, 347L, 348L, 349L, 350L, 351L, 355L, 356L, 357L, 358L, 359L, 360L, 361L, 362L, 363L, 364L, 365L, 366L, 367L, 368L, 369L, 370L, 371L, 372L, 373L, 374L, 377L, 381L, 382L, 384L, 385L, 386L, 387L, 388L, 389L, 390L, 391L, 392L, 393L, 394L, 395L, 396L, 397L, 398L, 399L, 400L, 401L, 402L, 403L, 404L, 405L, 406L, 407L, 408L, 409L, 410L, 417L, 420L, 421L, 422L, 423L, 424L, 425L, 426L, 427L, 428L, 429L, 430L, 431L, 432L, 433L, 434L, 435L, 436L, 437L, 438L, 439L, 440L, 441L, 442L, 443L, 444L, 446L, 447L, 448L, 449L, 450L, 451L, 452L, 454L, 455L, 456L, 457L, 458L, 459L, 460L, 461L, 462L, 469L, 470L, 471L, 472L, 473L, 474L, 475L, 476L, 477L, 478L, 479L, 480L, 481L, 482L, 483L, 484L, 485L, 486L, 487L, 488L, 489L, 490L, 491L, 492L, 493L, 494L, 495L, 496L, 498L, 499L, 500L, 501L, 502L, 503L, 505L, 506L, 507L, 508L, 509L, 510L, 511L, 512L, 513L, 514L, 515L, 516L, 517L, 518L, 519L, 520L, 521L, 522L, 523L, 524L, 525L, 526L, 527L, 528L, 529L, 530L, 531L, 532L, 533L, 534L, 535L, 536L, 537L, 538L, 539L, 540L, 541L, 542L, 543L, 544L, 545L, 546L, 550L, 551L, 552L, 553L, 554L, 555L, 556L, 557L, 558L, 559L, 560L, 561L, 562L, 563L, 564L, 565L, 566L, 567L, 568L, 569L, 570L, 571L, 572L, 573L, 574L, 575L, 576L, 577L, 578L, 579L, 580L, 581L, 582L, 583L, 584L, 585L, 587L, 588L, 590L, 591L, 592L, 593L, 594L, 595L, 596L, 597L, 598L, 599L, 600L, 601L, 602L, 603L, 604L, 605L, 606L, 607L, 608L, 609L, 610L, 611L, 612L, 613L, 614L, 615L, 616L, 617L, 618L, 619L, 620L, 621L, 622L, 623L, 624L, 625L, 626L, 627L, 628L, 629L, 637L, 638L, 639L, 640L, 641L, 642L, 643L, 652L, 653L, 654L, 655L, 656L, 657L, 658L, 659L, 660L, 661L, 662L, 663L, 664L, 665L, 666L, 667L, 668L, 669L, 670L, 671L, 672L, 673L, 679L, 680L, 681L, 683L, 684L, 685L, 686L, 688L, 689L, 692L, 693L, 694L, 695L, 696L, 697L, 698L, 699L, 700L, 701L, 702L, 703L, 704L, 705L, 706L, 707L, 708L, 709L, 710L, 711L, 712L, 713L, 714L, 715L, 716L, 717L, 718L, 719L, 720L, 721L, 722L, 723L, 737L, 738L, 743L, 744L, 745L, 746L, 747L, 748L, 749L, 750L, 751L, 752L, 753L, 754L, 755L, 756L, 757L, 758L, 759L, 760L, 761L, 762L, 763L, 764L, 765L, 766L, 767L, 769L, 770L, 780L, 781L, 782L, 783L, 784L, 785L, 786L, 787L, 788L, 789L)), row.names = c(NA, 6L), class = "data.frame") Preciso criar um novo data.frame com a seguinte estrutura: 1) manter a coluna de data intocável, 2) se a coluna 'ToKeep' for igual a FALSE, as colunas 'value1' e 'value2' recebem NA, 3) se a coluna 'ToKeep' for igual a TRUE, as colunas 'value1' e 'value2' nao mudam. Estou tentando usar ifelse, mas por enquanto não consegui encontrar o caminho certo para fazer a indexação. Alguém tem alguma sugestão? Saudações, -- Thiago V. dos Santos PhD student Land and Atmospheric Science University of Minnesota http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/ind… Phone: (612) 323 9898

2 1

função, if dentro fo loop
by Daniela Recchia 19 Nov '14

19 Nov '14

Boa tarde pessoal, tenho certeza que meu problema é bem simples mas não estou conseguindo visualizar a solução, criei a função: conv<- function(n){u<-matrix(nrow=n+1,ncol=1) for (i in 0:n) { if(i==0) f<-function(i) {pnorm(i,m,dp)} else f<- function(i) {pnorm(i*1/n,m,dp)} if(i==0) p<- function(i) {f(i)} else p<- function(i) {f(i)-f(i-1)} I<- function (i) {-p(i)*log2(p(i))} u[i,]<- I(i)} y<-u[0:n,]; y} Mas para o caso onde i==0 não está sendo feito corretamente, e preciso idendificar se não tenho problemas na condição f(i)-f(i-1), o que ele retorna no lugar do i=0 é o valor do i=1. Alguém saberia me dar uma dica? Abraços, Daniela -- Daniela Rodrigues Recchia “The standard of success in life is absolutely the amount of joy you feel” Esther & Jerry Hick

2 8

Não consigo estimar por máxima verossimilhança (problemas com logaritmo de argumento negativo)
by Romero Luiz M. Sales Filho 18 Nov '14

18 Nov '14

Caros colegas, estou com um problema que não consigo resolver... Na verdade acredito que já identifiquei o problemas, mas não sei como contornar esse tipo de situação no R. É o seguinte estou maximizando a função de log-verossimilhança de uma dist. de probabilidade, com o intuito de estimar os parâmetros da dist. Entretanto a otimização não está convergindo para o parâmetro conhecido. Essa dist. possui dois parâmetros, onde um deles o "z" assume valores menores que 2, e o outro parâmetro "eta" assume valores maiores que 0. Consegui simular situações para os casos que tenho 1<z<2 e para 0<z<1, entretanto para valores negativos de z, não estou conseguindo convergência para o parâmetro. Estou usando a função optim() do R, selecionando nessa função o método de Nelder-Mead (já usei outros métodos e nenhum convergiu). O problema que ocorre é que para z negativo, durante as iterações ocorrem o cálculo de logaritmos cujos argumentos estão com valores negativos, tenho quase certeza que o problema da não convergência está aí, aí eu gostaria de saber se tem como colocar alguma restrição durante o processo de iteração para que os argumentos das funções logarítmicas nunca sejam valores negativos. A seguir coloco o código que estou usando: ############################################## ##### Fornecendo os parâmetros ############################## n<-50 z<--1.95 eta<-0.1 z.est=vector() eta.est=vector() p=1 while (p<1001) { ###### Gerando valores da varíavel cujos parâmetros são z e eta u<-runif(n) y = eta*((1-(u)^((1-z)/(2-z)))/(1-z)) #### Obtendo as estimativas para z e eta ######################################## vero <- function(par,x){ z = par[1] eta = par[2] saida<-((sum(log(1 - ((1 - z)*x*(1/eta)))))/(1 - z)) + ((n)*log(2 - z)) + (n*log(1/eta)) return(-saida) } saida1<-optim(par=c(-1,4),fn=vero, method="Nelder-Mead",x=y ) z.est[p]<-saida1[1]$par[1] eta.est[p]<-saida1[1]$par[2] p=p+1 } mean(z.est) mean(eta.est) ############################################## Como estou replicando o procedimento mil vezes, irei utilizar a média das 1000 estimativas encontradas para z e para eta... Então fazendo isso tenho a seguinte saída: > mean(z.est) [1] -30177882 > mean(eta.est) [1] 964132.7 Conforme observamos as estimativas estão muito distantes dos parâmetros reais definidos no início do código (z = -1.95 e eta = 0.1). Dessa forma, desconfio que o problema seja decorrente dos vários erros de logaritmos, que se repetem durante a interação e são da seguinte forma: 1: In log(1 - ((1 - z) * x * (1/eta))) : NaNs produzidos Então tem alguma forma de controlar esse problema que ocorre durante a interação??? Tenho que fazer de alguma forma com que os argumentos fiquem sempre positivos, ou seja: (1-((1-z)*x*(1/eta))) tem que ser sempre positivo... Por favor, caso vcs possam me ajudar, eu agradeço imensamente. Att, Romero.

3 3

Não consigo estimar por máxima verossimilhança (problemas com logaritmo de argumento negativo)
by Romero Luiz M. Sales Filho 18 Nov '14

18 Nov '14

Caros, perdoem mas não estou conseguindo responder a mensagem diretamente pelo tópico já existente... Estou simplesmente dando o reply na mensagem que chega na minha conta e a mensagem retorna... Estou fazendo algum procedimento errado??? Bom, continuando com o tópico comentado por Walmes e Vinícius, faço os meus comentários: Olá Walmes e Vinícius, obrigado pelas dicas, pelo menos já estou buscando encontrar alguma forma de resolver, pois confesso, que já não estava vendo mais como sair do dilema. Vinicius, já tinha tentado fazer a otimização da função de verossimilhança (sem o log), entretanto, não sei o que acontece, mas as estimativas ficam todas iguais aos pontos iniciais, ou seja, para as mil replicações que estou fazendo, todos os resultados são iguais aos pontos iniciais... Com certeza, fazendo dessa forma, não deve estar chegando nem na segunda iteração (não sei o motivo). Walmes, a sua dica é boa, alguém já tinha falado em algo do tipo, mas eu não tinha entendido bem... acho que agora com sua explicação ficou mais claro, resumindo, por favor veja se é isso: preciso encontrar funções cujos possíveis resultado sejam pertencentes ao domínio dos meus dois parâmetros de interesse, com a exigência de que o argumento dessas funções sejam sempre valores pertencentes aos reais, certo?? Pensando assim, como o parâmetro z é sempre menor que 2 e o parâmetro eta é sempre um número positivo, então teoricamente posso usar a seguinte transformação: z= 2-exp(x) (argumento x real e resultado da função sempre menor que 2) e eta=exp (k) (argumento k real e resultado da função sempre maior que zero). É isso mesmo, vcs concordam??? Obrigado pela atençã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>

2 1

Re: [R-br] Problemas lendo .dbf no R
by Paulo Pimenta 16 Nov '14

16 Nov '14

Paulo Henrique de A. S. Pimenta. Graduando em Meteorologia (Bacharelado) - IAG/USP. Fone: +5511981318435. Obrigado pela ajuda!! Suas dicas foram de grande valor!! > ------------------------------ > > Message: 2 > Date: Thu, 13 Nov 2014 10:53:21 -0300 > From: Augusto Ribas <ribas.aca(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: Re: [R-br] Problemas lendo .dbf no R > Message-ID: > <CACMkfRxzSNWYuQcfjAExR5DA5Z8XdW1EZE+n19+YO08H+Rv1gg(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Opa, talvez você possa tentar o seguinte: > > > Um exemplo, suponto que vc tenha 20 arquivos: > > Quando voce da arquivos <- list.files(pattern="\\.dbf$"), você gostaria de > ver algo como isso que fiz no vetor1 certo? > > > vetor1<-paste("corte_estrada_aneis",1:20,".dbf",sep="") > > vetor1 > [1] "corte_estrada_aneis1.dbf" "corte_estrada_aneis2.dbf" > "corte_estrada_aneis3.dbf" > [4] "corte_estrada_aneis4.dbf" "corte_estrada_aneis5.dbf" > "corte_estrada_aneis6.dbf" > [7] "corte_estrada_aneis7.dbf" "corte_estrada_aneis8.dbf" > "corte_estrada_aneis9.dbf" > [10] "corte_estrada_aneis10.dbf" "corte_estrada_aneis11.dbf" > "corte_estrada_aneis12.dbf" > [13] "corte_estrada_aneis13.dbf" "corte_estrada_aneis14.dbf" > "corte_estrada_aneis15.dbf" > [16] "corte_estrada_aneis16.dbf" "corte_estrada_aneis17.dbf" > "corte_estrada_aneis18.dbf" > [19] "corte_estrada_aneis19.dbf" "corte_estrada_aneis20.dbf" > > > Mas o que você ve na verdade é o que está no vetor 2, certo, porque ta tudo > na ordem lexicografia,(Essa palavra é bonita né? Nem sei se escrevi certo, > desculpa) > > > vetor2<-sort(vetor1) > > vetor2 > [1] "corte_estrada_aneis10.dbf" "corte_estrada_aneis11.dbf" > "corte_estrada_aneis12.dbf" > [4] "corte_estrada_aneis13.dbf" "corte_estrada_aneis14.dbf" > "corte_estrada_aneis15.dbf" > [7] "corte_estrada_aneis16.dbf" "corte_estrada_aneis17.dbf" > "corte_estrada_aneis18.dbf" > [10] "corte_estrada_aneis19.dbf" "corte_estrada_aneis1.dbf" > "corte_estrada_aneis20.dbf" > [13] "corte_estrada_aneis2.dbf" "corte_estrada_aneis3.dbf" > "corte_estrada_aneis4.dbf" > [16] "corte_estrada_aneis5.dbf" "corte_estrada_aneis6.dbf" > "corte_estrada_aneis7.dbf" > [19] "corte_estrada_aneis8.dbf" "corte_estrada_aneis9.dbf" > > Mas beleza, você ta vendo o vetor 2 na saida do arquivos <- > list.files(pattern="\\.dbf$"), quando da arquivos. > Uma possibilidade seria reordenar o seu arquivos, antes de fazer a leitura, > da seguinte formula. > > Crie um vetor de numeros, usando uma expressão regular para pegar somente > os números, das string do arquivo, depois transforme em numero pro R > entender: > > > > numeros<-as.numeric(gsub("[a-z_.]", "\\1", vetor2)) > > numeros > [1] 10 11 12 13 14 15 16 17 18 19 1 20 2 3 4 5 6 7 8 9 > > vetor2 > [1] "corte_estrada_aneis10.dbf" "corte_estrada_aneis11.dbf" > "corte_estrada_aneis12.dbf" > [4] "corte_estrada_aneis13.dbf" "corte_estrada_aneis14.dbf" > "corte_estrada_aneis15.dbf" > [7] "corte_estrada_aneis16.dbf" "corte_estrada_aneis17.dbf" > "corte_estrada_aneis18.dbf" > [10] "corte_estrada_aneis19.dbf" "corte_estrada_aneis1.dbf" > "corte_estrada_aneis20.dbf" > [13] "corte_estrada_aneis2.dbf" "corte_estrada_aneis3.dbf" > "corte_estrada_aneis4.dbf" > [16] "corte_estrada_aneis5.dbf" "corte_estrada_aneis6.dbf" > "corte_estrada_aneis7.dbf" > [19] "corte_estrada_aneis8.dbf" "corte_estrada_aneis9.dbf" > > > Veja que a numeração ta na ordem dos arquivos, agora ordene pelos numeros > usando order, que ele vai te dar a ordem que você quer pros indices. > > order(numeros) > [1] 11 13 14 15 16 17 18 19 20 1 2 3 4 5 6 7 8 9 10 12 > > numeros[order(numeros)] > [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 > > > > Agora você aplica a ordem dos numeros no seu arquivos, que você pegou o ls > do diretorio. > > vetor2[order(numeros)] > [1] "corte_estrada_aneis1.dbf" "corte_estrada_aneis2.dbf" > "corte_estrada_aneis3.dbf" > [4] "corte_estrada_aneis4.dbf" "corte_estrada_aneis5.dbf" > "corte_estrada_aneis6.dbf" > [7] "corte_estrada_aneis7.dbf" "corte_estrada_aneis8.dbf" > "corte_estrada_aneis9.dbf" > [10] "corte_estrada_aneis10.dbf" "corte_estrada_aneis11.dbf" > "corte_estrada_aneis12.dbf" > [13] "corte_estrada_aneis13.dbf" "corte_estrada_aneis14.dbf" > "corte_estrada_aneis15.dbf" > [16] "corte_estrada_aneis16.dbf" "corte_estrada_aneis17.dbf" > "corte_estrada_aneis18.dbf" > [19] "corte_estrada_aneis19.dbf" "corte_estrada_aneis20.dbf" > > > > Ai você le tudo usando essa ordem > > for (nome_arquivo in arquivos){ > > g[i] <- read.dbf(paste("I:/qgis_1/exercicio_buffer_L1/testes/teste_1/", > nome_arquivo,sep="")) > > } > > Eu acho que isso deve funcionar. > Bem espero ter mais ajudado que atrabalhado :) > > > > Em 12 de novembro de 2014 03:08, Paulo Pimenta <paulopimenta6(a)hotmail.com> > escreveu: > > > Prezados, > > > > Estou com um problema ao ler arquivos .dbf no R. Abaixo segue a lógica que > > usei em meu código para fazer a leitura: > > > > > > ########################################################################################## > > > > > > > > ########################################################################################### > > ##comando para ler um arquivo ".dbf" > > > > ## > > ## > > > > ## > > ##library("foreign") > > > > ## > > ## > > > > ## > > ##read.dbf("I:/qgis_1/exercicio_buffer_L1/testes/teste_1/corte_estrada_aneis.dbf") > > ## > > ## > > > > ## > > > > ########################################################################################### > > ## > > ##Salvando tudo em forma de lista > > > > library("foreign") > > > > setwd("I:/qgis_1/exercicio_buffer_L1/testes/teste_1") > > > > arquivos <- list.files(pattern="\\.dbf$") > > > > x <- lapply(arquivos, read.dbf) > > > > ##Exportando os comprimentos de cada estrada cortada por um anel > > > > for (i in 1:91){ > > > > y[[i]] <- x[[i]] > > vector_length[i] <- sum(y[[i]][,2]) > > > > } > > > > > > ########################################################################################## > > > > Assim como mostrado no exemplo acima estou lendo um arquivo chamado "corte_estrada_aneis.dbf", > > todavia possui um banco com 91 desses arquivos numerados de forma que seus > > nomes ficam, por exemplo: corte_estrada_aneis.1dbf, > > corte_estrada_aneis2.dbf, corte_estrada_aneis3.dbf, ... até finalmente corte_estrada_aneis91.dbf. > > O problema na leitura dos .dbf consiste que a lógica acima usada por mim > > faz com que a ordenação saia errada, ou seja, ele faz a leitura dos > > arquivos corte_estrada_aneis1.dbf, corte_estrada_aneis12.dbf, corte_estrada_aneis13.dbf > > e assim por diante. Uma tentativa minha, porém não muito bem sucedida foi a > > seguinte: > > > > > > ######################################################################################### > > > > library("foreign") > > > > for (i in 1:9){ > > > > g[i] <- > > read.dbf("I:/qgis_1/exercicio_buffer_L1/testes/teste_1/corte_estrada_aneis[i].dbf") > > > > } > > > > g > > > > > > ######################################################################################### > > > > Mas como tudo esta contido em " " o índice [i] será lido literalmente. > > Peço ajuda aos colegas para poder solucinar este problema, seja com a minha > > lógica ou com alguma idéia diferente. > > > > Abraços a todos! > > > > Paulo Henrique de A. S. Pimenta. > > > > Graduando em Meteorologia (Bacharelado) - IAG/USP. > > Fone: +5511981318435. > > > > _______________________________________________ > > 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. > > > > > > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://recologia.com.br/ <http://augustoribas.heliohost.org> > Github: https://github.com/Squiercg > Lattes: http://lattes.cnpq.br/7355685961127056 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141113/eaa31741/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 47, assunto 14 > *****************************************

2 2

Permutação ou aleatorização
by Luciano F. Sgarbi 14 Nov '14

14 Nov '14

Bom dia Pessoal. Qual a diferença entre permutação e aleatorização? Se é que essas coisas são realmente diferentes... ou permutação é um tipo de aleatorização Em uma busca rápida no google encontrei (de modo bem simplificado) que permutação seria um rearranjo (ou recombinação) dos dados observados. No exemplo abaixo eu fiz uma aleatorização ou permutação? x<-1:50 y<-(.2*x)+rnorm(length(x)) plot(x,y,pch='*',cex=2);abline(lm(y~x),lwd=2) Ay<-sample(y) plot(x,Ay,pch='*',cex=2);abline(lm(Ay~x),lwd=2) Obrigado pela atenção! Luciano -- *Luciano F. Sgarbi* * Currículo Lattes <http://lattes.cnpq.br/9788198135883354>* *Laboratório de Ecologia Teórica e Síntese**Departamento de Ecologia, Instituto de Ciências Biológicas V.* *Universidade Federal de Goiás, campus II, Goiânia-GO, Brasil*

2 2

Semente
by Emerson Cotta Bodevan 13 Nov '14

13 Nov '14

Prezados, bom dia. Quero gerar uma sequência de números aleatórios. Mas quero guardar a semente para poder replicar a mesma sequência em outras simulações. Como faço para guardar esta semente e como faço para usá-la. Na verdade estarei gerando uma sequência de triplas ordenadas. Onde cada número da tripla representa uma localização. Agradeço qualquer ajuda. Att., *Emerson*

3 4