Outubro 2013 - R-br

[OFF -Tópic] Cotas em concursos para negros e índios!
by Andre Oliveira 12 Jun '49

12 Jun '49

Pessoal, tem um colegar meu que vou fazer um concurso ele é descendente de negro ele é pardo. Ele veio com a pergunta como provar que sou negro, pois, tem vagas reservada para negros e índios? Eu sinceramente não sei e não tenho ideia pois nunca tive esta curiosidade. O que de fato comprova sua raça? Algum documento? O que são usados em concursos para provar? obrigado André Oliveira Souza

5 6

Legenda eixo X
by Luciano F. Sgarbi 10 Jun '49

10 Jun '49

Prezados, Estou tentando centralizar a legenda do eixo X em uma figura com vários gráficos. Estou conseguindo centralizar a legenda da linha inferior como no exemplo, mas as linhas do meio não. par(mfrow=c(2,2)) plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 1",outer=TRUE) ### o problema está aqui, pois a legenda desce para a parte inferior da figura. plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 2",outer=TRUE)## Esse dá certo Alguém podia dar uma luz de como eu posso fazer isso? Abraços -- *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

2 2

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

04 Dez '14

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

4 3

gráfico linhas
by Maurício Lordêlo 14 Nov '13

14 Nov '13

Saudações a todos, Observando este exemplo (help da função "ts") z <- ts(matrix(rnorm(300), 100, 3), start = c(1961, 1), frequency = 12) plot(z) tem-se três gráficos empilhados, com uma série em cada um. Gostaria de saber se há possibilidade de que, em cada um dos gráficos, acrescentar outras linhas para efeito de comparação. Por exemplo, suponha que eu construa um único gráfico com quatro linhas da seguinte forma: time<-c(5,7,9,11,13) m1<-c(94.1,91.9,92.1,92.8,92.3) m2<-c(96.3,95.7,95.8,95.1,95.6) m3<-c(91.6,91.0,91.1,90.0,88.6) m4<-c(95.9,92.9,100.0,90.9,91.1) plot(tempo, m1,type="b",main="Percentual de IC´s que contém o parâmetro para X",xlab="tempo",ylab="%",col="blue",ylim=c(70,100)) lines(tempo, m2,col="red",type="b") lines(tempo,m3,col="green",type="b") lines(tempo, m4,col="black",type="b") legend(12.5,85,c("m1","m2","m3","m4"),col=c("blue","red","green","black"),pch=rep(20,2)) No exemplo do help, são três séries em três gráficos empilhados. Há possibilidade de acrescentar linhas a cada uma das séries, ou seja, em cada um dos três gráficos empilhados ter quatro linhas? Grato, Maurício

3 7

Manual
by Maicon J. Galiazzi 13 Nov '13

13 Nov '13

Olá, estou desenvolvendo um manual de fácil utilização ao usuário, e no momento estou descrevendo a interfacedo software R Stuio 0.97.551, assim como a barra de ferramentas do mesmo. Podem me indicar algum material no qual posso conferir? Desde já agradeço! Att Maicon J. Galiazzi

2 2

retidada de tendencia sem UTM
by Hélio Gallo Rocha 04 Nov '13

04 Nov '13

Caros da lista, Estudando outros dados, não houve problema na validação cruzada para retirada de tendencia polinomial de primeiro grau, achei que fosse pela quantidade de pontos, mais de 300. Como minhas coordenadas estão em UTM, subtrai todas elas pelo menor valor. Rodei o CRM... sem problemas. Para quem está em fase de estudos da geoestatistica, ainda faltam pontos de dúvida, favor darem suas opiniões, listo aqui minhas dúvidas e acrescentem outras: 1. Sempre temos de retirar tendencia? 2. Poderia retirar tendencia com uma covariavel fenológica da planta? penso que sim, já que a retirada polinomial de primeiro grau utiliza as coordenadas. 3. Alguém teria algum artigo para embasar este tipo de situação? Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 3

package ‘divisors’ is not available (for R version 3.0.2)
by Andre Oliveira 02 Nov '13

02 Nov '13

