performance icon indicating copy to clipboard operation
performance copied to clipboard

Published 20 hours ago • verified publisher icon easystats

R2 for spaMM

Open sskik opened this issue 3 years ago • 2 comments

Any chance you might add an r2 function for models produced using the spaMM package?

sskik avatar Apr 19 '22 21:04 sskik

Can you give a reproducible example of a spaMM model?

bwiernik avatar Apr 19 '22 21:04 bwiernik

g1<- data.frame( y= c(0.02120366393, 0.08525051821, -0.3703122665, -0.202579954, 0.2218212876, 0.2135410943, 0, 0.3289244879, -0.0146003462, -0.493927186, -0.3290933416, -0.2993991624, 0.5144957554, -0.05736474925, 0.2233112488, -0.1245225503, -0.1579905011, -0.1779939163, -0.2126805699, 0, -0.3538660498, 0.02826167095, 0.1237514472, 0.152956058, -0.01709264742, -0.2626128657, 0.07158710943, -0.02952134684, 0.0179120072, 0.2037445373, -0.22325522, -0.2492759943, 0.1031827878, -0.09729355931, -0.1052309677, 0.02120366393, 0.03887920042, 0.06925847715, 0.05018529447, 0.1431880214, 0.5354681677, 0.1410801544, 0.04205208877, 0.01691842996, 0.01106621469, -0.1789703248, -0.3335735585, 0.2993991624, -0.07731539427, 0, 0, 0.02826167095, 0, 0, 0.5125379437, 0.3361510643, 0, 0.195369854, 0.2604256263, 0.01182113897, -0.2443016845, 0.01411133998, 0.1214384639, 0.06496094919, -0.3091884982, 0.2130209406, 0.3773077409, -0.2933396643, -0.1947309068, -0.2019040317, -0.1284666006, 0.2416127038, 0.1284666006, 0.2453723334, 0.3103276864, -0.106345061, -0.1501215125, -0.3269394685, 0.1569564821, 0.2376142253, 0.1196614184, 0.6150317987, 0.3069942211, 0.2788124253, 0.007070714285, 0.007070714285, 0.007070714285, 0.007070714285, 0.007070714285), x = c(-0.300210406617295, 0.840223028276454, 0.114492660616796, 0.114492660616796, 0.114492660616796, 0.114492660616796, 0.351465841893419, 0.0404385414678508, -0.0632372253406717, 0.603249846999831, 0.514384904021097, 0.514384904021097, -0.0928588730002498, -0.0928588730002498, -0.0928588730002498, -0.240967111298139, -0.240967111298139, -0.240967111298139, -0.463129468744974, -0.463129468744974, -0.477940292574762, -0.596426883213074, -0.492751116404551, -0.492751116404551, -0.492751116404551, -0.566805235553496, -0.581616059383285, -0.581616059383285, -0.581616059383285, -0.359453701936451, -0.0336155776810939, -0.0336155776810939, -0.0188047538513051, -0.0188047538513051, -0.0188047538513051, -0.0188047538513051, -0.403886173425818, -0.403886173425818, -0.403886173425818, -0.403886173425818, -0.418696997255607, -0.418696997255607, -0.418696997255607, -0.418696997255607, 0.736547261467931, 0.736547261467931, 0.692114789978565, 0.603249846999831, -0.0484264015108828, -0.0484264015108828, 0.010816893808273, 0.144114308276374, 0.144114308276374, 0.144114308276374, -0.848210888319486, -0.848210888319486, 0.114492660616796, 0.395898313382786, 0.425519961042364, 0.395898313382786, 0.395898313382786, 0.395898313382786, 0.395898313382786, 0.425519961042364, 0.440330784872152, 0.33665501806363, 0.351465841893419, 0.484763256361519, 0.484763256361519, 0.514384904021097, 0.514384904021097, 0.514384904021097, 0.247790075084896, 0.46995243253173, 0.46995243253173, 0.46995243253173, 0.395898313382786, 0.61806067082962, 0.46995243253173, 0.395898313382786, 0.395898313382786, 0.677303966148776, 0.381087489552997, 0.351465841893419, -1.15923818874505, -1.15923818874505, -1.15923818874505, -1.15923818874505, -1.08518406959611), region = c('A','J','I','I','I','I','H','I','I','H','H','H','I','I','I','I','I','I','G','G','G','A','G','G','G','A','A','A','A','G','I','I','I','I','I','I','A','A','A','A','I','I','I','I','F','F','H','H','E','E','E','E','E','E','A','A','E','E','E','E','E','E','E','E','E','E','D','D','D','H','H','H','E','E','E','E','E','F','E','E','E','F','C','C','B','B','B','B','B'), LAT = c(42.153189, 37.297799, 40.191259, 40.191259, 40.191259, 40.191259, 39.073655, 40.391729, 41.333094, 38.720123, 38.739944, 38.739944, 41.050461, 41.050461, 41.050461, 41.73519, 41.73519, 41.319117, 42.143647, 42.143647, 42.10483, 42.726538, 42.254366, 42.254366, 42.254366, 42.682125, 42.645461, 42.645461, 42.645461, 41.807836, 40.806629, 40.806629, 40.780999, 40.780999, 40.780999, 40.780999, 42.265165, 42.265165, 42.265165, 42.265165, 42.08242, 42.08242, 42.08242, 42.08242, 37.007328, 37.007328, 38.367151, 38.519429, 43.702597, 43.702597, 43.224929, 42.836492, 42.836492, 42.836492, 43.810882, 43.810882, 42.82175, 39.655053, 39.738175, 39.788689, 39.788689, 39.788689, 39.788689, 39.775821, 39.872721, 40.371103, 39.505718, 39.282921, 39.282921, 39.169881, 39.169881, 39.169881, 40.407923, 40.398578, 40.398578, 40.398578, 39.97389, 37.735524, 38.898013, 39.776267, 39.776267, 37.625934, 38.543684, 39.444084, 45.296331, 45.296331, 45.296331, 45.296331, 44.807061), LON = c(-89.22871, -88.609806, -90.653806, -90.653806, -90.653806, -90.653806, -88.905039, -90.382165, -88.674397, -85.455282, -86.377949, -86.377949, -91.735715, -91.735715, -91.735715, -91.52571, -91.52571, -94.032623, -93.974468, -93.974468, -93.923812, -92.445311, -93.983794, -93.983794, -93.983794, -91.277902, -91.309445, -91.309445, -91.309445, -93.830591, -91.349861, -91.349861, -91.290163, -91.290163, -91.290163, -91.290163, -90.5374, -90.5374, -90.5374, -90.5374, -90.738654, -90.738654, -90.738654, -90.738654, -82.886132, -82.886132, -85.57063, -84.839868, -84.429404, -84.429404, -84.097254, -82.715595, -82.715595, -82.715595, -91.762898, -91.762898, -76.649504, -81.949676, -81.781398, -81.870826, -81.870826, -81.870826, -81.870826, -81.339806, -80.83136, -81.708198, -84.729517, -84.582095, -84.582095, -84.690264, -84.690264, -84.690264, -82.188372, -80.117196, -80.117196, -80.117196, -80.25197, -81.641589, -80.5776, -80.751696, -80.751696, -81.945147, -79.346929, -79.654829, -91.055685, -91.055685, -91.055685, -91.055685, -90.776971))

f<- y ~ x + (1 | region) + Matern(1|LON + LAT) library(spaMM) m <- fitme(formula = f, data = g1, family = gaussian())

sskik avatar Apr 19 '22 22:04 sskik