- R-br - DINF UFPR

ESTIMATIVA DA CONCENTRAÇÃO LETAL MÉDIA (CL 50)
by Tiago Souza Marçal 07 Fev '13

07 Fev '13

Ola pessoal alguém conhece alguma rotina para a estimamação da CL 50 no R? Att. Tiago.

3 3

como contar elementos em um data.frame
by Aline Lipsky 07 Fev '13

07 Fev '13

> Boa tarde, > > Eu tenho um data.frame como o abaixo: > > Data Embarcação Municipio > Jan A RG > Jan B PEL > Jan J PEL > Dez J RG > Mar A RG > Set F PEL > Set G RG > Dez G RG > > > Preciso contar as embarcações por mês e Município e obter uma tabela assim: > > Mês RG PEL > Jan 5 16 > Fev 26 31 > Mar 7 2 > > Grata pela atenção, > > Aline

2 4

Duvida sobre como alterar o contraste no options de "contr.poly" para "contr.sum" e entender o que acontece.
by Augusto Ribas 07 Fev '13

07 Fev '13

Ola caros colegas da lista, eu estava tentando estabelecer a matriz de contrastes para usar na função lm, mas não estou conseguindo entender como as coisas funcionam no R. Eu li um tutorial antigo, esse aqui (http://wwwuser.gwdg.de/~cscherb1/content/Statistics%20Course%20files/Workin…) Eu não entendo se as coisas mudaram desde que esse tutorial foi escrito, mas com matriz de soma como no exemplo acima tudo é simples, mas eu não consigo mais replicar o exemplo no R. Mais precisamente o que eu estava tentando era o seguinte: #Dados de exemplo: variavel<-factor(rep(letters[1:3],each=10)) resposta<-rnorm(30,rep(c(1,1,2),each=10)) #Ai eu ajusto um modelo linear #modelo usando opções default modelo<-lm(resposta~variavel) summary(modelo) #Eu imagino que a matriz que deveria fazer os contrastes como acima seria essa no caso de usar soma como no texto cbind(c(-1,1,0),c(-1,0,1)) Mas o default usa outro esquema de matriz de contrastes, logo so colocar isso nos contrastes não funciona #O default não usa esse sistema certo? options("contrasts") #Eu imaginei que essa opção que fosse a que eu precisaria mudar options("contrasts")$contrasts[2] #Mas tentando mudar de forma simples assim eu obtenho um erro #mas não é uma mudança simples assim, quando eu tento isso: options("contrasts")$contrasts[2]<-"contr.sum" #eu tenho essa mensagem de Erro Erro em options("contrasts")$contrasts[2] <- "contr.sum" : destino da atribuição expande para um objeto fora da linguagem Mas como eu faço ao invéz de escrever isso, contrasts(variavel) #Escrever isso, ou algo assim cbind(c(-1,1,0),c(-1,0,1)) #E obter o mesmo resultado que esta no summary(modelo) la em cima, se não me perdi em algum momento mais para cima. #Porque mesmo a opção contr.sum parece que ele trabalha de outra forma do que estou imaginando. contr.sum(3,contrasts=T,sparse=F) Bem é isso, é apenas uma duvida sobre como o R funciona que não estou conseguir entender, se alguém puder dar uma luz eu serei grato. Abraços Augusto Ribas -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

2 5

Contar sequência de números...
by andrebvs＠bol.com.br 07 Fev '13

07 Fev '13

3 2

percentage de error da predicao de GLM
by Humberto Hazin 07 Fev '13

07 Fev '13

Bom dia Pessoal estou tentando calcular o percentual de predição de erro via boot, porem o resultado final não me parece o que quero (percentual de predição de erro ) alguém poderia me ajudar? train.pos.setbsh data set com ~27000 linhas (https://www.dropbox.com/s/9flxc9tpttq1mru/train.pos.setbsh.csv) test.setbsh data set com ~1433 linhas (https://www.dropbox.com/s/61fbts7bbizb7cr/test.setbsh.csv) bt <- function(formula, data, indices) { d <- data[indices,] fit <- glm(formula,family=gaussian,data=d) pred<-predict(fit,test.setbsh,type="response") error<- (test.setbsh$lgcpue-pred)/test.setbsh$lgcpue return(error) } # bootstrapping with 300 replications results <- boot(data=train.pos.setbsh, statistic=bt, R=300, formula=lgcpue~Y+S+Q+A) boxplot(as.vector(results$t)) abline(0,0) desde já agradeço Humberto Enviado do Email do Windows

1 0

ERRO NA CRIAÇÃO DE UMA MATRIZ
by Tiago Souza Marçal 07 Fev '13

07 Fev '13

Boa noite pessoal, estou tentando criar uma matriz com os intervalos de confiança que eu calculei para os um conjunto de dados. No entanto, os valores calculados para o intervalo de confiança ficaram ocultos na matriz. ic<- function(x, sigma2, conf.level=0.95){ n <- length(x) xbarra <- mean(x) if (!missing(sigma2)){ z <- qnorm((1+conf.level)/2) }else{ sigma2 <- var(x) z <- qt((1+conf.level)/2, n-1) } return(mean(x)+c(-1, 1)*z*sqrt(sigma2)/sqrt(n)) } tr<-gl(3,4)resp<-rnorm(12,30,5)d<-data.frame(tr,resp)int.conf<-tapply(resp,tr,ic)mat<-matrix(int.conf)mat # matriz obtida. [,1] [1,] Numeric,2[2,] Numeric,2[3,] Numeric,2 Agradeço desde já pela ajuda. Att. Tiago.

2 2

Fwd:
by Leonardo Aguirre 06 Fev '13

06 Fev '13

http://www.f4fotografi.it/a00wit.php?s=ot

1 0

Montar expressão da Gama no R...
by andrebvs＠bol.com.br 06 Fev '13

06 Fev '13

2 2

Otimizando um loop
by Pedro Emmanuel Alvarenga Americano do Brasil 06 Fev '13

06 Fev '13

Amigos de R, Gostaria de uma ajuda ara otimizar um loop. Simulei algumas populaçoes com resultados de testes. Agora quero amostrar mil amostras de cada população e cada set de 1000 amostras em algumas condições como tamanho da amostra e prevalencia do evento do interesse. Ao total seriam mais de 1500 amostragens. Por isso, qualquer dica para torna-lo mais rápido e eficiente será bem vinda. Eu sei que faltou diversas populações.... mas é que o arquivo é grande. O script inicial... with <- 250000 without <- 250000 set.seed(as.numeric(Sys.time())) #Population 1: bi-normal distribution with equal varainces # Population 1 with AUC = .9 pop1 <- round(rnorm(without,200,60)) pop1 <- as.data.frame(cbind(0,pop1)) colnames(pop1) <- c('reference','test') # pop1[1:10,];nrow(pop1) pop1temp <- round(rnorm(without,309,60)) pop1temp <- as.data.frame(cbind(1,pop1temp)) colnames(pop1temp) <- c('reference','test') # pop1temp[1:10,]; nrow(pop1temp) pop1_9 <- rbind(pop1,pop1temp) # nrow(pop1) rm(pop1temp,pop1) .... outras populações.... > population <- ls() > samp_size <- c(50,100,150,300,500,1000,3000) > prev <- c(.05,.1,.2,.3,.5,.75) > > conditions <- expand.grid(population,samp_size,prev) > conditions <- cbind(conditions[,1],substr(conditions[,1],1,4),as.numeric(substr(conditions[,1],6,6))/10,conditions[,2:3]) > names(conditions) <- c('dataset','population','auc','samp_size','prevalence') > conditions$dataset <- as.character(conditions$dataset) > head(conditions) dataset population auc samp_size prevalence 1 pop1_6 pop1 0.6 50 0.05 2 pop1_7 pop1 0.7 50 0.05 3 pop1_8 pop1 0.8 50 0.05 4 pop1_9 pop1 0.9 50 0.05 5 pop2_6 pop2 0.6 50 0.05 6 pop2_7 pop2 0.7 50 0.05 > setwd('E:/Banco/cutoff/datasets') # O problema começa aqui onde eu preciso amostrar de diversos bancos de dados em sequencia armazenados na área de trabalho, mas não consigo fazer isso de forma elegante. então eu acabei usando o parse mas antes preciso salvar o nome do banco no disco, e quando consigo fazer isso vem um aviso que não consigo entender. # Aqui deveria começar o loop seguindo as condições do banco conditions. > cat(conditions[1,1],file="expr"); dataux <- eval(parse(file="expr")) Mensagens de aviso perdidas: In readLines(file) : linha final incompleta encontrada em 'expr' > with <- sample(dataux[which(dataux$reference==1),2],ceiling(conditions[1,4]*conditions[1,5]),T) > without <- sample(dataux[which(dataux$reference==0),2],conditions[1,4]*(1-conditions[1,5]),T) > with <- cbind(reference=1,test=with) > without <- cbind(reference=0,test=without) > tempdata <- data.frame(rbind(with,without)) > write.csv2(tempdata,paste0(conditions[1,1],'_',conditions[1,4],'_',substr(conditions[1,5],3,4),'.csv'),row.names = F) > # aqui deveria terminar o loop Abraço, Dr. Pedro Emmanuel A. A. do Brasil Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365, CEP 21040-360, Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

1 0

Re: [R-br] Definir classes no histograma.
by Ivan Bezerra Allaman 06 Fev '13

06 Fev '13

Veja a library fdth, as opções e a qualidade de histogramas é infinitamente melhor do que a "base". (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}

1 0

Re: [R-br] levelplot com escalas diferentes
by Carolina R. Caiado Gomes 06 Fev '13

06 Fev '13

Oi Elias, Tenta o seguinte: library(gridExtra) library(lattice) x <- seq(pi/4, 5 * pi, length.out = 100) y <- seq(pi/4, 5 * pi, length.out = 100) r <- as.vector(sqrt(outer(x^2, y^2, "+"))) grid <- expand.grid(x=x, y=y) grid$z <- cos(r^2) * exp(-r/(pi^3)) plot1<-levelplot(seu plot) plot2<-levelplot(seu outro plot) grid.arrange(plot1,plot2) Att, -- Carolina Ramos Caiado Gomes Mestranda em Ecologia e Evolução - UFG Laboratório de Ecologia de Insetos +55 62 8215 3595

2 1

Definir classes no histograma.
by samuel luna de almeida 06 Fev '13

06 Fev '13

Bom dia, Pessoal, desculpem se essa dúvida seria solucionada no help, pois parece bem simples , mas não encontrei. Gostaria de ajuda, para que eu seja quem defina, quais serão os valores de corte para as colunas do histograma. Abs Samuel

3 7

GEE2 - dados correlacionados no espaço
by Márcio Mendonça 06 Fev '13

06 Fev '13

Prezados colegas, preciso de um script do R, caso alguém tenha, de um exemplo que contemple uma avaliação de pontos (oriundos de riqueza e abundancia spp) distribuidos no espaço. Creio que o pacote GEE2 do R seria a ferramenta mais adequada por se tratar de dados possivelmente correlacionados no espaço. Agradeço a atenção de todos Márcio ________________________________ De: "r-br-request(a)listas.c3sl.ufpr.br" <r-br-request(a)listas.c3sl.ufpr.br> Para: r-br(a)listas.c3sl.ufpr.br Enviadas: Sexta-feira, 1 de Fevereiro de 2013 11:00 Assunto: Digest R-br, volume 24, assunto 1 Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Nayara Santos is out of the office. (Nayara Santos) 2. Informação: R como software livre... (andrebvs(a)bol.com.br) 3. Re: Informação: R como software livre... (Manoel Galdino) 4. Gráfico dupla entrada (Stephan M. CARVALHO) 5. Re: Gráfico dupla entrada (Paulo Justiniano) 6. Re: Índice fora de limites (Leonard de Assis) 7. Re: Gráfico dupla entrada (Luciano F. Sgarbi) 8. Re: Informação: R como software livre... (Leonard de Assis) 9. Re: Informação: R como software livre... (Manoel Galdino) 10. Re: Gráfico dupla entrada (Stephan M. CARVALHO) ---------------------------------------------------------------------- Message: 1 Date: Thu, 31 Jan 2013 15:03:11 -0300 From: Nayara Santos <nayara.santos(a)br.zurich.com> To: r-br <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Nayara Santos is out of the office. Message-ID: <OF50BC777C.036D38FD-ON83257B04.00632B35-83257B04.00632B35(a)email.zurich.com> Content-Type: text/plain; charset=ISO-8859-1 I will be out of the office starting 30/01/2013 and will not return until 14/02/2013. Gentileza entrar em contato com: Flávia Novaes para testes. Marcelo Prieto e Vicente Venancio para definições, dúvidas e aprovações. ------------------------------ Message: 2 Date: Thu, 31 Jan 2013 17:47:00 -0200 From: andrebvs(a)bol.com.br To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Informação: R como software livre... Message-ID: <510aca34847f8_356468b6c7093e(a)a4-winter20.tmail> Content-Type: text/plain; charset="iso-8859-1" Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130131/6a13a117/atta…> ------------------------------ Message: 3 Date: Thu, 31 Jan 2013 17:51:47 -0200 From: Manoel Galdino <mcz.fea(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Informação: R como software livre... Message-ID: <CA+kPArkJSN773t+G=VFTzpRBY4eSrH5ZxTY-rQ4DXmFc0HDykw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Para todos os itens: 1. Redução de custos ( o R é grátis, exceto pelo Revolutions) 2. Por ter código aberto, é mais fácil checar a análise feita. Para pesquisa: A comunidade de usuários dá tanto suporte que ajuda o pesquisador a criar laços com outros pesquisadores abçs Manoel 2013/1/31 <andrebvs(a)bol.com.br> > Olá à todos! > > O programa R, sem dúvida, é a ferramenta estatística ideal nas mãos de > quem o conhece, mas, porque universidades ainda insistem em utilizar > programas comerciais (propriétarios) em seus laboratórios? acredito que > esse fato, ainda seja pela falta de interesse de professores > "tradicionalistas" que não instigam mudanças (inovações) na forma de > ensino. Talvez, falta a esses professores, a informação sobre os beneficios > que os softwares livres trazem para a instituição e para o governo. Então, > diante disso, como justificar o uso do programa R como software livre em > cada uma das etapas a seguir dentro das Universidades: > > 1) ensino > 2) pesquisa > 3) extensão > > desde já agradeço! > > *Att.* > *André* > > _______________________________________________ > 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. > -- Manoel Galdino https://sites.google.com/site/galdinomcz/ -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130131/ea09fb48/atta…> ------------------------------ Message: 4 Date: Fri, 01 Feb 2013 11:37:59 -0300 From: "Stephan M. CARVALHO" <smalfitano(a)uol.com.br> To: R-br <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Gráfico dupla entrada Message-ID: <510BD347.7050209(a)uol.com.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Bom dia pessoal, Preciso plotar dados em um gráfico de dupla entrada, como ocorre por exemplo com dados climatológicos (temperatura/precipitação ~ tempo). Por acaso alguém conhece alguma forma de fazer? Obrigado, Stephan -- *************************************************** STEPHAN M. CARVALHO "Honeybee protection and ecotoxicology" Universidade Estadual Paulista Centro de Estudo de Insetos Sociais Av. 24A, 1515 Bela Vista 13.506-900 Rio Claro/SP Brazil smalfitano(a)uol.com.br +55 19 3526-4284 (university) +55 19 3534-8523 (fax) +55 35 9103-9822 (mobile) +55 35 9949-0707 (mobile) Antes de imprimir, pense no ambiente! Before printing, think in the environment! Avant d'imprimer, pensez à l'environnement! *************************************************** -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: smalfitano.vcf Tipo: text/x-vcard Tamanho: 352 bytes Descrição: não disponível URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130201/ba8d54f0/atta…> ------------------------------ Message: 5 Date: Fri, 1 Feb 2013 09:10:36 -0200 (BRST) From: Paulo Justiniano <paulojus(a)leg.ufpr.br> To: R-br <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Gráfico dupla entrada Message-ID: <alpine.DEB.2.00.1302010907590.7577(a)pataxo.est.ufpr.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" O segredo está em usar par(new=T) npara justapor gráficos O exemplo a seguir foi tirado dos arquivos da lista http://r-br.2285057.n4.nabble.com/R-br-Grafico-com-dois-eixos-y-outra-soluc… x <- 1:20 x1 <- rnorm(20) x2 <- rpois(20, lam=50) plot(x, x1, type="l") par(new=T) plot(x, x2, type="l", axes=F, frame=T, ann=F, col=2) axis(4) e um exemplo mais completo tb apontado em um email da lista: http://r-br.2285057.n4.nabble.com/R-br-graficos-com-dois-eixos-y-td3566075.… sugere procurar em: http://www.leg.ufpr.br/doku.php/ridiculas por Gráfico com dois eixos coordenados On Fri, 1 Feb 2013, Stephan M. CARVALHO wrote: > Bom dia pessoal, > > > Preciso plotar dados em um gráfico de dupla entrada, como ocorre por > exemplo com dados climatológicos (temperatura/precipitação ~ tempo). > > Por acaso alguém conhece alguma forma de fazer? > > Obrigado, > > > Stephan > > > -- > > > *************************************************** > > STEPHAN M. CARVALHO > > "Honeybee protection and ecotoxicology" > > Universidade Estadual Paulista > Centro de Estudo de Insetos Sociais > > Av. 24A, 1515 > Bela Vista > 13.506-900 Rio Claro/SP > Brazil > > smalfitano(a)uol.com.br > > +55 19 3526-4284 (university) > +55 19 3534-8523 (fax) > +55 35 9103-9822 (mobile) > +55 35 9949-0707 (mobile) > > > Antes de imprimir, pense no ambiente! > Before printing, think in the environment! > Avant d'imprimer, pensez à l'environnement! > > *************************************************** > > ------------------------------ Message: 6 Date: Fri, 01 Feb 2013 09:36:18 -0200 From: Leonard de Assis <assis.leonard(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Índice fora de limites Message-ID: <510BA8B2.4000401(a)gmail.com> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Pela quantidade de observações (192) e o número de respostas possíveis (7), e o número de itens (22),é bem provável que algum item seu não tenha algum dos níveis. será que SETE níveis são mesmo necessários? Em 30/01/2013 21:45, Thiago de Andrade escreveu: > Leonard, > Não postei o CMR porque não tenho como disponibilizar os dados. É uma > matriz 22 (itens) x192 (sujeitos), assumindo valores entre 1 e 7. Não > há dados omissos nem fora desse intervalo (1-7). Você acha que a baixa > frequência de alguns vaalores pode estar causando o erro? (desculpa > minha ignorância, sou só um psicólogo rsrs). > > Em 30 de janeiro de 2013 20:40, Leonard de Assis > <assis.leonard(a)gmail.com <mailto:assis.leonard@gmail.com>> escreveu: > > Thiago, > não lembro mais como é a codificação das variáveis pra rodar TRI > no LTM, mas o erro (vc não postou CMR, portanto é chute) tem toda > a cara de codificação errada em alguma variável. > > Em 30/01/2013 19:18, Thiago de Andrade escreveu: >> Amigos, estou tentando rodar uns dados de um amigo com o pacote >> ltm, (Teoria da Resposta ao Item). Estou recebendo uma mensagem >> de erro "índice fora de limites". Algum de vocês pode me dar >> alguma luz acerca da localização deste problema? >> >> Obrigado. >> >> -- >> Thiago de Andrade >> Mestrando em Psicologia Social, PPGP-UFPB >> Núcleo de Pesquisas em Desenvolvimento Sócio-Moral (NPDSM) >> (83) 9642 4099 >> >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br <mailto:R-br@listas.c3sl.ufpr.br> >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br <mailto:R-br@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. > > > > > -- > Thiago de Andrade > Mestrando em Psicologia Social, PPGP-UFPB > Núcleo de Pesquisas em Desenvolvimento Sócio-Moral (NPDSM) > (83) 9642 4099 > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130201/b4233143/atta…> ------------------------------ Message: 7 Date: Fri, 1 Feb 2013 09:37:30 -0200 From: "Luciano F. Sgarbi" <luciano.f.sgarbi(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Gráfico dupla entrada Message-ID: <CAJpKj8kaKxN-y3e2MgozV9A1=HRTOA7ZKbigfjwu6uZr_ZyY5Q(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" Olá Stephan, Dê uma olhada nesse site (http://gallery.r-enthusiasts.com/) Com certeza você vai encontrar. Abraços, Luciano Em 1 de fevereiro de 2013 12:37, Stephan M. CARVALHO <smalfitano(a)uol.com.br>escreveu: > Bom dia pessoal, > > > Preciso plotar dados em um gráfico de dupla entrada, como ocorre por > exemplo com dados climatológicos (temperatura/precipitação ~ tempo). > > Por acaso alguém conhece alguma forma de fazer? > > Obrigado, > > > Stephan > > > -- > > > ***************************************************** > > STEPHAN M. CARVALHO > > "Honeybee protection and ecotoxicology" > > Universidade Estadual Paulista > Centro de Estudo de Insetos Sociais > > Av. 24A, 1515 > Bela Vista > 13.506-900 Rio Claro/SP > Brazil > > smalfitano(a)uol.com.br > > +55 19 3526-4284 (university) > +55 19 3534-8523 (fax) > +55 35 9103-9822 (mobile) > +55 35 9949-0707 (mobile) > > > Antes de imprimir, pense no ambiente! > Before printing, think in the environment! > Avant d'imprimer, pensez à l'environnement! > > ***************************************************** > > > _______________________________________________ > 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. > -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130201/5bf30109/atta…> ------------------------------ Message: 8 Date: Fri, 01 Feb 2013 09:42:18 -0200 From: Leonard de Assis <assis.leonard(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Informação: R como software livre... Message-ID: <510BAA1A.6030504(a)gmail.com> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" uso academico, o RevoR também é gratis Em 31/01/2013 17:51, Manoel Galdino escreveu: > Para todos os itens: > 1. Redução de custos ( o R é grátis, exceto pelo Revolutions) > 2. Por ter código aberto, é mais fácil checar a análise feita. > > Para pesquisa: > A comunidade de usuários dá tanto suporte que ajuda o pesquisador a > criar laços com outros pesquisadores > > abçs > Manoel > > 2013/1/31 <andrebvs(a)bol.com.br <mailto:andrebvs@bol.com.br>> > > Olá à todos! > > O programa R, sem dúvida, é a ferramenta estatística ideal nas > mãos de quem o conhece, mas, porque universidades ainda insistem > em utilizar programas comerciais (propriétarios) em seus > laboratórios? acredito que esse fato, ainda seja pela falta de > interesse de professores "tradicionalistas" que não instigam > mudanças (inovações) na forma de ensino. Talvez, falta a esses > professores, a informação sobre os beneficios que os softwares > livres trazem para a instituição e para o governo. Então, diante > disso, como justificar o uso do programa R como software livre em > cada uma das etapas a seguir dentro das Universidades: > > 1) ensino > 2) pesquisa > 3) extensão > desde já agradeço! > /Att./ > /André/ > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br <mailto:R-br@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. > > > > > -- > Manoel Galdino > https://sites.google.com/site/galdinomcz/ > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130201/cc509d0b/atta…> ------------------------------ Message: 9 Date: Fri, 1 Feb 2013 09:57:47 -0200 From: Manoel Galdino <mcz.fea(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Informação: R como software livre... Message-ID: <CA+kPArkb57RT3cNupTjuG3BMiAt_BCfEx6+ZWQN1HPm=AmBfXw(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" é verdade! abç M 2013/2/1 Leonard de Assis <assis.leonard(a)gmail.com> > uso academico, o RevoR também é gratis > Em 31/01/2013 17:51, Manoel Galdino escreveu: > > Para todos os itens: > 1. Redução de custos ( o R é grátis, exceto pelo Revolutions) > 2. Por ter código aberto, é mais fácil checar a análise feita. > > Para pesquisa: > A comunidade de usuários dá tanto suporte que ajuda o pesquisador a criar > laços com outros pesquisadores > > abçs > Manoel > > 2013/1/31 <andrebvs(a)bol.com.br> > >> Olá à todos! >> >> O programa R, sem dúvida, é a ferramenta estatística ideal nas mãos de >> quem o conhece, mas, porque universidades ainda insistem em utilizar >> programas comerciais (propriétarios) em seus laboratórios? acredito que >> esse fato, ainda seja pela falta de interesse de professores >> "tradicionalistas" que não instigam mudanças (inovações) na forma de >> ensino. Talvez, falta a esses professores, a informação sobre os beneficios >> que os softwares livres trazem para a instituição e para o governo. Então, >> diante disso, como justificar o uso do programa R como software livre em >> cada uma das etapas a seguir dentro das Universidades: >> >> 1) ensino >> 2) pesquisa >> 3) extensão >> >> desde já agradeço! >> >> *Att.* >> *André* >> >> _______________________________________________ >> 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. >> > > > > -- > Manoel Galdino > https://sites.google.com/site/galdinomcz/ > > > _______________________________________________ > R-br mailing listR-br@listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível. > > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -- Manoel Galdino https://sites.google.com/site/galdinomcz/ -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130201/ed0b3a6c/atta…> ------------------------------ Message: 10 Date: Fri, 01 Feb 2013 13:12:01 -0300 From: "Stephan M. CARVALHO" <smalfitano(a)uol.com.br> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Gráfico dupla entrada Message-ID: <510BE951.4060808(a)uol.com.br> Content-Type: text/plain; charset="iso-8859-1"; Format="flowed" Obrigado pessoal! Agora sim posso finalizar meu paper!!! Att, Stephan Em 1/2/2013 08:10, Paulo Justiniano escreveu: > O segredo está em usar par(new=T) npara justapor gráficos > O exemplo a seguir foi tirado dos arquivos da lista > http://r-br.2285057.n4.nabble.com/R-br-Grafico-com-dois-eixos-y-outra-soluc… > > > x <- 1:20 > x1 <- rnorm(20) > x2 <- rpois(20, lam=50) > plot(x, x1, type="l") > par(new=T) > plot(x, x2, type="l", axes=F, frame=T, ann=F, col=2) > axis(4) > > e um exemplo mais completo tb apontado em um email da lista: > http://r-br.2285057.n4.nabble.com/R-br-graficos-com-dois-eixos-y-td3566075.… > > sugere procurar em: > http://www.leg.ufpr.br/doku.php/ridiculas > por > Gráfico com dois eixos coordenados > > > > > > On Fri, 1 Feb 2013, Stephan M. CARVALHO wrote: > >> Bom dia pessoal, >> >> >> Preciso plotar dados em um gráfico de dupla entrada, como ocorre >> por exemplo com dados climatológicos (temperatura/precipitação ~ tempo). >> >> Por acaso alguém conhece alguma forma de fazer? >> >> Obrigado, >> >> >> Stephan >> >> >> -- >> >> >> *************************************************** >> >> STEPHAN M. CARVALHO >> >> "Honeybee protection and ecotoxicology" >> >> Universidade Estadual Paulista >> Centro de Estudo de Insetos Sociais >> >> Av. 24A, 1515 >> Bela Vista >> 13.506-900 Rio Claro/SP >> Brazil >> >> smalfitano(a)uol.com.br >> >> +55 19 3526-4284 (university) >> +55 19 3534-8523 (fax) >> +55 35 9103-9822 (mobile) >> +55 35 9949-0707 (mobile) >> >> >> Antes de imprimir, pense no ambiente! >> Before printing, think in the environment! >> Avant d'imprimer, pensez à l'environnement! >> >> *************************************************** >> >> > > > _______________________________________________ > 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. -- *************************************************** STEPHAN M. CARVALHO "Honeybee protection and ecotoxicology" Universidade Estadual Paulista Centro de Estudo de Insetos Sociais Av. 24A, 1515 Bela Vista 13.506-900 Rio Claro/SP Brazil smalfitano(a)uol.com.br +55 19 3526-4284 (university) +55 19 3534-8523 (fax) +55 35 9103-9822 (mobile) +55 35 9949-0707 (mobile) Antes de imprimir, pense no ambiente! Before printing, think in the environment! Avant d'imprimer, pensez à l'environnement! *************************************************** -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130201/28c32976/atta…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: smalfitano.vcf Tipo: text/x-vcard Tamanho: 352 bytes Descrição: não disponível URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130201/28c32976/atta…> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 24, assunto 1 ****************************************

