Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/69.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 如何删除图形上的额外点_R - Fatal编程技术网

R 如何删除图形上的额外点

R 如何删除图形上的额外点,r,R,我编写了以下R函数: #initialprob <- c(0.4,0.5,0.1) f1 <- function(n,m,priceinitial,delta,mean, sd, ninterval){ initialprob <- c(1/3, 1/3, 1/3) traders <- vector(mode="character", length=n) traderscurrent <- vector(mode="character", length

我编写了以下R函数:

#initialprob <- c(0.4,0.5,0.1)
f1 <- function(n,m,priceinitial,delta,mean, sd, ninterval){
  initialprob <- c(1/3, 1/3, 1/3)
  traders <- vector(mode="character", length=n)
  traderscurrent <- vector(mode="character", length=n)
  price <- vector(mode="numeric")
  pricecurrent <- vector(mode="numeric")

  for(nint in 1:ninterval)
  {
    print(initialprob)
    L = floor(rnorm(1,mean,sd))
    print(L)
    for(i in 1:n)
    {
      traders[i] = sample(c("B", "S", "N"), size=1, prob=initialprob)  
    }
    print(table(traders))
    for(step in 1:L)
    {

      for(i in 1:n)
      {
        b <- sample(traders[-i], m)
        #print(b)
        #table(b)
        traderscurrent[i] <- sample(b,1)

      }
      print(table(traderscurrent))
      buy = length(which(traderscurrent == "B"))
      sell = length(which(traderscurrent == "S"))
      neutral = length(which(traderscurrent == "N"))
      total = buy+neutral+sell
      buyprop = buy/total
      sellprop = sell/total
      neutralprop = neutral/total
      pricecurrent[step] = priceinitial+buy*delta-sell*delta
      priceinitial = pricecurrent[step]
      traders <- traderscurrent
      #print(nint)
      #print
    initialprob <- c(sellprop,buyprop,neutralprop)
    print(initialprob)
    }
    a <- runif(1,0,1)
    b <- runif(1,0,1)
    c <- runif(1,0,1)
    total = a+b+c
    initialprob <- c(a/total, b/total, c/total)
    print(initialprob)
    price <- append(price,pricecurrent, after=length(price))
    #price <- price[-step]
    plot(price)

  }
  print(price)
  #plot(price)

}

#initialprob问题似乎是每次迭代后都没有重置pricecurrent向量。尝试添加
pricecurrent
for(nint in 1:ninterval)
{
  pricecurrent <- vector(mode="numeric") ## <- added
  print(initialprob)
  L = floor(rnorm(1,mean,sd))
  print(L)
  for(i in 1:n)
  ....