8 votes

Exploratory Spatial Regression: How to limit 60 variable candidates?

It seems you are doing research. Here are some advices: Ask the domain experts: if you ask from an expert about the exploratory variables, he/she can easily filter many of the variables due to their ...
Farid Cheraghi's user avatar
7 votes

How to conduct logistic regression between two rasters

To eleborate on Spacedman's answer (and to make it more clear where the problem may reside: Can you combine the data? s <- stack(sinks, tpi) If so, you should be able to do: v <- data.frame(...
Robert Hijmans's user avatar
7 votes
Accepted

Linear regresion at each raster pixel's value to predict future value (in R language)

Starting with your raster stack s: > s class : RasterStack dimensions : 15, 10, 150, 6 (nrow, ncol, ncell, nlayers) I'll show how to fit and predict in various ways. I'm going to try to ...
Spacedman's user avatar
  • 63.9k
6 votes

Difficulties applying zyp.sen() to a RasterBrick object in R

The problem is not entirely yours. zyp.sen is not well written: it fails unless you guess correctly concerning the syntax it is hoping for. (Whether that is its only problem is a question I have not ...
whuber's user avatar
  • 69.8k
6 votes

Exploratory Spatial Regression: How to limit 60 variable candidates?

Stepwise regression is generally an unwise approach, even though commonly used and published. For supported inference, evaluating model support and model fit, the canned ArcGIS tools are not a good ...
Jeffrey Evans's user avatar
6 votes

Calculating R-Squared Using the 'Residuals' Band from Google Earth Engine's linearRegression() Reducer

Thought I would add the answer to my question here for the next person who has this problem. GEE developer Nick Clinton presented his solution to it in the Google Earth Developer Forum (https://groups....
Elijah Tangenberg's user avatar
5 votes

Access to attribute in rasterlayer in R

You can access this table via levels: > levels(r) [[1]] ID OID Value Count LC_Code Land_Cover Area_Ha 1 0 0 0 25615512 255 No Data 2305400.00 2 1 1 ...
Spacedman's user avatar
  • 63.9k
5 votes
Accepted

Interpreting result of Geographically Weighted Regression (GWR)?

The parameter estimates are essentially the slopes of regression lines that fit your data averaged at that location. So if the parameter is positive then there are more auto thefts in places with more ...
Spacedman's user avatar
  • 63.9k
5 votes

Understanding regression function output in ArcGIS Kriging?

You are correct about value you referred to, the residuals: An excerpt from Esri's Regression analysis basics Residuals: These are the unexplained portion of the dependent variable, represented in ...
whyzar's user avatar
  • 12.1k
4 votes
Accepted

Getting Negative values in Regression based Biomass model!

The crux is most likely in the fact that you are using a MLR, as this method can indeed result in values below 0, while a negative biomass obviously does not exist in real life. As for the ...
LMB's user avatar
  • 1,166
4 votes

How to conduct logistic regression between two rasters

If your rasters have the same grid size and position, and the sinks raster is 0s and 1s (or TRUEs and FALSEs) then its just like doing a logistic regression any other way in R, but getting the values ...
Spacedman's user avatar
  • 63.9k
4 votes
Accepted

Linear regression between every 3×3 pixels between two rasters using R

I assume that by "lapse rate" you mean the regression slope. You are essentially describing a focal regression and since there is no multivariate version of raster::focal, there is no canned way to do ...
Jeffrey Evans's user avatar
3 votes
Accepted

How do you find the appropriate spatial weight to use for a spatial analysis?

I am going to tackle this because it is likely going to be closed as "too broad" and I doubt that you will get a satisfactory answer if migrated to something like Cross Validated. Although, there is ...
Jeffrey Evans's user avatar
3 votes
Accepted

Moran's I and GWR in QGIS

Well actually ArcMap uses points as well as the tool calculates the centroids first and then performs the Test. (http://edndoc.esri.com/arcobjects/9.2/net/shared/geoprocessing/...
Matte's user avatar
  • 6,235
3 votes
Accepted

Performing per-pixel linear regression over multiple rasters with NoData?

Using technique described here I populated the table of integer grid by values from 3 rasters of interest. I used IS NULL query to populate records with no match (NO DATA) by -999 during data ...
FelixIP's user avatar
  • 23k
3 votes

QGIS Regression Tools using only feature data

have you tried this one? Multiple Linear Regression(Shapes). http://www.saga-gis.org/saga_tool_doc/7.8.0/statistics_regression_12.html Interestingly the documentation refers to something called Trend ...
nr_aus's user avatar
  • 3,535
3 votes
Accepted

Surf3D, example from ggRandomForests: Random Forests for Regression

You need to create a bivariate partial dependency plot first. I believe the function you are after to create the object to pass to plot3D::surf3D is ggRandomForests::partial.rfsrc, ggRandomForests::...
Jeffrey Evans's user avatar
2 votes

Ignoring multi-collinearity in spatial stepwise regression modeling?

I think you must review a multicollinearity before you uses model regression. It is part of assumptions of multiple Linear Regression. And if you dont accomplish with all parameter, you have a bad ...
Dinosca's user avatar
  • 117
2 votes

Google Earth Engine: Evapotranspiration Linear Regression MODIS

The main issue is that you should not use arrays when you don't need to. If I'm interpreting your intention correctly, here it is the better way: var start = '2000-01-01'; var end = '2014-12-31'; ...
Nicholas Clinton's user avatar
2 votes

Linear Regression using Google Earth Engine

you can use ee.Reducer.linearRegression(numX, numY) API,the details search in GEE document.
Mass-Kanf's user avatar
2 votes
Accepted

Apply mann-kendall analysis on stacked raster in R as a function

Always check that your function for calc works before using it: > tsfun(runif(3)) Error in names(r) = c("Z", "Sen", "Old_Sen", "p", "s", "var(s)", "Tau") : 'names' attribute [7] must be the ...
Spacedman's user avatar
  • 63.9k
2 votes

Calculating Bare Soil Line using R?

Here's your model: NIR = αRed + β where NIR and Red correspond to the near-infrared and red bands of the satellite sensor, alpha is the slope, and beta the y intercept. An ordinary least ...
Spacedman's user avatar
  • 63.9k
2 votes

Display R2 in GEE chart

The R2 value is not shown on the chart by default (I don't know if there is a way to change this behaviour). It is displayed if you hover the mouse over the legend field for the line, on the displayed ...
Joris's user avatar
  • 321
2 votes

modifiedmk Mann-Kendall Test Over NC /Raster stack

First, you are using a double bracket where you do not need to and second, it is often better to use a numeric index. library(modifiedmk) x <- c(Nile) Here, if you use a single bracket it ...
Jeffrey Evans's user avatar
2 votes

When to use spatial regression for the proximity analysis

Sounds like an interesting piece of work. A simple statistical solution may be to map linear regression residuals of each house, based on distance vs property value change. (There are a million other ...
nr_aus's user avatar
  • 3,535
2 votes

Interpreting results of GWR - what do the values mean?

I had the same doubt myself, and contrary to what some other people suggested you to do in the answers above, from the source code it is not straightforward what all those concepts mean. However, ...
Perla Rivadeneyra's user avatar
2 votes

Regression in QGIS

In the processing toolbox (Processing > Toolbox) you can find a pile of suitable tools (search for "regression" and "sample"). You could first create a raster from the ...
Jochen Schwarze's user avatar
2 votes

How to calculate the distance of multiple points to a line in QGIS

Your measurements are a point layer, if I understood correctly. Create a new attribute on this layer with Field Calculator with this expression and replace line_layer with the name of the line layer. ...
Babel's user avatar
  • 71.2k
1 vote

How to conduct logistic regression between two rasters

To supplement the answers you already got, I believe that one viable strategy would entail: drawing a number of random sampling points within your study area; drawing them with a subjective minimum ...
NewAtGis's user avatar
  • 1,536
1 vote

Using R to do logistic geographically weighted regression(GWR) prediction

Is the gwmodel$SDF object not populated? The model object should contain a SpatialPixelsDataFrame object with the model estimates. It is very difficult to provide advice when you do not show us what ...
Jeffrey Evans's user avatar

Only top scored, non community-wiki answers of a minimum length are eligible