<div dir="ltr"><div>Hi Mitchell,</div><div><br></div><div>looking at your command there are several things to fix. Parameter columns should have only the predictors, so for example remove county_rast and urban_change_01_16. Parameter developed_column should be the response variable, so I suppose urban_change_01_16 is the one to use here. Lastly, without the dredge mode only use combination of predictors that are not correlated and don't represent the same info, for example, do not include dist_toh20_2001_km,dist_toh20_2004_km,dist_toh20_2006,dist_toh20_2008_km,dist_toh20_2011_km,dist_toh20_2013_km,dist_toh20_2016_km, since these all represent the same variable (I suppose). Without the dredge mode, all provided predictors will be used in the model, with the dredge mode, it will try different combinations and return the best model.</div><div><br></div><div>Note r.futures.potential is basically just a wrapper around R glmer (lme4) function, so if you are familiar with R, you can experiment with finding a model directly in R if you need more flexibility etc.</div><div><br></div><div>Please keep the conversation on this mailing list which is archived, so that others can add to it or benefit from it later.</div><div><br></div><div>Best,</div><div>Anna</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Thu, Feb 11, 2021 at 2:12 PM Mitchell Meads <<a href="mailto:mitchell.meads@tamu.edu" target="_blank">mitchell.meads@tamu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi Anna,<div><br></div><div>I've revisited the r.futures.potential submodel for the futures modeling and I still can't get the r.futures.potential submodel to run. I've rescaled my variables so they all should be on a similar scale but the program is telling me that the variables are still not in appropriate scale along with some other error messages that I have copied below. Please let me know what you think I'm doing wrong.</div><div><br></div><div>(Thu Feb 11 13:11:21 2021) <br>r.futures.potential --overwrite input=sampling_f_km@PERMANENT output=potential_f_km.csv columns=devclip_2016,block_rast,county_rast,devpress_0_5_2016,dist_to_evac_km,dist_to_pp_km,dist_toh20_2001_km,dist_toh20_2004_km,dist_toh20_2006,dist_toh20_2008_km,dist_toh20_2011_km,dist_toh20_2013_km,dist_toh20_2016_km,per_forest_2001,per_forest_2004,per_forest_2006,per_forest_2008,per_forest_2011,per_forest_2013,per_forest_2016,urban_change_01_16 developed_column=devclip_2016 subregions_column=county_rast min_variables=5 max_variables=10<br>Computing model...<br>WARNING: fixed-effect model matrix is rank deficient so dropping 2 columns / coefficients<br>Error in pwrssUpdate(pp, resp, tol = tolPwrss, GQmat = GQmat, compDev = compDev, :<br>pwrssUpdate did not converge in (maxit) iterations<br>Calls: glmer ... <Anonymous> -> <Anonymous> -> stopifnot -> fn -> pwrssUpdate<br>In addition: Warning message:<br>Some predictor variables are on very different scales: consider rescaling<br>Execution halted<br>WARNING: fixed-effect model matrix is rank deficient so dropping 2 columns / coefficients<br>Error in pwrssUpdate(pp, resp, tol = tolPwrss, GQmat = GQmat, compDev = compDev, :<br>pwrssUpdate did not converge in (maxit) iterations<br>Calls: glmer ... <Anonymous> -> <Anonymous> -> stopifnot -> fn -> pwrssUpdate<br>In addition: Warning message:<br>Some predictor variables are on very different scales: consider rescaling<br>Execution halted<br>ERROR: Running R script failed, check messages above<br>(Thu Feb 11 13:11:25 2021) Command finished (3 sec) <br></div><div><br clear="all"><div><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><font face="arial, sans-serif">Cheers,</font></div><div><font face="arial, sans-serif">Mitchell Meads, PhD Student</font></div><div><font face="arial, sans-serif">Texas A&M University at Galveston</font></div><div><font face="arial, sans-serif">Marine & Coastal Management Sciences; <a href="http://www.tamug.edu/ctbs/" target="_blank">CTBS</a></font></div><div><font face="arial, sans-serif">Email: <a href="mailto:mitchell.meads@tamu.edu" target="_blank">mitchell.meads@tamu.edu</a></font></div><div style="font-family:"trebuchet ms",sans-serif"><br></div><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><span style="font-size:14px"><font face="trebuchet ms, sans-serif"><br></font></span></div><div></div></div></div></div></div></div></div></div></div></div></div><div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Jan 23, 2021 at 9:54 PM Anna Petrášová <<a href="mailto:kratochanna@gmail.com" target="_blank">kratochanna@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Hi Mitchell,</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Fri, Jan 22, 2021 at 3:23 PM Mitchell Meads <<a href="mailto:mitchell.meads@tamu.edu" target="_blank">mitchell.meads@tamu.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hi all,<div><br></div><div>I'm getting an error when running r.futures.demand that says "Number of development raster maps doesn't not correspond to the number of observed times". I thought that a fix for this would be to make sure my input population csv had years that matched the input raster maps in the first column but this didn't change the error. If you could let me know what you think I'm doing wrong, that would be greatly appreciated!</div></div></blockquote><div><br></div><div>The error means the number of raster maps in the development parameter does not match the number of years in the observed_population CSV. So for example, if you specify development=urban_2001,urban_2006,urban_2011 (3 maps) then the observed_population file should look like this (header + 3 lines):</div><div><br></div><div>year,37037,37063,...<br>2001,19860,10980,...<br>2006,20760,12660,...<br>2011,21070,13090,...<br></div><div><br></div><div>If you can't see any obvious mistake with your data, I would have to see the r.futures.demand command you are trying to execute and the input CSV file.</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><br></div><div><br clear="all"><div><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><font face="arial, sans-serif">Cheers,</font></div><div><font face="arial, sans-serif">Mitchell</font></div><div style="font-family:"trebuchet ms",sans-serif"><br></div><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div><span style="font-size:14px"><font face="trebuchet ms, sans-serif"><br></font></span></div><div></div></div></div></div></div></div></div></div></div></div></div><div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://urldefense.com/v3/__https://lists.osgeo.org/mailman/listinfo/grass-user__;!!KwNVnqRv!S-XTr0kWrbJ704p6CIighh6PUpDBJchKEjSQMopVh0bXjHxBYSqNzjaTNlgopVQ1ulKp$" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote></div></div>
</blockquote></div>
</blockquote></div></div>