length(table(PRT332$PUU)) #Loetlen puude arvu PRT-l > #Mitu puud on sinu proovitükil? > PRT332=subset(puud2015,PRT=="332") #Teen eraldi andmestiku PRT332-st > length(table(PRT332$PUU)) #Loetlen puude arvu PRT-l [1] 121 VASTUS: Minu proovitükil on 121 puud Mitu 1. rinde puud on puuliikide kaupa sinu proovitükil? Esita tabel vaid proovitükil esinevate liikide kohta. #Mitu 1. rinde puud on puuliikide kaupa proovitükil. PUUD1=droplevels(subset(PRT332,RIN=="1")) #Teen 1. rinde puude andmestiku, viskan välja mitte vajaliku ehk nullid. table(PUUD1$PL) #Vaatan puid puuliikide kaupa #Ühendan saadud andmed ühte tabelisse P1=with(PUUD1,tapply(D,PL,length)) P1 PRT332RIN=data.frame(PL=names(P1),Rinne1=P1) PRT332RIN write.table(PRT332RIN,"clipboard",sep="t",dec=",",row.names=FALSE) > #Mitu 1. rinde puud on puuliikide kaupa proovitükil. > PUUD1=droplevels(subset(PRT332,RIN=="1")) #Teen 1. rinde puude andmestiku,
Kodutöö: Lineaarne regressioonanalüüs PD <- read.csv("puud15.CSV") # parameeter sep="," ja dec="." PD$d_k<-with(PD, ifelse(d2>0,(d1+d2)/2, d1)) PD.<-subset(PD, prt==642 & aasta==2001) PD.<-droplevels(PD.) plot(h~d_k,data=PD.) PD.H <- subset(PD., h>0 & hv>0) table(PD.H$pl) PD.KU<-subset(PD.H, pl=="KU") par(mar=c(4.5,4.5,1,1)) plot(NULL,xlim=c(0,40),ylim=c(0,25),xlab="diameeter, cm", ylab="kõrgus, m") abline(v=seq(0,40,10),lty=3,col="grey75") abline(h=seq(0,25,5),lty=3,col="grey75") # abijooned points(h~d_k,data=subset(PD.KU),lwd=1) with(subset(PD., pl=="KU"),rug(d_k)) 1. Sirge h=a+b*d M1 <- lm(h~d_k, data=PD.KU) summary(M1) D<-0:40 M1
annaabi.ee 
