Guest User

Untitled

a guest
May 21st, 2018
81
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.38 KB | None | 0 0
  1. Error in pic(independent, tree) :
  2. missing data in 'x': you may consider removing the species with missing data from your tree with the function 'drop.tip'.
  3.  
  4. library(ape)
  5. library(geiger)
  6. library(nlme)
  7. library(phytools)
  8.  
  9. data <- read.csv("App and ADW-Museum Ecology Data.csv", row.names=1)
  10. tree <- read.nexus('c://users/user/documents/formated zelditch.nex')
  11.  
  12. td = treedata(tree, data, warnings=FALSE)
  13. tree = td$phy
  14. data = td$data
  15.  
  16. dependent <- data[,"chroma"]
  17. independent <- data[,"mass..Pantheria."]
  18.  
  19. names(independent) <- names(dependent) <- rownames(data)
  20. dependent <- pic(dependent, tree)
  21. independent <- pic(independent, tree)
  22.  
  23. model <- lm(dependent ~ independent - 1)
  24. print(summary(model))
  25.  
  26. library(ape)
  27. library(geiger)
  28. library(nlme)
  29. library(phytools)
  30.  
  31. data <- read.csv("App and ADW-Museum Ecology Data.csv", row.names=1)
  32. tree <- read.nexus('c://users/alec/documents/formated zelditch.nex')
  33.  
  34. td = treedata(tree, data, warnings=FALSE)
  35. tree = td$phy
  36. data = td$data
  37.  
  38. dependent <- data[,"chroma"]
  39. independent <- data[,"mass..Pantheria."]
  40.  
  41. names(independent) <- names(dependent) <- rownames(data)
  42.  
  43. for (name in rownames(data)) {
  44.  
  45. if(is.na(dependent[name])) {
  46. drop.tip(tree, name)
  47. }
  48.  
  49. if(is.na(independent[name])) {
  50. drop.tip(tree, name)
  51. }
  52. }
  53.  
  54.  
  55. dependent <- pic(dependent, tree)
  56. independent <- pic(independent, tree)
  57.  
  58. model <- lm(dependent ~ independent - 1)
  59. print(summary(model))
Add Comment
Please, Sign In to add comment