Extracting confidence intervals from lme model












3















I have to make some transformations on the confidence intervals of multiple large models made with lme(). I'm using the intervals() function to get the intervals, however it is not possible to turn it into a dataframe.Is there any way I could turn this into accessable numbers?



library(nlme)

data(mtcars)

model = lme(mpg ~ disp + cyl, random=~1|gear,  data=mtcars)

CI<-intervals(model, which = "fixed")

Approximate 95% confidence intervals



 Fixed effects:

              lower        est.         upper

(Intercept) 29.43497446 34.66099475 39.8870150353

disp        -0.04163023 -0.02058363  0.0004629675

cyl         -3.04786061 -1.58727681 -0.1266930105

attr(,"label")

[1] "Fixed effects:"



data.frame(CI)

Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = 

stringsAsFactors) : 

  cannot coerce class ‘"intervals.lme"’ to a data.frame





r confidence-interval nlme







share|improve this question























  • have a look at str(CI), you will see it is a list. Try data.frame(CI$fixed)

    – user20650
    Nov 21 '18 at 9:53











  • Try to look here if you like the tidy approach github.com/bbolker/broom.mixed

    – DJV
    Nov 21 '18 at 10:27
















3















I have to make some transformations on the confidence intervals of multiple large models made with lme(). I'm using the intervals() function to get the intervals, however it is not possible to turn it into a dataframe.Is there any way I could turn this into accessable numbers?



library(nlme)

data(mtcars)

model = lme(mpg ~ disp + cyl, random=~1|gear,  data=mtcars)

CI<-intervals(model, which = "fixed")

Approximate 95% confidence intervals



 Fixed effects:

              lower        est.         upper

(Intercept) 29.43497446 34.66099475 39.8870150353

disp        -0.04163023 -0.02058363  0.0004629675

cyl         -3.04786061 -1.58727681 -0.1266930105

attr(,"label")

[1] "Fixed effects:"



data.frame(CI)

Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = 

stringsAsFactors) : 

  cannot coerce class ‘"intervals.lme"’ to a data.frame





r confidence-interval nlme







share|improve this question























  • have a look at str(CI), you will see it is a list. Try data.frame(CI$fixed)

    – user20650
    Nov 21 '18 at 9:53











  • Try to look here if you like the tidy approach github.com/bbolker/broom.mixed

    – DJV
    Nov 21 '18 at 10:27














3












3








3








I have to make some transformations on the confidence intervals of multiple large models made with lme(). I'm using the intervals() function to get the intervals, however it is not possible to turn it into a dataframe.Is there any way I could turn this into accessable numbers?



library(nlme)

data(mtcars)

model = lme(mpg ~ disp + cyl, random=~1|gear,  data=mtcars)

CI<-intervals(model, which = "fixed")

Approximate 95% confidence intervals



 Fixed effects:

              lower        est.         upper

(Intercept) 29.43497446 34.66099475 39.8870150353

disp        -0.04163023 -0.02058363  0.0004629675

cyl         -3.04786061 -1.58727681 -0.1266930105

attr(,"label")

[1] "Fixed effects:"



data.frame(CI)

Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = 

stringsAsFactors) : 

  cannot coerce class ‘"intervals.lme"’ to a data.frame





r confidence-interval nlme







share|improve this question














I have to make some transformations on the confidence intervals of multiple large models made with lme(). I'm using the intervals() function to get the intervals, however it is not possible to turn it into a dataframe.Is there any way I could turn this into accessable numbers?



library(nlme)

data(mtcars)

model = lme(mpg ~ disp + cyl, random=~1|gear,  data=mtcars)

CI<-intervals(model, which = "fixed")

Approximate 95% confidence intervals



 Fixed effects:

              lower        est.         upper

(Intercept) 29.43497446 34.66099475 39.8870150353

disp        -0.04163023 -0.02058363  0.0004629675

cyl         -3.04786061 -1.58727681 -0.1266930105

attr(,"label")

[1] "Fixed effects:"



data.frame(CI)

Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = 

stringsAsFactors) : 

  cannot coerce class ‘"intervals.lme"’ to a data.frame





r confidence-interval nlme




r confidence-interval nlme






share|improve this question













share|improve this question











share|improve this question




share|improve this question










asked Nov 21 '18 at 9:51









Fanni Daniella SzakálFanni Daniella Szakál

185




185













  • have a look at str(CI), you will see it is a list. Try data.frame(CI$fixed)

    – user20650
    Nov 21 '18 at 9:53











  • Try to look here if you like the tidy approach github.com/bbolker/broom.mixed

    – DJV
    Nov 21 '18 at 10:27



















  • have a look at str(CI), you will see it is a list. Try data.frame(CI$fixed)

    – user20650
    Nov 21 '18 at 9:53











  • Try to look here if you like the tidy approach github.com/bbolker/broom.mixed

    – DJV
    Nov 21 '18 at 10:27

















