Slope by row r

DOY<-c(102,102,102,102,102,102,102,102,102,102,212,212,212,212,212,212, 212,212,212,212)
LOCATION <- c(1,1,1,1,1,2,2,2,2,2,1,1,1,1,1,3,3,3,3,3)
response <-c(NA,NA,NA,NA,NA,7,10,15,20,30,2,4,6,NA,8,10,15,20,30,NA) 
ts <- c(0,10,20,30, 40,0,10,20,30,40,0,10,20,30,40,0,10,20,30,40)
test.data <- data.frame(cbind(DOY, LOCATION, response, ts))


library(tidyverse)
library(broom)
pos1 <-  possibly(lm, otherwise = NULL)
prsq <- possibly(pull, otherwise = NA)
test.data %>%
     group_by(DOY, LOCATION) %>%
     nest %>%
     mutate(model = map(data, ~ pos1(response~ ts, data = .x)),
            slope = map_dbl(model, ~ 
                            .x %>% 
                                tidy %>%
                                filter(term == 'ts') %>%
                                prsq(estimate)),
            R2 = map_dbl(model, ~ 
                             .x %>%
                                   glance %>%
                                   prsq(r.squared))) %>%
      select(-data, -model)
# A tibble: 4 x 4
#    DOY LOCATION  slope     R2
#  <dbl>    <dbl>  <dbl>  <dbl>
#1   102        1 NA     NA    
#2   102        2  0.56   0.953
#3   212        1  0.149  0.966
#4   212        3  0.650  0.966

Source: stackoverflow.com

Similar codes

2 votes141 views1 answer
R convert accented characters
R
0 vote139 views1 answer
For R
R
0 vote145 views2 answers
Remove column from matrix r
R
1 vote147 views1 answer
How to re-arrange weekdays in r
R
0 vote146 views1 answer
Save data frames in a loop r
R
0 vote140 views1 answer
Next element in a loop if error in r
R
0 vote144 views1 answer
Na by column r
R
0 vote142 views1 answer
Rnorm in r
R
1 vote153 views1 answer
Find length of a list or vector in r
R
0 vote141 views1 answer
R - create a new empty variable in a dataset
R