OLDMONK🍺

 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)

 values<-c(10,20,30,40,50,60,70,80) groups<-factor(c("A","B","A","B","A","B","A","B")) mean_tapply<-tapply(values,groups,mean) print("Mean values using tapply:") print(mean_tapply) split_values <- split(values, groups) mean_sapply<-sapply(split_values,mean) print("Mean values using sapply:") print(mean_sapply)

 group1<-as.numeric(strsplit(readline(prompt="enter values for group1 separated by space:"),"")[[1]]) group2<-as.numeric(strsplit(readline(prompt="enter values for group2 separated by space:"),"")[[1]]) result<-t.test(group1,group2) print(result)

 factorialrecursive <-function(n) { if(n==0||n==1) { return(1) }else{ return(n*factorialrecursive(n-1)) } } num<-as.integer(readline(prompt = "entre a number")) result<-factorialrecursive(num) print(paste("the factorial of",num,"is",result))

 temperatures=c(15,20,19,-5,18,20) for(i in 1:length(temperatures)){ if(temperatures[i]<0){ cat("Error:Negative temperature detected:",temperatures[i],"at index",i,"\n") break }else{ cat("temperature at index",i,"is",temperatures[i],"degrees\n") } }

  numbers=c(3,10,15,7,20,6) result=vector("character",length(numbers)) for(i in 1:length(numbers)){ if(numbers[i]>10){ result[i]="greater than 10" } else{ result[i]="10 or less" } } print(result)

  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))