R 用の indent 変換例 1

# 同一染色体上にある多型の数Ns
Ns<-4
# Ns個のマーカーのそれぞれのアレルの数をNaS
NaS<-c(3,4,5,2)
# アレル頻度をPaS
PaS<-list()
for(i in 1:Ns){
  PaS[[i]]<-c(rdirichlet(1,rep(1,NaS[i])))
}
NaS
PaS
 
hap1<-hap2<-rep(0,Ns)
for(i in 1:Ns){
  hap1[i]<-sample(1:NaS[i],1,prob=PaS[[i]])
  hap2[i]<-sample(1:NaS[i],1,prob=PaS[[i]])
}
hap1
hap2

これを以下のように変換

# 同一染色体上にある多型の数Ns
Ns <- 4
# Ns個のマーカーのそれぞれのアレルの数をNaS
NaS <- c(3, 4, 5, 2)
# アレル頻度をPaS
PaS <- list()
for (i in 1:Ns) {
    PaS[[i]] <- c(rdirichlet(1, rep(1, NaS[i])))
}
NaS
PaS
 
hap1 <- hap2 <- rep(0, Ns)
for (i in 1:Ns) {
    hap1[i] <- sample(1:NaS[i], 1, prob = PaS[[i]])
    hap2[i] <- sample(1:NaS[i], 1, prob = PaS[[i]])
}
hap1
hap2