Dezembro 2013 - R-br

Summary GLM
by ana paula coelho madeira 16 Dez '13

16 Dez '13

Boa tarde pessoal. Tenho dados de contagem de ácaros em função do tempo(8 tempos) e da cultivar(4 cultivares). Quando declaro o modelo NP1 <- glm (Y ~ T+C, family = poisson, data=dados) tenho a seguinte saída: Deviance Residuals: Min 1Q Median 3Q Max -5.7588 -1.0631 -0.0002 0.7864 4.7855 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -1.718e+01 9.994e+02 -0.017 0.98628 T15 3.932e-09 1.413e+03 0.000 1.00000 T29 1.951e+01 9.994e+02 0.020 0.98442 T43 2.088e+01 9.994e+02 0.021 0.98333 T57 2.139e+01 9.994e+02 0.021 0.98292 T71 2.140e+01 9.994e+02 0.021 0.98291 T85 2.081e+01 9.994e+02 0.021 0.98339 T99 1.745e+01 9.994e+02 0.017 0.98607 C2 -1.468e-01 5.659e-02 -2.594 0.00950 ** C3 -1.433e-01 5.653e-02 -2.535 0.01123 * C4 -1.800e-01 5.710e-02 -3.152 0.00162 ** Pegunto: Não aparece o efeito da cultivar 1 (C1) e do tempo 0 (T0), por que isso ocorre? Como falar sobre o efeito desses dois fatores? Desde já agradeço. Ana Paula

2 1

Re: [R-br] Digest R-br, volume 36, assunto 15
by Gilenio Borges Fernandes 16 Dez '13

16 Dez '13

