Improving my OLS Rolling Regression with lapply










0















This is the code of my rolling regression:



# betas from rolling regression will be in betas_rolling_250 matrix
fixed.window<-30
models<-c("A~B+E","A~B+C","A~B+D","A~B+C","A~B+F","A~B+B","A~B+F","A~B+C")
betas_rolling_250<-matrix(,nrow=length(models),ncol=250-fixed.window)
rownames(betas_rolling_250)<-models;

#Sample Matrix

sample_matrix_250<-matrix(rexp(10, rate=.1),nrow=250, ncol=6)
colnames(sample_matrix_250)<-c("A","B","C","D","E","F")

for(i in 1:(length(sample_matrix_250[,1])-fixed.window))
for(k in 1:length(rownames(betas_rolling_250)))
coefs<-lm(rownames(betas_rolling_250)[k],data=as.data.frame(sample_matrix_250[(221-i):(251-i),]))$coefficients[3]
betas_rolling_250[k,i]<-coefs



First question:



Is it possible to implement this by using lapply function?



I am asking this beacause I will have to do this 4 times. Because I have more 4 differents size dataframes: (sample_matrix_220, sample_matrix_200and sample_matrix_180). I dont want my code to be huge.



Ive been thinking in using list. Thats why I though about lapply function.



Is it possible?










share|improve this question






















  • zoo::rollapply is probably convenient

    – Richard Telford
    Nov 11 '18 at 15:32















0















This is the code of my rolling regression:



# betas from rolling regression will be in betas_rolling_250 matrix
fixed.window<-30
models<-c("A~B+E","A~B+C","A~B+D","A~B+C","A~B+F","A~B+B","A~B+F","A~B+C")
betas_rolling_250<-matrix(,nrow=length(models),ncol=250-fixed.window)
rownames(betas_rolling_250)<-models;

#Sample Matrix

sample_matrix_250<-matrix(rexp(10, rate=.1),nrow=250, ncol=6)
colnames(sample_matrix_250)<-c("A","B","C","D","E","F")

for(i in 1:(length(sample_matrix_250[,1])-fixed.window))
for(k in 1:length(rownames(betas_rolling_250)))
coefs<-lm(rownames(betas_rolling_250)[k],data=as.data.frame(sample_matrix_250[(221-i):(251-i),]))$coefficients[3]
betas_rolling_250[k,i]<-coefs



First question:



Is it possible to implement this by using lapply function?



I am asking this beacause I will have to do this 4 times. Because I have more 4 differents size dataframes: (sample_matrix_220, sample_matrix_200and sample_matrix_180). I dont want my code to be huge.



Ive been thinking in using list. Thats why I though about lapply function.



Is it possible?










share|improve this question






















  • zoo::rollapply is probably convenient

    – Richard Telford
    Nov 11 '18 at 15:32













0












0








0








This is the code of my rolling regression:



# betas from rolling regression will be in betas_rolling_250 matrix
fixed.window<-30
models<-c("A~B+E","A~B+C","A~B+D","A~B+C","A~B+F","A~B+B","A~B+F","A~B+C")
betas_rolling_250<-matrix(,nrow=length(models),ncol=250-fixed.window)
rownames(betas_rolling_250)<-models;

#Sample Matrix

sample_matrix_250<-matrix(rexp(10, rate=.1),nrow=250, ncol=6)
colnames(sample_matrix_250)<-c("A","B","C","D","E","F")

for(i in 1:(length(sample_matrix_250[,1])-fixed.window))
for(k in 1:length(rownames(betas_rolling_250)))
coefs<-lm(rownames(betas_rolling_250)[k],data=as.data.frame(sample_matrix_250[(221-i):(251-i),]))$coefficients[3]
betas_rolling_250[k,i]<-coefs



First question:



Is it possible to implement this by using lapply function?



I am asking this beacause I will have to do this 4 times. Because I have more 4 differents size dataframes: (sample_matrix_220, sample_matrix_200and sample_matrix_180). I dont want my code to be huge.



Ive been thinking in using list. Thats why I though about lapply function.



Is it possible?










share|improve this question














This is the code of my rolling regression:



# betas from rolling regression will be in betas_rolling_250 matrix
fixed.window<-30
models<-c("A~B+E","A~B+C","A~B+D","A~B+C","A~B+F","A~B+B","A~B+F","A~B+C")
betas_rolling_250<-matrix(,nrow=length(models),ncol=250-fixed.window)
rownames(betas_rolling_250)<-models;

