Diplomová práce – Prostorová neurčitost geodat v ... - kvhem

Příloha 3Skript pro výpočet pěti modelů logistické regrese z připravených dat jednégeodatabáze.___________________________________________________________________data = read.table ('CLC90.txt', dec = '.')datadata = data.frame (data)names (data)= c("kvadrat", "spol_hran", "plocha_les", "obvod_les", "plocha_pole", "obvod_pole", "vyskyt")data#generování matice pro jednotlivé prediktory a výpočty regrese#prediktor Apočet=628pm1=c(rep(0,počet))pm2=c(rep(0,počet))for (i in 1:počet){pm1[i]=data[i,2]pm2[i]=data[i,7]}pm1pm2prediktor_A=matrix(c(pm1,pm2),ncol=2,byrow=F)prediktor_A#pojmenování sloupcůprediktor_A = data.frame (prediktor_A)names (prediktor_A)= c("spol_hran", "vyskyt")#regresevysledek_A = glm (vyskyt ~ spol_hran, family = binomial, data = prediktor_A)summary (vysledek_A)confint (vysledek_A)#prediktor Bpočet=628pm1=c(rep(0,počet))pm2=c(rep(0,počet))for (i in 1:počet){pm1[i]=data[i,3]pm2[i]=data[i,7]}pm1pm2prediktor_B=matrix(c(pm1,pm2),ncol=2,byrow=F)prediktor_B#pojmenování sloupcůprediktor_B = data.frame (prediktor_B)names (prediktor_B)= c("plocha_les", "vyskyt")#regresevysledek_B = glm (vyskyt ~ plocha_les, family = binomial, data = prediktor_B)summary (vysledek_B)confint (vysledek_B)#prediktor Cpočet=628pm1=c(rep(0,počet))64