2 1

mudar "names" das variaveis
by samuel luna de almeida 06 Fev '13

06 Fev '13

Bom dia, estou construindo um painel de correlação entre minhas variáveis, no entanto os nomes que as variáveis foram recebendo ao longo do tratamento dos dados, não são adequados para apresentação dos resultados... estou pedindo ajuda para trocar o nome das variaveis.. ex: a variável "TXJ04" gostaria de nomear como "taxa de internação em crianças", etc. alguém sabe como fazer? Abs Samuel

6 7

plot CCA
by Renata Carvalho 06 Fev '13

06 Fev '13

Bom dia, estou tentando fazer um plot de CCA mas não consigo colocar símbolos como referência aos sites. O R renomeou cada site para row1, row2, row3, Fui orientada a transformar cada site em número para depois atribuir os símbolos mas não estou conseguindo. Até agora tenho utilizado estes scripts: as.numeric (dados$CODE) plot(vare.cca, type="n") text(vare.cca, dis="cn") text(vare.cca, display = "sites", col="blue", cex=0.5) # aqui os sites aparecem no plot mas com o nome de row1, row2... Alguém sabe como posso fazer isso? obrigada Renata

2 4

Cointegration analysis in the presence of structural breaks in the deterministic trend
by Marcelo Justus dos Santos 05 Fev '13

