- R-br - DINF UFPR

Regex em uma string de busca
by Fernando Gama 16 Abr '15

16 Abr '15

Bom dia amigos, estou tentando criar um padrão com as strings de busca do meu dataset. Essa é minha string: "municipio == Limeira" "municipio == Limeira & diasemanaincident == terca" "municipio == Limeira & diasemanaincident == terca & *trechoklmetros == ]12.200, 245.000]*" "municipio == Manaus" "municipio == Manaus & diasemanaincident == sexta" "municipio == Manaus & diasemanaincident == sexta & turnoincidente == manha" o resultado que pretendo alcançar é: municipio == ' Limeira ' municipio == ' Limeira ' & diasemanaincident == ' terca ' municipio == ' Limeira ' & diasemanaincident == ' terca ' & ** *trechoklmetros > 12.200 & * *trechoklmetros <= **245.000* ) municipio == 'Manaus' municipio == 'Manaus' & diasemanaincident == 'sexta' municipio == 'Manaus' & diasemanaincident == 'sexta' & turnoincidente =='manha' Ou seja, em termos gerais, inserir *aspas simples (ou duplas*) depois do sinal de atribuição == e a segunda condicao é *ajustar o intervalo ]12.200, 245.000]* (em negrito). Tentei gsub, sub, package stringr, porém sem sucesso, nao sou muito familiarizado com regex :/ Quem puder me ajudar... -- Att, Fernando Gama da Mata

2 2

Importar arquivos de imagem para o R
by Pedro Emmanuel Alvarenga Americano do Brasil 15 Abr '15

15 Abr '15

Amigos R, Eu estou fazendo uma analise com um pacote cujas funçoes graficas utilizam o pacote grid. O que eu gostaria de fazer é de salvar dois graficos no mesmo arquivo lado a lado. Eu não tenho familiaridade com grid então estou batendo cabeça. POstei uma mensagem há dois dias e ninguem respondeu então acho que estou num beco sem saida. Me ocorreu que os amigos que trabalham com mapas podem ter a solução. Seria possivel salvar os graficos separadamente, e depois importar ambos os graficos em diferentes janelas do mesmo dispositivo? Se isso for possível, alguem pode me indicar algum exemplo que eu pudesse seguir? Abraço forte, Pedro Brasil

2 1

Re: [R-br] [Regressao logística] - comparacao de coeficientes beta
by Tiago Souza Marçal 15 Abr '15

15 Abr '15

Acho que você deseja fazer a analise de identidade de modelos. O procedimento foi descrito por Regazzi (2004, 2010) para o delineamento inteiramente ao acaso e blocos ao acaso respectivamente. Caso tenha interesse nesta metodologia em R é só entrar em contato. Att. Tiago. Manuella Lech Cantuaria <manuellalech(a)gmail.com> escreveu: Olá a todos, Alguém de vcs conhece algum teste estatístico para comparacao de coeficientes beta obtidos por regressao logistica? Estou analisando diferentes modelos nos quais vario uma das variáveis independentes e minhas amostras sao dependentes. Desde já, muito obrigada :-) Manuella _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

3 2

Análise de Fatores PCA e Dendograma
by Andre Oliveira 15 Abr '15

15 Abr '15

Pessoal boa tarde,estou com dois bancos de dados com a estrutura do USArrests que vem no R só que repetido no tempo. USArrests2007, USArrests2009, USArrests2010. Gostaria de orientação de como anexar estes dados para em um banco de dados único e rodar PCA, Plotar, dendograma e a análise de fatores. Pelo que percebi em meus devo fazer o seguinte. Montar um banco de dados com os três bancos e tirar a soma ou média. Não sei se fere princípios da análise. dados= data.frame(USArrests2007, USArrests2009, USArrests2010) # Não sei se é data.frame aqui mesmo! V1=tapply(V1,grupos,mean) .....Vn=tapply(Vn,grupos,mean) dadosmedio=cbind(V1, ........................,Vn)hclust(dist(dadosmedio)prcomp(dadosmédio) Ou usar a estrutura com Vn=tapply(Vn,grupos,sum) em vez de Vn=tapply(Vn,grupos,mean). Ou usar o próprio banco formado por USArrests2007, USArrests2009, USArrests2010. Caso falei besteira, desconsidere. Ficarei agradecido por sugestões. obrigado pela atenção 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

5 11

Coluna categorizando variável por intervalo
by Fátima Lima Paula 15 Abr '15

15 Abr '15

Prezados, tenho uma coluna com variáveis com uma letra e três números, conforme exemplo abaixo. a<-c("B095","B096","B784","B561","N590","N980","N067","C678","C234","C009","B095", "B561","N980","N980","C234","C009") b<-c("Maria","Paola","José","Pedro","Kleber","Henrique","Marcelo","João","Lia", "Ingrid","Renato","Paulo","Mariana","Mario","Luiza","Luiz") banco<-data.frame(a,b) banco Quero criar uma coluna de forma que as linhas que tiverem de C000 a C300 recebam a palavra"verde", as que tiverem de C301 a C700, recebam "amarelo". As que tiverem de B000 a B550 recebam vermelho e de B 551 a B999 recebam "branco". O restante ficaria NA. Alguém poderia me ajudar? Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

3 4

[Regressao logística] - comparacao de coeficientes beta
by Manuella Lech Cantuaria 14 Abr '15

14 Abr '15

Olá a todos, Alguém de vcs conhece algum teste estatístico para comparacao de coeficientes beta obtidos por regressao logistica? Estou analisando diferentes modelos nos quais vario uma das variáveis independentes e minhas amostras sao dependentes. Desde já, muito obrigada :-) Manuella

1 0

Re: [R-br] Gráfico/CL50
by Tiago Souza Marçal 13 Abr '15

13 Abr '15

Use a função points() Att. Tiago. Flavia Batista Gomes - Embrapa Amazonia Ocidental - CPAA <flavia.b.gomes(a)embrapa.br> escreveu: Olá, Tenho um gráfico dose-resposta e preciso marcar um ponto específico, o ponto que representa a CL50. Como insiro essa informação no gráfico? Obrigada, Flávia ----- Mensagem original ----- De: "walmes ." <walmeszeviani(a)gmail.com> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Segunda-feira, 19 de janeiro de 2015 22:49:32 Assunto: Re: [R-br] Probit/CL50 O pacote MASS tem funções para cálculo de dl50. Como a DL50 é uma função (não linear) dos parâmetros do preditor linear de um modelo glm, sua estimação intervalar geralmente se dá por meio da aplicação do método delta. O pacote car pode ser considerado par isso. Se não estiver enganado, a função chama se deltaMethod(). Consulte a documentação para exemplos de uso. À disposição. Walmes. -- ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ========================================================================== _______________________________________________ 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. ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure. _______________________________________________ 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.

2 1

Residuos quantile randomizados (qresiduals{stamod})
by Fernando Antonio de souza 13 Abr '15

13 Abr '15

Alguém pode me explicar os residuo quantile aleatorizados? Como posso usá-lo para avaliar ajuste do modelo de distibuição de meu modelo através de quantil-quantil plot, uma vez que o mesmo utiliza aleatorizações para gerar residuos de distribuição contínua, ou seja a cada rodada ele dará valores diferentes. Tem sentido utilizá-lo com a função shapiro.test()? -- ======================================================================= Fernando Souza Zootecnista, DSc. Produção Animal celular: (+55) 82 8113-8781 e-mail:nandodesouza@gmail.com https://producaoanimalcomr.wordpress.com/ ========================================================================

2 3

Sobre R Markdow
by Marcio Leles 12 Abr '15

12 Abr '15

Boa tarde. Estou fazendo um script referente a uma prova e queria que somente saísse os resultados e não os comandos. Estou usando o Rmarkdow no Rstudio. Gostaria de saber se alguém sabe o comando para que ficasse uma questão a cada página, ou seja o comando para quebrar página. Como abaixo eu gostaria que a quatro 6 ficasse em uma página e a 8 na outra. O pedaço do script é. **Resolução Para a Questão 6** ```{r echo=FALSE} dic(DIC1$trat, DIC1$htponto, quali = TRUE, mcomp = "tukey", sigT = 0.05, sigF = 0.05) ``` ```{r include=FALSE} DBC1<-read.table("/Users/Marcio/Documents/_R/rdbc_prova1_2015.txt",header=TRUE, dec = ",") ``` **Resolução Para a Questão 8** ```{r echo=FALSE} dbc(DBC1$trat, DBC1$bloco,DBC1$tabponto, quali = TRUE, mcomp = "tukey", sigT = 0.05, sigF = 0.05) ``` -- Atenciosamente; Prof. Marcio Leles R. de Oliveira http://marcioromarco.webnode.com/ Mensuração e Manejo Florestal Universidade Federal dos Vales do Jequitinhonha e Mucuri - UFVJM Departamento de Engenharia Florestal Campus JK - Rodovia MGT367-KM 583, 5000 - Alto da Jacuba Diamantina - MG, CEP: 39100-000 Tel: (38) 9158-1435

5 4

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 11 Abr '15

11 Abr '15

Oi Wagner, isso mesmo... eu havia conseguido apenas resolver a integral mais interna... a outra realmente não consigo... Mesmo assim, muito obrigado pela ajuda! Vou ter que resolver numericamente mesmo, e isso com certeza vai consumir muito tempo computacional, mas é o jeito... Muito Obrigado! Em 10 de abril de 2015 05:37, Wagner Bonat <wbonat(a)gmail.com> escreveu: > Ola Romero, > > Eu tentei resolver a sua integral a interna tem solução analitica que > segue a abaixo, mas a externa de 0 to Infinity me parece que não tem > solução analitica. A não ser que tenha mais algum tipo de restrição nos > parâmetros. > > Agora vc decide o que fica melhor. > Flw > > (eta^(beta/q-1+1)*q*%e^(beta*log(x)-log((q-1)*%e^(beta*log(x))+eta^beta)/q-1))/(beta*eta^beta*q-2*beta*eta^beta)-(eta^(beta/q-1+1)*%e^(beta*log(x)-log((q-1)*%e^(beta*log(x))+eta^beta)/q-1))/(beta*eta^beta*q-2*beta*eta^beta)+(eta^(beta/q-1+1)*%e^(-log((q-1)*%e^(beta*log(x))+eta^beta)/q-1))/(beta*q-2*beta)-eta^(beta/q-1+1)/(beta*(eta^beta)^(1/q-1)*q-2*beta*(eta^beta)^(1/q-1)) > > > Em 8 de abril de 2015 23:25, Romero Luiz M. Sales Filho < > romero.sfilho(a)gmail.com> escreveu: > >> Oi Wagner, desculpa... realmente esqueci de adicionar esse detalhe... >> Segue: >> >> q < 2, >> beta, eta e teta > 0. >> >> Att, >> >> Romero. >> >> _______________________________________________ >> 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. >> > > > > -- > Wagner Hugo Bonat > > ---------------------------------------------------------------------------------------------- > Department of Mathematics and Computer Science (IMADA) > University of Southern Denmark (SDU) and > Laboratório de Estatística e Geoinformação (LEG) > Universidade Federal do Paraná (UFPR) > > -- *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>

