Stationery distribution

 p=matrix(c(0.5, 0.5,0.2, 0.8),nrow = 2,byrow = TRUE)

Find_stationary_distribution=function(p)

{

n=nrow(p)

a=t(p)-diag(1,n)

a=rbind(a,rep(1,n))

b=c(rep(0,n),1)

stationary_distribution=qr.solve(a, b)

return(stationary_distribution)

}

stationary=Find_stationary_distribution(p)

cat("stationary Distribution:","\n",stationary)

Comments

Post a Comment

Popular posts from this blog

Reading and writing files

 numbers <- c(1, 2, 3, 4, 5) write(numbers, file = "numbers.txt") read_numbers <- read.table("numbers.txt") print(read_numbers) data<-data.frame(name=c("john","mary","david"),age=c(25,31, 42)) write.table(data,file="data.csv",sep=",",row.names=FALSE,col.names=TRUE) read_data <- read.csv("data.csv") print(read_data)
Read more

Mean mode median:

  dataset <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10)   mean_value <- mean(dataset)   print(paste("Mean:", mean_value))   mode_value <- names(which.max(table(dataset)))   print(paste("Mode:", mode_value))   median_value <- median(dataset)   print(paste("Median:", median_value))
Read more