Extract random effects from CARBayes output Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?Extract values from raster into points: data loss using extract function from raster packageSpatial weight matrix with boundary effects in R?What to do with the intercept output from regression analsyis?Is there really random distribution of points?Using Random Forest for big data?Predicting Enhanced Vegetation Index using Random Forest in RHow to extract specific values with point coordinates from Kriging interpolations made in R?Plotting/Exporting the Spatial Random Effects in BYM model - CARBayes - in REstimate spatial kernel density from raster and output kernel values in a gridExtract random coordinate points using classes (ranges) of values from a raster image

Is there a kind of relay only consumes power when switching?

Where are Serre’s lectures at Collège de France to be found?

Is it fair for a professor to grade us on the possession of past papers?

What does "lightly crushed" mean for cardamon pods?

Denied boarding although I have proper visa and documentation. To whom should I make a complaint?

On SQL Server, is it possible to restrict certain users from using certain functions, operators or statements?

Why are there no cargo aircraft with "flying wing" design?

Do jazz musicians improvise on the parent scale in addition to the chord-scales?

How do I make this wiring inside cabinet safer? (Pic)

What is the longest distance a player character can jump in one leap?

When a candle burns, why does the top of wick glow if bottom of flame is hottest?

また usage in a dictionary

When coming out of haste, do attackers have advantage on you?

How can I use the Python library networkx from Mathematica?

Is there any way for the UK Prime Minister to make a motion directly dependent on Government confidence?

Why wasn't DOSKEY integrated with COMMAND.COM?

Is safe to use va_start macro with this as parameter?

Why aren't air breathing engines used as small first stages

Wu formula for manifolds with boundary

What are the out-of-universe reasons for the references to Toby Maguire-era Spider-Man in ITSV

How come Sam didn't become Lord of Horn Hill?

8 Prisoners wearing hats

Find the length x such that the two distances in the triangle are the same

Do I really need recursive chmod to restrict access to a folder?



Extract random effects from CARBayes output



Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)
Announcing the arrival of Valued Associate #679: Cesar Manara
Unicorn Meta Zoo #1: Why another podcast?Extract values from raster into points: data loss using extract function from raster packageSpatial weight matrix with boundary effects in R?What to do with the intercept output from regression analsyis?Is there really random distribution of points?Using Random Forest for big data?Predicting Enhanced Vegetation Index using Random Forest in RHow to extract specific values with point coordinates from Kriging interpolations made in R?Plotting/Exporting the Spatial Random Effects in BYM model - CARBayes - in REstimate spatial kernel density from raster and output kernel values in a gridExtract random coordinate points using classes (ranges) of values from a raster image



.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;








0















We conducted a study at the U.S. census tract level looking at the relationship between environmental factors and social mobility in the ten largest U.S. cities. We used CARBayes::S.CARmultilevel in R to control for spatial autocorrelation between neighboring tracts and higher-order nested effects of cities (the city was included as a random effect variable).



We want to know how each city influenced the results in the fully adjusted model: In other words, what is the median coefficient for each city, relative to the reference group. We were able to calculate the median random effect coefficients as well as standard deviations and confidence intervals using the following code:



CAR_nwide_m4_re <- as.data.frame(CAR_nwide_m4[["samples"]][["zeta"]])
CAR_nwide_m4_re_median <- as.data.frame(apply(CAR_nwide_m4_re,2,FUN=median))
CAR_nwide_m4_re_ci <- as.data.frame(apply(as.matrix(CAR_nwide_m4_re), 2, function(x)mean(x)+c(-1.96,1.96)*sd(x)/sqrt(length(x))))
CAR_nwide_m4_re_sd <- as.data.frame(apply(CAR_nwide_m4_re,2,FUN=sd))


Our questions relate to which city is "V1" "V2" etc in the output?



Also, there are values for all 10 variables, so we are unclear what city was left out of the model and is therefore the reference group.










share|improve this question
























  • Welcome to GIS SE. As a new user, please take the Tour, which emphasizes the importance of asking One question per Question. Additional questions can be asked in other Questions, but you should probably go through the process once before adding new questions. Please remember to always specify the exact version of software in the body of the question, and to format code with the button.

    – Vince
    Apr 10 at 11:05

















0















We conducted a study at the U.S. census tract level looking at the relationship between environmental factors and social mobility in the ten largest U.S. cities. We used CARBayes::S.CARmultilevel in R to control for spatial autocorrelation between neighboring tracts and higher-order nested effects of cities (the city was included as a random effect variable).



We want to know how each city influenced the results in the fully adjusted model: In other words, what is the median coefficient for each city, relative to the reference group. We were able to calculate the median random effect coefficients as well as standard deviations and confidence intervals using the following code:



CAR_nwide_m4_re <- as.data.frame(CAR_nwide_m4[["samples"]][["zeta"]])
CAR_nwide_m4_re_median <- as.data.frame(apply(CAR_nwide_m4_re,2,FUN=median))
CAR_nwide_m4_re_ci <- as.data.frame(apply(as.matrix(CAR_nwide_m4_re), 2, function(x)mean(x)+c(-1.96,1.96)*sd(x)/sqrt(length(x))))
CAR_nwide_m4_re_sd <- as.data.frame(apply(CAR_nwide_m4_re,2,FUN=sd))


Our questions relate to which city is "V1" "V2" etc in the output?



Also, there are values for all 10 variables, so we are unclear what city was left out of the model and is therefore the reference group.










