【订正版】研一R语言选修课程第二次课后习题分享（含源代码）

怼烎@ 2023-09-27 12:31 30阅读 0赞

#### R语言作业分享 ####

*  一、写在前面的话
 *  二、题目分享
 *  三、源代码分享
 *  四、结果展示

## 一、写在前面的话 ##

昨天下午发的版本不够完善，老师希望我们也能用`outer`函数解决问题，但是我对R语言了解不深，在一番学习之后，简单实现了老师的要求，这篇博客将补全另外六种情况。

有更好的方法的同学可以私聊我，给我一个求学的机会。

## 二、题目分享 ##

![在这里插入图片描述][424401c1e73b417abeca5ba9558f94d6.png]

## 三、源代码分享 ##

sink("./23.R-Project/output2.doc", append = TRUE, split = TRUE)
    # R语言第二次作业，环境：R，编译器：VSCode('PPT第一题,未要求写，可以不做！')
    # 恰有一次出现正面
    p_1 <- dbinom(1, size = 3, prob = 0.5)
    print(p_1)
    
    # 至少有一次出现正面
    p_2 <- 1 - dbinom(0, size = 3, prob = 0.5)
    print((p_2))
    
    
    print('-------------------放回的情况-------------------')
    x <- c('白1', '白2', '白3', '白4', '红1', '红2')
    
    # 取到两只都是白球的概率
    p_3 <- dbinom(2, size = 2, prob = 2 / 3) 
    print(p_3)
    
    # 取到两只球颜色相同的概率
    p_4 <- dbinom(2, size = 2, prob = 1 / 3) + dbinom(2, size = 2, prob = 2 / 3)
    print(p_4)
    
    # 取到两只球中至少有一只是白球的概率
    p_5 <- 1 - dbinom(0, size = 2, prob = 2 / 3)
    print(p_5)
    
    print('-------------------不放回的情况-------------------')
    # 考虑不放回的情况
    # 取到两只都是白球的概率
    sam <- choose(4, 1) * choose(3, 1)
    all <- choose(6, 1) * choose(5, 1)
    p_6 <- sam / all
    print(p_6)
    
    # 取到两只球颜色相同的概率
    p_7 <- p_6 + (choose(2 , 1) / all)
    print(p_7)
    
    # 取到两只球中至少有一只是白球的概率
    p_8 <- 1 - (choose(2, 1) / all)
    print(p_8)
    
    print('-------------------利用outer函数的放回情况-------------------')
    output <- outer(x, x, FUN = 'paste')
    print(output)
    counts_3 <- 0
    counts_4 <- 0
    for (x_1 in x){
        
        for (x_2 in x){
        
            if (substr(x_1, 1, 1) == '白' & substr(x_2, 1, 1) == '白'){
        
                counts_3 <- counts_3 + 1
            }
            if (substr(x_1, 1, 1) == '红' & substr(x_2, 1, 1) == '红'){
        
                counts_4 <- counts_4 + 1
            }
        }
    }
    all_2 <- length(x) **2
    a_4 <- counts_3 / all_2
    print(a_4)
    
    a_5 <- (counts_3 + counts_4) / all_2
    print(a_5)
    
    a_6 <- 1 - counts_4 / all_2
    print(a_6)
    
    
    print('-------------------利用outer函数的不放回情况-------------------')
    counts_1 <- 0
    counts_2 <- 0
    x <- c('白1', '白2', '白3', '白4', '红1', '红2')
    for (x_1 in x){
        
        y <- x[-which(x == x_1)]
        out <- outer(x_1, y, FUN = "paste")
        print(out)
        for (y_1 in y){
        
            if (substr(x_1, 1, 1) == '白' & substr(y_1, 1, 1) == '白'){
        
                counts_1 <- counts_1 + 1
            }
            if (substr(x_1, 1, 1) == '红' & substr(y_1, 1, 1) == '红'){
        
                counts_2 <- counts_2 + 1
            }
        }
    }
    all <- length(x) * length(y)
    a_1 <- counts_1 / all
    print(a_1)
    
    a_2 <- a_1 + (counts_2 / all)
    print(a_2)
    
    a_3 <- 1 - (counts_2 / all)
    print(a_3)
    
    sink()

## 四、结果展示 ##

> source("d:\\Code Project\\23.R-Project\\homework_2.r", encoding = "UTF-8")
    [1] 0.375
    [1] 0.875
    [1] "-------------------放回的情况-------------------"
    [1] 0.4444444
    [1] 0.5555556
    [1] 0.8888889
    [1] "-------------------不放回的情况-------------------"
    [1] 0.4
    [1] 0.4666667
    [1] 0.9333333
    [1] "-------------------利用outer函数的放回情况-------------------"
         [,1]      [,2]      [,3]      [,4]      [,5]      [,6]
    [1,] "白1 白1" "白1 白2" "白1 白3" "白1 白4" "白1 红1" "白1 红2"
    [2,] "白2 白1" "白2 白2" "白2 白3" "白2 白4" "白2 红1" "白2 红2"
    [3,] "白3 白1" "白3 白2" "白3 白3" "白3 白4" "白3 红1" "白3 红2"
    [4,] "白4 白1" "白4 白2" "白4 白3" "白4 白4" "白4 红1" "白4 红2"
    [5,] "红1 白1" "红1 白2" "红1 白3" "红1 白4" "红1 红1" "红1 红2"
    [6,] "红2 白1" "红2 白2" "红2 白3" "红2 白4" "红2 红1" "红2 红2"
    [1] 0.4444444
    [1] 0.5555556
    [1] 0.8888889
    [1] "-------------------利用outer函数的不放回情况-------------------"
         [,1]      [,2]      [,3]      [,4]      [,5]
    [1,] "白1 白2" "白1 白3" "白1 白4" "白1 红1" "白1 红2"
         [,1]      [,2]      [,3]      [,4]      [,5]
    [1,] "白2 白1" "白2 白3" "白2 白4" "白2 红1" "白2 红2"
         [,1]      [,2]      [,3]      [,4]      [,5]
    [1,] "白3 白1" "白3 白2" "白3 白4" "白3 红1" "白3 红2"
         [,1]      [,2]      [,3]      [,4]      [,5]
    [1,] "白4 白1" "白4 白2" "白4 白3" "白4 红1" "白4 红2"
         [,1]      [,2]      [,3]      [,4]      [,5]
    [1,] "红1 白1" "红1 白2" "红1 白3" "红1 白4" "红1 红2"
         [,1]      [,2]      [,3]      [,4]      [,5]
    [1,] "红2 白1" "红2 白2" "红2 白3" "红2 白4" "红2 红1"
    [1] 0.4
    [1] 0.4666667
    [1] 0.9333333
    >

[424401c1e73b417abeca5ba9558f94d6.png]: https://img-blog.csdnimg.cn/424401c1e73b417abeca5ba9558f94d6.png