Opened 10 years ago

Last modified 5 years ago

#2376 new enhancement

Add p-value for regression slope in r series

Reported by: mastho Owned by: grass-dev@…
Priority: normal Milestone: 7.6.2
Component: Raster Version: svn-releasebranch70
Keywords: r.series, p-value regression Cc: pvanbosgeo
CPU: Unspecified Platform: Unspecified

Description

The calculation of t-values has been added to r.series in r53742 (method=tvalue). It would be nice to have also the corresponding p-values for the regression as method in r.series.

Also discussed here: http://osgeo-org.1560.x6.nabble.com/Calculate-p-value-for-regression-slope-in-r-series-td5014228.html

Change History (12)

comment:1 by neteler, 10 years ago

Keywords: r.series p-value regression → r.series, p-value regression
Milestone: 7.0.0
Version: unspecifiedsvn-releasebranch70

comment:2 by pvanbosgeo, 10 years ago

Cc: pvanbosgeo added

You can calculate the p-value using the scipy library (from http://stackoverflow.com/questions/17559897/python-p-value-from-t-statistic). For example, if tt is the t value

import scipy.stats as stats
import numpy as np
tt = 2.4
n = 10
pval = stats.t.sf(np.abs(tt), n-1)*2  (two-sided p-value)

Perhaps this can be used either to create a script or in r.series?

in reply to:  2 comment:3 by mmetz, 10 years ago

Replying to pvanbosgeo:

You can calculate the p-value using the [python] scipy library

Considering that you want to calculate a p value for each single pixel, the chances that you get randomly some very small p values increase with the number of pixels. You would need to correct for that, with Bonferroni or some method that escapes the iron claw of Bonferroni. In any case, all your p values will most probably become not significant after any method of correction for multiple testing.

Therefore it is more appropriate to make prior assumptions about slope, intercept, and effect size, then judge the results according to these prior assumptions.

comment:4 by martinl, 8 years ago

Milestone: 7.0.07.0.5

comment:5 by martinl, 8 years ago

Milestone: 7.0.57.3.0

comment:6 by martinl, 8 years ago

Milestone: 7.3.07.4.0

Milestone renamed

comment:7 by dvictori, 6 years ago

You can calculate regression p-value using R and the raster package. Here is a StackOverflow post about it:

https://stackoverflow.com/questions/20262999/how-to-output-regression-summarye-g-p-value-and-coeff-into-a-rasterbrick

comment:8 by neteler, 6 years ago

Milestone: 7.4.07.4.1

Ticket retargeted after milestone closed

comment:9 by neteler, 6 years ago

Milestone: 7.4.17.4.2

comment:10 by martinl, 6 years ago

Milestone: 7.4.27.6.0

All enhancement tickets should be assigned to 7.6 milestone.

comment:11 by martinl, 5 years ago

Milestone: 7.6.07.6.1

Ticket retargeted after milestone closed

comment:12 by martinl, 5 years ago

Milestone: 7.6.17.6.2

Ticket retargeted after milestone closed

Note: See TracTickets for help on using tickets.