Boxplot e um histograma juntos O Pacote que substituiu o simple é UsingR. Abs Gilenio Fernandes Em 15 de dezembro de 2013 12:00, <r-br-request(a)listas.c3sl.ufpr.br>escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Ajuda com gráficos no R - Diagrama volumétrico dos espaços no > solo (Wenceslau Teixeira) > 2. Histograma e boxplot juntos (Andre Oliveira) > 3. Re: Histograma e boxplot juntos (Tiago Souza Marçal) > 4. Re: Histograma e boxplot juntos (Andre Oliveira) > 5. Re: Ajuda com gráficos no R - Diagrama volumétrico dos > espaços no solo (Benilton Carvalho) > 6. Gráfico de regressão linear (Andre Oliveira) > 7. Re: Gráfico de regressão linear (Rodrigo Coster) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sat, 14 Dec 2013 12:31:22 -0200 > From: Wenceslau Teixeira <wgt007(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Ajuda com gráficos no R - Diagrama volumétrico dos > espaços no solo > Message-ID: > <CACqRER=2kyBPdwUcX0V0QzTa=2c1WXPx= > 4FMY6sNT08qyLB2Zw(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Caro Walnes, > > > Agradeço a atenção,excelente, esta isto mesmo que precisávamos. > Agora só me falta descobrir como fazer uma legenda usando os "símbolos" > usados para diferenciar os poligonos. Qualquer dica é bem vinda > > Sds gráficas > > Wenceslau > > > Message: 11 > Date: Wed, 11 Dec 2013 19:56:22 -0200 > From: "walmes ." <walmeszeviani(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Ajuda com gráficos no R - Diagrama volumétrico dos > espaços no solo > Message-ID: > <CAFU=EkY= > nDFvVPz3vAc1txSYxDURzv1xCPodDjLHNpXDZ9kU7Q(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > A partir dessa tabela que você forneceu, você pode criar as listas que > demarcam os vértices dos polígonos que representam cada porção do solo. > Algo como > > plot(x=NULL, y=NULL, xlim=c(0,1), ylim=c(1.2,0), > xlab="Umidade (%)", ylab="Profundidade (m)") > > ## points(expand.grid(seq(0,1,0. > 1), seq(0,1.2,0.2))) > ## l <- locator() > ## lines(l) > ## l <- lapply(l, round, digits=1) > ## dput(l) > > l1 <- list(x = c(0, 0.2, 0.2, 0.4, 0.4, 0.5, 0.5, 0), > y = c(1.2, 1.2, 0.8, 0.8, 0.4, 0.4, 0, 0)) > l2 <- list(x = c(0.2, 0.7, 0.7, 0.6, 0.6, 0.6, 0.5, > 0.5, 0.4, 0.4, 0.2, 0.2), > y = c(1.2, 1.2, 0.8, 0.8, 0.6, 0, 0, 0.4, > 0.4, 0.8, 0.8, 1.2)) > l3 <- list(x = c(0.7, 1, 1, 0.6, 0.6, 0.7), > y = c(1.2, 1.2, 0, 0, 0.8, 0.8)) > > polygon(l1$x, l1$y, col="gray90") > polygon(l2$x, l2$y, density=3, angle=33) > polygon(l3$x, l3$y, density=12, angle=55, col="blue") > > > À disposição. > Walmes. > > > -- > Wenceslau Teixeira > > wgt007(a)gmail.com > skype: wgt007 > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131214/10382342/atta… > > > > ------------------------------ > > Message: 2 > Date: Sat, 14 Dec 2013 10:15:27 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Histograma e boxplot juntos > Message-ID: > <1387044927.86401.YahooMailNeo(a)web125905.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Pessoal, > estava precisando de um pacote que plota o boxplot e um histograma > juntos. Dei uma olhando a internet achei, mas não deu certo parece que o > pacote Simple saiu dos repositórios. Tem um pacote que plota y~x e na > margem de x e y plota o histograma. Estou precisando desta duas coisas eu > sei que existe, mas não me lembro qual pacote faz estas dois gráficos. > > library("Simple") > data(emissions) > attach(emissions) > simple.scatterplot(perCapita,CO2) > title("GDP/capita vs. CO2 emissions 1999") > > > > obrigado > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131214/0dcd7061/atta… > > > > ------------------------------ > > Message: 3 > Date: Sat, 14 Dec 2013 23:22:03 +0300 > From: Tiago Souza Marçal <tiagosouzamarcal(a)hotmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Histograma e boxplot juntos > Message-ID: <BLU170-W38E483DCD7AC28E10DBD5D9DE0(a)phx.gbl> > Content-Type: text/plain; charset="utf-8" > > Você pode fazer isso usando o comando par() para dividir a janela gráfica > e depois usar as funções boxplot() e hist(). > > Att. > > Tiago. > > > ################################################################# Tiago de > Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação > Científica da área de Genética e Melhoramento de Plantas > ################################################################# > > Date: Sat, 14 Dec 2013 10:15:27 -0800 > From: andreolsouza(a)yahoo.com.br > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Histograma e boxplot juntos > > Pessoal,estava precisando de um pacote que plota o boxplot e um histograma > juntos. Dei uma olhando a internet achei, mas não deu certo parece que o > pacote Simple saiu dos repositórios. Tem um pacote que plota y~x e na > margem de x e y plota o histograma. Estou precisando desta duas coisas eu > sei que existe, mas não me lembro qual pacote faz > estas dois gráficos. > library("Simple")data(emissions)attach(emissions)simple.scatterplot(perCapita,CO2)title("GDP/capita > vs. CO2 emissions 1999") > > > obrigado > > > > _______________________________________________ > 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/20131214/7a59de80/atta… > > > > ------------------------------ > > Message: 4 > Date: Sat, 14 Dec 2013 14:07:26 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Histograma e boxplot juntos > Message-ID: > <1387058846.73250.YahooMailNeo(a)web125906.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > Eu fiz isto Tiago. Queria a outra função só para colocar em uma > apresentação sobre R apenas para fim didático! > > > > > > > > André Oliveira Souza > > > > Em Sábado, 14 de Dezembro de 2013 18:22, Tiago Souza Marçal < > tiagosouzamarcal(a)hotmail.com> escreveu: > > Você pode fazer isso usando o comando par() para dividir a janela gráfica > e depois usar as funções boxplot() e hist(). > > Att. > > Tiago. > > > > ################################################################# > > Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES > > Bolsista de Iniciação Científica da área de Genética e Melhoramento de > Plantas > > ################################################################# > > > > ________________________________ > Date: Sat, 14 Dec 2013 10:15:27 -0800 > From: andreolsouza(a)yahoo.com.br > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Histograma e boxplot juntos > > > Pessoal, > estava precisando de um pacote que plota o boxplot e um histograma > juntos. Dei uma olhando a internet achei, mas não deu certo parece que o > pacote Simple saiu dos repositórios. Tem um pacote que plota y~x e na > margem de x e y plota o histograma. Estou precisando desta duas coisas eu > sei que existe, mas não me lembro qual pacote faz estas dois gráficos. > > library("Simple") > data(emissions) > attach(emissions) > simple.scatterplot(perCapita,CO2) > title("GDP/capita vs. CO2 emissions 1999") > > > > obrigado > > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne?a > c?digo m?nimo reproduz?vel. > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça > código mínimo reproduzível. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131214/f1130aeb/atta… > > > > ------------------------------ > > Message: 5 > Date: Sat, 14 Dec 2013 21:16:16 -0200 > From: Benilton Carvalho <beniltoncarvalho(a)gmail.com> > To: r-br <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Ajuda com gráficos no R - Diagrama volumétrico dos > espaços no solo > Message-ID: > < > CAO-arWMqcX481cpDvopuU-sAj6y9UdwukvKu0BBneCLHariJJQ(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > ?legend > > usando os comandos q o Walmes te enviou, uma tentativa natural e': > > legend('right', c('txtA', 'txtB', 'txtC'), fill=c('gray90', 'black', > 'blue'), angle=c(0, 33, 55), density=c(100, 3, 12)) > > > Em 14 de dezembro de 2013 12:31, Wenceslau Teixeira <wgt007(a)gmail.com > >escreveu: > > > Caro Walnes, > > > > > > Agradeço a atenção,excelente, esta isto mesmo que precisávamos. > > Agora só me falta descobrir como fazer uma legenda usando os "símbolos" > > usados para diferenciar os poligonos. Qualquer dica é bem vinda > > > > Sds gráficas > > > > Wenceslau > > > > > > Message: 11 > > Date: Wed, 11 Dec 2013 19:56:22 -0200 > > From: "walmes ." <walmeszeviani(a)gmail.com> > > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > > Subject: Re: [R-br] Ajuda com gráficos no R - Diagrama volumétrico dos > > espaços no solo > > Message-ID: > > <CAFU=EkY= > > nDFvVPz3vAc1txSYxDURzv1xCPodDjLHNpXDZ9kU7Q(a)mail.gmail.com> > > Content-Type: text/plain; charset="iso-8859-1" > > > > A partir dessa tabela que você forneceu, você pode criar as listas que > > demarcam os vértices dos polígonos que representam cada porção do solo. > > Algo como > > > > plot(x=NULL, y=NULL, xlim=c(0,1), ylim=c(1.2,0), > > xlab="Umidade (%)", ylab="Profundidade (m)") > > > > ## points(expand.grid(seq(0,1,0. > > 1), seq(0,1.2,0.2))) > > ## l <- locator() > > ## lines(l) > > ## l <- lapply(l, round, digits=1) > > ## dput(l) > > > > l1 <- list(x = c(0, 0.2, 0.2, 0.4, 0.4, 0.5, 0.5, 0), > > y = c(1.2, 1.2, 0.8, 0.8, 0.4, 0.4, 0, 0)) > > l2 <- list(x = c(0.2, 0.7, 0.7, 0.6, 0.6, 0.6, 0.5, > > 0.5, 0.4, 0.4, 0.2, 0.2), > > y = c(1.2, 1.2, 0.8, 0.8, 0.6, 0, 0, 0.4, > > 0.4, 0.8, 0.8, 1.2)) > > l3 <- list(x = c(0.7, 1, 1, 0.6, 0.6, 0.7), > > y = c(1.2, 1.2, 0, 0, 0.8, 0.8)) > > > > polygon(l1$x, l1$y, col="gray90") > > polygon(l2$x, l2$y, density=3, angle=33) > > polygon(l3$x, l3$y, density=12, angle=55, col="blue") > > > > > > À disposição. > > Walmes. > > > > > > -- > > Wenceslau Teixeira > > > > wgt007(a)gmail.com > > skype: wgt007 > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > > _______________________________________________ > > 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/20131214/558be745/atta… > > > > ------------------------------ > > Message: 6 > Date: Sun, 15 Dec 2013 04:33:51 -0800 (PST) > From: Andre Oliveira <andreolsouza(a)yahoo.com.br> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Gráfico de regressão linear > Message-ID: > <1387110831.59894.YahooMailNeo(a)web125903.mail.ne1.yahoo.com> > Content-Type: text/plain; charset="iso-8859-1" > > Bom dia, > alguem da lsta sabe qual função do R faz este gráfico. > https://www.dropbox.com/s/iyfef9sy4tfv4r4/teste.png > > obrigado > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131215/39a0e6fa/atta… > > > > ------------------------------ > > Message: 7 > Date: Sun, 15 Dec 2013 11:08:40 -0200 > From: Rodrigo Coster <rcoster(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, Andre > Oliveira <andreolsouza(a)yahoo.com.br> > Subject: Re: [R-br] Gráfico de regressão linear > Message-ID: > < > CAKU4wov0PYNJXtYhCT-4HKikOHzJ0z7tPF3tnn76yoQUezruVA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Uma simples busca no google resulta em diversas soluções: > > https://www.google.com/search?q=scatter+plot+with+histogram+r< > https://www.google.com/search?q=scatter+plot+with+histogram+r&oq=scatter+pl… > > > > > O 2o link é mais amigável para quem não tem 'intimidade' com o ggplot2... > > > 2013/12/15 Andre Oliveira <andreolsouza(a)yahoo.com.br> > > > Bom dia, > > alguem da lsta sabe qual função do R faz este gráfico. > > https://www.dropbox.com/s/iyfef9sy4tfv4r4/teste.png > > > > obrigado > > > > > > > > > > _______________________________________________ > > 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/20131215/50469f52/atta… > > > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > > ------------------------------ > > Fim da Digest R-br, volume 36, assunto 15 > ***************************************** > -- Gilenio Borges Fernandes Professor Associado IV Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

1 0

Re: [R-br] Gráfico de regressão linear
by Andre Oliveira 15 Dez '13

15 Dez '13

require(UsingR) x=rnorm(999) y=x+rnorm(999) simple.scatterplot(x,y) André Oliveira Souza Em Domingo, 15 de Dezembro de 2013 15:50, Andre Oliveira <andreolsouza(a)yahoo.com.br> escreveu: André Oliveira Souza Em Domingo, 15 de Dezembro de 2013 11:35, Rodrigo Coster <rcoster(a)gmail.com> escreveu: Acho que tu nao entrou no link que eu falei... Enfim, segue o texto do link (http://www.r-bloggers.com/example-10-3-enhanced-scatterplot-with-marginal-h…) scatterhist <- function(x, y, xlab = "", ylab = "", plottitle="", xsize=1, cleanup=TRUE,...){ # save the old graphics settings-- they may be needed def.par <- par(no.readonly = TRUE) zones <- matrix(c(1,1,1, 0,5,0, 2,6,4, 0,3,0), ncol = 3, byrow = TRUE) layout(zones, widths=c(0.3,4,1), heights = c(1,3,10,.75)) # tuning to plot histograms nicely xhist <- hist(x, plot = FALSE) yhist <- hist(y, plot = FALSE) top <- max(c(xhist$counts, yhist$counts)) # for all three titles: # drop the axis titles and omit boxes, set up margins par(xaxt="n", yaxt="n",bty="n", mar = c(.3,2,.3,0) +.05) # fig 1 from the layout plot(x=1,y=1,type="n",ylim=c(-1,1), xlim=c(-1,1)) text(0,0,paste(plottitle), cex=2) # fig 2 plot(x=1,y=1,type="n",ylim=c(-1,1), xlim=c(-1,1)) text(0,0,paste(ylab), cex=1.5, srt=90) # fig 3 plot(x=1,y=1,type="n",ylim=c(-1,1), xlim=c(-1,1)) text(0,0,paste(xlab), cex=1.5) # fig 4, the first histogram, needs different margins # no margin on the left par(mar = c(2,0,1,1)) barplot(yhist$counts, axes = FALSE, xlim = c(0, top), space = 0, horiz = TRUE) # fig 5, other histogram needs no margin on the bottom par(mar = c(0,2,1,1)) barplot(xhist$counts, axes = FALSE, ylim = c(0, top), space = 0) # fig 6, finally, the scatterplot-- needs regular axes, different margins par(mar = c(2,2,.5,.5), xaxt="s", yaxt="s", bty="n") # this color allows traparency & overplotting-- useful if a lot of points plot(x, y , pch=19, col="#00000022", cex=xsize, ...) # reset the graphics, if desired if(cleanup) {par(def.par)} } scatterhist(rnorm(1000), rchisq(1000, 7)) 2013/12/15 Andre Oliveira <andreolsouza(a)yahoo.com.br> Já havia buscado em meus arquivos pessoais e na internet com construir este gráfico com uma rodada só de camandos. Algo deste tipo: > > >library("Simple") >data(emissions) >attach(emissions) >simple.scatterplot(perCapita,CO2) >title("GDP/capita vs. CO2 emissions 1999") > > >Só que não deu certo aqui, parece que esta biblioteca saiu dos repositórios. > > > > > > > >Em Domingo, 15 de Dezembro de 2013 11:08, Rodrigo Coster <rcoster(a)gmail.com> escreveu: > >Uma simples busca no google resulta em diversas soluções: > > >https://www.google.com/search?q=scatter+plot+with+histogram+r > > > > > >O 2o link é mais amigável para quem não tem 'intimidade' com o ggplot2... > > > >2013/12/15 Andre Oliveira <andreolsouza(a)yahoo.com.br> > >Bom dia, >>alguem da lsta sabe qual função do R faz este gráfico. >>https://www.dropbox.com/s/iyfef9sy4tfv4r4/teste.png >> >> obrigado >> >> >> >> >>_______________________________________________ >>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

Gráfico de regressão linear
by Andre Oliveira 15 Dez '13

15 Dez '13

Bom dia, alguem da lsta sabe qual função do R faz este gráfico. https://www.dropbox.com/s/iyfef9sy4tfv4r4/teste.png obrigado

2 1

Ajuda com gráficos no R - Diagrama volumétrico dos espaços no solo
by Wenceslau Teixeira 14 Dez '13

14 Dez '13

Caro Walnes, Agradeço a atenção,excelente, esta isto mesmo que precisávamos. Agora só me falta descobrir como fazer uma legenda usando os "símbolos" usados para diferenciar os poligonos. Qualquer dica é bem vinda Sds gráficas Wenceslau Message: 11 Date: Wed, 11 Dec 2013 19:56:22 -0200 From: "walmes ." <walmeszeviani(a)gmail.com> To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Ajuda com gráficos no R - Diagrama volumétrico dos espaços no solo Message-ID: <CAFU=EkY=nDFvVPz3vAc1txSYxDURzv1xCPodDjLHNpXDZ9kU7Q(a)mail.gmail.com> Content-Type: text/plain; charset="iso-8859-1" A partir dessa tabela que você forneceu, você pode criar as listas que demarcam os vértices dos polígonos que representam cada porção do solo. Algo como plot(x=NULL, y=NULL, xlim=c(0,1), ylim=c(1.2,0), xlab="Umidade (%)", ylab="Profundidade (m)") ## points(expand.grid(seq(0,1,0. 1), seq(0,1.2,0.2))) ## l <- locator() ## lines(l) ## l <- lapply(l, round, digits=1) ## dput(l) l1 <- list(x = c(0, 0.2, 0.2, 0.4, 0.4, 0.5, 0.5, 0), y = c(1.2, 1.2, 0.8, 0.8, 0.4, 0.4, 0, 0)) l2 <- list(x = c(0.2, 0.7, 0.7, 0.6, 0.6, 0.6, 0.5, 0.5, 0.4, 0.4, 0.2, 0.2), y = c(1.2, 1.2, 0.8, 0.8, 0.6, 0, 0, 0.4, 0.4, 0.8, 0.8, 1.2)) l3 <- list(x = c(0.7, 1, 1, 0.6, 0.6, 0.7), y = c(1.2, 1.2, 0, 0, 0.8, 0.8)) polygon(l1$x, l1$y, col="gray90") polygon(l2$x, l2$y, density=3, angle=33) polygon(l3$x, l3$y, density=12, angle=55, col="blue") À disposição. Walmes. -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

2 1

Histograma e boxplot juntos
by Andre Oliveira 14 Dez '13

14 Dez '13

Pessoal, estava precisando de um pacote que plota o boxplot e um histograma juntos. Dei uma olhando a internet achei, mas não deu certo parece que o pacote Simple saiu dos repositórios. Tem um pacote que plota y~x e na margem de x e y plota o histograma. Estou precisando desta duas coisas eu sei que existe, mas não me lembro qual pacote faz estas dois gráficos. library("Simple") data(emissions) attach(emissions) simple.scatterplot(perCapita,CO2) title("GDP/capita vs. CO2 emissions 1999") obrigado

2 2

Re: [R-br] Otimização (boost + simplex)
by Ivan Bezerra Allaman 13 Dez '13

13 Dez '13

Debug a função para que tenha um diagnóstico melhor. library(debug) mtrace(simplex) mtrace(simplex1) O erro está acontecendo dentro a função simplex1 que é interna da função simplex. O erro ocorre em um while. Então, só você mesmo pra saber o porque do erro. (s,f,p) Allaman

1 0

Otimização (boost + simplex)
by Rodolfo Timoteo da Silva 13 Dez '13

13 Dez '13

Senhores, ao tentar otimizar a partir do arquivo abaixo: /library(boot)// //z = c(1,1,1,1,1,1,1)// //res2 = matrix(c(1,0,0,1,1,1,1,1,1,0,0,1,1,1,1,1,1,0,0,1,1,1,1,1,1,0,0,1,1,1,1,1,1,0,0,0,1,1,1,1,1,0,0,0,1,1,1,1,1),7,7)// //simplex(a=z,A1=NULL,b1=NULL,A2=rbind(res2),b2=c(18,12,15,19,14,16,11),A3=NULL,b3=NULL,maxi=FALSE)// / Recebo a mensagem: /Erro em pivot(tableau, prow, pcol) : // // NAs não são permitidos em atribuições por subscritos// //Calls: simplex -> simplex1 -> pivot// //Execução interrompida// / A título de elucidação, abaixo o arquivo utilizado pelo gurobi (www.gurobi.com) ao qual realizo a tarefa normalmente. /Minimize// // x1 + x2 + x3 + x4 + x5 + x6 + x7// //Subject To// // x1 + x2 + x3 + x4 + x5 >= 18// // x2 + x3 + x4 + x5 + x6 >= 12// // x3 + x4 + x5 + x6 + x7 >= 15// // x4 + x5 + x6 + x7 + x1 >= 19// // x5 + x6 + x7 + x1 + x2 >= 14// // x6 + x7 + x1 + x2 + x3 >= 16// // x7 + x1 + x2 + x3 + x4 >= 11// //Bound// // x1 >= 0// // x2 >= 0// // x3 >= 0// // x4 >= 0// // x5 >= 0// // x6 >= 0// // x7 >= 0// //Integers// // x1 x2 x3 x4 x5 x6 x7// //End/ Alias, como garantir lá no R que os resultados sejam números inteiros apenas? Obrigado Att

1 0

[OFF] Banco de dados
by Andre Oliveira 12 Dez '13

12 Dez '13

Boa noite, estou para fazer uma apresentação sobre análise descritiva de dados no R e não tenho um banco de dados com descrição das variáveis. Algo semelhante aquele conjunto de dados milsa que circula na internet. Alguém teria algum site para indicar onde eu possa conseguir alguns conjuntos de dados? obrigado

5 5

Re: [R-br] Ajuda com gráficos no R - Diagrama volumétrico dos espaços no solo
by Wenceslau Teixeira 12 Dez '13

12 Dez '13

Caros colegas do R, > > Um amigo me pediu ajuda para fazer um diagrama volumétrico de algumas > propriedades do solo (agua retida a diferentes potenciais, volume de > sólidos, etc). > > Este é uma figura que ilustra todo a distribuição do espaço poroso de um > perfil de solo. Mas não consigo montar este diagrama no R. Talvez alguém > aqui da lista possa me dar uma ajuda ou já tenha visto algo semelhante. Ela > era usualmente feita como figura no Corel. > > A figura/diagrama que comento está neste endereço > > https://www.dropbox.com/l/zPSzOqAKPVSEhHFopbu3Xc > > No diagrama no eixo x fica o volume do solo (variando de 0,0 - 1,0 m3/m3. > No eixo y fica os horizontes/profundidades em cm. > E o volume fica dividido entre os valores de i) CAD1 capacidade de água > disponível – [0,010 – 1,5 MPA], ii) CAD2 Capacidade de água disponível [ > 0,033 – 1,5 MPA], iii) Macroporos; iv) Água fortemente retida; v) água > efêmera (conteúdo de água retida entre 0,01 e 0,033 MPA), vi) Matéria > sólida. Alguns deste parâmetros se sobrepõe. > > Sei que a regra e tentar enviar um CMR. Mas não sei por onde começar. > > > Cordialmente, > > > Wenceslau Teixeira > > Embrapa Solos > > > Horiz. > > Profund. > > CAD1 (0,01) > > CAD2 (0,03) > > Macro > > Água fortemente retida (PM) > > Água retida entre (0,01 e 0,03) > > Matéria sólida > > A > > 0-18 > > 0,079 > > 0,062 > > 0,311 > > 0,206 > > 0,08 > > 0,396 > > AB > > 18-37 > > 0,085 > > 0,057 > > 0,252 > > 0,260 > > 0,016 > > 0,387 > > Bw1 > > 37-90 > > 0,079 > > 0,051 > > 0,251 > > 0,265 > > 0,017 > > 0,388 > > Bw2 > > 90-150 > > 0,089 > > 0,054 > > 0,215 > > 0,263 > > 0,019 > > 0,414 > > > > > > > > > > > > > > > > > > -- > Wenceslau Teixeira > > wgt007(a)gmail.com > skype: wgt007 > > > > > > > > > > > > > > > > > > > > > > > > > > > > -- Wenceslau Teixeira wgt007(a)gmail.com skype: wgt007

3 2