Pessoas, preciso da ajuda de vocês estou tentado instalar estes pacotes e obtenho a aviso. install.packages("divisors",dep=T) package ‘divisors’ is not available (for R version 3.0.2) Minha versão consta a informação. R version 3.0.2 (2013-09-25) -- "Frisbee Sailing" Copyright (C) 2013 The R Foundation for Statistical Computing Platform: i686-pc-linux-gnu (32-bit) R é um software livre e vem sem GARANTIA ALGUMA. Você pode redistribuí-lo sob certas circunstâncias. Digite 'license()' ou 'licence()' para detalhes de distribuição. R é um projeto colaborativo com muitos contribuidores. Digite 'contributors()' para obter mais informações e 'citation()' para saber como citar o R ou pacotes do R em publicações. Digite 'demo()' para demonstrações, 'help()' para o sistema on-line de ajuda, ou 'help.start()' para abrir o sistema de ajuda em HTML no seu navegador. Digite 'q()' para sair do R. O que pode estar ocorrendo? obrigado André Oliveira Souza

4 8

Re: [R-br] erro na validação cruzada com pacote geoR
by Hélio Gallo Rocha 31 Out '13

31 Out '13

Caro Eder, Já havia feito a retirada de tendencia com o atributo altura da planta, que está ligado com a produção. O variograma gerado ficou diferente se usar polinomial de 1st, mas não sei como validar isso depois... Abraço Em 31 de outubro de 2013 13:15, Éder Comunello <comunello.eder(a)gmail.com>escreveu: > Boa tarde, Helio! > > Vi seu email, mas vou precisar reservar um tempo pra tratar dele. Você > chegou a usar as glebas como covariável na modelagem? Para isso cada ponto > deveria ter uma var a mais identificando a gleba origem. > > Atte., > > > > Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> > Dourados, MS - [22 16.5'S, 54 49'W] > > > Em 29 de outubro de 2013 22:25, Hélio Gallo Rocha < > heliogallorocha(a)gmail.com> escreveu: > >> Caro Eder, >> >> Só hoje vi seu post com CRM anexo. >> >> Muito obrigado pela revisão do CRM. >> >> Como também está a procura de uma maior entendimento da geoestatística, >> vamos seguindo. >> >> >> Vou numerar para separar as observações: >> >> 1. Pra entender melhor a área, realmente é café, o recorte esta meio >> estranho, mas na realidade são três talhões de 0.5 ha cada. >> Cada um tem uma característica bem diferenciada com relação a >> produção:alta, baixa e média. >> >> Segundo li, numa situação em que a média é flutuante, como parece ser o >> caso, teríamos de usar krigagem ordinária. >> >> >> A maioria das observações tinha notado, a linha: >> plot(dados, lowess=T, trend='1st') >> achei interessante, mostrando o resultado da retirada de tendencia. >> >> 2. no final : >> >> ### Lembra do resíduo que separei anteriormente? >> ### Se fosse usar... >> va.r = variog(dados.res, uvec=uvec, pairs.min=pairs, max.dist=hmax, >> estimator="classical") >> va.r.env <- variog.mc.env(dados.res, obj.var = va.r) >> plot(va.r, xlab="Dist. (m)", ylim=c(0,summary(va.r$v)[6]*1.1), >> env=va.r.env, main="Residual") >> >> Foi possível gerar o semivariograma sem tendencia, com os resíduos, >> fazer o ajuste assentimento usando: eyefit(va.r) >> >> Dai fazer a validação, fazendo a validação ou não e gerar o mapa sem >> tendencia. >> >> 3. Quanto ao envelopamento, não sei se estou certo, mas ele simula >> valores máximos e minimos para dependencia, mas como é uma simulação, a >> cada vez que vc. roda, pode ter pontos com dependencia, nem que seje um, e >> as vezes não dá nenhum, mas de qualquer forma indica que o valor do alcance >> é pequeno em relação ao todo e baixa dependencia espacial, como nos meus >> dados. >> >> Esse item do envelopamento vale discussão, principalmente em situações >> que a variancia sobe além do limite superior e desce, isso significa algo >> importante, mas o que? >> >> 4. Quanto aos links do Landim já tinha visto, usando o surfer, também >> vale atenção. >> Fiz o recomendado por ele, mas dobrou os valores do mapa., vc. tentou? >> >> 5. Quanto ao ajuste por verosimilhança colocada pelo Prof. Paulo: >> Na primeira vez que havia rodado os dados achei muito estranho, pois o >> semivariograma subia indefinidamente e o semivariograma teórico por >> verossimilhança, usando likfit, baixou muito o patamar. >> Interessante se usar likfit com trend="1st", o semivariograma teórico >> baixa mais ainda, das dai não consigo krigar >> >> 6. Não sei se posso fazer o ajuste assentimento usando o semivariograma >> com tendencia mas usar os parametros sem tendencia. >> >> 7. Tenho de retirar a tendencia sempre? >> Tem casos que há uma patamar claro, mas se rodar o variog com trend, fica >> efeito pepita puro...e ai? >> >> 8. Outra coisa que não sei se estou entendendo: >> após a retirada da tendencia, fazer o mapa com o residuo, mas os valores >> não tem nada a ver com os reais, onde estou entendendo errado? >> >> Abraço, >> >> Hélio >> >> >> >> >> >> Em 26 de outubro de 2013 11:05, Eder Comunello [via R-br] < >> ml-node+s2285057n4660724h12(a)n4.nabble.com> escreveu: >> >>> Senhores, bom dia! >>> >>> Agradeço novamente ao Prof. Paulo pela atenção dispensada. >>> >>> Hélio, quanto ao script que eu tinha ficado de olhar, vou colocar minhas >>> considerações abaixo. >>> >>> Procurei não mudar muito, mas acabei alterando um pouco. Os comentários >>> estão no corpo do script. >>> >>> Lembrando sempre que também estou aprendendo sobre o assunto, então pode >>> ter muita coisa errada ou mal interpretada aí. Mas acho que pode te ajudar. >>> >>> Outra coisa é que não tenho muito ideia do que os dados são, como foram >>> obtidos e as condições da área de onde vieram, o que pra mim é fundamental >>> pra análise. Pelos indícios nos dados (localização, estrutura), vou >>> arriscar que são dados da produção de café ou citros e que o terreno está >>> numa pendente da direita pra esquerda, com maior produção na parte baixa. O >>> que chamou a atenção é o formato da área que é bem particular e não parece >>> muito 'natural'. >>> >>> Mas como o foco principal é o uso do R/geoR, vou tentar me deter nisso. >>> >>> Use por sua própria conta e risco! :D >>> >>> >>> ### <BEGIN> >>> >>> require(geoR) #;require(tcltk2) >>> setwd("C:/LAB/RGIS/trend"); dir() >>> par.ori <- par(no.readonly = TRUE) >>> >>> url1 <- 'https://dl.dropboxusercontent.com/u/117618178/helio/dados.txt' >>> url2 <- 'https://dl.dropboxusercontent.com/u/117618178/helio/recorte.txt >>> ' >>> >>> ### Como criar uma pasta pública no Dropbox... >>> browseURL('https://www.dropbox.com/help/16/pt_BR') >>> >>> ### 1. Carregando dados >>> ------------------------------------------------------ >>> dados <- read.geodata(url1, coords.col=1:2, data.col=3, sep="", header=T) >>> dados$borders <- read.table(url2, header=T) >>> summary(dados); length(dados$data) >>> ls(dados);ls.str(dados) >>> >>> ### 2. Visualizando dados >>> ---------------------------------------------------- >>> >>> ### Visão geral >>> par(mfrow=c(1,1)); par()$mfrow >>> points(dados) >>> points(dados, pt.div="quart") >>> points(dados, pt.div="quint") >>> ### tem pelos três sub-áreas bem marcadas... >>> >>> ### Verificando tendências com o plot.geodata >>> -------------------------------- >>> ### Antes de partir pros variogramas, use o plot.geodata para >>> pré-avaliar tendências >>> ### Bom para esclarecer opções de trend, pois notei certa confusão >>> (cte,1st,2nd) >>> ### Recomendo ler o item 'trend' na ajuda da função plot.geodata >>> ?plot.geodata >>> >>> plot(dados, lowess=T) >>> plot(dados, lowess=T, trend='cte') ### mesmo que anterior >>> ### linhas de comando anteriores são equivalentes: trend='cte' é >>> 'default' da geoR! >>> ### nos gráf. 2 e 3, plota-se as coord vs. dados. Observar as distorções >>> na linha 'lowess' >>> >>> plot(dados, lowess=T, trend=~coords) >>> plot(dados, lowess=T, trend='1st') >>> plot(dados, lowess=T, trend=~dados$coords[,1]+dados$coords[,2]) >>> ### as três linhas anteriores são equivalentes >>> ### o uso da polinomial de primeira ordem melhora o comportamento do >>> histograma >>> ### mas a linha de lowess continua com comportamento anômalo >>> >>> plot(dados, lowess=T, trend='2nd') >>> plot(dados, lowess=T, >>> trend=~I(coords)+I(coords^2)+I(coords[,1]*coords[,2])) >>> ### as duas linhas anteriores são equivalentes (polinomal de segunda >>> ordem)! >>> ### não aparenta agregar nenhuma melhoria >>> >>> ### Mais sobre modelos de tendência com trend.spatial >>> ------------------------ >>> ### Pra melhorar o entendimento dos modelos, dá pra dar uma olhada nas >>> matrizes >>> ?trend.spatial >>> head(unclass(trend.spatial('cte', dados))) ### constante >>> head(unclass(trend.spatial(~coords, dados))) ### polinomal 1a ordem, >>> usando coords >>> head(unclass(trend.spatial('1st', dados))) ### mesmo que anterior >>> head(unclass(trend.spatial('2nd', dados))) ### polinomal 2a ordem, >>> usando coords >>> head(unclass(trend.spatial(~coords[,1], dados))) ### em função de X >>> head(unclass(trend.spatial(~coords[,2], dados))) ### em função de Y >>> head(unclass(trend.spatial(~coords+I(coords[,1]^2), dados))) ### coords >>> + X**2 >>> head(unclass(trend.spatial(~I(coords)+I(coords^2), dados))) >>> head(unclass(trend.spatial(~I(coords)+I(coords^2)+I(coords[,1]*coords[,2]), >>> dados))) ### mesmo que 2nd >>> head(unclass(trend.spatial(~I(coords)+I(coords[,1]^2), dados))) >>> head(unclass(trend.spatial(~I(coords^2), dados))) >>> head(unclass(trend.spatial(~coords[,1]*coords[,2], dados))) >>> head(unclass(trend.spatial(~I(coords[,1]*coords[,2]), dados))) >>> head(unclass(trend.spatial(~coords[,1]+I(coords[,1]^2)+coords[,2], >>> dados))) >>> >>> ### Separando o drift com base na matriz trend.spatial >>> ----------------------- >>> trend.matrix <- unclass(trend.spatial(trend='1st', geodata=dados)); >>> nrow(trend.matrix) >>> trend.lm <- lm(dados$data ~ trend.matrix) >>> trend.res <- lm(dados$data ~ trend.matrix)$residuals >>> trend.fit <- lm(dados$data ~ trend.matrix)$fitted.values >>> >>> ### 'Empacotando' os resíduos em um >>> geodata----------------------------------- >>> dados.res <- dados >>> names(trend.res) <- NULL; trend.res >>> dados.res$data <- trend.res; dados.res$data >>> str(dados.res) >>> ### Pode usar dados.res pra fazer a análise variográfica no lugar de >>> dados, >>> ### caso tenha problemas com krige.conv()!!! >>> >>> ### Embora o uso da polinomial de 1a ordem tenha representado alguma >>> melhoria, >>> ### eu não fiquei muito entusiamado por conta do 'lowess' apresentado >>> ### Talvez fosse o caso de testar covariáveis. >>> ### Uma fácil de obter e que pode ter relação é a altitude... >>> ### Mas vamos investigar o efeito nos variogramas >>> >>> ### 3. Distancia máxima e parâmetros para semivariograma >>> --------------------- >>> hmax <- summary(dados)[[3]][[2]]*.8 ; hmax ### 80% da máxima distância >>> uvec <- 12; pairs <- 20 ### pairs.min >>> >>> ### 4. Variogramas >>> ----------------------------------------------------------- >>> ### Não precisaria rodar todos, afinal como visto, tá repetindo >>> modelos!!! >>> ### v0:média constante; v1:1st; v2:2nd; v3:cte; v4:coords >>> v0 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, >>> estimator="classical") >>> v1 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, >>> estimator="classical", trend="1st") >>> v2 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, >>> estimator="classical", trend="2nd") >>> v3 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, >>> estimator="classical", trend="cte") >>> v4 <- variog(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, >>> estimator="classical", trend=~coords) >>> >>> ### envelopes >>> v0.env <- variog.mc.env(dados, obj.var = v0) >>> v1.env <- variog.mc.env(dados, obj.var = v1) >>> v2.env <- variog.mc.env(dados, obj.var = v2) >>> v3.env <- variog.mc.env(dados, obj.var = v3) >>> v4.env <- variog.mc.env(dados, obj.var = v4) >>> >>> ### 5. Visualizando variogramas >>> ---------------------------------------------- >>> par(par.ori) >>> par(mfrow=c(2,3)) >>> >>> points(dados, pt.div="quartiles",xlab="Leste", ylab="Norte") >>> >>> plot(v0, xlab="Dist. (m)", ylim=c(0,summary(v0$v)[6]*1.1), env=v0.env, >>> main="Isotrópico com tendência") >>> >>> text(v0$u,v0$v,round(v0$n,1),pos=1);text(v0$u,v0$v,round(v0$u,1),pos=3);lines(v0$u,v0$v);abline(h=mean(v0$v)) >>> >>> plot(v1, xlab="Dist. (m)", ylim=c(0,summary(v1$v)[6]*1.1), env=v1.env, >>> main="Isotrópico trend 1st") >>> >>> text(v1$u,v1$v,round(v1$n,1),pos=1);text(v1$u,v1$v,round(v1$u,1),pos=3);lines(v1$u,v1$v);abline(h=mean(v1$v)) >>> >>> plot(v2, xlab="Dist. (m)", ylim=c(0,summary(v2$v)[6]*1.1), env=v2.env, >>> main="Isotrópico trend 2nd") >>> >>> text(v2$u,v2$v,round(v2$n,1),pos=1);text(v2$u,v2$v,round(v2$u,1),pos=3);lines(v2$u,v2$v);abline(h=mean(v2$v)) >>> >>> plot(v3, xlab="Dist. (m)", ylim=c(0,summary(v3$v)[6]*1.1), env=v3.env, >>> main="Isotrópico trend cte") >>> >>> text(v3$u,v3$v,round(v3$n,1),pos=1);text(v3$u,v3$v,round(v3$u,1),pos=3);lines(v3$u,v3$v);abline(h=mean(v3$v)) >>> >>> plot(v4, xlab="Dist. (m)", ylim=c(0,summary(v4$v)[6]*1.1), env=v4.env, >>> main="Isotrópico trend coords") >>> >>> text(v4$u,v4$v,round(v4$n,1),pos=1);text(v4$u,v4$v,round(v4$u,1),pos=3);lines(v4$u,v4$v);abline(h=mean(v4$v)) >>> >>> ### Analisando os variogramas >>> ------------------------------------------------ >>> ### Realmente o comportamento do histograma de trend='cte' parece conter >>> uma tendência, >>> ### mas os modelos de 'retirada' usados deixam os dados sem uma >>> estrutura de >>> ### variação. Dificilmente vai conseguir ajustar um bom modelo nesses >>> casos, pois >>> ### o comportamento é muito errático. É possível que O erro que você >>> estava obtendo >>> ### seja decorrente da fragilidade do modelo ajustado. >>> >>> ### Realmente não ficou bom! Mas e aí? O que dá pra fazer? Tenta ssim >>> mesmo? >>> ### Ainda tem que testar a anisotropia, que não tem no seu script >>> original. >>> ?variog4 >>> >>> va <- variog4(dados, uvec=uvec, pairs.min=pairs, max.dist=hmax, >>> estimator="classical") >>> par(par.ori) >>> plot(va) >>> >>> ### Pronto! Apareceu anisotropia... >>> ### Tem duas sugestões de leitura abaixo pra dar uma ajuda... >>> browseURL(' >>> http://www.rc.unesp.br/igce/aplicada/DIDATICOS/LANDIM/tkrigagem.pdf') >>> browseURL(' >>> https://www.soils.org/publications/sssaj/abstracts/61/1/SS0610010298') >>> >>> ### Você vai ter que estudar adequadamente a questão da anisotropia em >>> seus dados >>> ### Isso não te isenta do efeito combinado anisotropia + tendência >>> ### Só pra testar as possibilidades, vamos fazer o variograma na direção >>> 135 graus >>> >>> ### Vai ter que entrar a direções, sendo mais fácil usar PI/n >>> ### Tabela das direções (PI/n) para variogramas >>> denom <- c(1:4,6,8,10,12,24,36,45,60); k <- (pi/180) ### Denominadores >>> data.frame(rad=paste0('pi/', denom), rad2=round(pi/denom,4), >>> deg=pi/denom/k) >>> ### 135 é 3*45 ou 3*pi/4 >>> >>> va.135 <- variog(dados, uvec=uvec, pairs.min=pairs, dir=3*pi/4, >>> max.dist=hmax, estimator="classical") >>> plot(va.135) >>> ### Aqui parece ser mais fácil ajustar. Talvez tenha que ajustar o 'lag' >>> ### Avaliar a questão da distância entre amostras... >>> >>> ### Lembra do resíduo que separei anteriormente? >>> ### Se fosse usar... >>> va.r = variog(dados.res, uvec=uvec, pairs.min=pairs, max.dist=hmax, >>> estimator="classical") >>> va.r.env <- variog.mc.env(dados.res, obj.var = va.r) >>> plot(va.r, xlab="Dist. (m)", ylim=c(0,summary(va.r$v)[6]*1.1), >>> env=va.r.env, main="Residual") >>> >>> ### Vou parar por aqui, porque antes de prosseguir vai precisar ver >>> ### adequadamente essa questão da anisotropia. >>> ### No caso do efeito ser só da anisotropia, vai simplificar a krigagem >>> e você >>> ### não deve mais ter o problema que vinha tendo. >>> >>> ### <END> >>> >>> Éder Comunello <[hidden email]<http://user/SendEmail.jtp?type=node&node=4660724&i=0>[hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4660724&i=1>> >>> Dourados, MS - [22 16.5'S, 54 49'W] >>> >>> >>> >>> >>> _______________________________________________ >>> R-br mailing list >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4660724&i=2> >>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br >>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça >>> código mínimo reproduzível. >>> >>> ------------------------------ >>> If you reply to this email, your message will be added to the >>> discussion below: >>> >>> http://r-br.2285057.n4.nabble.com/Re-R-br-erro-na-validacao-cruzada-com-pac… >>> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >>> . >>> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >>> >> >> >> >> -- >> Hélio Gallo Rocha >> IFSULDEMINAS - Câmpus Muzambinho >> >> _______________________________________________ >> 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. >> > > -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

1 0

Legenda em gráfico de regressão
by Edson Lira 31 Out '13

31 Out '13

Caros amigos, tenho a rotina abaixo: require(PASWR) data(vit2005) attach(vit2005) names(vit2005) head(vit2005) Elevator<-as.factor(elevator) contrasts(Elevator) modSimpl<-lm(totalprice~area) summary(modSimpl) plot(totalprice~area) Este banco tem uma variável que é o Elevator (presença ou ausência do elevador), gostaria de plotar o gráfico da regressão totalprice~area, identificando os pontos aonde tem a presença ou ausência do elevador. Como faço? plot(totalprice~area) Edson Lira Estatístico Manaus-Amazonas

3 2

ggplot2
by geovane barbosa 31 Out '13

31 Out '13

OLá pessoal tudo bem, estou usando o seguinte código e esta dando o seguinte erro Error: Aesthetics must either be length one, or the same length as the dataProblems:IDADEA, factor(Cor2) Alguém poderia me dar uma dica aguardo #################################### #################################### Prof. Geovane Carlos Barbosa UCL - Faculdade do Centro Leste ################################### ###################################

3 2