Abril 2015 - R-br - DINF UFPR

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