1 0

Dados com muitos zeros (sugestão de como analisar)
by Fernando Antonio de souza 11 Abr '15

11 Abr '15

Caros amigos Tenho os seguintes dados abaixo. Desejo avaliar a resposta NEFA em função dos fatores TRATAMENTO E TEMPO. Acontece que a variável NEFA possui muitos zeros, os quais são também respostas. Gostaria de sugestões de como posso analisar tal tipo de dado. att dados<-structure(list(ANIMAL = structure(c(13L, 13L, 13L, 13L, 13L, 13L, 13L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 17L, 17L, 17L, 17L, 17L, 17L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("20", "18", "3", "2", "5", "6", "13", "9", "14", "4", "17", "7", "1", "12", "11", "8", "19", "15", "16", "10"), class = c("ordered", "factor")), TRATAMENTO = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("0.146", "0.154", "0.182", "0.214"), class = "factor"), TEMPO = structure(c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 7L, 1L, 2L, 3L, 4L, 5L, 6L, 1L, 2L, 3L, 4L, 5L, 6L), .Label = c("0", "6", "12", "24", "36", "60", "84" ), class = "factor"), NEFA = c(0.28762, 0, 0, 0, 0.0344, 0, 0, 0.12422, 0, 0, 0, 0.00382, 0, 0, 0.08313, 0, 0, 0, 0.00192, 0, 0, 0.2064, 0, 0, 0, 0.02007, 0, 0, 0.16244, 0, 0, 0, 0, 0, 0, 0.16626, 0, 0, 0, 0.08982, 0, 0, 0.23602, 0.01242, 0, 0, 0, 0, 0, 0.5179, 0, 0, 0, 0, 0, 0, 0.18728, 0, 0, 0, 0.0086, 0, 0, 0.71283, 0, 0, 0, 0.00764, 0, 0, 0.30386, 0, 0, 0, 0.0086, 0, 0, 0.2943, 0, 0, 0, 0.01911, 0, 0, 0.17009, 0.0258, 0, 0, 0, 0, 0, 0.20066, 0, 0, 0, 0, 0, 0, 0.5781, 0, 0, 0, 0.01338, 0, 0, 0.5867, 0, 0, 0, 0.01624, 0, 0, 0.23315, 0, 0, 0, 0, 0, 0, 0.06211, 0, 0, 0, 0, 0, 0, 0.54083, 0, 0, 0, 0.01816, 0, 0.05638, 0.05829, 0, 0, 0.02484, 0)), .Names = c("ANIMAL", "TRATAMENTO", "TEMPO", "NEFA"), class = c("nffGroupedData", "nfGroupedData", "groupedData", "data.frame"), row.names = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 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, 52L, 53L, 54L, 55L, 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, 85L, 86L, 87L, 88L, 89L, 90L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 98L, 99L, 100L, 101L, 102L, 103L, 104L, 105L, 106L, 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, 134L, 135L, 136L, 137L, 138L, 139L), formula = NEFA ~ TEMPO | ANIMAL, FUN = function (x) max(x, na.rm = TRUE), order.groups = TRUE, na.action = structure(c(133L, 140L), .Names = c("133", "140"), class = "omit")) -- ======================================================================= Fernando Souza Zootecnista, DSc. Produção Animal celular: (+55) 82 8113-8781 e-mail:nandodesouza@gmail.com https://producaoanimalcomr.wordpress.com/ ========================================================================

2 1

desassociar arquivo
by Hélio Gallo Rocha 10 Abr '15

10 Abr '15

Boa tarde Listeiros... No pacote EBImage tem o comando display. Na primeira execução do comando, não havia nenhum programa associado ao "display" e selecionei um visualizador de imagem, mas na realidade deveria associar um navegador, pois o display geral um html. Como faço para associar o display a um navegador? Desde já agradeço -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Ler arquivo .zip
by Luciano F. Sgarbi 10 Abr '15

10 Abr '15

Bom dia pessoal. Gostaria de saber se algum de vocês tem a solução para o seguinte problema: Estou baixando uns arquivos raster (pelo R) e gostaria de salva-los em uma pasta para posteriormente usá-los. O problema é que para um formato de arquivo e site eu estou conseguindo (Exemplo 1 no código abaixo), entretanto, para o outro conjunto eu não consigo ler os arquivos de jeito nenhum. Qualquer ajuda sera super bem vinda.. ### ##Inicio: Exemplo 1 require(raster) url1<-"ftp://srtm.csi.cgiar.org/SRTM_v41/SRTM_Data_ArcASCII/srtm_03_08.zip" file1<-'srtm_03_08'##Tamanho do arquivo: 395k utils::download.file(url1,destfile=paste(file1,".zip",sep=""),mode='b',method='wget') r1<-raster::raster(unzip(paste(file1,".zip",sep=""),paste(file1,".asc",sep=""))) r1 ##Fim: Exemplo 1 ##Aqui está o problema ##Tentativa 1 url2<-" ftp://ftp.glcf.umd.edu/glcf/SRTM/media/GLCF.TSM.B3-010.00.Finished_B-01sec-… " file2<-'SRTM_fB01_n000w177'##Tamanho do arquivo: 66k utils::download.file(url2,destfile=paste(file2,".zip",sep=""),mode='b',method='wget') r2<-raster::raster(unzip(paste(file2,".zip",sep=""),paste(file2,".tif",sep=""))) zip.file.extract(paste(file2,".zip",sep=""),paste(file2,".tif",sep="")) ##Tentativa 2 temp<-tempfile() utils::download.file(url2,destfile=temp,mode='b',method='wget') raster::raster(temp) raster(paste(temp,'/',file2,".tif",sep="")) ###Informações do OS sysname release "Linux" "3.2.0-39-generic-pae" version nodename "#62-Ubuntu SMP Wed Feb 27 22:25:11 UTC 2013" "aa" machine login "i686" "unknown" user effective_user "aa" "aa" > sessionInfo() R version 3.1.3 (2015-03-09) Platform: i686-pc-linux-gnu (32-bit) Running under: Ubuntu precise (12.04.5 LTS) locale: [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8 LC_PAPER=pt_BR.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] RCurl_1.95-4.1 bitops_1.0-6 rgdal_0.9-2 raster_2.2-31 sp_1.0-15 downloader_0.3 loaded via a namespace (and not attached): [1] digest_0.6.4 grid_3.1.3 lattice_0.20-30 tools_3.1.3 -- *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*

3 4

Probit/CL50
by Flavia Batista Gomes - Embrapa Amazonia Ocidental - CPAA 10 Abr '15

10 Abr '15

Olá pessoal, Preciso de orientação para fazer a análise de Probit, estimar a CL50 e fazer o gráfico. Meu experimento avaliou o efeito inseticida de doses de uma substância. Se puderem...será muito válido! Flávia Batista ____________________________________________________________________________ Aviso de confidencialidade Esta mensagem da Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), empresa publica federal regida pelo disposto na Lei Federal no. 5.851, de 7 de dezembro de 1972, e enviada exclusivamente a seu destinatario e pode conter informacoes confidenciais, protegidas por sigilo profissional. Sua utilizacao desautorizada e ilegal e sujeita o infrator as penas da lei. Se voce a recebeu indevidamente, queira, por gentileza, reenvia-la ao emitente, esclarecendo o equivoco. Confidentiality note This message from Empresa Brasileira de Pesquisa Agropecuaria (Embrapa), a government company established under Brazilian law (5.851/72), is directed exclusively to its addressee and may contain confidential data, protected under professional secrecy rules. Its unauthorized use is illegal and may subject the transgressor to the law's penalties. If you are not the addressee, please send it back, elucidating the failure.

5 11

Medidas tomadas em diferentes tempos
by Fernando Souza 10 Abr '15

10 Abr '15

