R语言教程课后习题答案（持续更新中~~）

R语言教程网址如下

https://www.math.pku.edu.cn/teachers/lidf/docs/Rbook/html/_Rbook/index.html

目录

source()函数可以运行保存在一个文本文件中的源程序

R向量下标和子集

数值型向量及其运算

日期功能

R因子类型

---

source()函数可以运行保存在一个文本文件中的源程序

如下内容保存在文件ssq.r中

sum.of.squares <- function(x){
  sum(x^2)
}

用source()命令运行

source("ssq.r")

可以运用sum.of.squares(1:5)

sum.of.squares(1:5)

R向量下标和子集

class.csv

name,sex,age,height,weight
Alice,F,13,56.5,84
Becka,F,13,65.3,98
Gail,F,14,64.3,90
Karen,F,12,56.3,77
Kathy,F,12,59.8,84.5
Mary,F,15,66.5,112
Sandy,F,11,51.3,50.5
Sharon,F,15,62.5,112.5
Tammy,F,14,62.8,102.5
Alfred,M,14,69,112.5
Duke,M,14,63.5,102.5
Guido,M,15,67,133
James,M,12,57.3,83
Jeffrey,M,13,62.5,84
John,M,12,59,99.5
Philip,M,16,72,150
Robert,M,12,64.8,128
Thomas,M,11,57.5,85
William,M,15,66.5,112

#取出其中的name和age列到变量name和age中
d.class <- read.csv("class.csv", header=TRUE, stringsAsFactors=FALSE)
name <- d.class[,"name"]
age <- d.class[,"age"]

1.求出age中第3, 5, 7号的值

age[c(3,5,7)]
输出 [1] 14 12 11

2.用变量age, 求出达到15岁及以上的那些值

age[age>=15]
输出 [1] 15 15 15 16 15

3.用变量name和age, 求出Mary与James的年龄

mary_age <- age[name == "Mary"]
james_age <- age[name == "James"]

> mary_age
[1] 15
> james_age
[1] 12

4.求age中除Mary与James这两人之外的那些人的年龄值，保存到变量age1中

age1<-age[name!="Mary" & name!="James"]
> age1
 [1] 13 13 14 12 12 11 15 14 14 14 15 13 12 16 12 11 15

5.假设向量x长度为n, 其元素是{1,2,...,n}的一个重排。 可以把x看成一个i到x[i]的映射(i在{1,2,...,n}中取值)。 求向量y, 保存了上述映射的逆映射，即： 如果x[i]=j, 则y[j]=i

n <- length(x)
y <- vector("integer", length = n)

for (i in 1:n) {
  j <- x[i]
  y[j] <- i
}

数值型向量及其运算

显示1到100的整数的平方根和立方根

for (i in 1:100) {
  sqrt_val <- sqrt(i)
  cubert_val <- i^(1/3)
  print(paste("整数:", i, " 平方根:", sqrt_val, " 立方根:", cubert_val))
}

生成[0,1]区间上等间隔的100个格子点存入变量x中

x <- seq(0, 1, length.out = 100)
print(x)

日期功能

#dates.csv中包含如下内容
"出生日期","发病日期"
"1941/3/8","2007/1/1"
"1972/1/24","2007/1/1"
"1932/6/1","2007/1/1"
"1947/5/17","2007/1/1"
"1943/3/10","2007/1/1"
"1940/1/8","2007/1/1"
"1947/8/5","2007/1/1"
"2005/4/14","2007/1/1"
"1961/6/23","2007/1/2"
"1949/1/10","2007/1/2"

date1 <- dates.tab[,'出生日期']
date2 <- dates.tab[,'发病日期']

1.把date1、date2转换为R的POSIXct日期型

dates_converted1 <- as.POSIXct(date1, format = "%Y/%m/%d")
> dates_converted
 [1] "1941-03-08 CST" "1972-01-24 CST" "1932-06-01 CST" "1947-05-17 CDT"
 [5] "1943-03-10 CDT" "1940-01-08 CST" "1947-08-05 CDT" "2005-04-14 CST"
 [9] "1961-06-23 CST" "1949-01-10 CST"
dates_converted2 <- as.POSIXct(date2, format = "%Y/%m/%d")
> dates_converted
 [1] "2007-01-01 CST" "2007-01-01 CST" "2007-01-01 CST" "2007-01-01 CST"
 [5] "2007-01-01 CST" "2007-01-01 CST" "2007-01-01 CST" "2007-01-01 CST"
 [9] "2007-01-02 CST" "2007-01-02 CST"

2.date1中的各个出生年

library(lubridate)

date1 <- c("1941/3/8", "1972/1/24", "1932/6/1", "1947/5/17", "1943/3/10",
           "1940/1/8", "1947/8/5", "2005/4/14", "1961/6/23", "1949/1/10")

birth_years1 <- year(as.Date(date1, format = "%Y/%m/%d"))

print(birth_years1)

birth_years2 <- year(as.Date(date2, format = "%Y/%m/%d"))

print(birth_years2)

3. 计算发病时的年龄，以周岁论（过生日才算）

ages <- year(date2) - year(date1)
for (i in seq_along(ages)) {
     if (month(date2[i]) < month(date1[i])) {
         ages[i] <- ages[i] - 1
     }
 }
print(ages)

#输出
 [1] 65 35 74 59 63 67 59  1 45 58

4.把date2中发病年月转换为'monyy'格式，这里mon是如FEB这样英文三字母缩写， yy是两数字的年份

> formatted_date <- format(as.Date(date2), "%b%y")
> 
> formatted_date
 [1] "1月07" "1月07" "1月07" "1月07" "1月07" "1月07" "1月07" "1月07" "1月07"
[10] "1月07"

R因子类型

1.把sex变量转换为一个新的因子，F显示成"Female"，M显示成"Male"

factor(sex, levels = c("F", "M"), labels = c("Female", "Male")