#Sample Matrix

sample_matrix_250<-matrix(rexp(10, rate=.1),nrow=250, ncol=6)
colnames(sample_matrix_250)<-c("A","B","C","D","E","F")

for(i in 1:(length(sample_matrix_250[,1])-fixed.window))
for(k in 1:length(rownames(betas_rolling_250)))
coefs<-lm(rownames(betas_rolling_250)[k],data=as.data.frame(sample_matrix_250[(221-i):(251-i),]))$coefficients[3]
betas_rolling_250[k,i]<-coefs



First question:



Is it possible to implement this by using lapply function?



I am asking this beacause I will have to do this 4 times. Because I have more 4 differents size dataframes: (sample_matrix_220, sample_matrix_200and sample_matrix_180). I dont want my code to be huge.



Ive been thinking in using list. Thats why I though about lapply function.



Is it possible?







r lapply






share|improve this question













share|improve this question











share|improve this question




share|improve this question










asked Nov 11 '18 at 14:42









LinkmanLinkman

441313




441313












  • zoo::rollapply is probably convenient

    – Richard Telford
    Nov 11 '18 at 15:32

















  • zoo::rollapply is probably convenient

    – Richard Telford
    Nov 11 '18 at 15:32
















zoo::rollapply is probably convenient

– Richard Telford
Nov 11 '18 at 15:32





zoo::rollapply is probably convenient

– Richard Telford
Nov 11 '18 at 15:32












1 Answer
1






active

oldest

votes


















1














Consider generalizing your setup by passing the only changed variable, matrix column size, as input parameter. Then, past list of numeric values into an lapply call.



Below replaces any 250 to a variable (adjust as needed). The static models and fixed.window should be assigned once outside the function.



get_coeffs <- function(mat_size) 
betas_rolling <- matrix(,nrow=length(models), ncol=mat_size-fixed.window)
rownames(betas_rolling) <- models

sample_matrix <- matrix(rexp(10, rate=.1), nrow=mat_size, ncol=6)
colnames(sample_matrix) <- LETTERS[1:6]

for(i in 1:(length(sample_matrix[,1])-fixed.window))
for(k in 1:length(rownames(betas_rolling)))
coefs <- lm(rownames(betas_rolling)[k],
as.data.frame(sample_matrix[(mat_size-fixed.window)-i):(mat_size+1-i),]))$coefficients[3]
betas_rolling[k,i] <- coefs



return(betas_rolling)