Caros amigos, tenho um banco de dados organizados da seguinte forma e gostaria de uma opiniao para confirmar se estou procedendo corretamente > head(dadosmetabolitos[,1:4]) ANIMAL TRATAMENTO TEMPO PROTEINA 1 1 0.146 0 5.70 2 1 0.146 6 3.37 3 1 0.146 12 4.41 4 1 0.146 24 5.68 5 1 0.146 36 5.37 6 1 0.146 60 6.63 Os tratamentos foram aleatoriamente distribuídos aos animais em um DIC. A resposta PROTEINA, foi medida em cada animal em tempos diferentes. Estou em dúvida quanto a forma de avaliar este trabalho. Entendo que as medidas tomadas em diferentes tempos no mesmo animal estão correlacionadas, por isso estou utilizando o seguinte modelo: #install.packages("nlme") library(nlme) modelo2 <- gls(PROTEINA~TRATAMENTO*as.factor(TEMPO),correlation=corAR1(form=~1+TEMPO|Animal),na.action=na.omit,data=dados) onde TRATAMENTO e TEMPO são fatores uma variável de númerica. Minha dúvida é a seguinte. Seria melhor eu tratar a variável tempo como uma variavel do tipo fator e analisar a interação, ou tratar o tempo como uma variável numerica (neste caso uma covariável). O objetivo do trabalho é avaliar se a resposta dos animais submetidos a cada tratamento foi influenciada pelo tempo que os animais permaneceram no tratamento. dados<-structure(list(ANIMAL = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, oi 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 11L, 11L, 11L, 11L, 11L, 11L, 11L, 12L, 12L, 12L, 12L, 12L, 12L, 12L, 13L, 13L, 13L, 13L, 13L, 13L, 13L, 14L, 14L, 14L, 14L, 14L, 14L, 14L, 15L, 15L, 15L, 15L, 15L, 15L, 15L, 16L, 16L, 16L, 16L, 16L, 16L, 16L, 17L, 17L, 17L, 17L, 17L, 17L, 17L, 18L, 18L, 18L, 18L, 18L, 18L, 18L, 19L, 19L, 19L, 19L, 19L, 19L, 19L, 20L, 20L, 20L, 20L, 20L, 20L, 20L), .Label = c("1", "2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", "15", "16", "17", "18", "19", "20"), class = "factor"), TRATAMENTO = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L), .Label = c("0.146", "0.154", "0.182", "0.214"), class = "factor"), TEMPO = c(0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L, 0L, 6L, 12L, 24L, 36L, 60L, 84L), PROTEINA = c(5.7, 3.37, 4.41, 5.68, 5.37, 6.63, 5.39, 7, 4.84, 4.83, 5.81, 5.99, 7.65, 6.16, 5.53, 5.1, 5.99, 5.29, 4.46, 7.31, 6.57, 5.44, 5.45, 3.86, 3.83, 5.15, 6.13, 4.91, 5.96, 5.45, 5.2, 5.54, 4.54, 6.75, 5.84, 5.47, 5.37, 4.67, 3.98, 5.01, 6, 4.05, 6.03, 4.35, 4.32, 4.9, 4.22, 6.27, 5.71, 6.15, 4.4, 4.86, 5.81, 4.91, 6.31, 4.21, 6.74, 4.82, 6.42, 5.69, 5.67, 6.2, 6.13, 5.71, 5.53, 3.75, 4.6, 4.74, 6.26, 5.53, 5.37, 3.54, 5.63, 5.1, 3.4, 5.63, 3.94, 6.62, 5.43, 4.76, 3.69, 3.75, 6.48, 4.71, 6.74, 4.25, 4.11, 4.84, 5.13, 6.74, 4.06, 7.06, 5.69, 4.3, 6.06, 5.47, 6.8, 5.44, 6.63, 5.92, 4.3, 6.31, 5.64, 6.57, 5.03, 5.07, 4.64, 4.64, 4.71, 4.38, 5.52, 4.14, 4.64, 5.75, 4.67, 5.45, 5.1, 6.08, 6.09, 6.32, 3.33, 4.49, 4.69, 5.97, 5.74, 5.25, 5.74, 5.67, 4.37, 5.03, 4.12, 6.02, NA, 6.67, 6.32, 4.64, 4.75, 5.83, 6.06, NA)), .Names = c("ANIMAL", "TRATAMENTO", "TEMPO", "PROTEINA"), class = "data.frame", row.names = c(NA, -140L)) -- Fernando Antônio de Souza Zootecnista,Dsc. Nutrição e alimentação animal (UFMG) Lattes: http://lattes.cnpq.br/6519538815038307 e-mail:nandodesouza@gmail.com (82)9618-6457 (tim) (82)8113-8781 (vivo) (82)3313-8781 (net) (31)9784-5273 (vivo-MG)

4 5

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 08 Abr '15

08 Abr '15

Oi Wagner, desculpa... realmente esqueci de adicionar esse detalhe... Segue: q < 2, beta, eta e teta > 0. Att, Romero.

1 0

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 08 Abr '15

08 Abr '15

Wagner, um caso particular que tbm é do meu interesse é quando o Beta (na integral acima) é igual a 1... Mesmo assim achei analiticamente complexo de se resolver... Mas caso vc tenha uma luz já ajudará bastante! Abraç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

Fwd: Plotar duas escalas no eixo x
by Eder Comunello 08 Abr '15

08 Abr '15

Segue uma ideia inicial, que você pode ajustar às suas necessidades. ### <code r> a <- c(3.4, 4.3) b <- c(.53, .56) aRange <- c(0, 6) bRange <- c(0, .6) a1 <- t(scale(a-aRange[1], center=F, scale=diff(aRange))) b1 <- t(scale(b-bRange[1], center=F, scale=diff(bRange))) plot1 <- barplot(-a1, col=8, bor="light gray", horiz=T, beside=F, xlim=c(-1,1), axes=F, space=NULL, names.arg=c("1985-1991","1991-1999"), cex.names=.75, ylab="Time Interval") plot2 <- barplot(b1, col=8, bor="light gray", horiz=T, beside=F, xlim=c(-1,1), axes=F, space=NULL, add=T) axis(1, cex.axis=.75, seq(-1,0,len=2), rev(aRange)) axis(3, cex.axis=.75, seq(0,1,len=2), bRange) abline(v=0, lty=1) abline(v=.54/.6, lwd=1, col="red", lty=2) text(x=-a1+.1, y=plot1, labels=round(a,2), cex=.6, col=1, xpd=TRUE) text(x=b1-.1, y=plot2, labels=round(b,2), cex=.6, col=1, xpd=TRUE) title(xlab="Interval change Area (percent of map)", cex.lab=.75, adj=0) title(xlab="Anual change Area (percent of map)", cex.lab=.75, adj=1) ### </code> [image: Imagem inline 1] Em 7 de abril de 2015 17:28, Rodrigo Muniz <uenf.rodrigo(a)gmail.com> escreveu: > Olá Eder, tentei uma estratégia semelhante a sua, um pouco menos elegante, > porém a maior dificuldade é que o zero teria de ficar no meio do gráfico e, > para direita teria uma escala que varia de 0 a 0.6, para esquerda de 0 a 5. > No seu exemplo a origem é na esquerda, desta forma todas as barras partem > do mesmo ponto, o que não é desejado. > Segue um exemplo do gráfico que eu gostaria de programar. > > Não consegui programar o gráfico como eu gostaria, então como uma solução > bastante grosseira e imediata eu "desenhei" o gráfico, como se segue: > Porém continua a dificuldade em faze-lo de uma outra forma. > > #******Criando o rótulo do eixo x > eixox=c("5","4","3","2","1","0", "0.1", "0.2","0.3","0.4","0.5","0.6") > #********Criando a sequensia compativel com o rótulo de x > eixoy=seq(1:12) > #*********Criando o dataframe > eixos=data.frame(eixox, eixoy) > > #*********Configurando parametros gráficos, o mgp é o parâmetro mais > importante, pois define > #********as ditâncias dos ticks em relação ao gráfico > par(mar = c(0.5, 10, 0, 0.3), mgp = c(0.5, -0.25, 0), > bg="gray94", bty="n", cex.axis=0.9, mar = c(4,8, 4, 4), ps=15) > > #***********Criando o gráfico de barras em branco com os rótulos de x > desejado, > #a escala de x é a relação entre número de rótulos dividido por 10 > barplot(with(eixos,(tapply(eixoy,eixox,mean))),axes=FALSE, > names.arg=as.factor(c("5","4","3","2","1","0", > "0.1", "0.2","0.3","0.4","0.5","0.6")), col="white", > border="white", > ylab="1985-1991 1991-1999", > xlab="Interval change Area (percent of map) Anual > change Area (percent of map)") > #*****Colocando as barras verticais sobre os rótulos de x, a sequência é > determinada pelo > #número de rótulos dividido por 10 > abline(v=c(0.7,1.9,3.1,4.3,5.5,6.7,7.9,9.1,10.3,11.5,12.7,13.9), lwd=0.9, > col="gray70") > > #****Criando as barras horizontais > rect(2.5,1, 13.5,4, col="gray", border=FALSE) > rect(1.3,7, 13,10, col="gray", border=FALSE) > > #****Criando as linhas verticais > abline(v=6.7, lwd=2) > abline(v=13.1, lwd=1, col="red", lty=2) > > #******Criando as setas > arrows(13.2, 10.5,15, 10.5, cex=0.5, col="red", length=0.1) > arrows(13, 10.5,11.6, 10.5, cex=0.5, col="red", length=0.1) > > #Definido o local do texto. Obs: os comando devem ser executados um de > cada vez. > text(locator(1), "Slow Fast") > text(locator(1), "0.54 = Uniform Intensity") > > >

3 2

ajuste Regressão logística com o uso do cluster
by Luciane Maria Pilotto 08 Abr '15

08 Abr '15

Bom dia, gostaria de ajuda para realizar o ajuste de um modelo de regressão logística usando cluster. Como estimar a variância robusta ou o odds proporcional? #Utilizei os seguintes comandos: dstrat<-svydesign(ids=~1, probs=cod_mun, strata = NULL, variable = NULL, fpc=NULL, data = id15, nest = FALSE, weights=NULL,pps=FALSE,method="logistic") m15<-svyolr(q13~q11,design=dstrat, method="logistic") #modelo bruto #m156<-svyolr(q13~q11+q06+q10+q08+q12+q14,design=dstrat, method="logistic") #modelo ajustado (ctable <- coef(summary(m15))) p <- pnorm(abs(ctable[, "t value"]), lower.tail = FALSE) * 2 (ctable <- cbind(ctable, "p value" = p)) (ci <- confint(m15)) # default method gives profiled CIs exp(coef(m15)) exp(cbind(OR = coef(m15), ci)) regTermTest(m15, ~q11, method="Wald") LogRegR2(m15) #erro Error in if (!(model$family$family == "binomial" && (model$family$link == : missing value where TRUE/FALSE needed Atenciosamente, Luciane Pilotto

1 0

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 08 Abr '15

08 Abr '15

Bom dia Wagner, na verdade coloquei essa integral apenas para ilustrar... A minha integral na verdade é bem mais complexa... já tentei resolve-la de forma analítica e não consegui... Em anexo envio a forma ad minha integral... se vc conseguir me ajudar a resolver analiticamente será realmente muito bom! [image: Imagem inline 1] Abraç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>

1 0

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 08 Abr '15

08 Abr '15

Wagner, isso mesmo! Vc conseguiu! Muito obrigado! Estava já perdendo as esperanças! Era isso que eu precisava! Abraç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

Organização de banco de dados
by Andre Oliveira 07 Abr '15

07 Abr '15

