2 votos

¿Cómo se llama exactamente el tipo de análisis de regresión en el que se intenta ver si el modelo es significativo sobre *múltiples* valores de inicio/final?

Un ejemplo es este: http://www.reddit.com/r/askscience/comments/ine4x/regarding_the_recent_lapse_of_global_warming_in/c2554al

Seguro que está relacionado con las estadísticas robustas. Pero estoy seguro de que hay una etiqueta más específica que eso.

\==

Bien, he intentado realizar análisis de regresión sobre los conjuntos de datos del calentamiento global de los últimos 10 años. Lo que quería demostrar era que no había una tendencia significativa al calentamiento en los últimos 10 años. Sin embargo, como hay ciertos años anómalos, quería que el argumento fuera más convincente, así que intenté diferentes valores del año de inicio. Mi intención era demostrar que la premisa del artículo (http://www.physorg.com/news/2011-07-global-linked-sulfur-china.html ) era básicamente correcta. No soy un negacionista del cambio climático ni mucho menos, ya que yo mismo trabajo con científicos del clima (odio a los negacionistas del cambio climático tanto como a cualquier otro científico), sin embargo, la cuestión es que China ha liberado tanto azufre a la atmósfera que el cambio climático prácticamente se ha estancado durante los últimos 10 años (dicho esto, no tengo ninguna duda de que se reanudará una vez que China reduzca la cantidad de carbón que contiene azufre que utiliza)

Utilizaré el conjunto de datos de http://www.ncdc.noaa.gov/cmb-faq/anomalies.php - y en particular - el ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/annual.land_ocean.90S.90N.df_1901-2000mean.dat de datos (las anomalías terrestres + oceánicas), que es más riguroso científicamente

\==

De todos modos, al realizar análisis de regresión sobre estos datos...

Para el periodo 1997-2010, para el coeficiente de x, obtengo un valor de 0,00722 grados por año. Esto no es estadísticamente significativo, ya que el valor p es de 0,13

1998-2010: 0,00515 grados por año/valor de 0,19. No es estadísticamente significativo.

Para el periodo 1999-2010, para el coeficiente de x, obtengo un valor de 0,01182 grados por año. Esto apenas alcanza la significación estadística, ya que el valor p es de 0,0457

2000-2010: 0,0089 grados por año/valor de 0,17. No es estadísticamente significativo

2001-2010: 0,0016 grados por año/valor de 0,759. No es estadísticamente significativo

2002-2010: -0,00167 grados por año/valor de 0,784. Pero aquí tenemos una pendiente negativa

2003-2010: -0,00141 grados por año/valor de 0,86. De nuevo, pendiente negativa

4voto

Boris Tsirelson Puntos 191

Sigo pensando que se busca una regresión rodante. Ejemplo:

Data <- read.table("ftp://ftp.ncdc.noaa.gov/pub/data/anomalies/annual.land_ocean.90S.90N.df_1901-2000mean.dat")
colnames(Data) <- c('Year','Anom')
Data <- Data[Data$Anom!=-999,]

require(xts)
Data <- xts(Data,ISOdate(Data$Year,1,1))
indexFormat(Data) <- "%Y"  

GlobalWarming <- rollapply(Data, width = 10,
  FUN = function(z) coef(lm(Anom ~ Year, data = as.data.frame(z))),
  by.column = FALSE, align = "right")
plot(GlobalWarming)

Rolling Regression Results

Como se puede ver, la solidez de la tendencia a 10 años ha disminuido en los últimos años (y no es la primera vez que ocurre)

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X