Diagramm der Ergebnisse in Intervall in R

Question

n=50 
p=0.32 
P=matrix(c(p, 1-p, 0, 0, 0, 0, 
p, 0, 1-p, 0, 0, 0, 
p, 0, 0, 1-p, 0, 0, 
0, p, 0, 0, 1-p, 0, 
0, 0, p, 0, 0, 1-p, 
0, 0, 0, p, 0, 1-p), 
ncol=6, nrow=6, byrow = T)  

X=2 
for(j in 1:n) 
{Y=runif(1)  
k=P[X[j],]  
k=cumsum(k) 
if(Y<=k[1])  
{X[j+1]=1} 
else if (Y<=k[2]) 
{X[j+1]=2} 
else if (Y<=k[3]) 
{X[j+1]=3} 
else if (Y<=k[4]) 
{X[j+1]=4} 
else if (Y<=k[5]) 
{X[j+1]=5} 
else {X[j+1]=6}} 

mean(X) 

x=c(1,2,3,4,5,6) 
y=c(0.1,0.15,0.22,0.29,0.38,0.45) 
approx(x,y,xout=mean(X)) 

Ich habe den obigen Code verwendet, um einen Mittelwert (y) durch lineare Interpolation mit einem festen p zu erhalten. Aber jetzt, wie man den Code ändert, um ein Diagramm von Mittelwert (y) gegen p [0: 1] zu zeichnen ??? Ich bekam immer nur einen (y), hilf mir bitte. Psi braucht nur ca. $ y, das ist, wo ich bin stecken :(

Answer 1

Ich denke, es wäre besser, Ihren Code in function() zu ändern und es mit sapply(p.vector, ...) verwenden.

die Funktion

func <- function(p) { 
    P = matrix(c(p, 1-p, 0, 0, 0, 0, 
       p, 0, 1-p, 0, 0, 0, 
       p, 0, 0, 1-p, 0, 0, 
       0, p, 0, 0, 1-p, 0, 
       0, 0, p, 0, 0, 1-p, 
       0, 0, 0, p, 0, 1-p), 
       ncol=6, nrow=6, byrow = T) 
    X = 2 
    for(j in 1:n) 
    {Y=runif(1)  
    k=P[X[j],]  
    k=cumsum(k) 
    if(Y<=k[1])  
    {X[j+1]=1} 
    else if (Y<=k[2]) 
    {X[j+1]=2} 
    else if (Y<=k[3]) 
    {X[j+1]=3} 
    else if (Y<=k[4]) 
    {X[j+1]=4} 
    else if (Y<=k[5]) 
    {X[j+1]=5} 
    else {X[j+1]=6}} 
    return(approx(x, y, xout = mean(X))$y) 
} 

die Verwendung Funktion mit p [0: 1]

p.vec <- seq(0, 1, 0.01)      # preparation of p[0:1] as a vector 
n = 50          # defining other paramaters 
x = c(1, 2, 3, 4, 5, 6) 
y = c(0.1, 0.15, 0.22, 0.29, 0.38, 0.45) 

y.vec <- sapply(p.vec, func)     # calculation of the y about p[0:1] 

plot(p.vec, y.vec, type="o") # for example