Boa tarde pessoas para eu fazer uma análise preciso de manipular os dados baixo. > teste = matrix(c(Altura, Peso, Alt_pai, Alt_mae, Peso_pai, Peso_mae),ncol = 6) > teste [,1] [,2] [,3] [,4] [,5] [,6] [1,] 156 49 165 163 65 59 [2,] 160 50 170 150 70 50 [3,] 165 67 180 170 90 60 [4,] 149 41 169 158 73 59 [5,] 170 61 160 162 73 68 [6,] 165 65 183 162 90 75 [7,] 156 49 166 160 71 64 [8,] 173 55 168 150 86 77 [9,] 162 52 171 156 89 80 [10,] 174 79 170 160 72 64 [11,] 172 75 170 165 70 60 [12,] 175 60 175 161 80 57 [13,] 176 65 176 152 95 48 [14,] 175 69 170 168 81 78 [15,] 174 68 172 168 98 70 [16,] 163 61 171 160 90 62 [17,] 176 83 160 158 80 61 [18,] 181 73 177 160 77 50 [19,] 175 74 165 160 62 60 [20,] 166 60 160 156 63 60 [21,] 170 63 167 160 63 60 [22,] 165 69 176 150 90 55 [23,] 170 71 168 160 85 70 [24,] 171 72 170 152 84 62 [25,] 179 76 170 170 80 88 [26,] 179 60 160 160 65 60 [27,] 174 70 160 165 75 70 [28,] 173 65 170 160 80 60 [29,] 151 48 175 170 80 78 [30,] 160 51 165 152 62 68 [31,] 155 55 158 156 85 65 [32,] 150 50 168 160 80 82 [33,] 158 51 155 153 60 53 [34,] 173 76 175 160 80 67 [35,] 170 65 171 158 72 55 [36,] 165 62 180 165 80 75 [37,] 167 67 185 162 86 76 [38,] 153 43 161 152 75 62 [39,] 177 69 176 153 72 80 [40,] 172 66 169 163 98 61 [41,] 168 62 165 162 75 62 [42,] 179 74 170 170 85 80 [43,] 154 52 164 153 88 62 [44,] 175 75 174 159 80 62 [45,] 183 96 190 165 86 76 [46,] 165 60 166 164 80 73 [47,] 180 60 165 160 62 58 [48,] 174 102 168 162 66 62 [49,] 157 50 164 150 96 54 [50,] 158 58 164 166 53 75 [51,] 180 74 165 160 70 56 [52,] 173 66 170 163 85 64 [53,] 171 76 169 161 73 62 [54,] 188 76 181 165 86 65 [55,] 170 70 175 165 85 68 [56,] 190 72 173 168 74 69 [57,] 156 56 170 168 85 68 [58,] 174 85 170 165 70 104 [59,] 182 74 168 165 86 75 [60,] 182 82 179 160 80 96 [61,] 162 53 179 1 59 80 [62,] 160 53 165 150 60 79 [63,] 175 90 175 165 100 60 [64,] 172 74 170 168 80 70 [65,] 155 51 170 153 80 50 [66,] 163 62 180 155 80 85 [67,] 177 64 191 169 110 75 [68,] 171 87 168 165 80 84 [69,] 180 66 182 157 70 65 [70,] 179 75 163 160 76 65 [71,] 175 104 168 156 64 56 [72,] 181 89 162 161 70 60 [73,] 170 74 190 170 80 95 [74,] 180 64 182 165 74 68 [75,] 173 89 185 165 98 62 [76,] 173 70 168 166 80 75 [77,] 182 83 165 157 69 57 colnames(teste)=c("X1", "X", "X2", "X3", "X4", "X5") Preciso de usar rownames(teste) para inserir os rótulos das linhas automaticamente. Alguém de vocês tem uma ideia de como fazer isto automaticamente? Os rótulos estão abaixo! Como este rótulos baixo estão sem padrão se puderem sugerir algo do tipo que envolva x <- gl(5, 2, labels=LETTERS[1:26]) [1] Ativ. Culturais Praia Balada Praia [5] Praia Praia Outros Praia [9] Praia Balada Ativ. Esportivas Balada [13] Balada Outros Ativ. Esportivas Outros [17] Praia Balada Balada Ativ. Esportivas[21] Balada Praia Ativ. Culturais Ativ. Esportivas[25] Ativ. Esportivas Ativ. Esportivas Praia Ativ. Culturais [29] Ativ. Esportivas Praia Praia Praia [33] Praia Ativ. Culturais Balada Outros [37] Outros Outros Ativ. Esportivas Ativ. Esportivas[41] Ativ. Esportivas Outros Ativ. Esportivas Ativ. Esportivas[45] Ativ. Esportivas Ativ. Esportivas Ativ. Esportivas Outros [49] Praia Ativ. Esportivas Ativ. Esportivas Outros [53] Ativ. Esportivas Ativ. Esportivas Outros Ativ. Esportivas[57] Outros Ativ. Esportivas Outros Ativ. Esportivas[61] Outros Outros Outros Ativ. Esportivas[65] Ativ. Culturais Ativ. Culturais Balada Ativ. Esportivas[69] Ativ. Esportivas Ativ. Esportivas Ativ. Culturais Ativ. Esportivas[73] Balada Ativ. Esportivas Ativ. Esportivas Ativ. Esportivas[77] Ativ. Esportivas 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

2 2

Plotar duas escalas no eixo x
by Rodrigo Muniz 07 Abr '15

07 Abr '15

Caros companheiros da lista, Tenho de fazer alguns gráficos de barra na posição horizontal, onde no eixo y serão as classes e no eixo x os valores, porém o eixo x tem duas escalas diferentes, por exempo: -10 a 0 e 100000 a 200000. Ou seja, o eixo x deve ter duas escalas diferentes. Fiz uma busca pela internet em português e em inglês, mas não obtive sucesso em encontrar um material que me ajude a iniciar a programação. Gostaria de saber se alguém teria alguma sugestão, alguma ideia ou material para eu começar a estudar sobre como fazer esses gráficos? Desde já agradeço a atenção. -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutor em Ciências (ESALQ/USP) E-mail - muniz.ra(a)usp.br Cel (19) 98300-4333 (Pessoal)

3 4

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 07 Abr '15

07 Abr '15

Oi, Mauro... Sim o Wolfram resolve... O problema é que como preciso resolver várias integrais mudando os seus parâmetros fica inviável resolver no wolfram... No caso o 2 e o 3 que aparecem na função são valores particulares de parâmetros da que irão variar.. E preciso resolver isso pelo menos pra umas mil possibilidades... Aí fazer um por um não dá... Tem que usar o R pra rodar tudo. Abraço, Romero.

3 2

Dummy
by Alexandre Loures 07 Abr '15

07 Abr '15

Boa noite! Quero construir dummies da seguinte forma: Tenho uma variável origem (os 27 estados), uma ano (1998 e 1999) e outra destino (os 27 estados mais 51 países), então gostaria de forma as seguintes dummies: 1º) 27 dummies para origem e ano 1998 e 78 dummies para destino e ano 1998; e 2º) 27 dummies para origem e ano 1999 e 78 dummies para destino e ano 1999. > dput(head(mydata)) structure(list(origin = structure(c(1L, 1L, 1L, 1L, 1L, 1L), .Label = c("ac", "al", "am", "ap", "ba", "ce", "df", "es", "go", "ma", "mg", "ms", "mt", "pa", "pb", "pe", "pi", "pr", "rj", "rn", "ro", "rr", "rs", "sc", "se", "sp", "to"), class = "factor"), destination = structure(c(3L, 3L, 5L, 5L, 7L, 7L), .Label = c("ac", "africa", "al", "algeria", "am", "angola", "ap", "arab", "arabia", "argentina", "australia", "austria", "ba", "belgium", "bolivia", "canada", "ce", "chile", "china", "colombia", "df", "ecuador", "egypt", "es", "finland", "france", "germany", "go", "greece", "hong_kong", "india", "indonesia", "iran", "iraq", "ireland", "israel", "italy", "japan", "kingdom", "korea", "ma", "malaysia", "mexico", "mg", "ms", "mt", "netherlands", "nigeria", "norway", "pa", "paraguay", "pb", "pe", "peru", "philippines", "pi", "poland", "portugal", "pr", "rj", "rn", "ro", "rr", "rs", "russian", "sc", "se", "singapore", "sp", "spain", "sweden", "switzerland", "taiwan", "thailand", "to", "turkey", "uruguay", "usa"), class = "factor"), ano = c(1998L, 1999L, 1999L, 1998L, 1999L, 1998L)), .Names = c("origin", "destination", "ano"), row.names = c(NA, 6L), class = "data.frame") Alguém poderia me ajudar? Desde já muito obrigado! -- *Alexandre Rodrigues Loures*

4 4

Re: [R-br] Sorting
by Robert Iquiapaza 07 Abr '15

07 Abr '15

