复制到本帖中一些简单的函数
感谢帖子:https://bbs.pinggu.org/thread-3136233-1-1.html
#本脚本中所有程序来自:《复杂数据统计方法》一书的附录
#####################################################
#实践1(最初几步)
x=1:100
(x=1:100)
sample(x,20)
set.seed(0);sample(1:10,3)
z=sample(1:200000,100000)
z[1:10
y=c(1,3,7,3,4,2)
z[y
(z=sample(x,100,rep=T))
(z1=unique(z))
length(z1)
xz=setdiff(x,z)
sort(union(xz,z))
setequal(union(xz,z),x)
intersect(1:10,7:50)
sample(1:100,20,prob=1:100)
#实践2(一些简单运算)
pi*10^2
"*"(pi,"^"(10,2))
pi*(1:10)^-2.3
x=pi*10^2
x
print(x)
(x=pi*10^2)
pi^(1:5)
print(x,digits=12)
#实践3(关于R对象的类型等)
x=pi*10^2
class(x)
typeof(x)
class(cars)
typeof(cars)
names(cars)
summary(cars)
head(cars)
tail(cars)
str(cars)
row.names(cars)
attributes(cars)
class(dist~speed)
plot(dist~speed,cars)
plot(cars$speed,cars$dist)
#实践4(包括简单自变量为定量变量及定性变量的回归)
ncol(cars)
nrow(cars)
dim(cars)
lm(dist~speed,data=cars)
cars$qspeed=cut(cars$speed,breaks=quantile(cars$speed),include.lowest=TRUE)
names(cars)
cars[3
table(cars[3])
is.factor(cars$qspeed)
plot(dist~qspeed,data=cars)
(a=lm(dist~qspeed,data=cars))
summary(a)
#实践5(简单样本描述统计量)
x
summary(x)
min(x)
max(x)
range(x)
median(x)