<div dir="ltr"><span style="font-family:arial,sans-serif;font-size:13px">Mauro, bom dia.</span><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><font face="arial, sans-serif">Não sei qual seu objetivo no uso desta função, assim vou te passar o histórico e se quiser dando idéias, será bem vindo</font></div><div style="font-family:arial,sans-serif;font-size:13px"><font face="arial, sans-serif">Estou estudando geoestatística e é interessante saber que regiões tem seus atributos "parecidos".</font></div><div style="font-family:arial,sans-serif;font-size:13px"><font face="arial, sans-serif">A idéia inicial era utilizar a função skater do pacote spdep.</font></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Como tenho pontos teria de criar poligonos, que foi realizado com a função deldir, e listar os poligonos.</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">A listagem era um problemão, pois teria de crira poligono por poligono, como no CRM inicial.</div><div style="font-family:arial,sans-serif;font-size:13px">Problema resolvido com a preciosa ajuda do Rafael. que sugeriu usar a função poly2nb do pacote sddep</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Com a função deldir apenas gera poligonos assim foi necessário inserir dados aos mesmo e torná-los um SpatialPolygonsDataFrame<br></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Dai pra frente o código está em estudo, veja:</div><div style="font-family:arial,sans-serif;font-size:13px"><div>require(deldir) ;require(spdep)</div><div>## 1. Coordenadas</div><div>n <- 40 </div><div>set.seed(47)</div><div>x <- runif(n)</div><div>y <- runif(n)</div><div>plot(x,y)</div><div>## 2. Dados</div><div> set.seed(1389)</div><div>dados <- data.frame(rnorm(n,10,sd=5), rnorm(n,15,sd=5),rnorm(n,20,sd=5))</div><div>colnames(dados)<- paste("var",1:3,sep="")</div><div>row.names(dados) <- paste("reg", 1:n, sep = "")</div><div>dados</div><div>class(dados)</div><div><br></div><div>## 3. gerando coordenadas dos poligonos com função deldir</div><div>d <- deldir(x,y)</div><div><br></div><div>ttt <- triang.list(d)</div><div>plot(ttt)</div><div>plot(ttt,border="red",showrect=TRUE)</div><div>sss <- tile.list(d)</div><div>plot(sss)</div><div>plot(ttt,add=TRUE,border="blue")</div><div><br></div><div>## 4. Gerando poligonos</div><div>### Rafael Tieppo---- G R A N D E A J U D A ! ! ! !</div><span class="im"><div><br></div><div><br></div><div>FUN_VOR_MAT <- function (A,TT) #TT é o arquivo, N o número do POL</div><div> {</div><div> AA <- 1</div><div> BB <- 2</div><div> NROW <- length (A[[TT]]$x) + 1 #add primeiro ponto</div><div><br></div><div><br></div><div> VORONOI_MATRIX <- matrix(0,NROW,2)</div><div> VORONOI_MATRIX[1:NROW,AA] <- A[[TT]]$x[1:NROW]</div><div> VORONOI_MATRIX[1:NROW,BB] <- A[[TT]]$y[1:NROW]</div><div> VORONOI_MATRIX[NROW,AA] <- A[[TT]]$x[1] #+ 1 add primeiro ponto</div><div> VORONOI_MATRIX[NROW,BB] <- A[[TT]]$y[1] #+ 1 add primeiro ponto</div><div><br></div><div> if(TT == NROW){</div><div> print("17/10/2014, [hidden email] ")}</div><div> return(VORONOI_MATRIX)</div><div> }</div><div><br></div></span><div>## Agora vc vai usar a funcao:</div><div>AREAS<- lapply(1:length(sss), function (x) ( FUN_VOR_MAT(sss,x)) )</div><div> summary(AREAS)</div><div>## 5. Formando Polygons</div><span class="im"><div>POLS<- lapply(1:length(AREAS), function(x) Polygons(list(Polygon(AREAS[[x]])),</div><div> ID = paste("reg", x, sep = "")))</div></span><div> summary(POLS)</div><div> class(POLS)</div><div> </div><div><br></div><div> SPol <- SpatialPolygons(POLS) </div><div> class(SPol) </div><div> plot(SPol)</div><div> bbox(SPol)</div><div> </div><div>## 6. Inserindo dados aos poligonos </div><div> SPolDF <- SpatialPolygonsDataFrame(SPol, dados) </div><div>plot(SPolDF)</div><div> </div><div> </div><div>## 7. Código para visualizar áreas "parecidas"</div><div><br></div><div>## em estudo</div><div>(coordenadas=cbind(x,y))</div><div>class(coordenadas)</div><div><br></div><div><br></div><div>## </div><div>rn <- sapply(slot(SPolDF, "polygons"), function(x) slot(x, "ID"))</div><div><br></div><div>knn1 <- knearneigh(coordenadas, 1)</div><div><br></div><div>knn2 <- knearneigh(coordenadas, 2)</div><div><br></div><div>nb1 <- knn2nb(knn1, row.names=rn)</div><div><br></div><div>nb2 <- knn2nb(knn2, row.names=rn)</div><div><br></div><div>##1</div><div>diffs <- diffnb(nb2, nb1)</div><div>plot(SPolDF, border="grey")</div><div>plot(nb1, coordenadas, add=TRUE)</div><div>plot(diffs, coordenadas, add=TRUE, col="red")</div><div>title(main=paste("Differences (red) in Columbus GAL weights (black)", "and polygon generated queen weights", sep="\n"))</div><div> </div><div> </div><div>##2 </div><div>xxx <- poly2nb(SPolDF, queen=FALSE)</div><div>dxxx <- diffnb(xxx, nb1)</div><div>plot(SPolDF, border = "grey")</div><div>plot(nb1, coordenadas, add = TRUE)</div><div>plot(dxxx, coordenadas, add = TRUE, col = "red")</div><div>title(main=paste("Differences (red) in Columbus GAL weights (black)",</div><div> "and polygon generated rook weightsvvv", sep="\n"))</div><div><br></div><div>## 3</div><div>cards <- card(nb2)</div><div>maxconts <- which(cards == max(cards))</div><div>if(length(maxconts) > 1) maxconts <- maxconts[1]</div><div>fg <- rep("grey", length(cards))</div><div>fg[maxconts] <- "red"</div><div>fg[nb2[[maxconts]]] <- "green"</div><div>plot(SPolDF, col=fg)</div><div>title(main="Region with largest number of contiguities")</div><div><br></div><div><br></div><div>plot(SPolDF, border="grey")</div><div>plot(nb1, coordenadas, add=TRUE)</div><div>plot(diffs, coordenadas, add=TRUE, col="red", lty=2)</div></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Abraço</div><div style="font-family:arial,sans-serif;font-size:13px"><br></div><div style="font-family:arial,sans-serif;font-size:13px">Hélio</div></div><div class="gmail_extra"><br><div class="gmail_quote">Em 21 de outubro de 2014 23:41, Mauro Sznelwar [via R-br] <span dir="ltr"><<a href="mailto:ml-node+s2285057n4663338h46@n4.nabble.com" target="_blank">ml-node+s2285057n4663338h46@n4.nabble.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<p class="MsoNormal"><span style="font-family:"Arial","sans-serif";color:#1f497d">Tentei
rodar vosso script, mas a partir daí não foi mais, qual o motivo?<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1f497d">> plot(sss)<u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1f497d">>
polygon(p1,border="green",lwd=3) <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1f497d">Error in xy.coords(x,
y) : <u></u><u></u></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1f497d"> 'x' is a list, but
does not have components 'x' and 'y'</span><span lang="EN-US"><u></u><u></u></span></p><div><div class="h5">
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";color:#1f497d"><u></u> <u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Caros Elias e Rafael.<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Rafael, obrigado pelo código, vou testar, mas depois de
conseguir separar os polígonos.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">Na realidade é essa a minha dificuldade.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Fiz uns testes que deram certo, mas estou separando um
polígono por vez, assim fica muito trabalhoso, de uma olhada no código que é a
partir do enviado no primeiro post:<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## objeto que possui as coordenadas geradas na função deldir
é o objeto sss<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## coordenadas do poligono 01<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">(x1 <- sss[[1]][2])<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">(y1 <- sss[[1]][3])<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## coordenadas do poligono 02<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">x2 <- sss[[2]][2]<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">y2 <- sss[[2]][3]<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">.<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">## coordenadas do poligono 20<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">x20 <- sss[[20]][2]<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">y20 <- sss[[20]][3]<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">## Agrupando as coordenadas num só objeto<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">(p1 <-c(x1,y1))<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"> <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">(p2 <-c(x2,y2))<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">(p20 <-c(x20,y20))<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## visualizando os poligonos<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">plot(sss)<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">polygon(p1,border="green",lwd=3) <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">polygon(p2,border="red",lwd=3) <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">polygon(p20,border="blue",lwd=3) <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## Os passos seguintes diferem do seu<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## tentativa de gerar um SpatialPolygons, que tem de
estar fechado, ou seja a primeira e a ultima linhas de coordenadas tem de ser a
mesma<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<div>
<p class="MsoNormal">(x11 <- c(p1$x, p1$x[1]))<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">(y11 <- c(p1$y, p1$y[1]))<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">(p1 <- cbind(x11,y11))<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## transformando em SpatialPolygons<u></u><u></u></p>
</div>
<div>
<div>
<p class="MsoNormal">### Passo 1<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">linep1=Polygon(p1);class(linep1) #Polygon<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">linep1<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## Passo 2<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">linep1a=Polygons(list(linep1),ID="a");class(linep1a)
# Polygons<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">linep1a<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## Passo 3<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">linesp=SpatialPolygons(list(linep1a));class(linesp)
#SpatialPolygons<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">linesp<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">## Passo 4<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">proj4string(linesp) <- CRS("+proj=utm + zone=23s
+datum=WGS84")<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">summary(linesp)<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal">plot(linesp)<u></u><u></u></p>
</div>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Se puder dar uma mão, agradeço desde já..<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Abraço<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Hélio<u></u><u></u></p>
</div>
</div>
</div></div><div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<p class="MsoNormal">Em 20 de outubro de 2014 13:22, Rafael Tieppo [via R-br]
<<a href="http://user/SendEmail.jtp?type=node&node=4663338&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>>
escreveu:<u></u><u></u></p><div><div class="h5">
<div>
<div>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">Hélio,<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">Após dividir sua área em polígonos, vc terá um conjunto de
coordenadas, onde cada conjunto será um polígono.<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">Com esse objeto pronto em forma de lista:<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">Passo 1: Formar "Polygons"<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">POLS <- lapply(1:length(AREAS), function(x)
Polygons(list(Polygon(AREAS[[x]])),<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"> ID = paste("reg", x, sep =
"")))<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">Passo 2:Converter em "SpatialPolygons"<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">POLS <- as.SpatialPolygons.PolygonsList(POLS,
proj4string=CRS(as.character(NA)))<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">Agora vc já pode usar a função poly2nb()<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"> poly2nb(POLS)<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black">Tem um tutorial bem detalhado do professor:<u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"><a href="http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase31.html" rel="nofollow" link="external" target="_blank">Classes para dados espaciais: o pacote sp</a><u></u><u></u></span></p>
</div>
<div style="margin-top:3.75pt;margin-bottom:3.75pt">
<table border="0" cellspacing="0" cellpadding="0" width="450" style="width:337.5pt;border-collapse:collapse">
<tbody><tr style="height:.75pt">
<td colspan="8" style="background:#e5e5e5;padding:0cm 0cm 0cm 0cm;height:.75pt">
<div>
<p class="MsoNormal" style="background:#e5e5e5"><span style="font-size:1.0pt"> <u></u><u></u></span></p>
</div>
</td>
</tr>
<tr style="height:4.5pt">
<td width="1" rowspan="5" style="width:.75pt;background:#e5e5e5;padding:0cm 0cm 0cm 0cm;height:4.5pt">
<div>
<p class="MsoNormal" style="background:#e5e5e5"><span style="font-size:1.0pt"> <u></u><u></u></span></p>
</div>
</td>
<td width="168" rowspan="5" style="width:126.0pt;background:black;padding:0cm 0cm 0cm 0cm;height:4.5pt">
<p class="MsoNormal" align="center" style="text-align:center"><a href="http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase31.html" rel="nofollow" link="external" target="_blank"><span style="text-decoration:none"><img border="0" width="168" height="168" src="https://ec.yimg.com/mail?url=http%3A%2F%2Fwww.leg.ufpr.br%2F%7Epaulojus%2Fembrapa%2FRembrapa%2FRembrapa202x.png&t=1413818683&sig=EnbxwD80kyOV.7nKWeJpzw--~B" alt="image"></span></a><u></u><u></u></p>
</td>
<td width="1" rowspan="5" style="width:.75pt;background:#e5e5e5;padding:0cm 0cm 0cm 0cm;height:4.5pt">
<div>
<p class="MsoNormal" style="background:#e5e5e5"><span style="font-size:1.0pt"> <u></u><u></u></span></p>
</div>
</td>
<td width="14" rowspan="5" style="width:10.5pt;background:white;padding:0cm 0cm 0cm 0cm;height:4.5pt">
<div>
<p class="MsoNormal" style="background:white"><span style="font-size:14.0pt"> <u></u><u></u></span></p>
</div>
</td>
<td colspan="2" style="background:white;padding:0cm 0cm 0cm 0cm;height:4.5pt">
<div>
<p class="MsoNormal" style="background:white"><span style="font-size:6.0pt"> <u></u><u></u></span></p>
</div>
</td>
<td width="20" rowspan="5" style="width:15.0pt;background:white;padding:0cm 0cm 0cm 0cm;height:4.5pt">
<div>
<p class="MsoNormal" style="background:white"><span style="font-size:20.0pt"> <u></u><u></u></span></p>
</div>
</td>
<td width="1" rowspan="5" style="width:.75pt;background:#e5e5e5;padding:0cm 0cm 0cm 0cm;height:4.5pt">
<div>
<p class="MsoNormal" style="background:#e5e5e5"><span style="font-size:1.0pt"> <u></u><u></u></span></p>
</div>
</td>
</tr>
<tr>
<td width="100%" colspan="2" style="width:100.0%;padding:0cm 0cm 0cm 0cm">
<div>
<div>
<p class="MsoNormal" style="line-height:12.4pt;background:white;word-break:break-all"><span style="font-family:"Georgia","serif""><a href="http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase31.html" rel="nofollow" link="external" target="_blank"><span style="font-size:13.5pt;color:black">Classes para dados
espaciais: o pacote sp</span></a><u></u><u></u></span></p>
<div>
<p class="MsoNormal" style="line-height:15.0pt;background:white;word-break:break-all"><span style="font-size:10.0pt;font-family:"Georgia","serif";color:#999999">[next] [prev] [prev-tail] [tail] [up] 31 Classes para dados
espaciais: o pacote sp <u></u><u></u></span></p>
</div>
</div>
</div>
</td>
</tr>
<tr style="height:3.0pt">
<td colspan="2" style="background:white;padding:0cm 0cm 0cm 0cm;height:3.0pt"></td>
</tr>
<tr>
<td style="padding:0cm 0cm 0cm 0cm">
<div>
<p class="MsoNormal"><span style="font-size:1.0pt;font-family:"Arial","sans-serif""><a href="http://www.leg.ufpr.br/~paulojus/embrapa/Rembrapa/Rembrapase31.html" rel="nofollow" link="external" target="_blank"><span style="font-size:7.0pt;color:#999999">View on <b>www.leg.ufpr.br</b></span></a><u></u><u></u></span></p>
</div>
</td>
<td width="100" style="width:75.0pt;padding:0cm 0cm 0cm 0cm">
<p class="MsoNormal" align="right" style="text-align:right;line-height:8.25pt"><span style="font-size:7.0pt;font-family:"Arial","sans-serif";color:#999999">Preview
by Yahoo</span><span style="font-size:1.0pt;font-family:"Arial","sans-serif""><u></u><u></u></span></p>
</td>
</tr>
<tr style="height:6.75pt">
<td colspan="2" style="background:white;padding:0cm 0cm 0cm 0cm;height:6.75pt"></td>
</tr>
<tr style="height:.75pt">
<td colspan="8" style="background:#e5e5e5;padding:0cm 0cm 0cm 0cm;height:.75pt">
<div>
<p class="MsoNormal" style="background:#e5e5e5"><span style="font-size:1.0pt"> <u></u><u></u></span></p>
</div>
</td>
</tr>
<tr height="0">
<td width="1" style="border:none"></td>
<td width="168" style="border:none"></td>
<td width="1" style="border:none"></td>
<td width="5" style="border:none"></td>
<td width="139" style="border:none"></td>
<td width="128" style="border:none"></td>
<td width="7" style="border:none"></td>
<td width="1" style="border:none"></td>
</tr>
</tbody></table>
</div>
<div>
<p class="MsoNormal"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"><u></u> <u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="background:white"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"> <u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt;background:white"><b><span style="font-size:7.5pt;font-family:"Verdana","sans-serif";color:black">Rafael
Tieppo</span></b><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"><u></u><u></u></span></p>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt;background:white"><span style="font-size:10.0pt;font-family:"Lucida Console";color:black"><u></u> <u></u></span></p>
</div>
<div>
<div>
<div>
<div>
<p class="MsoNormal" style="background:white"><span style="font-family:"Arial","sans-serif";color:black">On Sunday, October 19, 2014 7:52 PM, "<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=0" rel="nofollow" link="external" target="_blank">[hidden email]</a>" <<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>> wrote:</span><span style="font-family:"Helvetica","sans-serif";color:black"><u></u><u></u></span></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt;background:white"><span style="font-family:"Helvetica","sans-serif";color:black"><u></u> <u></u></span></p>
<div>
<p class="MsoNormal" style="background:white"><span style="font-family:"Helvetica","sans-serif";color:black">Enviar submissões para a lista de discussão R-br para <br>
<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=2" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
<br>
Para se cadastrar ou descadastrar via WWW, visite o endereço<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
ou, via email, envie uma mensagem com a palavra 'help' no assunto ou<br>
corpo da mensagem para <br>
<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=3" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
<br>
Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br>
endereço<br>
<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=4" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
<br>
Quando responder, por favor edite sua linha Assunto assim ela será<br>
mais específica que "Re: Contents of R-br digest..."<br>
<br>
<br>
Tópicos de Hoje:<br>
<br>
1. Integral (<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=5" rel="nofollow" link="external" target="_blank">[hidden email]</a>)<br>
2. skater e deldir (Hélio Gallo Rocha)<br>
3. Re: Integral (walmes .)<br>
4. Re: Integral (Diogo Ferrari)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Sun, 19 Oct 2014 17:23:19 -0200<br>
From: <a href="http://user/SendEmail.jtp?type=node&node=4663312&i=6" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
To: "<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=7" rel="nofollow" link="external" target="_blank">[hidden email]</a>" <<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=8" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
Subject: [R-br] Integral<br>
Message-ID: <<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=9" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/d7d3470d/attachment-0001.html" rel="nofollow" link="external" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/d7d3470d/attachment-0001.html</a>><br>
-------------- Próxima Parte ----------<br>
Um anexo não-texto foi limpo...<br>
Nome: Integral.jpg<br>
Tipo: image/jpeg<br>
Tamanho: 38510 bytes<br>
Descrição: não disponível<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/d7d3470d/attachment-0001.jpg" rel="nofollow" link="external" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/d7d3470d/attachment-0001.jpg</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Sun, 19 Oct 2014 18:29:49 -0200<br>
From: Hélio Gallo Rocha <<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=10" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
To: r-br <<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=11" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
Subject: [R-br] skater e deldir<br>
Message-ID:<br>
<CA+<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=12" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
Content-Type: text/plain; charset="utf-8"<u></u><u></u></span></p>
<div>
<div>
<p class="MsoNormal" style="background:white"><span style="font-family:"Helvetica","sans-serif";color:black"><br>
<br>
Caros da lista,<br>
<br>
Estou tentando implementar o uso de SKATER do pacote spdep.<br>
Gostaria de contar com a boa vontade dos listeiros para resolver esta<br>
dificuldade.<br>
<br>
Para iniciar o skater tenho de dividir minha área em poligonos.<br>
Para tal estou usando o pacote deldir. Abaixo um exemplo da própria página<br>
CRAN.<br>
<br>
<br>
require(deldir)<br>
set.seed(42)<br>
x <- runif(20)<br>
y <- runif(20)<br>
d <- deldir(x,y)<br>
(ttt <- triang.list(d) )<br>
plot(ttt,border="red",showrect=TRUE)<br>
(sss <- tile.list(d))<br>
plot(sss)<br>
plot(ttt,add=TRUE,border="blue")<br>
<br>
<br>
O objeto ttt lista os conjuntos de três pontos que formam o poligono<br>
<br>
Teria uma forma automatizada para gerar estes poligonos no formato shape,<br>
pois o CRM do skater que tenho utiliza o comando "readShapePoly" do
spdep.<br>
<br>
Tentei exportar para o formato XLS, mas não foi possível pois o objeto ttt<br>
não é um data.frame<br>
<br>
<br>
<br>
Agradeço desde já<br>
-- <br>
Hélio Gallo Rocha<br>
IFSULDEMINAS - Câmpus Muzambinho<u></u><u></u></span></p>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt;background:white"><span style="font-family:"Helvetica","sans-serif";color:black">-------------- Próxima
Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/e658abd9/attachment-0001.html" rel="nofollow" link="external" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/e658abd9/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Sun, 19 Oct 2014 19:40:44 -0200<br>
From: "walmes ." <<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=13" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
To: <a href="http://user/SendEmail.jtp?type=node&node=4663312&i=14" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
Subject: Re: [R-br] Integral<br>
Message-ID:<br>
<CAFU=EkZ_E8hLAgXoUk79T7u6qpi4LLtbqrNAvjCXx+<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=15" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
?<br>
Eu gosto muito de desafios de cálculo mas aprecio muito mais o desafio<br>
quando sei do que se trata, quando tem contexto. Você não forneceu nenhum.<br>
Meu julgamento adivinhatório de imprecisão acha que tem a ver com<br>
estatísticas de ordem. Essa lista tem a prerrogativa de não resolver<br>
problemas de sala de aula para não tirar do aluno o privilégio do<br>
conhecimento a se adquirir em resolvê-lo. Espero que não seja esse o caso.<br>
A lista também pede para não encaminhar conteúdo em anexo. Como se trata de<br>
uma sequência de desenvolvimento matemático, um opção seria fazer no<br>
eqneditor (<a href="http://www.codecogs.com/latex/eqneditor.php" rel="nofollow" link="external" target="_blank">http://www.codecogs.com/latex/eqneditor.php</a>)
e colar o código<br>
latex correspondente, como por exemplo esse código<br>
<br>
\begin{align*}<br>
f(x) &= \theta_y+\theta_2 (x-\theta_x)^2 \\<br>
&= \theta_y+\theta_2 (x^2-2\theta_x x+\theta_x^2)\\<br>
&=<br>
\underbrace{(\theta_y+\theta_2\theta_x^2)}_{\beta_0}-\underbrace{2\theta_2\theta_x}_{-\beta_1}<br>
x+\underbrace{\theta_2}_{\beta_2} x^2\\<br>
&= \beta_0+\beta_1 x+\beta_2 x^2<br>
\end{align*}<br>
<br>
Aí a coisa funciona como CMR, onde o código é latex. É só recortar esse<br>
código e colar no eqneditor. O mais experts em latex podem editar sem<br>
precisar renderizar, inclusive.<br>
<br>
À disposição.<br>
Walmes.<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/587e8e00/attachment-0001.html" rel="nofollow" link="external" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/587e8e00/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Sun, 19 Oct 2014 17:52:04 -0400<br>
From: Diogo Ferrari <<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=16" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
To: <a href="http://user/SendEmail.jtp?type=node&node=4663312&i=17" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
Subject: Re: [R-br] Integral<br>
Message-ID: <<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=18" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
Content-Type: text/plain; charset="iso-8859-1"<br>
<br>
Caro Andre<br>
<br>
Segue abaixo a resposta. Acho porem que existem alguns fóruns mais apropriados
para perguntas desse tipo. Procure fóruns sobre matemática, ou calculo
especificamente. Acho que essa lista e mais adequada para questões envolvendo
programação estatística.<br>
<br>
Use integração por partes (no exemplo q estou enviando, não inclui a constante
na integral indefinida pq vc tem uma integral definida. Portanto, apenas inclua
os limites de integração no exemplo).<br>
<br>
<br>
<br>
<br>
Abracos<br>
<br>
<br>
-- <br>
Diogo Ferrari<br>
PhD Student<br>
University of Michigan, Ann Arbor<br>
Department of Political Science<br>
E-mail: <a href="http://user/SendEmail.jtp?type=node&node=4663312&i=19" rel="nofollow" link="external" target="_blank">[hidden email]</a> <mailto:<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=20" rel="nofollow" link="external" target="_blank">[hidden email]</a>><br>
<br>
Open Source! Use R! Use Linux!<br>
<br>
"A vida é a arte de tirar conclusões suficientes de dados
insuficientes"<br>
<br>
<br>
<u></u><u></u></span></p>
<div>
<p class="MsoNormal" style="background:white"><span style="font-family:"Helvetica","sans-serif";color:black"><br>
> On Oct 19, 2014, at 3:23 PM, <a href="http://user/SendEmail.jtp?type=node&node=4663312&i=21" rel="nofollow" link="external" target="_blank">[hidden email]</a> wrote:<br>
> <br>
> Aonde estou errando nessa integral, o caminho não seria esse?<br>
> <br>
> Att.<br>
> André<br>
> <Integral.jpg>_______________________________________________<br>
> R-br mailing list<br>
> <a href="http://user/SendEmail.jtp?type=node&node=4663312&i=22" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
> <a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
> Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo
reproduzível.<u></u><u></u></span></p>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt;background:white"><span style="font-family:"Helvetica","sans-serif";color:black"><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/f9144d27/attachment.html" rel="nofollow" link="external" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/f9144d27/attachment.html</a>><br>
-------------- Próxima Parte ----------<br>
Um anexo não-texto foi limpo...<br>
Nome: latexit-drag.pdf<br>
Tipo: application/pdf<br>
Tamanho: 59430 bytes<br>
Descrição: não disponível<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/f9144d27/attachment.pdf" rel="nofollow" link="external" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/f9144d27/attachment.pdf</a>><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/f9144d27/attachment-0001.html" rel="nofollow" link="external" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141019/f9144d27/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Legenda do Digest<br>
<br>
_______________________________________________<br>
R-br mailing list<br>
<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=23" rel="nofollow" link="external" target="_blank">[hidden email]</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
------------------------------<br>
<br>
Fim da Digest R-br, volume 46, assunto 20<br>
*****************************************<u></u><u></u></span></p>
</div>
</div>
</div>
</div>
</div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
_______________________________________________ <br>
R-br mailing list <br>
<a href="http://user/SendEmail.jtp?type=node&node=4663312&i=24" rel="nofollow" link="external" target="_blank">[hidden email]</a> <br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo
reproduzível. <u></u><u></u></p>
<div class="MsoNormal" align="center" style="text-align:center">
<hr size="1" width="100%" noshade style="color:#cccccc" align="center">
</div>
<div>
<div>
<p class="MsoNormal"><b><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif";color:#444444">If you reply to this email, your message will be added to the
discussion below:<u></u><u></u></span></b></p>
</div>
<p class="MsoNormal"><span style="font-size:9.0pt;font-family:"Tahoma","sans-serif";color:#444444"><a href="http://r-br.2285057.n4.nabble.com/R-br-skater-e-deldir-tp4663303p4663312.html" rel="nofollow" link="external" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-skater-e-deldir-tp4663303p4663312.html</a>
<u></u><u></u></span></p>
</div>
<div>
<div>
<div style="margin-top:4.8pt">
<p class="MsoNormal" style="line-height:18.0pt"><span style="font-size:8.5pt;font-family:"Tahoma","sans-serif";color:#666666">To unsubscribe from R-br, <a rel="nofollow" link="external">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" link="external" target="_blank"><span style="font-size:7.0pt;font-family:"Times New Roman","serif"">NAML</span></a>
<u></u><u></u></span></p>
</div>
</div>
</div>
</div></div></div><div><div class="h5">
<p class="MsoNormal"><br>
<br clear="all">
<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal">-- <br>
Hélio Gallo Rocha<br>
IFSULDEMINAS - Câmpus Muzambinho<u></u><u></u></p>
</div></div></div>
</div>
<br><br>
<hr style="border:none;color:#909090;background-color:#b0b0b0;min-height:1px;width:99%">
<table style="border-collapse:collapse;border:none">
<tbody><tr>
<td style="border:none;padding:0px 15px 0px 8px">
<a href="http://www.avast.com/" rel="nofollow" link="external" target="_blank">
<img border="0" src="http://static.avast.com/emails/avast-mail-stamp.png">
</a>
</td>
<td>
<p style="color:#3d4d5a;font-family:"Calibri","Verdana","Arial","Helvetica";font-size:12pt">
Este email está limpo de vírus e malwares porque a proteção do <a href="http://www.avast.com/" rel="nofollow" link="external" target="_blank">avast! Antivírus</a> está ativa.
</p>
</td>
</tr>
</tbody></table><span class="">
<br>
<br>_______________________________________________
<br>R-br mailing list
<br><a href="http://user/SendEmail.jtp?type=node&node=4663338&i=1" rel="nofollow" link="external" target="_blank">[hidden email]</a>
<br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="nofollow" link="external" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" rel="nofollow" link="external" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.
<br>
<br>
<hr noshade size="1" color="#cccccc">
</span><div style="color:#444;font:12px tahoma,geneva,helvetica,arial,sans-serif"><span class="">
<div style="font-weight:bold">If you reply to this email, your message will be added to the discussion below:</div>
</span><a href="http://r-br.2285057.n4.nabble.com/R-br-skater-e-deldir-tp4663303p4663338.html" target="_blank">http://r-br.2285057.n4.nabble.com/R-br-skater-e-deldir-tp4663303p4663338.html</a>
</div><div class="HOEnZb"><div class="h5">
<div style="color:#666;font:11px tahoma,geneva,helvetica,arial,sans-serif;margin-top:.4em;line-height:1.5em">
To unsubscribe from R-br, <a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=3357982&code=aGVsaW9nYWxsb3JvY2hhQGdtYWlsLmNvbXwzMzU3OTgyfC0xMzQ3NTkwMDY4" target="_blank">click here</a>.<br>
<a href="http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml" rel="nofollow" style="font:9px serif" target="_blank">NAML</a>
</div></div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>Hélio Gallo Rocha<br>IFSULDEMINAS - Câmpus Muzambinho<br>
</div>