O que vc considera "tamanho" de "m"? já que "m" é um objeto de duas dimensões. Se considerar o produto das dimensões (número de elementos) poderia fazer algo assim: um = list(id=1, m=matrix(1:9, nrow=3)) dois = list(id=2, m=matrix(1:4, nrow=2)) tres = list(id=3, m=matrix(1, nrow=1)) quatro = list(id=4, m=matrix(1:6, nrow=3)) lista = list(um, dois, tres,quatro) (dims=sapply(lista,function(z)prod(dim(z$m)))) (indxs=order(dims)) lapply(indxs,function(z)lista[[z]]) Robert Em 6 de abril de 2015 22:04, David Feitosa <davidfeitosa(a)gmail.com> escreveu: > Caros, boa noite. > > Uma dúvida de ordenação. > Suponha que eu tenha 3 list compostas: > > um = list(id=1, m=matrix(1:9, nrow=3)) > dois = list(id=2, m=matrix(1:4, nrow=2)) > tres = list(id=3, m=matrix(1, nrow=1)) > > e crie um vetor com as 3 list: > > lista = c(um, dois, tres) > > Haveria alguma forma de ordená-las pelo tamanho do elemento m, > ou seja, de forma que a lista final fosse equivalente a: > > c(tres, dois, um) ? > > > Atenciosamente, > > David Feitosa > > (\_(\ > (=°;°) > (("")("") > > _______________________________________________ > 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. >

2 1

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 06 Abr '15

06 Abr '15

Muito obrigado, David, mas infelizmente não sei se essa função resolve o problema, uma vez que o pacote que tal função está contida, não está mais disponível para instalar no R. Consegui baixar uma versão do pacote que estava na internet mas não consegui instalar, acho que pelo fato do arquivo ser bastante antigo (versão de 1999)... Continuo tentando resolver o problema... Se alguém puder ajudar agradeço bastante! Romero.

3 2

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 06 Abr '15

06 Abr '15

Ah entendi agora o que vc quis falar... Mas o meu problema é justamente esse... Existe sim possibilidades de integrais em que o limite de uma variável dependa da outra variável... O problema é que no R, só vejo possibilidade de integrar como vc falou com limites numéricos... mas o problema que tenho pra resolver tem um limite como variável... Mesmo assim valeu pela ajuda!! Vamos ver se mais alguma boa dica aparece! Abraç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>

1 0

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 06 Abr '15

06 Abr '15

Oi David, não entendi o que vc questionou no último e-mail... Mas simplificando apenas escrevi a função considerando as duas variáveis X e Y da seguinte forma: f <- function(x,y) (2*exp(-2*x))*(3*exp(-3*y)) e em seguida, pedi para integrar usando a função do pacote rmutil, ou seja: int2(f, a=c(0,x), b=c(0,Inf)) Vc acha que precisa de mais alguma complementação?? Assumi que quando escrevi a função as variáveis já teriam sido definidas... X e Y... Valeu!

2 1

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 06 Abr '15

06 Abr '15

Oi David, consegui instalar o rmutil... baixei a versão que esá nesse link aqui: http://www.commanster.eu/rcode.html Infelizmente não funcionou... Ainda não é essa função que resolve o meu problema... Pior é que esse tipo de função é tão comum em outros programas... mas no R não encontrei nada que me ajudasse... Preciso resolver no R, pois vou colocar uma rotina pra ele resolver uma grande quantidade de integrais... Olha aí o resultado usando a função rmutil : > f <- function(x,y) (2*exp(-2*x))*(3*exp(-3*y)) > int2(f, a=c(0,x), b=c(0,Inf)) Erro em is.vector(a, mode = "numeric") : objeto 'x' não encontrado Continuo tentando aqui... qualquer coisa me avisa! Obrigado, 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

Sorting
by David Feitosa 06 Abr '15

06 Abr '15

Caros, boa noite. Uma dúvida de ordenação. Suponha que eu tenha 3 list compostas: um = list(id=1, m=matrix(1:9, nrow=3)) dois = list(id=2, m=matrix(1:4, nrow=2)) tres = list(id=3, m=matrix(1, nrow=1)) e crie um vetor com as 3 list: lista = c(um, dois, tres) Haveria alguma forma de ordená-las pelo tamanho do elemento m, ou seja, de forma que a lista final fosse equivalente a: c(tres, dois, um) ? Atenciosamente, David Feitosa (\_(\ (=°;°) (("")("")

1 0

Re: [R-br] Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 06 Abr '15

06 Abr '15

Oi David, cheguei a tentar instalar por fora mas não consegui... como te disse pensei que seria talvez incompatibilidade com a versão mai recente do R... Vc lembra em qual versão do R usou o pacote? Tem como disponibilizar o link de onde baixou o arquivo zip com o pacote?? Obrigado!

2 1

Como resolver uma integral dupla definida quando um dos limites dependem de uma das variáveis???
by Romero Luiz M. Sales Filho 06 Abr '15

06 Abr '15

Caros amigos, estou com uma dúvida simples entretanto não consigo solucionar o problema... Já procurei solucionar usando os pacotes R2Cuba e o cubature, entretanto sem sucesso... Quero resolver uma integral dupla entretanto o limite superior de uma das variáveis depende da outra variável, e nos pacotes listados acima consigo apenas colocar limites numéricos... Existe algo que eu possa fazer para resolver esse tipo de problema no R?? A função que desejo integrar é a seguinte: f=(2*exp(-2*x))*(3*exp(-3*y)) Resumindo quero resolver o seguinte: Inf x ∫ ∫ (2*exp(-2*x))*(3*exp(-3*y)) dy dx 0 0 Agradeço desde Já!!! Romero.

2 1

Como fazer esse grafico no ggplot2
by Pedro Emmanuel Alvarenga Americano do Brasil 06 Abr '15

06 Abr '15

Amigos de R, Essa semana eu decidi que iria aprender a usar o ggplot2. Logo na primeira tarefa empaquei uns 40 min. Como fazer esse grafico que faço no graphics com o ggplot2? > rec.cum b2 r3 r24 r25 1 1 2014-04-09 Sim 2014-04-09 10 2 2014-04-15 Não 2014-04-15 12 3 2014-04-16 Sim 2014-04-16 21 4 2014-05-16 Sim 2014-05-16 28 5 2014-05-28 Sim 2014-05-28 35 6 2014-07-29 Sim 2014-07-29 42 7 2014-08-27 Sim 2014-08-27 49 8 2014-10-07 Sim 2014-10-08 55 9 2014-10-22 Sim 2014-10-22 61 10 2014-10-27 Não 2014-10-27 78 11 2014-10-29 Sim 2014-10-29 84 12 2014-11-10 Sim 2014-11-10 89 13 2014-11-10 Sim 2014-11-10 93 14 2014-11-19 Não 2014-11-19 94 15 2014-12-02 Sim 2014-12-02 95 16 2014-12-08 Sim 2014-12-08 99 17 2015-02-25 Sim 2015-02-25 100 18 2015-03-03 Sim 2015-03-03 101 19 2015-03-10 Sim 2015-03-10 102 20 2015-03-17 Sim 2015-03-17 plot.recrut <- function(){ plot(rec.cum$r3,1:nrow(rec.cum),type='b',xlab='',ylab='Número de pacientes rastreados/recrutados',lty=2,col=2,pch=6,cex=1,lwd=3,xlim=c(min(rec.cum$r3,na.rm=T),max(rec.cum$r25,na.rm=T)),axes=F);grid();box();axis(2) tmp <- rec.cum[order(rec.cum$r25),] tmp <- tmp[which(tmp$r24=='Sim'),] lines(tmp$r25,1:nrow(tmp),type='b',lty=3,col=4,pch=2,cex=1,lwd=3) rm(tmp) axis.Date(1,rec.cum$r3,format='%d') # axis.Date(1,rec.cum$r3,format='%b',line=1.5,tick=F) axis.Date(1,rec.cum$r3,format='%Y',line=2.5,tick=F) legend('bottomrigh',c(paste('Primeiro rastreado:',min(rec.cum$r3,na.rm=T)),paste('Último rastreado:',max(rec.cum$r3,na.rm=T))),bty='n') legend('topleft',c(paste('Primeiro recrutado:',min(rec.cum$r25,na.rm=T)),paste('Último recrutado:',max(rec.cum$r25,na.rm=T))),bty='n') legend('top',c('Rastreamento','Recrutamento'),lty=c(2,3),col=c(2,4),pch=c(6,2),pt.cex=1.5,lwd=3,bty='n',xpd=T,inset=-.15,horiz=T) } plot.recrut() Abraço forte, Pedro Brasil

3 5

successive projection algorithm (SPA)
by Ari Clecius 01 Abr '15

01 Abr '15

Boa noite estou realizando uma seleção de variáveis para análise discriminante linear e gostaria de usar o algoritmo SPA, alguém saberia me dizer se existe implementação em R? Att, Ari Clecius Alves de Lima Engenheiro Químico Me. Engenharia Civil (085)84174333 (085)33669042

3 4

Mudar formatação de uma timestamp separando em colunas.
by Michelle Bau Graczyk 31 Mar '15

31 Mar '15

Boa tarde, Eu tenho um arquivo nesse formato e queria que a coluna de timestamp fosse separada em uma de data e outra de hora pois preciso calcular a media do volume por horário. Para isso eu tentei usar o programa abaixo mas a saida foi NA: > INTC<-read.table("/Users/bau/Documents/Mestrado/Dadosbaixadosdothebonnotgang/INTC_1m.txt", header=TRUE, sep=";", dec=",") > INTC[1:25,] symbol nr timestamp open high low close volume 1 INTC 1 2011-09-09 09:30:59 19.7500 19.86 19.75 19.8400 615100 2 INTC 2 2011-09-09 09:31:59 19.8350 19.86 19.79 19.8600 242900 3 INTC 3 2011-09-09 09:32:58 19.8501 19.86 19.80 19.8001 314400 4 INTC 4 2011-09-09 09:33:56 19.8068 19.81 19.77 19.7800 183500 5 INTC 5 2011-09-09 09:34:59 19.7900 19.81 19.78 19.7901 150900 6 INTC 6 2011-09-09 09:35:59 19.8000 19.81 19.74 19.7700 264600 7 INTC 7 2011-09-09 09:36:57 19.7700 19.81 19.77 19.8001 174300 8 INTC 8 2011-09-09 09:37:57 19.8100 19.81 19.76 19.7700 146200 9 INTC 9 2011-09-09 09:38:59 19.7600 19.83 19.76 19.8300 192100 10 INTC 10 2011-09-09 09:39:59 19.8400 19.87 19.83 19.8600 172500 11 INTC 11 2011-09-09 09:40:59 19.8700 19.87 19.76 19.7700 275000 12 INTC 12 2011-09-09 09:41:58 19.7700 19.77 19.72 19.7500 269200 13 INTC 13 2011-09-09 09:42:59 19.7500 19.82 19.75 19.8100 209600 14 INTC 14 2011-09-09 09:43:57 19.8197 19.82 19.79 19.7900 303200 15 INTC 15 2011-09-09 09:44:57 19.7900 19.82 19.78 19.7900 148000 16 INTC 16 2011-09-09 09:45:58 19.7900 19.79 19.75 19.7550 168700 17 INTC 17 2011-09-09 09:46:59 19.7600 19.85 19.75 19.8500 224200 18 INTC 18 2011-09-09 09:47:57 19.8594 19.88 19.85 19.8700 158600 19 INTC 19 2011-09-09 09:48:01 19.8800 19.91 19.87 19.9000 219500 20 INTC 20 2011-09-09 09:49:59 19.8900 19.92 19.89 19.9200 121000 21 INTC 21 2011-09-09 09:50:59 19.9200 19.93 19.90 19.9200 241400 22 INTC 22 2011-09-09 09:51:59 19.9100 19.95 19.90 19.9300 214000 23 INTC 23 2011-09-09 09:52:59 19.9300 19.94 19.92 19.9300 74200 24 INTC 24 2011-09-09 09:53:59 19.9292 19.93 19.88 19.9000 274400 25 INTC 25 2011-09-09 09:54:59 19.9000 19.92 19.87 19.9200 75800 > strptime(INTC$timestamp," %H:%M:%S") [1] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [31] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [61] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [91] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [121] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [151] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [181] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [211] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [241] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [271] NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA NA [ reached getOption("max.print") -- omitted 335637 entries ] Muito obrigada!!!

2 3

Manipulação matriz com laço
by Michelle Bau Graczyk 31 Mar '15

31 Mar '15

Boa noite, Eu estou normalizando autovetores de uma matriz e colocando eles em outra. Para isso eu montei o programa a seguir. A primeira forma, ele só calcula o ultimo elemento da coluna 4 da matriz, e da segunda forma os valores do resultado estão errados. Primeira forma: > a<-read.table("/Users/bau/Documents/Mestrado/teste/eigen1.txt", header=TRUE) > b<-as.matrix(a) > matriznorm<-matrix(0,30,4) > matriznorm[,1]<-b[,1]# só vou usar este valores do arquivo acima > show(matriznorm[,1]) [1] 10.373173554 6.125265914 3.148949951 2.722401273 2.045041946 1.498078581 1.350122780 [8] 1.335457110 1.190051106 0.999610112 0.850212256 0.716720583 0.600704889 0.474433269 [15] 0.431826856 0.388576986 0.266709936 0.176960029 0.045935755 0.020126910 -0.003897135 [22] -0.030136936 -0.074815367 -0.145410630 -0.294553778 -0.453450244 -0.551343912 -0.710794973 [29] -1.222537293 -1.273419528 > for( i in 1:30){ + matriznorm[i,2]<-matriznorm[i,1]*matriznorm[i,1] + } > > somaQuadrados<-sum(matriznorm[,2]) > raiz<-sqrt(somaQuadrados) > matriznorm[,3]<-raiz > > for(j in i:30){ + matriznorm[j,4]<-matriznorm[j,1]/raiz + } > > write.table(matriznorm,"/Users/bau/Documents/Mestrado/teste/eigenValuesNorm1.txt") > show(matriznorm) [,1] [,2] [,3] [,4] [1,] 10.373173554 1.076027e+02 13.46876 0.00000000 [2,] 6.125265914 3.751888e+01 13.46876 0.00000000 [3,] 3.148949951 9.915886e+00 13.46876 0.00000000 [4,] 2.722401273 7.411469e+00 13.46876 0.00000000 [5,] 2.045041946 4.182197e+00 13.46876 0.00000000 [6,] 1.498078581 2.244239e+00 13.46876 0.00000000 [7,] 1.350122780 1.822832e+00 13.46876 0.00000000 [8,] 1.335457110 1.783446e+00 13.46876 0.00000000 [9,] 1.190051106 1.416222e+00 13.46876 0.00000000 [10,] 0.999610112 9.992204e-01 13.46876 0.00000000 [11,] 0.850212256 7.228609e-01 13.46876 0.00000000 [12,] 0.716720583 5.136884e-01 13.46876 0.00000000 [13,] 0.600704889 3.608464e-01 13.46876 0.00000000 [14,] 0.474433269 2.250869e-01 13.46876 0.00000000 [15,] 0.431826856 1.864744e-01 13.46876 0.00000000 [16,] 0.388576986 1.509921e-01 13.46876 0.00000000 [17,] 0.266709936 7.113419e-02 13.46876 0.00000000 [18,] 0.176960029 3.131485e-02 13.46876 0.00000000 [19,] 0.045935755 2.110094e-03 13.46876 0.00000000 [20,] 0.020126910 4.050925e-04 13.46876 0.00000000 [21,] -0.003897135 1.518766e-05 13.46876 0.00000000 [22,] -0.030136936 9.082349e-04 13.46876 0.00000000 [23,] -0.074815367 5.597339e-03 13.46876 0.00000000 [24,] -0.145410630 2.114425e-02 13.46876 0.00000000 [25,] -0.294553778 8.676193e-02 13.46876 0.00000000 [26,] -0.453450244 2.056171e-01 13.46876 0.00000000 [27,] -0.551343912 3.039801e-01 13.46876 0.00000000 [28,] -0.710794973 5.052295e-01 13.46876 0.00000000 [29,] -1.222537293 1.494597e+00 13.46876 0.00000000 [30,] -1.273419528 1.621597e+00 13.46876 -0.09454616 > Da segunda forma: > a<-read.table("/Users/bau/Documents/Mestrado/teste/eigen1.txt", header=TRUE) > b<-as.matrix(a) > matriznorm<-matrix(0,30,4) > matriznorm[,1]<-b[,1] > show(matriznorm[,1]) [1] 10.373173554 6.125265914 3.148949951 2.722401273 2.045041946 1.498078581 1.350122780 [8] 1.335457110 1.190051106 0.999610112 0.850212256 0.716720583 0.600704889 0.474433269 [15] 0.431826856 0.388576986 0.266709936 0.176960029 0.045935755 0.020126910 -0.003897135 [22] -0.030136936 -0.074815367 -0.145410630 -0.294553778 -0.453450244 -0.551343912 -0.710794973 [29] -1.222537293 -1.273419528 > for( i in 1:30){ + matriznorm[i,2]<-matriznorm[i,1]*matriznorm[i,1] + + somaQuadrados<-sum(matriznorm[,2]) + raiz<-sqrt(somaQuadrados) + matriznorm[,3]<-raiz + + for(j in i:30){ + matriznorm[j,4]<-matriznorm[j,1]/raiz + } + } > write.table(matriznorm,"/Users/bau/Documents/Mestrado/teste/eigenValuesNorm1.txt") > show(matriznorm) [,1] [,2] [,3] [,4] [1,] 10.373173554 1.076027e+02 13.46876 1.0000000000 [2,] 6.125265914 3.751888e+01 13.46876 0.5084624670 [3,] 3.148949951 9.915886e+00 13.46876 0.2528991576 [4,] 2.722401273 7.411469e+00 13.46876 0.2135962692 [5,] 2.045041946 4.182197e+00 13.46876 0.1584251428 [6,] 1.498078581 2.244239e+00 13.46876 0.1152793173 [7,] 1.350122780 1.822832e+00 13.46876 0.1033376904 [8,] 1.335457110 1.783446e+00 13.46876 0.1016853664 [9,] 1.190051106 1.416222e+00 13.46876 0.0902440218 [10,] 0.999610112 9.992204e-01 13.46876 0.0755856426 [11,] 0.850212256 7.228609e-01 13.46876 0.0641564605 [12,] 0.716720583 5.136884e-01 13.46876 0.0540043425 [13,] 0.600704889 3.608464e-01 13.46876 0.0452163558 [14,] 0.474433269 2.250869e-01 13.46876 0.0356888678 [15,] 0.431826856 1.864744e-01 13.46876 0.0324667092 [16,] 0.388576986 1.509921e-01 13.46876 0.0292025272 [17,] 0.266709936 7.113419e-02 13.46876 0.0200398900 [18,] 0.176960029 3.131485e-02 13.46876 0.0132951406 [19,] 0.045935755 2.110094e-03 13.46876 0.0034511674 [20,] 0.020126910 4.050925e-04 13.46876 0.0015121392 [21,] -0.003897135 1.518766e-05 13.46876 -0.0002927926 [22,] -0.030136936 9.082349e-04 13.46876 -0.0022641887 [23,] -0.074815367 5.597339e-03 13.46876 -0.0056207917 [24,] -0.145410630 2.114425e-02 13.46876 -0.0109238800 [25,] -0.294553778 8.676193e-02 13.46876 -0.0221227473 [26,] -0.453450244 2.056171e-01 13.46876 -0.0340370871 [27,] -0.551343912 3.039801e-01 13.46876 -0.0413498303 [28,] -0.710794973 5.052295e-01 13.46876 -0.0532327969 [29,] -1.222537293 1.494597e+00 13.46876 -0.0911767942 [30,] -1.273419528 1.621597e+00 13.46876 -0.0945461625 Muito obrigada!

2 4

Sub-matrizes em R
by Tiago Souza Marçal 31 Mar '15

31 Mar '15

Boa noite, alguém sabe como selecionar 2 ou mais sub-matrizes de uma matriz maior de forma direta sem usar um loop???? CMR: mat <- matrix(1:25,nc=5) z = c(1,4)lz = c(2,5)k <- numeric(0) for(j in 1:2) print( k <- mat[z[j]:lz[j],z[j]:lz[j]] ) # Gostaria de chegar a este resultado sem o loop Estou tentando evitar este loop, pois, ele esta dentro de while. Att. Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

2 2

Analise fatorial
by Bruce Kelly 30 Mar '15

30 Mar '15

Boa tarde colegas, Estou tentando usar a função 'factanal' para fazer uma analise fatorial ja utilizei mas a base de dados era um pouco, a mensagem de erro é essa Erro em cov.wt(z) : 'x' deve conter apenas valores finitos Alguém já se deparou com esse problema, ou sabe como resolver? Obrigado Bruce Kelly

2 2

Látice quadrado parcialmente balanceado
by Rodrigo Lima 30 Mar '15

30 Mar '15

Prezados listeiros Bom dia. Eu tenho um experimento em látice 8x8 com duas repetições. Fiz a análise de variância utilizando a função "PBIB" - Analysis of the Partially Balanced Incomplete Block Design. Script: exp01.gy<-PBIB.test(bloco,gen,rep,GY,k=8,group=TRUE, console=TRUE, method="REML") Minha dúvida é como fazer para estimar os resíduos desse modelo, uma vez que não estou conseguindo com a função "residuals". O que você me sugerem fazer? Desde já agradeço. Rodrigo. -- Prof. Rodrigo Oliveira de Lima Federal University of Viçosa Department of Plant Sciences Viçosa, MG, Brazil, 36570-900, Phone: +55(31)3899-2781 Fax: +55(31)3899-2614 Email: rodrigoodelima(a)ufv.br CV: http://lattes.cnpq.br/2761299616351806

2 4

problema GLMM negative binomial
by Ludmila Rattis 30 Mar '15

30 Mar '15

Ola, Tenho uma serie de modelos para rodar e quando tento colocar a funcao glmmadmb, pertencente ao pacote glmmADMB, num looping, retorna o seguinte erro: Error in if (rchar == "") rchar <- "1" : missing value where TRUE/FALSE needed Procurei por esse erro no google e nao encontrei nada semelhante. Quando tento fora do looping, a funcao roda normalmente. Eis um codigo de exemplo: install.packages("glmmADMB", repos="http://r-forge.r-project.org", type="source") data<-matrix(0,180,43);for(i in 1:43){data[,i]<-rnorm(180,mean=0,sd=1) } landscape<-factor(sort(rep(1:30,6))) resp<-paste("resp_",1:37,sep='') pred<-paste("pred_",1:6,sep='') data<-cbind(data,landscape) colnames(data)<-c(resp,pred,"landscape") formulas<-paste("resp",1:37,"~pred",1:6,sep='') fitted.glmm<-list() for (i in 1:length(formulas)){ fitted.glmm[[]]<-glmmadmb(formulas[i],data=data,random=1|~landscape,family="nbinom",zeroInflation=FALSE)} Alguem poderia me dizer o que esta acontecendo? Muito obrigada, Lud -- *Ludmila Rattis* Programa de Pós-graduação em Ecologia/UNICAMP Geomatics and Landscape Ecology Research Laboratory (Carleton University - Canada) http://www.glel.carleton.ca Conservation Biogeography Lab (Goias Federal University - Brazil) rdloyola.wix.com/cblab <http://sites.ffclrp.usp.br/ficus>

2 2

Re: [R-br] Organização de banco de dados!
by Fernando Mayer 30 Mar '15

30 Mar '15

Andre, (re-postando para a lista) no seu exemplo nao aparece nenhuma UF repetida, por isso essa abordagem funcionaria. Se no seu outro arquivo de dados tiver um nome repetido isso já não vai funcionar mesmo, pois não é sintaticamente permitido no R. Uma opção (não muito ideal) nesse caso seria importar os dados com a primeira coluna e colocar um índice: dados <- read.csv("total.txt") rownames(dados) <- make.names(dados$UF, unique = T) --- Fernando de Pol Mayer Laboratório de Estatística Ambiental - LEA Instituto de Matemática, Estatística e Física - IMEF Universidade Federal do Rio Grande - FURG URL: http://fernandomayer.github.io e-mail: fernando.mayer [@] {gmail.com, furg.br} 2015-03-29 18:32 GMT-03:00 Andre Oliveira <andreolsouza(a)yahoo.com.br>: > Pois é Fernando já havia olhado ajuda! > > dados<-read.table("Jaq4.txt", h=T, row.names = 1) > > Erro em read.table("Jaq4.txt", h = T, row.names = 1) : > duplicate 'row.names' are not allowed > > > E além de tudo cria uma outra coluna que sempre não esta no arquivo txt. > Tenho que remover as duas colunas. > 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 > > > > > Em Domingo, 29 de Março de 2015 18:23, Fernando Mayer < > fernandomayer(a)gmail.com> escreveu: > > > ?read.table > dados <- read.table("total.txt", row.names = 1) # assumindo que a coluna > UF seja a 1 > > > > > --- > Fernando de Pol Mayer > Laboratório de Estatística Ambiental - LEA > Instituto de Matemática, Estatística e Física - IMEF > Universidade Federal do Rio Grande - FURG > URL: http://fernandomayer.github.io > e-mail: fernando.mayer [@] {gmail.com, furg.br} > > 2015-03-29 17:38 GMT-03:00 Andre Oliveira <andreolsouza(a)yahoo.com.br>: > > Estou tentando arrumar um banco de dados como sugerido no link > http://lincolntneves.weebly.com/blog/anlise-de-cluster-no-r. No caso > preciso que a coluna UF sem sem rótulo de UF como mostro abaixo. Se eu > deixar sem o UF o R não lê o banco de dados. Como eu poderia resolver isto? > Alem do mais sempre que importo com o rótulo da coluna UF ainda cria uma > outra coluna com as numerações 1,2,3,...27 antes da coluna UF. Enfim eu > preciso do data frame exatamente como mostro abaixo. Sem rótulo UF e sem > coluna antes dos estados. Obrigado por qualquer contribuição! > > > dados=read.table("total.txt")dados A B B C > AC 1 38 242 983 > AL 10 65 281 2194 > AM 32 108 319 3462 > AP 0 33 242 297 > BA 135 214 480 12068 > CE 178 301 617 13156 > DF 169 313 585 5655 > ES 63 189 456 6995 > GO 40 137 391 8365 > MA 8 63 271 3868 > MG 1027 2022 5663 77191 > MS 95 194 397 4916 > MT 19 75 267 4140 > PA 19 76 295 4858 > PB 25 59 266 2883 > PE 133 269 477 7418 > PI 42 95 283 3830 > PR 284 542 1360 25057 > RJ 542 827 1368 16173 > RN 41 110 352 3055 > RO 9 52 243 2842 > RR 2 38 266 492 > RS 226 393 1028 20525 > SC 113 269 833 14969 > SE 10 50 252 886 > SP 772 1713 4842 86742 > TO 11 68 244 2495 > > > > > > > > > > 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 > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > > >

3 2

Re: [R-br] Organização de banco de dados!
by Fernando Mayer 29 Mar '15

29 Mar '15

Andre (favor responder para o email da lista, assim a conversa fica arquivada) o USArrests está exatamene como o primeiro exemplo que você mandou, nenhum nome de linha se repete. Se tiver pelo menos um nome de linha repetido então não tem como transformar isso em rownames. --- Fernando de Pol Mayer Laboratório de Estatística Ambiental - LEA Instituto de Matemática, Estatística e Física - IMEF Universidade Federal do Rio Grande - FURG URL: http://fernandomayer.github.io e-mail: fernando.mayer [@] {gmail.com, furg.br} 2015-03-29 18:50 GMT-03:00 Andre Oliveira <andreolsouza(a)yahoo.com.br>: > Este banco este como eu gostaria USArrests! > > > > > 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 > > > > > Em Domingo, 29 de Março de 2015 18:32, Andre Oliveira < > andreolsouza(a)yahoo.com.br> escreveu: > > > Pois é Fernando já havia olhado ajuda! > > dados<-read.table("Jaq4.txt", h=T, row.names = 1) > > Erro em read.table("Jaq4.txt", h = T, row.names = 1) : > duplicate 'row.names' are not allowed > > > E além de tudo cria uma outra coluna que sempre não esta no arquivo txt. > Tenho que remover as duas colunas. > 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 > > > > > Em Domingo, 29 de Março de 2015 18:23, Fernando Mayer < > fernandomayer(a)gmail.com> escreveu: > > > ?read.table > dados <- read.table("total.txt", row.names = 1) # assumindo que a coluna > UF seja a 1 > > > > > --- > Fernando de Pol Mayer > Laboratório de Estatística Ambiental - LEA > Instituto de Matemática, Estatística e Física - IMEF > Universidade Federal do Rio Grande - FURG > URL: http://fernandomayer.github.io > e-mail: fernando.mayer [@] {gmail.com, furg.br} > > 2015-03-29 17:38 GMT-03:00 Andre Oliveira <andreolsouza(a)yahoo.com.br>: > > Estou tentando arrumar um banco de dados como sugerido no link > http://lincolntneves.weebly.com/blog/anlise-de-cluster-no-r. No caso > preciso que a coluna UF sem sem rótulo de UF como mostro abaixo. Se eu > deixar sem o UF o R não lê o banco de dados. Como eu poderia resolver isto? > Alem do mais sempre que importo com o rótulo da coluna UF ainda cria uma > outra coluna com as numerações 1,2,3,...27 antes da coluna UF. Enfim eu > preciso do data frame exatamente como mostro abaixo. Sem rótulo UF e sem > coluna antes dos estados. Obrigado por qualquer contribuição! > > > dados=read.table("total.txt")dados A B B C > AC 1 38 242 983 > AL 10 65 281 2194 > AM 32 108 319 3462 > AP 0 33 242 297 > BA 135 214 480 12068 > CE 178 301 617 13156 > DF 169 313 585 5655 > ES 63 189 456 6995 > GO 40 137 391 8365 > MA 8 63 271 3868 > MG 1027 2022 5663 77191 > MS 95 194 397 4916 > MT 19 75 267 4140 > PA 19 76 295 4858 > PB 25 59 266 2883 > PE 133 269 477 7418 > PI 42 95 283 3830 > PR 284 542 1360 25057 > RJ 542 827 1368 16173 > RN 41 110 352 3055 > RO 9 52 243 2842 > RR 2 38 266 492 > RS 226 393 1028 20525 > SC 113 269 833 14969 > SE 10 50 252 886 > SP 772 1713 4842 86742 > TO 11 68 244 2495 > > > > > > > > > > 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 > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > > > > > > >

1 0

Organização de banco de dados!
by Andre Oliveira 29 Mar '15

29 Mar '15

Estou tentando arrumar um banco de dados como sugerido no link http://lincolntneves.weebly.com/blog/anlise-de-cluster-no-r. No caso preciso que a coluna UF sem sem rótulo de UF como mostro abaixo. Se eu deixar sem o UF o R não lê o banco de dados. Como eu poderia resolver isto? Alem do mais sempre que importo com o rótulo da coluna UF ainda cria uma outra coluna com as numerações 1,2,3,...27 antes da coluna UF. Enfim eu preciso do data frame exatamente como mostro abaixo. Sem rótulo UF e sem coluna antes dos estados. Obrigado por qualquer contribuição! dados=read.table("total.txt") dados A B B C AC 1 38 242 983 AL 10 65 281 2194 AM 32 108 319 3462 AP 0 33 242 297 BA 135 214 480 12068 CE 178 301 617 13156 DF 169 313 585 5655 ES 63 189 456 6995 GO 40 137 391 8365 MA 8 63 271 3868 MG 1027 2022 5663 77191 MS 95 194 397 4916 MT 19 75 267 4140 PA 19 76 295 4858 PB 25 59 266 2883 PE 133 269 477 7418 PI 42 95 283 3830 PR 284 542 1360 25057 RJ 542 827 1368 16173 RN 41 110 352 3055 RO 9 52 243 2842 RR 2 38 266 492 RS 226 393 1028 20525 SC 113 269 833 14969 SE 10 50 252 886 SP 772 1713 4842 86742 TO 11 68 244 2495 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 1

for dentro de while
by Tiago Souza Marçal 29 Mar '15

29 Mar '15

Boa noite, estou com um código que necessita o uso de um for dentro de um while. Entretanto, o uso do for aumentou o tempo de execução. Gostaria de saber se alguém tem alguma sugestão para contornar este problema???? Acredito que a questão seja organizacional, portanto, não estou enviando um CMR mas caso seja necessário será encaminhado!!!!!!! dif <- 1 tol <- 1e-05 i <- 0 neff <- c(2) while(dif > tol){ i <- i + 1 LHS <- rBind( cBind(XpX, XpZ), cBind(ZpX, (ZpZ + Ginv)) ) invLHS <- solve(LHS) est <- solve(LHS,RHS) b <- est[1:ncol(X.f),] u <- est[(ncol(X.f) + 1):nrow(est),] sig2E <- (ypy - t(b)%*%Xpy - t(u)%*%Zpy)/(N - rankX) sig2U <- numeric(0) for(j in 1:length(neff)){ sig2U <- c(sig2U, ( t(est[z[j]:lz[j],])%*%est[z[j]:lz[j],] + sig2E*sum(diag( invLHS[z[j]:lz[j], z[j]:lz[j]] )) )/neff[j]) } sig2U <- do.call(rBind,lapply(sig2U,as.matrix)) rel <- as.numeric(sig2E)/sig2U if(!is.null(var)){ lamb <- rel ratio <- rep(lamb,neff) newvar <- rBind(sig2U,sig2E) cds <- sig2U/sum(newvar) dif <- as.numeric(max(abs(newvar - ovar))) } if(is.null(var)){ lamb <- rel ratio <- rep(lamb,neff) newvar <- rBind(sig2U,sig2E) cds <- sig2U/sum(newvar) dif <- as.numeric(max(abs(cds - ocds))) } ovar <- newvar ocds <- cds oratio <- ratio olamb <- lamb Ginv <- Diagonal(sum(neff),oratio) }Desde já agradeço por qualquer colaboração. Atenciosamente, Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

2 1

Manipulação de banco de dados
by Andre Oliveira 29 Mar '15

29 Mar '15

Pessoal tenho estes dados dentro do R, gostaria de criar um outro banco com a soma de cada item dos estados ao longo destes 10 anos. str(dados) data.frame': 270 obs. of 7 variables: $ ANO : int 2005 2005 2005 2005 2005 2005 2005 2005 2005 2005 ... $ UF : Factor w/ 27 levels "AC","AL","AM",..: 1 2 3 4 5 6 7 8 9 10 ... $ OURO : int 0 1 2 0 15 14 7 7 2 0 ... $ PRATA : int 15 15 15 15 15 15 15 15 15 15 ... $ BRONZE: int 15 15 15 15 15 15 15 15 15 15 ... $ MENCAO: int 136 543 718 98 2976 2242 306 377 859 830 ... $ TOTAL : int 166 574 750 128 3021 2286 343 414 891 860 ...head(dados) | ANO UF OURO PRATA BRONZE MENCAO TOTAL 1 2005 AC 0 15 15 136 166 2 2005 AL 1 15 15 543 574 3 2005 AM 2 15 15 718 750 4 2005 AP 0 15 15 98 128 5 2005 BA 15 15 15 2976 3021 6 2005 CE 14 15 15 2242 2286 | | | | | > | | 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

3 2

Leitura de dados
by Andre Oliveira 28 Mar '15

28 Mar '15

Pessoal,como posso remover as posições da coluna 1 que não está no banco de dados? Este exemplo é da lista mesmo, usei só para ilustração! require(gdata) # contém a função read.xls() para importar planilhasteste <- read.xls("http://www.leg.ufpr.br/~walmes/foo/teste.xls", sheet=1, method="tab") teste$Municipio <- iconv(teste$Muni, "latin1", "UTF-8") head(teste) Sigla Codigo Municipio a2000 1 AC 1200013 Acrelândia 1042.38 2 AC 1200054 Assis Brasil 429.22 3 AC 1200104 Brasiléia 2243.25 4 AC 1200138 Bujari 786.68 5 AC 1200179 Capixaba 603.84 6 AC 1200203 Cruzeiro do Sul 9363.32 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 1

Função identify
by ana paula coelho madeira 27 Mar '15

27 Mar '15

Caros colegas, boa noite! Ao fazer o diagnóstico de um GLM tento identificar os pontos de alavanca e pontos influentes utilizando a função identify( ), porém o R fia executando e nunca libera o resultado, daí tenho que parar a computação. Alguém tem ideia do pode estar acontecendo? Desde já agradeço. Ana Paula

4 5

Diferença entre dados coletados do yahoo usando diferentes funções.
by Michelle Bau Graczyk 26 Mar '15

26 Mar '15

Boa tarde, Eu estava verificando os diferentes modos de coletar dados intraday do yahoo e percebi que para cada função que eu uso, os dados obtidos são diferentes. Por exemplo, se eu uso: dados <- readLines(' http://chartapi.finance.yahoo.com/instrument/1.0/AAPL/chartdata;type=quote;…',sep = ",", dec = ".") eu obtenho: "18" "1427203800,127.1800,127.2600,127.0000,127.2300,631900" "19" "1427203868,126.9200,127.1900,126.8750,127.1700,202000" "20" "1427203920,126.9600,127.0600,126.8700,126.9100,189500" "21" "1427203981,126.9500,126.9900,126.8800,126.8900,198200" "22" "1427204099,126.7799,126.9900,126.7700,126.9500,217500" "23" "1427204100,126.6700,126.9000,126.6401,126.7799,286300" "24" "1427204219,126.7100,126.8000,126.6500,126.6800,222800" "25" "1427204221,126.8300,126.8700,126.5600,126.7000,238500" "26" "1427204280,126.7100,126.9200,126.6700,126.8200,167500" "27" "1427204340,126.8300,126.9300,126.7000,126.7166,114900" "28" "1427204457,126.6400,126.8300,126.6000,126.8100,156100" "29" "1427204460,126.6500,126.7900,126.6385,126.6900,102800" se eu uso: dados <- getURL(' http://chartapi.finance.yahoo.com/instrument/1.0/AAPL/chartdata;type=quote;…',sep = ",", dec = ".") obtenho: 1427203805,127.1900,127.2600,127.0796,127.2000,626100 1427203862,126.9200,127.1900,126.8750,127.1800,207800 1427203979,126.9000,127.0600,126.8700,126.9100,198000 1427203980,126.9500,126.9900,126.8800,126.9000,172200 1427204099,126.7799,126.9900,126.7700,126.9500,235000 1427204100,126.6700,126.9000,126.6401,126.7799,286300 1427204219,126.7100,126.8000,126.6500,126.6800,222800 1427204221,126.8300,126.8700,126.5600,126.7000,238500 1427204280,126.7100,126.9200,126.6700,126.8200,167500 1427204340,126.8300,126.9300,126.7000,126.7166,114900 1427204459,126.6915,126.8300,126.6100,126.8100,160900 1427204460,126.6500,126.7900,126.6385,126.6900,98000 Usando read.csv também dá dados diferentes de começo. Alguem saberia me dizer o por que? E em qual devo confiar? Porque nenhuma é igualzinha ao que aparece quando você coloca esse endereço http://chartapi.finance.yahoo.com/instrument/1.0/AAPL/chartdata;type=quote;… no navegador. O resultado deste é: 1427203800,127.1200,127.2300,127.0000,127.2300,631900 1427203900,127.0300,127.1900,127.0000,127.1700,87400 1427203978,126.9200,127.0400,126.8700,126.9050,309600 1427204038,126.9600,126.9900,126.8800,126.9399,174900 1427204099,126.7799,126.9900,126.7700,126.9600,235300 1427204100,126.6700,126.9000,126.6401,126.7799,286300 1427204219,126.7100,126.8000,126.6500,126.6800,222800 1427204221,126.8300,126.8700,126.5600,126.7000,238500 1427204280,126.7100,126.9200,126.6700,126.8200,167500 1427204340,126.8300,126.9300,126.7000,126.7166,114900 1427204459,126.6915,126.8300,126.6000,126.8100,160900 1427204460,126.6500,126.7900,126.6385,126.6900,98000 1427204520,126.7100,126.7200,126.6000,126.6500,108700 1427204581,126.7100,126.8400,126.6800,126.7100,129600 1427204641,126.7400,126.7700,126.6600,126.7300,98700 Obrigada, Michelle

3 3

Pacote QCC
by Edson Lira 26 Mar '15

26 Mar '15

Caros amigos, estou tentando replicar um exemplo do Livro do Montgomery, Introduction to Stastistical Quality Control - 7 edição, pg 240. Estou usando o pacote QCC do R. Abaixo o banco e o script Download iqcpg240.csv | | | | | | | | | | | Download iqcpg240.csvDownload iqcpg240.csv free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Download exemplo iqc wiley pg240.R | | | | | | | | | | | Download exemplo iqc wiley pg240.RDownload exemplo iqc wiley pg240.R free. File hosted free by DataFileHost. | | | | Visualizar em www.datafilehost... | Visualizado por Yahoo | | | | | Tenho algumas dúvidas. Para os limites x de barra usa-seXbarra+-A2*Rbarra, aonde o A2 é tabelado de acordo com o tamanho da amostra. No pacote QCC, fiz algumas tentativas e não vi como ou aonde ele usa o A2, visto que quando plota os lintel ele só soma ou diminui o desvio padrão. Alguém poderia me ajudar?. No livro do Montgomery os limites para o barra são (25 amostra de tamanho 5)1,5056+-(0,577)*(0,32521) Não consegui replicar esses resultados no gráfico. Usei as 25 primeiras amostras do banco, que é o que o livro usa no exemplo. Alguém poderia me ajudar? [ ]'s.Edson Lira Estatístico Manaus-Amazonas

1 0