matrix_list <- lapply(c(180, 200, 220, 250), get_coeffs)





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%2f53249806%2fimproving-my-ols-rolling-regression-with-lapply%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









    1














    Consider generalizing your setup by passing the only changed variable, matrix column size, as input parameter. Then, past list of numeric values into an lapply call.



    Below replaces any 250 to a variable (adjust as needed). The static models and fixed.window should be assigned once outside the function.



    get_coeffs <- function(mat_size) 
    betas_rolling <- matrix(,nrow=length(models), ncol=mat_size-fixed.window)
    rownames(betas_rolling) <- models

    sample_matrix <- matrix(rexp(10, rate=.1), nrow=mat_size, ncol=6)
    colnames(sample_matrix) <- LETTERS[1:6]

    for(i in 1:(length(sample_matrix[,1])-fixed.window))
    for(k in 1:length(rownames(betas_rolling)))
    coefs <- lm(rownames(betas_rolling)[k],
    as.data.frame(sample_matrix[(mat_size-fixed.window)-i):(mat_size+1-i),]))$coefficients[3]
    betas_rolling[k,i] <- coefs



    return(betas_rolling)


    matrix_list <- lapply(c(180, 200, 220, 250), get_coeffs)





    share|improve this answer



























      1














      Consider generalizing your setup by passing the only changed variable, matrix column size, as input parameter. Then, past list of numeric values into an lapply call.



      Below replaces any 250 to a variable (adjust as needed). The static models and fixed.window should be assigned once outside the function.



      get_coeffs <- function(mat_size) 
      betas_rolling <- matrix(,nrow=length(models), ncol=mat_size-fixed.window)
      rownames(betas_rolling) <- models

      sample_matrix <- matrix(rexp(10, rate=.1), nrow=mat_size, ncol=6)
      colnames(sample_matrix) <- LETTERS[1:6]

      for(i in 1:(length(sample_matrix[,1])-fixed.window))
      for(k in 1:length(rownames(betas_rolling)))
      coefs <- lm(rownames(betas_rolling)[k],
      as.data.frame(sample_matrix[(mat_size-fixed.window)-i):(mat_size+1-i),]))$coefficients[3]
      betas_rolling[k,i] <- coefs



      return(betas_rolling)


      matrix_list <- lapply(c(180, 200, 220, 250), get_coeffs)





      share|improve this answer

























        1












        1








        1







        Consider generalizing your setup by passing the only changed variable, matrix column size, as input parameter. Then, past list of numeric values into an lapply call.



        Below replaces any 250 to a variable (adjust as needed). The static models and fixed.window should be assigned once outside the function.



        get_coeffs <- function(mat_size) 
        betas_rolling <- matrix(,nrow=length(models), ncol=mat_size-fixed.window)
        rownames(betas_rolling) <- models

        sample_matrix <- matrix(rexp(10, rate=.1), nrow=mat_size, ncol=6)
        colnames(sample_matrix) <- LETTERS[1:6]

        for(i in 1:(length(sample_matrix[,1])-fixed.window))
        for(k in 1:length(rownames(betas_rolling)))
        coefs <- lm(rownames(betas_rolling)[k],
        as.data.frame(sample_matrix[(mat_size-fixed.window)-i):(mat_size+1-i),]))$coefficients[3]
        betas_rolling[k,i] <- coefs



        return(betas_rolling)


        matrix_list <- lapply(c(180, 200, 220, 250), get_coeffs)





        share|improve this answer













        Consider generalizing your setup by passing the only changed variable, matrix column size, as input parameter. Then, past list of numeric values into an lapply call.



        Below replaces any 250 to a variable (adjust as needed). The static models and fixed.window should be assigned once outside the function.



        get_coeffs <- function(mat_size) 
        betas_rolling <- matrix(,nrow=length(models), ncol=mat_size-fixed.window)
        rownames(betas_rolling) <- models

        sample_matrix <- matrix(rexp(10, rate=.1), nrow=mat_size, ncol=6)
        colnames(sample_matrix) <- LETTERS[1:6]

        for(i in 1:(length(sample_matrix[,1])-fixed.window))
        for(k in 1:length(rownames(betas_rolling)))
        coefs <- lm(rownames(betas_rolling)[k],
        as.data.frame(sample_matrix[(mat_size-fixed.window)-i):(mat_size+1-i),]))$coefficients[3]
        betas_rolling[k,i] <- coefs



        return(betas_rolling)


        matrix_list <- lapply(c(180, 200, 220, 250), get_coeffs)






        share|improve this answer












        share|improve this answer



        share|improve this answer










        answered Nov 11 '18 at 15:50









        ParfaitParfait

        51.1k84370




        51.1k84370



























            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%2f53249806%2fimproving-my-ols-rolling-regression-with-lapply%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

            𛂒𛀶,𛀽𛀑𛂀𛃧𛂓𛀙𛃆𛃑𛃷𛂟𛁡𛀢𛀟𛁤𛂽𛁕𛁪𛂟𛂯,𛁞𛂧𛀴𛁄𛁠𛁼𛂿𛀤 𛂘,𛁺𛂾𛃭𛃭𛃵𛀺,𛂣𛃍𛂖𛃶 𛀸𛃀𛂖𛁶𛁏𛁚 𛂢𛂞 𛁰𛂆𛀔,𛁸𛀽𛁓𛃋𛂇𛃧𛀧𛃣𛂐𛃇,𛂂𛃻𛃲𛁬𛃞𛀧𛃃𛀅 𛂭𛁠𛁡𛃇𛀷𛃓𛁥,𛁙𛁘𛁞𛃸𛁸𛃣𛁜,𛂛,𛃿,𛁯𛂘𛂌𛃛𛁱𛃌𛂈𛂇 𛁊𛃲,𛀕𛃴𛀜 𛀶𛂆𛀶𛃟𛂉𛀣,𛂐𛁞𛁾 𛁷𛂑𛁳𛂯𛀬𛃅,𛃶𛁼

            Edmonton

            Crossroads (UK TV series)