Guest User

Untitled

a guest
Nov 22nd, 2017
79
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.82 KB | None | 0 0
  1. #These data are the same for CaseA and CaseB:
  2. df$data1 <- log(explanatory variable)
  3. df$data2 <- explanatory variable
  4. df$data3 <- explanatory variable #(between 0 and 1)
  5.  
  6. #These data are related to the Case, IE CaseA or CaseB.
  7. df$data4a <- measured value #(between 0 and 1)
  8. df$data4b <- 1 - df$data4a #(between 0 and 1)
  9. df$data5a <- measured value #(between 0 and 1)
  10. df$data5b <- 1 - df$data5a #(between 0 and 1)
  11.  
  12.  
  13. df$perc_a <- measured value #(between 0 and 1)
  14. dfb$perc_b <- 1 - df$perc_a #(between 0 and 1)
  15.  
  16. ModelA <- lm(perc_a~df$data4a+df$data5a+df$data1+df$data2+df$data3,data=df)
  17. ModelB <- lm(perc_b~df$data4b+df$data5b+df$data1+df$data2+df$data3,data=df)
  18.  
  19. stepA <- step(ModelA)
  20. stepB <- step(ModelB)
  21.  
  22. stepA$call$formula
  23. RETURNS:
  24. perc_a~df$data4a+df$data1+df$data3
  25.  
  26. stepB$call$formula
  27. RETURNS:
  28. perc_b~df$data4b+df$data5b`
Add Comment
Please, Sign In to add comment