Clase 3

Patron de puntos

set.seed(164)
x=runif(100,0,1)
y=runif(100,0,1)
plot(x,y,col="darkblue", pch=19)

colores=ifelse(x^0.5+y/(x+y)<1, "red", "darkgreen")
par(mar=c(5.1,4.1,4.1,8.1), xpd=T)
plot(x,y,col=colores, pch=19)
grid(10,10, col="lightblue")
legend("topright", legend=c("sana", "enferma"), col=colores, title="Estado", text.font = 4, bg="lightblue"
       ,inset=c(-0.4,0), fill=colores)

#################
###Lote 1 tenia 450 plantas y el viento tumbo 40
##Lote 2 tenia 420 plantas y el viento tumbo 30
## Lote 3 tenia 400 plantas y el viento tumbo 25
(p1=(40/450)*100)

## [1] 8.888889

(p2=(30/420)*100)

## [1] 7.142857

(p3=(25/400)*100)

## [1] 6.25

#Porporcion promedio bien 
(40+30+25)/(450+420+400)*100

## [1] 7.480315

#Proporcion promedio mal 
mean(c(p1,p2,p3))

## [1] 7.427249

Libreria Spatstat

library(digest)
library(spatstat)

## Loading required package: spatstat.data

## Loading required package: nlme

## Loading required package: rpart

## 
## spatstat 1.59-0       (nickname: 'J'ai omis les oeufs de caille') 
## For an introduction to spatstat, type 'beginner'

set.seed(1964)
nclust<- function(x0,y0,radius, n){
  return(runifdisc(n,radius, centre=c(x0,y0)))
}
plot(rPoissonCluster(5,0.2,nclust, radius=0.2, n=100))

plot(rPoissonCluster(7,0.2,nclust, radius=0.1, n=100))

dat<- matrix(stats::rnorm(2000), ncol=2)
ch<- chull(dat)
coords <- dat[c(ch, ch[1]),]
plot(dat, pch=19, cex=0.5)
lines(coords, col="red")

library(spatstat)
data("swedishpines")
(x<- swedishpines)

## Planar point pattern: 71 points
## window: rectangle = [0, 96] x [0, 100] units (one unit = 0.1 metres)

# windows(7,7) # Para sacar una ventana afuera
plot(x)

class(x)

## [1] "ppp"

###

set.seed(1964)
x<- runif(40,0,1)
y<- runif(40,0,1)
plot(x,y,pch=19, cex= 0.5)

dataxy=ppp(x,y)
Q=t(quadratcount(dataxy, nx=22, ny=22))
vc=data.frame(Q)
colores=c(ifelse(vc$Freq==0, "darkred", "darkblue"))
plot(t(Q), col=colores, cex=0.8)

cua=1


#X<- clickppp(15)
#plot(X)
#####