Commit e703bd34 authored by Boulangeat Isabelle's avatar Boulangeat Isabelle
Browse files

hauteur log

parent 3a946e7e
......@@ -200,7 +200,7 @@ vars = c("alti", "slope", "lf", "northing", "easting", "snowfree_time", "cumgdd0
dat = na.omit(releves_all[-which(releves_all$ref_typoveg %in% c("HUM", "LAND", "NA", "EBOU")),c(vars, "hmean")])
dat$lf = as.factor(dat$lf)
dim(dat)
mod = randomForest(hmean ~ ., data = dat, ntree = 100, mtry=16)
mod = randomForest(log(hmean+1) ~ ., data = dat, ntree = 100, mtry=16)
mod
varImpPlot(mod)
......@@ -249,10 +249,10 @@ vars = c("alti", "slope", "lf", "northing", "easting", "snowfree_time", "cumgdd0
dat = na.omit(releves_all[which(releves_all$ref_typoveg =="PROD"),c(vars, "hmean")])
# dat$lf = as.factor(dat$lf)
dim(dat)
mod = randomForest(hmean ~., data = dat, ntree = 100, mtry=16)
mod = randomForest(log(hmean+1) ~., data = dat, ntree = 100, mtry=16)
mod
varImpPlot(mod)
varImpPlot(mod, main="PROD")
```
```{r ,fig=TRUE}
......@@ -267,10 +267,10 @@ vars = c("alti", "slope", "lf", "northing", "easting", "snowfree_time", "cumgdd0
dat = na.omit(releves_all[which(releves_all$ref_typoveg =="ALP"),c(vars, "hmean")])
# dat$lf = as.factor(dat$lf)
dim(dat)
mod = randomForest(hmean ~., data = dat, ntree = 100, mtry=16)
mod = randomForest(log(hmean+1) ~., data = dat, ntree = 100, mtry=16)
mod
varImpPlot(mod)
varImpPlot(mod, main="ALP")
```
```{r ,fig=TRUE}
......@@ -285,10 +285,10 @@ vars = c("alti", "slope", "lf", "northing", "easting", "snowfree_time", "cumgdd0
dat = na.omit(releves_all[which(releves_all$ref_typoveg =="SUB"),c(vars, "hmean")])
# dat$lf = as.factor(dat$lf)
dim(dat)
mod = randomForest(hmean ~., data = dat, ntree = 100, mtry=16)
mod = randomForest(log(hmean+1) ~., data = dat, ntree = 100, mtry=16)
mod
varImpPlot(mod)
varImpPlot(mod, main="SUB")
```
```{r ,fig=TRUE}
......@@ -306,28 +306,25 @@ dim(dat)
mod = randomForest(hmean ~., data = dat, ntree = 100, mtry=16)
mod
varImpPlot(mod)
varImpPlot(mod, main="NIV")
```
### Interannuel
```{r ,fig=TRUE}
vars = c("alti", "slope", "lf", "northing", "easting", "snowfree_time", "cumgdd0", "cumgdd60d", "albedo.gdd300",
"albedo.gdd600", "albedo.gdd900","temperature.gdd300", "frost.gdd300", "frost.gdd600", "frost.gdd900",
vars = c("id_site", "cumgdd0", "albedo.gdd300","temperature.gdd300", "frost.gdd300",
"gddspeed.gdd300", "gddspeed.gdd600", "gddspeed.gdd900", "radiations.gdd300",
"radiations.gdd600", "radiations.gdd900", "rainfall.gdd300", "rainfall.gdd600",
"rainfall.gdd900", "snowdays.gdd300", "snowdays.gdd600", "snowdays.gdd900")
#vars = c("slope", "lf", "northing", "easting", "snowfree_time", "cumgdd60d", "ref_typoveg")
"rainfall.gdd300")
# dat = na.omit(releves_all[-which(releves_all$ref_typoveg %in% c("HUM", "LAND", "NA", "EBOU")),c(vars, "hmean")])
sort(table(releves_all$id_site), dec=TRUE)
dat = na.omit(releves_all[-which(releves_all$ref_typoveg %in% c("HUM", "LAND", "NA", "EBOU")),c(vars, "hmean")])
# dat$lf = as.factor(dat$lf)
dim(dat)
dat
mod = randomForest(hmean ~., data = dat, ntree = 100, mtry=16)
mod = randomForest(log(hmean+1) ~., data = dat, ntree = 100, mtry=16)
mod
varImpPlot(mod)
varImpPlot(mod, main="interannuel")
```
This diff is collapsed.
This diff is collapsed.
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment