0

我是一名植物生态学家，我想设计一个现场操作实验。为了在实验地块内的植物布置中实现良好程度的随机化，我想使用 R（我很熟悉）。

图片是设计的示意图：总共将有 6 个地块，每个地块包含 24 种不同的植物物种。我总共有 36 种植物，我想从每个地块中随机抽取 24 种植物（地块之间的差异很大）。每个物种应出现在 6 个地块中的 4 个，但位置不同，沿边界两次（两种暖色）和两次在核心区域（每种冷色一次）；在图片中，我使用 4 个不同的物种（A、B、C、D）直观地解释了我的意思。

任何人都可以提出一种方法或一些见解来编写函数吗？

情节设计

4

1 回答 1

1

我从这个线程和这个线程中借鉴了一些见解。基本上，您制作一个圆形矩阵并混洗组。

coords <- list(c(1,1),c(1,2),c(2,1),c(3,1),
c(1,3),c(1,4),c(2,4),c(3,4),
c(2,2),c(3,2),c(2,3),c(3,3),
c(4,1),c(5,1),c(6,1),c(6,2),
c(4,2),c(5,2),c(4,3),c(5,3),
c(4,4),c(5,4),c(6,3),c(6,4))

Matrix <- matrix(c(LETTERS,0:9)[1:36][matrix(1:36,36+1,36+1,byrow=T)[c(1,36:2),1:36]],36,36)
PlotLayouts <- Matrix[(1:6*6),1:24][,unlist(lapply(split(1:24,rep(1:6,each=4)),sample,4))]
PlotLayouts <- split(PlotLayouts,sample(1:6,6))

Result <- lapply(PlotLayouts,function(Vector){
  Layout <- matrix(NA,nrow=6,ncol=4)
  for(i in 1:24){
    Layout[coords[[i]][1],coords[[i]][2]] <- Vector[i]
  }
  Layout
})

#Species Counts
table(unlist(Result))
0 1 2 3 4 5 6 7 8 9 A B C D E F G H I J K L M N O P Q R S T U V W X Y Z 
4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 

Result
$`1`
     [,1] [,2] [,3] [,4]
[1,] "I"  "H"  "M"  "L" 
[2,] "K"  "R"  "S"  "N" 
[3,] "J"  "P"  "Q"  "O" 
[4,] "U"  "Z"  "0"  "2" 
[5,] "V"  "Y"  "X"  "1" 
[6,] "T"  "W"  "4"  "3" 

$`2`
     [,1] [,2] [,3] [,4]
[1,] "0"  "Z"  "4"  "3" 
[2,] "2"  "9"  "A"  "5" 
[3,] "1"  "7"  "8"  "6" 
[4,] "C"  "H"  "I"  "K" 
[5,] "D"  "G"  "F"  "J" 
[6,] "B"  "E"  "M"  "L" 
...
于 2020-03-26T16:52:06.163 回答