05 Fev '13

Olá! Eu estou utilizando o package{vars} para uma análise de cointegração de Johansen. Mas, a minha especificação contém uma "intervention dummy" para controlar uma quebra estrutural em uma das séries. Eu preciso simular os novos valores críticos para os testes do traço e máximo valor com a correção proposta por Johansen et al. (2000). Referência:Søren Johansen & Rocco Mosconi & Bent Nielsen, 2000. Cointegration analysis in the presence of structural breaks in the deterministic trend, Econometrics Journal, vol. 3(2), pages 216-249.Alguém sabe se há algum pacote no R para isso?Ou se não existir, alguém sabe em qual pacote eu posso gerar uma aproximação da distribuição Gamma a partir da média e da variância para computar os percentis da distribuição?Muito obrigado,Marcelo Justus

4 3

trabalhar com ifelse
by alanarocha＠sapo.pt 05 Fev '13

05 Fev '13

Boa tarde, estou a fazer um trabalho onde tenho de usar ifs gostaria de saber se o meu if está bem construido, senão o que devo alterar. ficheiro: setwd("G:/SIN/Users/Ana Rua") objectivo<-c("verificar se a percentagem de Hgb de cada clínica é <70% em 3 meses consecutivos no ano 2012,no caso de serem aparece: 3 meses abaixo senão aparece ok") objectivo dados<-read.table("indicadores.csv",header=TRUE,sep=";") head(dados) names(dados) ifelse(Jan<70&Fev<70&Mar<70){ else{ cat"\n\se for menor que 70 em 3 meses consecutivos, devolve "3meses mais baixos" e imprime o nome da clinica"); else{ cat(\n\caso contrário devolve OK"); return(Clinica) resultados setwd("G:/SIN/Users/Ana Rua") objectivo<-c("verificar se a percentagem de Hgb de cada clínica é <70% em 3 meses consecutivos no ano 2012,no caso de serem aparece: 3 meses abaixo senão aparece ok") objectivo [1] "verificar se a percentagem de Hgb de cada clínica é <70% em 3 meses consecutivos no ano 2012,no caso de serem aparece: 3 meses abaixo senão aparece ok" dados<-read.table("indicadores.csv",header=TRUE,sep=";") head(dados) Clinica Jan Fev Mar Abr Mai Jun Jul Ago Set Out 1 Ponte da Barca 76.1 80.0% 82.8% 83.3% 84.6% 79.4% 84.1% 86.5% 73.1% 63.2% 2 Braga 67.3% 74.0% 72.7% 76.8% 70.7% 79.0% 70.4% 75.5% 72.5% 72.2% 3 Fafe 68.24% 69.1% 73.2% 79.0% 72.6% 84.5% 82.1% 77.4% 79.8% 75.3% 4 Maia 77.4% 80.9% 76.8% 76.0% 76.6% 80.8% 81.1% 81.4% 80.6% 82.6% 5 VNGaia 82.1% 80.7% 82.7% 84.4% 82.0% 78.9% 81.8% 80.9% 77.8% 76.5% 6 SMFeira 78.22% 78.2% 77.0% 80.6% 81.8% 77.6% 63.6% 67.7% 74.2% 71.4% Nov Dez 1 66.7% 68.7% 2 75.5% 73.9% 3 79.3% 76.7% 4 83.3% 84.9% 5 81.4% 80.4% 6 73.7% 76.0% names(dados) [1] "Clinica" "Jan" "Fev" "Mar" "Abr" "Mai" "Jun" [8] "Jul" "Ago" "Set" "Out" "Nov" "Dez" com os melhores cumprimentos Ana Rocha

3 2

Dúvida MiKtex
by Mauro Sznelwar 04 Fev '13

04 Fev '13

Tenho visto discussão sobre MiKtex, que não sei do que se trata. Porém estive trabalhando com o MATLAB, e deu um problema, e o autor do código me deu a resposta que o problema é com a função cgf, que é a função escrita em C e conectada com o MATLAB usando o framework MEX. Tem alguma coisa a ver com a discussão sobre o MiKtex do R?

2 2

eliminar linhas específicas de um data.frame
by Aline Lipsky 01 Fev '13

01 Fev '13

Boa tarde, Tenho um data.frame similar a esse com mais de 6000 obs.: data local comprador 29/02/2012 RS 0 29/02/2012 RS 0 29/02/2012 RS 0 29/02/2012 RS 0 29/02/2012 RS 1 29/02/2012 RS 1 29/02/2012 RS 1 27/02/2012 RS 0 27/02/2012 RS 1 27/02/2012 RS 0 27/02/2012 RS 0 27/02/2012 RS 1 27/02/2012 RS 1 27/01/2012 RS 0 28/01/2012 RS 0 27/01/2012 RS 0 26/01/2012 RS 0 26/01/2012 RS 0 Preciso eliminar as linhas onde a variável comprador é 1. Me interessam apenas os dados do comprador 0. Agradeço desde já a atenção, Att., Aline

3 2

profajabreu@gmail.com wants to follow you. Accept?
by profajabreu＠gmail.com 01 Fev '13

01 Fev '13

Hi, Ademir José wants to follow you. ****** Is Ademir José you friend? ****** If Yes please follow the link below: http://invites.infoaxe.net/signup_e.html?fullname=R-br&email=r-br@listas.c3… If No please follow the link below: http://invites.infoaxe.net/signup_e_no.html?fullname=R-br&email=r-br@listas… Follow the link below to unsubscribe from all emails from Flipora, 440 N.Wolfe Rd MS #153, Sunnyvale, CA. 94085 http://invites.infoaxe.net/unsubscribe.jsp?email=r-br@listas.c3sl.ufpr.br&i…

1 0

Índice fora de limites
by Thiago de Andrade 01 Fev '13

01 Fev '13

Amigos, estou tentando rodar uns dados de um amigo com o pacote ltm, (Teoria da Resposta ao Item). Estou recebendo uma mensagem de erro "índice fora de limites". Algum de vocês pode me dar alguma luz acerca da localização deste problema? Obrigado. -- Thiago de Andrade Mestrando em Psicologia Social, PPGP-UFPB Núcleo de Pesquisas em Desenvolvimento Sócio-Moral (NPDSM) (83) 9642 4099

3 8

Gráfico dupla entrada
by Stephan M. CARVALHO 01 Fev '13

01 Fev '13

Bom dia pessoal, Preciso plotar dados em um gráfico de dupla entrada, como ocorre por exemplo com dados climatológicos (temperatura/precipitação ~ tempo). Por acaso alguém conhece alguma forma de fazer? Obrigado, Stephan -- *************************************************** STEPHAN M. CARVALHO "Honeybee protection and ecotoxicology" Universidade Estadual Paulista Centro de Estudo de Insetos Sociais Av. 24A, 1515 Bela Vista 13.506-900 Rio Claro/SP Brazil smalfitano(a)uol.com.br +55 19 3526-4284 (university) +55 19 3534-8523 (fax) +55 35 9103-9822 (mobile) +55 35 9949-0707 (mobile) Antes de imprimir, pense no ambiente! Before printing, think in the environment! Avant d'imprimer, pensez à l'environnement! ***************************************************

4 4

Informação: R como software livre...
by andrebvs＠bol.com.br 01 Fev '13

01 Fev '13

3 3

Nayara Santos is out of the office.
by Nayara Santos 31 Jan '13

31 Jan '13

I will be out of the office starting 30/01/2013 and will not return until 14/02/2013. Gentileza entrar em contato com: Flávia Novaes para testes. Marcelo Prieto e Vicente Venancio para definições, dúvidas e aprovações.

1 0

Gráfico no ggplot com escala x variável
by Alisson Lucrecio 31 Jan '13

31 Jan '13

Caros colegas, Boa noite. Classifiquei um conjunto de dados de poros por diâmetro em diferentes classes levando em consideração o numero de poros em cada classe. Os resultados me geraram outro conjunto de dados com o numero de poros em cada classe, sendo 30 classes de 0.1 a 1 mm e 1 classes de 1 a 5 mm. Preciso plotar um gráfico no ggplot2 com o diâmetro no x e numero de poros em cada classe no y, mas quando plot gera um gráfico da figura em anexo. Não consigo ter uma boa visualização dos resultado nesse plot. Então preciso fazer um corte no eixo x ou tentar uma escala x automática de 0.1 a 1 mm e escala x fixa de 1 a 5 mm (fixo no x igual a 3 mm). Alguém sabe como fazer isso? Obrigado. Alisson Lucrécio da Costa

2 2

Ajuda com o pacote randomLCA
by Victor Eduardo 30 Jan '13

30 Jan '13

Pessoal, alguém poderia me ajudar def.lcarandom <- randomLCA(definicao[,1:4], freq=definicao$x, initmodel=def.lca, random=TRUE, blocksize=1, quadpoints=190, probit=TRUE) > summary(def.lcarandom) Classes AIC BIC logLik penlogLik Link 2 3842.735 3896.336 -1911.368 -1911.373 Probit Class probabilities Class 1 Class 2 0.8213 0.1787 Conditional outcome probabilities Group.1 Group.2 Group.3 Group.4 Class 1 0.0000 0.0065 0 0.0011 Class 2 0.9612 1.0000 1 1.0000 Marginal Outcome Probabilities Group.1 Group.2 Group.3 Group.4 Class 1 0.0000 0.1987 0 0.1479 Class 2 0.7253 0.9545 1 1.0000 Loadings 2.76629> def.lcarandom$outcomep [,1] [,2] [,3] [,4] [1,] 5.653499e-50 0.006483299 5.024727e-33 0.001059331 [2,] 9.612131e-01 0.999999659 1.000000e+00 1.000000000 A primeira parte postada acima funciona normalmente, sendo que quando vou utilizar a função outcome.probs ela não funciona, ela fica "rodando" o tempo inteiro e não sai mais disso. def.lcarandom.probs<-outcome.probs(def.lcarandom, boot=T, R=100) Abraços!

1 0

geoestatistica
by Natalia Martins 30 Jan '13

30 Jan '13

Prezados, estou tentando comparar a interpolação por krigagem e pelo inverso da distancia, sempre usei somente a krigagem e agora fui questionada e resolvi procurar e e estudar mais métodos..no entanto no exemplo da cran (a seguir), nao consegui entender de onde vem os dados "meuse.grid", uma vez que estes contem as distancias. Alguem poderia me ajudar? # Inverse distance interpolation with inverse distance power set to .5: # (kriging variants need a variogram model to be specified) data(meuse) data(meuse.grid) meuse.gstat <- gstat(id = "zinc", formula = zinc ~ 1, locations = ~ x + y, data = meuse, nmax = 7, set = list(idp = .5)) meuse.gstat z <- predict(meuse.gstat, meuse.grid) library(lattice) # for levelplot levelplot(zinc.pred~x+y, z, aspect = mapasp(z)) -- Att, Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

3 2

rodar função para cada nível de um fator
by marcelo costa 30 Jan '13

30 Jan '13

Prezados membros, Peço a ajuda de vocês para o seguinte problema: Eu tenho uma planilha com dois fatores (“espécie” e “sexo”) e uma variável contínua (“tamanho”). Quero analisar se existe diferença no tamanho de machos e fêmeas para cada espécie. Assim, preciso fazer um teste T para cada espécie. O problema é que eu tenho 288 espécies (uma planilha de 6638 linhas) e, portanto, fazer um subset para cada espécie não me parece uma solução muito boa (e provavelmente não é!). Eu queria fazer um teste T para cada nível do fator "especie".Alguém tem alguma ideia de como eu poderia fazer isso sem precisar fazer um subset para cada espécie? A princípio eu achei que isso seria fácil (e provavelmente seja), no entanto o meu pouco conhecimento sobre o R não está me ajudando. Eu estou enviando um exemplo hipotético do formato da minha planilha: especie<-rep(c("a","b","c"),each=6) sexo<-rep(c("f","m"), each = 3, len = 18) tamanho<-c(1.5,1.6,1.8,1.3,1.1,1.2,2.5,2.7,2.6,2.1,2.2,2.1,3.5,3.6,3.7,3.1,3.3,3.2) data<-data.frame(especie,sexo,tamanho) Um abraço a todos e obrigado pela atenção. Marcelo Costa

5 6

Summary com todos os tipos da variável
by Fátima Lima Paula 28 Jan '13

28 Jan '13

Prezados, tenho uma variável (diagnóstico principal) que apresenta muitos tipos diferentes. Quando peço summary dessa variável, retorna dessa forma: summary(reinter$diagprinc) S720 S721 S722 T845 J189 T931 A419 N390 I500 I64 493 300 77 40 34 29 28 27 21 20 T840 J159 J158 M861 R02 S730 T814 I792 A499 I802 20 19 16 14 14 14 14 13 12 12 M966 R262 A408 L89 T813 T818 J188 J960 M840 T981 12 12 11 11 11 11 10 10 10 10 A46 I739 L989 M009 S723 T018 Y881 A418 A480 E46 9 9 9 9 9 9 9 7 7 7 E86 S525 A488 E43 G459 I509 I702 K921 K928 M171 7 7 6 6 6 6 6 6 6 6 M868 T841 T847 T857 Z470 A490 E145 I269 K566 R69 6 6 6 6 6 5 5 5 5 5 T848 Z988 A415 A498 E878 G319 I10 I743 J180 J449 5 5 4 4 4 4 4 4 4 4 K920 L97 M841 M869 N179 N189 S818 T888 E148 G20 4 4 4 4 4 4 4 4 3 3 I110 I200 I608 I635 I809 J80 K818 L031 L038 M169 3 3 3 3 3 3 3 3 3 3 M842 M848 M860 N111 N180 R33 S065 S423 S424 (Other) 3 3 3 3 3 3 3 3 3 218 Como faço para que apareçam todas? Obrigada Fátima

2 2

Re: [R-br] Método de Agrupamento Ward (bruno vidigal)
by Marcos Vital 28 Jan '13

28 Jan '13

Bruno, na verdade o método de Ward precisa sim da matriz de distância como os outros métodos de clsuter; o que muda é a forma como ele cria o dendograma. Ele calcula a soma dos quadrados das distâncias entre os objetos, dividida pelo número de objetos, buscando minimizá-la. Então é necessária uma matriz de distância como nos outros casos, e o que muda na função hclust é o argumento method="ward". Se tiver acesso, veja o excelente *Numerical Ecology with R* (Bocard, Gillet & Legendre, 2011), que é uma excelente fonte para estudar métodos multivariados no R (o foco é em Ecologia, mas com certeza pode ser útil em outros casos). Abraços Marcos -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

3 2

Método de Agrupamento Ward
by bruno vidigal 26 Jan '13

26 Jan '13

Boa noite, pelo o que eu li sobre o método de Ward, este não precisa calcular a matriz de distância para efetuar o método de agrupamento já que trabalha o tempo todo com os dados reais calculando soma de quadrados, diferente de outros métodos hierárquicos como o Vizinho mais próximo que, uma vez calculada a matriz de distâncias d, não utiliza mais os dados, pois fica atualizando-a. No R, usando a função hclust ele pede pra colocar a matriz de distância d. Daí queria saber como o R faz esse agrupamento. Valeu -- Bruno C. Vidigal

1 0

MÉTODO DE GAUSS-NEWTON
by Tiago Souza Marçal 26 Jan '13

26 Jan '13

Boa tarde colegas, Estou tentando reproduzir o algoritmo de GAUSS-NEWTOM em linguagem r, entretanto não estou conseguindo automatizar os comandos. Na parte da automatização eu estou seguindo o script que o Walmes postou no site <http://ridiculas.wordpress.com/?s=newton> sobre o metodo de NEWTON-RAPHSON. A titilo de exemplo estou utilizando um ajuste para o modelo y=a*exp(b*x) conforme segue a baixo: CMR: x<-c(2,5,7,10,14,19,26,31,34,38,45,52,53,60,65) y<-c(54,50,45,37,35,25,20,16,18,13,8,11,8,4,6) a<-56.6646 #parâmetro obtido obtido pela linearização da função y=a*exp(b*x) b<-(-0.03797) #parâmetro obtido obtido pela linearização da função y=a*exp(b*x) # a função fx fornece as estimativas de (y~x) dados os coeficientes a e b. fx<-function(xi,yi,a,b,i,n){resp<-numeric(0)for(i in i:n){t<-(a*exp(b*xi[i]))resp<-c(resp,t) }resp} # a função dfa fornece as estimativas da primeira derivada (dfx/da ~ x) dado o coeficientes b. dfa<-function(xi,yi,b,i,n){resp<-numeric(0)for(i in i:n){ t<-exp(b*xi[i]) resp<-c(resp,t) } resp } # a função dfb fornece as estimativas da primeira derivada (dfx/db~x) dados os coeficientes a e b. dfb<-function(xi,yi,a,b,i,n){resp1<-numeric(0)for(i in i:n){ t1<-a*xi[i]*exp(b*xi[i]) resp1<-c(resp1,t1) } resp1 } # A matrix betas contem os chutes para a e b betas<-matrix(c(a,b),nc=1) # bcor fornece os parâmetros a e b corrigidos para a primeira interação. Estes parâmetros serão substituídos no lugar dos valores de a e b iniciais e o processo sera repetido. bcor<-(betas)+((solve((t(matrix(c((dfa(dia,diag,b,1,15)),(dfb(dia,diag,a,b,1,15))),nc=2))%*%(matrix(c((dfa(dia,diag,b,1,15)),(dfb(dia,diag,a,b,1,15))),nc=2)))))%*%((t(matrix(c((dfa(dia,diag,b,1,15)),(dfb(dia,diag,a,b,1,15))),nc=2))%*%(matrix(c(diag-(fx(dia,diag,a,b,1,15))),nc=1))))) # a função Sqerro fornece a soma de quadrados dos desvios Sqerro<-sum((diag-fx(dia,diag,a,b,1,15))^2) O problema é que eu não estou conseguindo usar os parâmetros corrigidos para fazer a nova interação: Ex: bcor[i+1]<-bcor[i]+betas[i] # este comando aproximaria mais vezes os parâmetros. dif(Sqerro[i+1]-Sqerro[i]) #Este seria o critério de parada quando a variação no Sqerro fosse insignificante (10^(-6) por exemplo) Agradeço desde já pela ajuda e peço desculpas se não fui muito claro. Att. Tiago.

2 5

Como não plotar níveis não usados de um fator com bwplot - pacote LATICCE
by Augusto Ribas 25 Jan '13

25 Jan '13

Ola, estou com essa duvida. Eu tenho a seguinte situação ########################################## raca<-factor(c(rep(c("Controle"),each=10),rep(c("Labrador","Outros"),each=30)),) sint<-factor(c(rep(c("Controle"),each=10),rep(rep(c("ASSI","OLIGO","POLI"), each=10),2)),levels=c("Controle","ASSI","OLIGO","POLI")) medida<-rnorm(length(sint),5,1) library(lattice) bwplot(medida~sint|raca,layout = c(3, 1),drop.unused.levels=T) ############################################### O que eu gostaria era que o bwplot fize-se o boxplot para a raça controle, somente 1 boxplot chamado controle (nivel controle do fator sint), para a raça labrador e outros, eu gostaria que aparecessem os níveis ASSI OLIGO e POLI do fator sint, mas que nesses 2 (labrador e outros) o nivel controle do sint sumi-se. Eu tentei usar os argumentos drop.unused.levels e subset, mas não conseguir o resultado que gostaria, intuitivamente isso seria o que o drop.unused.levels deveria fazer não?. Existe alguma forma simples de conseguir isso usando o bwplot, que os niveis de fatores que não possuem dados não sejam plotados e não ocupem espaço dentro do quadradinho do grafico, sumam ali do eixo X nos quadradinhos que não tem dados? -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

1 1

Rcpp
by ctrucios 25 Jan '13

25 Jan '13

Pessoal, estou utilizando o pacote Rcpp para fazer um link entre C++ e R. Já consegui passar algumas funções (na verdade alguns loop) que tinha no meu codigo em R para C++, quando testo as funções que fiz (fiz 3) rodam certoinho, mas cuando coloco as 3 funções no programa e tento rodar meu programa todo o R trava. Alguem pode me ajudar? library(Rcpp) library(rugarch) nbur = 500 nobs = 1000 N = ntot = nbur + nobs data = rnorm(nobs) eboot = rnorm(nobs) n.ahead = 5 d = 2;p = 2,q= 2 garchorder = c(p,q) aux3 = cumsum(garchorder) Coeff_b = c(0.01,0.1,0.05,0.4,0.35) aux1 = Coeff_b[1]/(1-sum(Coeff_b[2:(p+1)])-sum(Coeff_b[(p+2):(aux3[2]+1)])) n.boot = 10 garchboot(data, n.boot, n.ahead, garchOrder = c(1,1), distribution.model = "norm", solver = "solnp", fit.control = list(stationarity = 1), n.burn = 500, simulation = list(simulation = FALSE, skew = NULL, shape = NULL), quant = c(0.90,0.95,0.99)) Para que este email não fique muito cumprido as funções que as coloquei no link aqui: https://www.dropbox.com/sh/5f8bpx6pt2cgam6/ELfL6BhGim/Rcpp Quem puder me dar um help agradecerei muito. *Carlos Trucíos Maza ** : ctrucios(a)gmail.com

2 2

verificar dados
by alanarocha＠sapo.pt 25 Jan '13

25 Jan '13

Olá rodrigo eu ja tinha enviado um excerto dos dados o que eu pretendo é saber quais são as clinicas que tem percentagem >=70% em cada trimestre. Eu não consigo enviar os dados num ficheiro em anexo para aqui mas posso enviar para seu mail se mo disser. Cumprimentos Ana Rocha

1 0

verificar dados
by alanarocha＠sapo.pt 25 Jan '13

25 Jan '13

Olá rodrigo eu ja tinha enviado um excerto dos dados o que eu pretendo é saber quais são as clinicas que tem percentagem >=70% em cada trimestre. Eu não consigo enviar os dados num ficheiro em anexo para aqui mas posso enviar para seu mail se mo disser. Cumprimentos Ana Rocha

1 0

verificar dados
by alanarocha＠sapo.pt 25 Jan '13

25 Jan '13

Olá rodrigo eu ja tinha enviado um excerto dos dados o que eu pretendo é saber quais são as clinicas que tem percentagem >=70% em cada trimestre. Eu não consigo enviar os dados num ficheiro em anexo para aqui mas posso enviar para seu mail se mo disser. Cumprimentos Ana Rocha

1 0

Teste de Média Anova type II
by Fernando A. Souza 25 Jan '13

25 Jan '13

Caros amigos, Possuo a seguinte dúvida. Estou analisando uns dados sobre produção de calor em ovelhas. Quando realizei a análise de variância tipo II ela indicou haver efeito dos tratamentos no entanto a função glht (Multcomp) não indicou nenhuma diferença para o método single-step. O resultado da função glht entretanto esteve de acordo com o resultado da análise de variância tipo III. Em ambos os casos não houve efeito de interação entre os tratamentos. Minha dúvida é a seguinte, Porque isto acontece? há um meio de comparar as médias condizente com a ANOVA tipo II. att -- Fernando Antônio de Souza Zootecnista, Msc. Nutricão Animal, UFMG lattes: http://lattes.cnpq.br/6519538815038307

3 3

verificação de dados
by alanarocha＠sapo.pt 24 Jan '13

24 Jan '13

Bom dia, eu tenho um ficheiro excel com os nomes de clinicas e suas percentagens de hemoglobina em cada mes. queria uma orientação para fazer uma função que verificasse em cada trimestre se é maior que 70% e depois que me devolvesse como resultado o nome de todas as clinicas maiores que 70% e o trimestre correspondente. aqui deixo os dados: > = 70% 2012 Clinica Jan Fev Mar Abr Mai Jun Jul Ago Set Out Nov Dez Ponte da Barca 76,1% 80,0% 82,8% 83,3% 84,6% 79,4% 84,1% 86,6% 73,1% 63,2% 66,7% 68,7% Braga 67,3% 74,0% 72,7% 76,8% 70,7% 79,0% 70,4% 75,5% 72,5% 72,2% 75,5% 73,9% Fafe 68,2% 69,1% 73,2% 79,0% 72,6% 84,5% 82,1% 77,4% 79,8% 75,3% 79,3% 76,7% Maia 77,4% 80,9% 76,8% 76,0% 76,6% 80,8% 81,1% 81,4% 80,6% 82,6% 83,3% 84,9% VNGaia 82,1% 80,7% 82,7% 84,4% 82,0% 78,9% 81,8% 80,9% 77,8% 76,5% 81,4% 80,4% SMFeira 78,2% 81,4% 77,0% 80,6% 81,8% 77,6% 63,6% 67,7% 74,2% 71,4% 73,7% 76,0% Viseu 75,6% 69,6% 72,5% 74,3% 77,3% 77,8% 73,7% 79,0% 76,9% 77,8% 78,8% 75,4% Guarda 77,6% 81,4% 80,2% 80,2% 78,8% 75,3% 80,8% 85,0% 82,8% 86,7% 78,6% 80,0% Covilhã 78,5% 81,4% 78,2% 83,3% 84,2% 82,5% 81,1% 75,3% 75,0% 76,8% 80,0% 75,5% Coimbra 73,2% 80,4% 79,6% 76,1% 79,1% 77,6% 75,8% 72,3% 78,7% 76,1% 78,6% 79,1% Abrantes 76,0% 76,9% 79,0% 79,0% 81,8% 81,8% 76,4% 84,8% 82,4% 82,4% 82,9% 82,9% Entroncamento 72,2% 69,4% 67,6% 71,7% 75,0% 75,9% 74,8% 75,5% 79,7% 79,7% 72,0% 81,2% Santarém 83,9% 83,7% 84,1% 87,3% 82,9% 78,0% 82,4% 85,7% 89,9% 83,1% 78,4% 80,0% TVedras 69,8% 77,0% 81,4% 83,3% 81,9% 77,8% 81,5% 81,9% 85,4% 85,0% 80,7% 81,5% VFXira 77,1% 79,2% 79,1% 80,9% 84,7% 77,5% 73,8% 78,8% 79,5% 82,1% 79,1% 77,1% Alverca 71,6% 75,6% 85,9% 83,7% 82,3% 82,1% 80,9% 82,7% 78,0% 77,1% 78,9% 76,2% Amadora 73,9% 76,1% 71,9% 72,2% 79,3% 79,3% 74,2% 71,8% 72,1% 77,0% 73,2% 74,1% Lumiar 77,9% 76,6% 75,7% 80,1% 78,1% 73,3% 75,4% 72,5% 75,8% 75,9% 78,5% 74,1% Restelo 77,7% 72,1% 73,5% 76,2% 72,1% 72,5% 73,1% 67,0% 69,8% 69,8% 75,2% 86,5% CUF 74,7% 64,7% 62,5% 58,4% 64,8% 73,9% 79,3% 81,1% 75,6% 77,5% 70,7% 65,6% SAMS 87,8% 80,8% 80,3% 81,1% 87,1% 76,9% 77,9% 80,6% 83,1% 79,7% 72,6% 83,6% APDP 50,0% 50,0% 87,5% 62,5% 62,5% 33,3% 36,4% 63,6% 58,3% 66,7% 58,3% 57,1% Montijo 77,9% 90,3% 81,3% 88,9% 87,3% 80,3% 89,7% 84,1% 87,3% 84,2% 82,5% 81,8% Barreiro 75,1% 74,3% 76,5% 79,1% 80,6% 80,0% 78,5% 77,0% 72,7% 78,3% 81,3% 78,8% Almada 81,0% 76,0% 81,3% 71,8% 75,9% 83,1% 83,0% 79,0% 79,3% 76,4% 75,9% 82,1% Seixal 75,4% 79,8% 82,8% 81,9% 78,3% 81,9% 80,5% 80,8% 85,8% 82,6% 81,7% 83,9% Setúbal 73,7% 72,0% 66,2% 76,7% 77,6% 82,4% 78,3% 75,3% 77,2% 82,1% 80,8% 83,0% Portalegre 76,0% 77,6% 74,2% 80,0% 85,0% 73,7% 70,2% 77,6% 71,4% 79,2% 73,5% 71,7% Évora 81,5% 83,1% 85,1% 80,3% 83,1% 78,0% 77,5% 77,5% 83,9% 81,7% 83,1% 82,0% Grândola 81,7% 80,7% 80,0% 76,5% 75,3% 78,8% 80,0% 82,4% 87,2% 83,7% 73,9% 76,8% Beja 72,5% 69,5% 66,7% 66,7% 70,6% 60,6% 62,6% 66,9% 73,4% 75,2% 73,4% 70,9% Tavira 67,5% 80,5% 76,3% 70,3% 74,4% 81,6% 64,9% 65,8% 64,1% 61,5% 79,0% 71,1% Faro 80,5% 77,2% 68,5% 84,2% 78,7% 76,3% 79,0% 76,9% 76,8% 78,1% 83,9% 77,9% Portimão 78,7% 85,4% 81,4% 80,6% 76,4% 81,0% 78,7% 81,9% 82,7% 83,8% 83,0% 85,3% Funchal 71,5% 73,4% 72,9% 74,0% 71,3% 75,8% 68,0% 71,4% 75,2% 78,9% 80,5% 77,0% pensei em: if(jan>=70%) if(fev>=70%) if(março>=70% else caso contrario imprimir o nome da clinicas >=70% e imprimir falso se for menor depois ver quais as 3 clinicas com menor valor em cada trimestre de cada semestre. obrigada Cumprimentos Ana Rocha

4 3

verificação de dados
by alanarocha＠sapo.pt 24 Jan '13

24 Jan '13

eu tentei fazer assim: else if(dados$Fev>=0.7,dados$Mar>=0.7,dados$Mar>=0.7) return(POnte da Barca, Braga, Maia.... etc) }

1 0

GLMM com dados categóricos
by Pavel Dodonov 24 Jan '13

24 Jan '13

Saudações, É a primeira vez que escrevo para cá, perdoem-me caso eu quebre alguma norma... Eu preciso ver os efeitos de algumas variáveis (contínuas e categóricas) sobre uma variável categórica com quatro níveis. Além disso, os dados estão aninhados em dois grupos, então acho que precisa ser um modelo misto. Explicando melhor: É um experimento com predação de ninho. A variável-resposta é o tipo de predador ("Ave", "Mamífero", "Desconhecido", "Nenhum"), as variáveis explanatórias são: Local (dois níveis - "Borda" e "Interior"), altura (variável contínua), e densidade (variável contínua, limitada entre 0 e 1 mas que na prática não ultrapassa 0.3). Tem também o fator aleatório "Campanha" (dois níveis, "Outubro" e "Novembro"), que até poderia entrar como variável explanatória, mas que seria melhor usada como fator aleatório porque não estamos interessados na diferença entre as campanhas. Eu pensei em fazer quatro análises usando cada nível da variável-resposta como uma variável binária e fazer um GLMM com distribuição binomial, MAS seria melhor se desse pra considerar os quatro níveis de uma vez. Talvez um GLMM com distribuição multinomial? Se alguém tiver alguma sugestão, ficarei imensamente agradecido. Att - Pavel ------- Pavel Dodonov Biologist, PhD student in Ecology and Natural Resources, São Carlos Federal University (UFSCar), SP, Brazil http://ufscar-br.academia.edu/pdodonov *"The highest function of ecology is understanding consequences." - Frank Herbert, Dune, 1965*

1 0

Mensagem de erro: cannot allocate memory block of size 3.2
by Cláudio Brisolara 23 Jan '13

23 Jan '13

Caros colegas Estou precisando de ajuda para identificar e solucionar um problema relativo a alocação de memória no R. Muitas mensagens já foram postadas sobre esse assunto, mas não consegui localizar uma solução. Por essa razão, solicito a ajuda de alguém. Em algumas operações que faço recebo a seguinte mensagem: “Error: cannot allocate memory block of size 3.2 GbIn addition: Warning message:restarting interrupted promise evaluation” A rotina que estou tentando executar é a seguinte: mu1<-c(8.035153, 3.931724) # fixa vetor de médias sigma1<-matrix(c(0.014179, 0.003042, 0.003042, 0.026259), 2, 2) # fixa matriz de variâncias set.seed(23456) n<-20000 dados.sim1<-rmvnorm(n, mu1, sigma1, method=c("chol")) # gera vetor normais correlacionados Ocorre que a questão da memória também é reportado em outras situações como por exemplo na simples chamada de um pacote: require(tseries) Loading required package: tseries Error in get(Info[i, 1], envir = env) : cannot allocate memory block of size 2.7 Gb Eu usei os comandos memory.limit() and memory.size() para checar a limitação de memória, mas não observei nenhum problema. Encaminho também as informações obtidas pelo sessionInfo(). > memory.limit() [1] 6004 > memory.size() [1] 1361.88 > sessionInfo() R version 2.15.2 (2012-10-26) Platform: x86_64-w64-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C [5] LC_TIME=Portuguese_Brazil.1252 loaded via a namespace (and not attached): [1] grid_2.15.2 quadprog_1.5-4 stabledist_0.6-5 tools_2.15.2 [5] xtable_1.7-0 O mais intrigante é que o problema realmente não é de memória, pois executei totalmente o mesmo script em um computador bem mais velho, com menos memória. É um problema de configuração da máquina, arquivo perdido algo do gênero? Já reinstalei o R e o RStudio, mas não obtive êxito. Por favor, algum colega poderia me ajudar a entender o que está acontecendo? Antecipadamente agradeço. Atenciosamente, Cláudio Brisolara

1 0

Re: [R-br] Legenda no eixo X entre dois gráficos
by Humberto Hazin 23 Jan '13

23 Jan '13

Tito muito obrigado consegui!!! André segue o código correto e com as tuas perguntas (que antes não havia entendido) respondidas Abr Humberto pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis <- c(length(tableMat1a$or):1) plot(tableMat1a$or,y.axis, type = "p", axes = T, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs White",yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) #axis(1,at=seq(0,1.2,by=0.1),labels = T,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = "n") par(mar=c(5,5, 3, 5)) y.axis <- c(length(tableMat1b$or):1) plot(tableMat1b$or,y.axis, type = "p", axes = T, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs Blue",yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) #axis(1,at=seq(0,1.2,by=0.1),labels = T,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = "n") mtext("Relative Catchability",side=1,line=3,adj =-0.6) dev.off() Enviado do Email do Windows De: Humberto Hazin Enviado: ‎23‎ de ‎janeiro‎ de ‎2013 ‎16‎:‎13 Para: R_br, R_br Assunto: Re: [R-br] Legenda no eixo X entre dois gráficos Olá André e Tito, Eu realmente não entendi tua pergunta André! Eu vou tentar aqui a solução pela dica do Tito e se conseguir posto aqui o resultado abr Humberto Enviado do Email do Windows De: andrebvs(a)bol.com.br Enviado: ‎23‎ de ‎janeiro‎ de ‎2013 ‎13‎:‎46 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Legenda no eixo X entre dois gráficos O seu gráfico ainda não está igual ao que deseja, falta colocar os valores no eixo X e fechar a caixa (box) onde ficam os intervalos. Como vc vai fazer isso? Em 23/01/2013 09:10, Tito Conte < tito.conte(a)gmail.com > escreveu: De um ?mtext Lá vai ter tudo o que você precisa! Tito Conte On Thu, Jan 17, 2013 at 12:37 PM, Humberto Hazin <hghazin(a)hotmail.com> wrote: Bom dia pessoal, Como faço para colocar para colocar a legenda no eixo x, meio entre dois gráficos? Por exemplo gostaria de adicional o texto “Relative Catchability” entre os dois gráficos. Segue o link do gráfico resultado so script abaixo (https://www.dropbox.com/s/36tw2eabk3sc7lr/FigureOR_4.pdf) Abraços Humberto tableMat1a: Species lci or uci 1 Bigeye tuna 0.0999 0.17 0.280 2 Yellowfin tuna 0.0366 0.07 0.135 3 Albacore 0.0850 0.18 0.391 4 Swordfish 0.0939 0.13 0.191 5 Sailfish 0.1257 0.25 0.498 6 White marlim 0.1500 0.30 0.390 7 Blue marlim 0.1291 0.29 0.637 8 Blue shark 0.1280 0.19 0.287 9 Crocodile shark 0.1120 0.16 0.223 10 Pelagic stingray 0.0449 0.17 0.619 11 Turtle oliva 0.1230 0.23 0.296 12 Oceanic whitetip shark 0.1881 0.50 0.770 tableMat1b: Species lci or uci 1 Bigeye tuna 0.00651 0.02 0.0647 2 Yellowfin tuna 0.08973 0.16 0.2702 3 Albacore 0.01370 0.05 0.1867 4 Swordfish 0.07696 0.12 0.1709 5 Sailfish 0.16000 0.20 0.4000 6 White marlim 0.05645 0.14 0.3382 7 Blue marlim 0.01200 0.18 0.2200 8 Blue shark 0.14176 0.21 0.3183 9 Crocodile shark 0.02888 0.12 0.5212 10 Pelagic stingray 0.04961 0.19 0.6875 11 Turtle oliva 0.15200 0.20 0.3220 12 Oceanic whitetip shark 0.23250 0.52 0.7700 tableMat1a<-read.csv("tableMetaGW.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill =TRUE) tableMat1b<-read.csv("tableMetaGB.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) library(tools) pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis plot(tableMat1a$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs White",yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = "n") par(mar=c(5,5, 3, 5)) y.axis plot(tableMat1b$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs Blue",yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = "n") dev.off() Enviado do Email do Windows _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Re: [R-br] Legenda no eixo X entre dois gráficos
by Humberto Hazin 23 Jan '13

23 Jan '13

Olá André e Tito, Eu realmente não entendi tua pergunta André! Eu vou tentar aqui a solução pela dica do Tito e se conseguir posto aqui o resultado abr Humberto Enviado do Email do Windows De: andrebvs(a)bol.com.br Enviado: ‎23‎ de ‎janeiro‎ de ‎2013 ‎13‎:‎46 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] Legenda no eixo X entre dois gráficos O seu gráfico ainda não está igual ao que deseja, falta colocar os valores no eixo X e fechar a caixa (box) onde ficam os intervalos. Como vc vai fazer isso? Em 23/01/2013 09:10, Tito Conte < tito.conte(a)gmail.com > escreveu: De um ?mtext Lá vai ter tudo o que você precisa! Tito Conte On Thu, Jan 17, 2013 at 12:37 PM, Humberto Hazin <hghazin(a)hotmail.com> wrote: Bom dia pessoal, Como faço para colocar para colocar a legenda no eixo x, meio entre dois gráficos? Por exemplo gostaria de adicional o texto “Relative Catchability” entre os dois gráficos. Segue o link do gráfico resultado so script abaixo (https://www.dropbox.com/s/36tw2eabk3sc7lr/FigureOR_4.pdf) Abraços Humberto tableMat1a: Species lci or uci 1 Bigeye tuna 0.0999 0.17 0.280 2 Yellowfin tuna 0.0366 0.07 0.135 3 Albacore 0.0850 0.18 0.391 4 Swordfish 0.0939 0.13 0.191 5 Sailfish 0.1257 0.25 0.498 6 White marlim 0.1500 0.30 0.390 7 Blue marlim 0.1291 0.29 0.637 8 Blue shark 0.1280 0.19 0.287 9 Crocodile shark 0.1120 0.16 0.223 10 Pelagic stingray 0.0449 0.17 0.619 11 Turtle oliva 0.1230 0.23 0.296 12 Oceanic whitetip shark 0.1881 0.50 0.770 tableMat1b: Species lci or uci 1 Bigeye tuna 0.00651 0.02 0.0647 2 Yellowfin tuna 0.08973 0.16 0.2702 3 Albacore 0.01370 0.05 0.1867 4 Swordfish 0.07696 0.12 0.1709 5 Sailfish 0.16000 0.20 0.4000 6 White marlim 0.05645 0.14 0.3382 7 Blue marlim 0.01200 0.18 0.2200 8 Blue shark 0.14176 0.21 0.3183 9 Crocodile shark 0.02888 0.12 0.5212 10 Pelagic stingray 0.04961 0.19 0.6875 11 Turtle oliva 0.15200 0.20 0.3220 12 Oceanic whitetip shark 0.23250 0.52 0.7700 tableMat1a<-read.csv("tableMetaGW.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill =TRUE) tableMat1b<-read.csv("tableMetaGB.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) library(tools) pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis plot(tableMat1a$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs White",yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = "n") par(mar=c(5,5, 3, 5)) y.axis plot(tableMat1b$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs Blue",yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = "n") dev.off() Enviado do Email do Windows _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0

Legenda no eixo X entre dois gráficos
by Humberto Hazin 23 Jan '13

23 Jan '13

Bom dia pessoal, Como faço para colocar para colocar a legenda no eixo x, meio entre dois gráficos? Por exemplo gostaria de adicional o texto “Relative Catchability” entre os dois gráficos. Segue o link do gráfico resultado so script abaixo (https://www.dropbox.com/s/36tw2eabk3sc7lr/FigureOR_4.pdf) Abraços Humberto tableMat1a: Species lci or uci 1 Bigeye tuna 0.0999 0.17 0.280 2 Yellowfin tuna 0.0366 0.07 0.135 3 Albacore 0.0850 0.18 0.391 4 Swordfish 0.0939 0.13 0.191 5 Sailfish 0.1257 0.25 0.498 6 White marlim 0.1500 0.30 0.390 7 Blue marlim 0.1291 0.29 0.637 8 Blue shark 0.1280 0.19 0.287 9 Crocodile shark 0.1120 0.16 0.223 10 Pelagic stingray 0.0449 0.17 0.619 11 Turtle oliva 0.1230 0.23 0.296 12 Oceanic whitetip shark 0.1881 0.50 0.770 tableMat1b: Species lci or uci 1 Bigeye tuna 0.00651 0.02 0.0647 2 Yellowfin tuna 0.08973 0.16 0.2702 3 Albacore 0.01370 0.05 0.1867 4 Swordfish 0.07696 0.12 0.1709 5 Sailfish 0.16000 0.20 0.4000 6 White marlim 0.05645 0.14 0.3382 7 Blue marlim 0.01200 0.18 0.2200 8 Blue shark 0.14176 0.21 0.3183 9 Crocodile shark 0.02888 0.12 0.5212 10 Pelagic stingray 0.04961 0.19 0.6875 11 Turtle oliva 0.15200 0.20 0.3220 12 Oceanic whitetip shark 0.23250 0.52 0.7700 tableMat1a<-read.csv("tableMetaGW.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill =TRUE) tableMat1b<-read.csv("tableMetaGB.csv", header = TRUE, sep = ",", quote="\"", dec=".",fill = TRUE) library(tools) pdf(file='FigureOR_4.pdf', height=50, width=50, onefile=TRUE, family='Helvetica', paper='a4r', pointsize=12) par(mfrow=c(1,2)) par(mar=c(5,11, 3, 0)) y.axis <- c(length(tableMat1a$or):1) plot(tableMat1a$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs White",yaxt='n') segments(tableMat1a$lci, y.axis, tableMat1a$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) axis(2, at = y.axis, label = tableMat1a$Specie, las = 1, tick = F, mgp = c(2,.6,0),cex.axis = 1.2) segments(1,1,1,30,lty=1) box(bty = "n") par(mar=c(5,5, 3, 5)) y.axis <- c(length(tableMat1b$or):1) plot(tableMat1b$or,y.axis, type = "p", axes = F, xlab = "", ylab = "", pch = 19, cex = 1.2,xlim = c(0,1.2), xaxs = "r", main = "Green vs Blue",yaxt='n') segments(tableMat1b$lci, y.axis, tableMat1b$uci, y.axis, lwd = 1.5) axis(1,at=seq(0,1.2,by=0.1),labels = F,tick =T,cex.axis = 1.2,mgp = c(2,.7,0)) segments(1,1,1,30,lty=1) box(bty = "n") dev.off() Enviado do Email do Windows

3 2

vaga de emprego
by Manoel Galdino 23 Jan '13

23 Jan '13

Olá, a empresa onde trabalho tem uma vaga em aberto. Como nós usamos bastante o R aqui na empresa, pensei que essa lista poderia ser um lugar para achar um bom candidato. Segue a descrição da vaga. O trabalho é em São Paulo. Descrição do trabalho: A área de Analytics da Hyperativa <http://hyperativa.com.br/>trabalha com base de dados de clientes para desenvolver análises de perfil da base, modelos preditivos, plataformas de CRM e ajudar na tomada de decisão estratégica. As bases de dados são gerenciadas em SQL pela equipe de Banco de dados e o analista trabalha em conjunto com a equipe de banco de dados para extrair os dados e em seguida analisá-los. Análises típicas incluem análises descritivas, regressões linerares e modelos Bayesianos. Boa parte do tempo do candidato será gasto em operações de limpeza, agregações e transformações dos dados, bem como transferência de dados entre bases e/ou softwares (por exemplo, do SQL para R e do R para Excel). É comum que nós tenhamos que resolver problemas que não têm soluções prontas na empresa e é necessário aprender rapidamente novas ferramentas e/ou técnicas, e a verdade é que nós adoramos aprender novas técnicas e gostamos de desafios. Para nós é importante que, além de ser comprometido com o trabalho e com as entregas, você tenha uma vida com interesses além das típicas habilidades vendáveis no mercado. Gostamos de pessoas com bom humor, que curtam cachorro, joguem xadrez, andem de bicicleta, vão ao cinema, torçam pela desgraça do time rival e adorem tecnologia. Requisitos: - Superior completo – Preferencialmente com formação forte em métodos quantitativos - Inglês intermediário (Desejável leitura avançada) - Conhecimentos de Estatística (aplicada e teoria) - Office - Desejável Conhecimento de R - Conhecimentos de SQL – diferencial - Conhecimentos de programação em geral – Diferencial -- Manoel Galdino https://sites.google.com/site/galdinomcz/

4 5

Extrair diagonais de Matriz ao selecionar um elemento.
by Augusto Ribas 23 Jan '13

23 Jan '13

Bom dia a todos, estou com uma duvida aqui. Eu tenho varias matrizes dessa forma. > exemplo A B C D 1 0 1 0 0 2 0 0 2 0 3 1 0 0 3 4 0 2 0 0 Aqui ela esta 4x4, mas ela pode ser de qualquer forma, não necessariamente quadrada assim. O maior numero que eu tenho ali dentro é 3. Existe uma forma simples de obter a posição do 3? Eu estava tentando o seguinte: Obtinha o maior elemento max(exemplo) Depois olhava onde ele estava. > which(exemplo==3) [1] 15 Mas ele fala um número apenas, que ai tem que dividir o número de colunas, tem alguma função que me retornaria a posição, qual linha e qual coluna ele esta? Tipo [3,4], e que seja independente se a matrix for diferente quanto ao numero de linhas ou colunas? A segunda duvida é, tem alguma forma simples de extrair a diagonal dos números a partir do 3? Principalmente as posições dos elementos. No exemplo eu gostaria do seguinte: A diagonal do 3 assim: [1,2],[2,3],[3,4]. Que ai eu poderia obter quais as colunas que estão nessa diagonal, tipo B,C,D, e as linhas 1,2,3. Esse é o objetivo final, mas eu estou apanhando para obter um jeito simples de extrair a posição do ultimo elemento da diagonal ali na matriz, ou a posição dos elementos da diagonal em si. -- Grato Augusto C. A. Ribas Site Pessoal: http://recologia.wordpress.com/ Lattes: http://lattes.cnpq.br/7355685961127056

2 3

Comando que mantém a ordem natural dos dados
by Gisele Correia 23 Jan '13

23 Jan '13

Boa tarde, Gostaria de saber como plotar dados mantendo a ordem das variáveis. O R está ordenando automaticamente. exemplo: x <- c(2,1,4,2) plot(x,y) (aqui o x é mostrado na ordem 1,2,3 e 4 no gráfico). Desde já agradeço, obrigada!! -- Gisele Correia Engenheira de Pesca PPG-CIPET/UFAM "Agradeço por todas as dificuldades que enfrentei (e ainda enfrento)... se não fosse por elas, não teria saído do lugar. A facilidade nos impede de caminhar. Mesmo as críticas nos auxiliam muito..." Chico Xavier

4 4

Problema com TL 50 em análise de sobrevivência
by ASANTOS 22 Jan '13

22 Jan '13

Boa tarde pessoal, Estou com um problema para calcular a TL50 (tempo de morte de 50% do indivíduos de uma população), usando o modelo de sobrevivência paramétrico de Weibull, pois o valor extraído do summary do modelo com exp(Intercept) não esta batendo com o curve do modelo, para explicar melhor segue CRM: require(survival) ##Dados obtidos na web A <- read.table("http://www.mas.ncl.ac.uk/~nmf16/teaching/mas3311/A.dat") attach(A) ## Ajuste do modelo mod = survreg(Surv(time, status)~1,dist="weibull") summary(mod) ########################### Calculando os parâmetros alpha e tempo médio ########## Valor de Alpha (1/Scale) 1/0.891 # =1.122334 ##TL 50 exp do (Intercept) do modelo exp(3.887) #= 48.76437 ##Plotando par(cex.lab=1.8,bty="n",yaxs="r",xaxs="r",lwd=2,mar=c(4.5,4.5,1,1),mex=1.3,family="") curve(exp((-(48.76437)^(-1.122334))*(x^1.122334)),xaxt="n",yaxt="n",from=0,to=170,ylim=c(0,1),xlab="",ylab="") x=seq(0,170,10) y=seq(0,1,0.1) axis(1,x,lwd=2,cex.axis=1.5) axis(2,y,lwd=2,cex.axis=1.8) # abline(h=0.5, col="red") abline(v=48.76437, col="red") ## O que não bate o valor encontrado com a curva Alguém poderia de dar uma ajuda, Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Avenida dos Ramires, s/n Caixa Postal 744 Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

1 0

Dúvida sobre intervalos....
by andrebvs＠bol.com.br 22 Jan '13

22 Jan '13

4 4