share|improve this question
























  • Welcome to GIS SE. As a new user, please take the Tour, which emphasizes the importance of asking One question per Question. Additional questions can be asked in other Questions, but you should probably go through the process once before adding new questions. Please remember to always specify the exact version of software in the body of the question, and to format code with the button.

    – Vince
    Apr 10 at 11:05













0












0








0








We conducted a study at the U.S. census tract level looking at the relationship between environmental factors and social mobility in the ten largest U.S. cities. We used CARBayes::S.CARmultilevel in R to control for spatial autocorrelation between neighboring tracts and higher-order nested effects of cities (the city was included as a random effect variable).



We want to know how each city influenced the results in the fully adjusted model: In other words, what is the median coefficient for each city, relative to the reference group. We were able to calculate the median random effect coefficients as well as standard deviations and confidence intervals using the following code:



CAR_nwide_m4_re <- as.data.frame(CAR_nwide_m4[["samples"]][["zeta"]])
CAR_nwide_m4_re_median <- as.data.frame(apply(CAR_nwide_m4_re,2,FUN=median))
CAR_nwide_m4_re_ci <- as.data.frame(apply(as.matrix(CAR_nwide_m4_re), 2, function(x)mean(x)+c(-1.96,1.96)*sd(x)/sqrt(length(x))))
CAR_nwide_m4_re_sd <- as.data.frame(apply(CAR_nwide_m4_re,2,FUN=sd))


Our questions relate to which city is "V1" "V2" etc in the output?



Also, there are values for all 10 variables, so we are unclear what city was left out of the model and is therefore the reference group.










share|improve this question
















We conducted a study at the U.S. census tract level looking at the relationship between environmental factors and social mobility in the ten largest U.S. cities. We used CARBayes::S.CARmultilevel in R to control for spatial autocorrelation between neighboring tracts and higher-order nested effects of cities (the city was included as a random effect variable).



We want to know how each city influenced the results in the fully adjusted model: In other words, what is the median coefficient for each city, relative to the reference group. We were able to calculate the median random effect coefficients as well as standard deviations and confidence intervals using the following code:



CAR_nwide_m4_re <- as.data.frame(CAR_nwide_m4[["samples"]][["zeta"]])
CAR_nwide_m4_re_median <- as.data.frame(apply(CAR_nwide_m4_re,2,FUN=median))
CAR_nwide_m4_re_ci <- as.data.frame(apply(as.matrix(CAR_nwide_m4_re), 2, function(x)mean(x)+c(-1.96,1.96)*sd(x)/sqrt(length(x))))
CAR_nwide_m4_re_sd <- as.data.frame(apply(CAR_nwide_m4_re,2,FUN=sd))


Our questions relate to which city is "V1" "V2" etc in the output?



Also, there are values for all 10 variables, so we are unclear what city was left out of the model and is therefore the reference group.







r spatial-statistics regression






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited Apr 10 at 11:04









TomazicM

1,8171516




1,8171516










asked Apr 10 at 10:11









mutelbrowningmutelbrowning

1




1












  • Welcome to GIS SE. As a new user, please take the Tour, which emphasizes the importance of asking One question per Question. Additional questions can be asked in other Questions, but you should probably go through the process once before adding new questions. Please remember to always specify the exact version of software in the body of the question, and to format code with the button.

    – Vince
    Apr 10 at 11:05

















  • Welcome to GIS SE. As a new user, please take the Tour, which emphasizes the importance of asking One question per Question. Additional questions can be asked in other Questions, but you should probably go through the process once before adding new questions. Please remember to always specify the exact version of software in the body of the question, and to format code with the button.

    – Vince
    Apr 10 at 11:05
















Welcome to GIS SE. As a new user, please take the Tour, which emphasizes the importance of asking One question per Question. Additional questions can be asked in other Questions, but you should probably go through the process once before adding new questions. Please remember to always specify the exact version of software in the body of the question, and to format code with the button.

– Vince
Apr 10 at 11:05





Welcome to GIS SE. As a new user, please take the Tour, which emphasizes the importance of asking One question per Question. Additional questions can be asked in other Questions, but you should probably go through the process once before adding new questions. Please remember to always specify the exact version of software in the body of the question, and to format code with the button.

– Vince
Apr 10 at 11:05










0






active

oldest

votes












Your Answer








StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "79"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);













draft saved

draft discarded


















StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f318347%2fextract-random-effects-from-carbayes-output%23new-answer', 'question_page');

);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes















draft saved

draft discarded
















































Thanks for contributing an answer to Geographic Information Systems Stack Exchange!


  • Please be sure to answer the question. Provide details and share your research!

But avoid


  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.

To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f318347%2fextract-random-effects-from-carbayes-output%23new-answer', 'question_page');

);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







Popular posts from this blog

Crop image to path created in TikZ? Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)Crop an inserted image?TikZ pictures does not appear in posterImage behind and beyond crop marks?Tikz picture as large as possible on A4 PageTransparency vs image compression dilemmaHow to crop background from image automatically?Image does not cropTikzexternal capturing crop marks when externalizing pgfplots?How to include image path that contains a dollar signCrop image with left size given

រឿង រ៉ូមេអូ និង ហ្ស៊ុយលីយេ សង្ខេបរឿង តួអង្គ បញ្ជីណែនាំ

Ромео және Джульетта Мазмұны Қысқаша сипаттамасы Кейіпкерлері Кино Дереккөздер Бағыттау мәзірі