have a look at str(CI), you will see it is a list. Try data.frame(CI$fixed)

– user20650
Nov 21 '18 at 9:53





have a look at str(CI), you will see it is a list. Try data.frame(CI$fixed)

– user20650
Nov 21 '18 at 9:53













Try to look here if you like the tidy approach github.com/bbolker/broom.mixed

– DJV
Nov 21 '18 at 10:27





Try to look here if you like the tidy approach github.com/bbolker/broom.mixed

– DJV
Nov 21 '18 at 10:27












1 Answer
1






active

oldest

votes


















2














You can use CI$fixed to get the values, then storing in a data frame. E.g.



df <- data.frame(CI$fixed)





share|improve this answer























    Your Answer






    StackExchange.ifUsing("editor", function () {
    StackExchange.using("externalEditor", function () {
    StackExchange.using("snippets", function () {
    StackExchange.snippets.init();
    });
    });
    }, "code-snippets");

    StackExchange.ready(function() {
    var channelOptions = {
    tags: "".split(" "),
    id: "1"
    };
    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: true,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: 10,
    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%2fstackoverflow.com%2fquestions%2f53409315%2fextracting-confidence-intervals-from-lme-model%23new-answer', 'question_page');
    }
    );

    Post as a guest















    Required, but never shown

























    1 Answer
    1






    active

    oldest

    votes








    1 Answer
    1






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes









    2














    You can use CI$fixed to get the values, then storing in a data frame. E.g.



    df <- data.frame(CI$fixed)





    share|improve this answer




























      2














      You can use CI$fixed to get the values, then storing in a data frame. E.g.



      df <- data.frame(CI$fixed)





      share|improve this answer


























        2












        2








        2







        You can use CI$fixed to get the values, then storing in a data frame. E.g.



        df <- data.frame(CI$fixed)





        share|improve this answer













        You can use CI$fixed to get the values, then storing in a data frame. E.g.



        df <- data.frame(CI$fixed)






        share|improve this answer












        share|improve this answer



        share|improve this answer










        answered Nov 21 '18 at 9:55









        CIAndrewsCIAndrews

        28817




        28817
































            draft saved

            draft discarded




















































            Thanks for contributing an answer to Stack Overflow!


            • 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%2fstackoverflow.com%2fquestions%2f53409315%2fextracting-confidence-intervals-from-lme-model%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

            MongoDB - Not Authorized To Execute Command

            Image
            .everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{ height:90px;width:728px;box-sizing:border-box; } 0 I have successfully enabled authorization on MongoDB and I have created an account on the admin database and then I created an account for my database called test. The following connection string to connect to my test database works successfully: mongo --host 192.168.17.52 --port 27017 -u user1 -p password --authenticationDatabase test Only problem I have now is, I cannot execute commands such as: show dbs. I get the following error when I try to do so: "errmsg" : "not authorized on admin to execute command { listDatabases: 1.0, lsid: { id: UUID("a1d5bc0d-bc58-485e-b232-270758a89455") }, $db: "admin...
            Read more

            How to fix TextFormField cause rebuild widget in Flutter

            Image
            .everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{ height:90px;width:728px;box-sizing:border-box; } 0 I've got the same issue at #9280. I tried all solutions from that topic with/without GlobalKey for Scaffold and Form but the problem wasn't fixed. I tried to build a authentication form with login and signup. But the keyboard doesn't show when I tap the TextFormField at the first time. Then I try to input some letter from physical keyboard, the widget is rebuilt. Like this: login screen Nothing to be log when I run flutter run --verbose, so I logged manually. Every tapping on field, the widget called dispose -> init -> build. [√] Flutter (Channel beta, v1.0.0, on Microsoft Windows [Version 10.0.17763.195], ...
            Read more

            in spring boot 2.1 many test slices are not allowed anymore due to multiple @BootstrapWith

            Image
            1 I tried to upgrade a yummy sandwich made of two test slices (@JsonTest and @JdbcTest in my case, crunchy test code in between) adding spring boot 2.1 flavour to it. But it seems it was not much of a success. I cannot annotate my tests with many @...Test since they are now each bringing their own XxxTestContextBootstrapper. It used to work when they all used same SpringBootTestContextBootstrapper. @RunWith(SpringRunner.class) @JdbcTest @JsonTest public class Test { @Test public void test() { System.out.printn("Hello, World !"); } } The error I get from BootstrapUtils is illegalStateException : Configuration error: found multiple declarations of @BootstrapWith for test class I understand I might be doing something wrong here but is there an easy way I could load both Json and Jdbc contex...
            Read more