| 共 146 篇文章 |
|
ANOVA.summary(fit) # display Type I ANOVA table drop1(fit,~.,test="F") # type III SS and F Tests.# Two-way Interaction Plot attach(mtcars) gears <- factor(gears) cyl <- factor(cyl) interaction.plot(cyl, gear, mpg, type="b", col=c(1:3), leg.bty="o", leg.bg="beige", lwd=... 阅148 转2 评0 公众公开 18-06-18 18:18 |
二、函数tapply（进行分组统计）： 使用格式为：tapply(X, INDEX, FUN = NULL, ..., simplify = TRUE)四、函数lapply与函数sapply： lapply的使用格式为：lapply(X, FUN, ...)lapply的返回值是和一个和X有相同的长度的list对象，这个list对象中的每个元素是将函数FUN应用到X的每一个元素。五、函数mapply： 函数mapply是函数sapply的变形版，ma... 阅595 转4 评0 公众公开 18-03-17 22:18 |
> lines(density(x))> lines(density(y))> lines(density(y),col="red")运行的结果想要得到多张图表，但总是前一张图被后一张图覆盖，如果是想在屏幕上显示多个图形，在Windows系统中只能在plot()后输入windows()重新生成一个图形窗口，然后就可以plot()了，然后再windows()，依次类推，例如。使用低级绘图函数，第一个用pl... 阅11062 转11 评0 公众公开 18-03-13 22:33 |
阅4586 转34 评0 公众公开 18-03-13 22:32 |
confint<-function(x,sigma=-1,alpha=0.05) { n<-length(x) xb<-mean(x) if(sigma>=0) { tmp<-sigma/sqrt(n)*qnorm(1-alpha/2);df<- n-1 } data.frame(mean=xb,df=df,a=xb-tmp,b=xb+tmp) }这个函数的使用： 如果不知道方差，则confint（x,alpha） 知道方差，则confint（x,sigma,alpha)... 阅2248 转9 评0 公众公开 18-03-13 11:38 |
mydata <- rbind(df.patients, df.population)mydata$Group <- as.logical(mydata$Sample == ''''''''Patients'''''''')mydata$Distress <- ifelse(mydata$Sex == ''''''''Male'''''''', age(nrow... 阅11090 转118 评2 公众公开 17-11-29 16:37 |
R语言绘制RSI交叉线 RSI= Relative Strength Index ,相对强弱指标其背后的基本逻辑是：股票买入的力量大于卖出的力量，则股票的价格会上涨。RSI=100*(up/(up+down))式中：up-表示一个特定时间间隔内的上涨幅度平均值down-表示一个特定时间间隔内的下跌幅度平均值从公式可以看到：RSI的取值范围是[0,100]当RSI=0时，说明up/(up+down)=0... 阅203 转2 评0 公众公开 17-10-31 22:36 |
从R中调用Cytoscape绘制复杂网络。最近在Cytoscape的mail-list上一个很热门的话题就是如何在R中调用Cytoscape绘制复杂的生物网络，从而将R强大的统计功能，尤其是igraph，sna等复杂网络分析包跟Cytoscape灵活的可视化复杂网络分析功能及众多插件结合在一起。昨天，Cytoscape的核心开发人员，UCSD的Keiichiro Ono在Cytoscape的wiki上发了一篇对... 阅1554 转6 评0 公众公开 17-10-31 22:36 |
R语言如何计算样本量？最近在学习统计学，并且试着用R语言解决习题，然而出师不利，在开始不久的一道样本量计算上就卡住了：一家公司调查员工的医疗费用，希望样本均值（这里我理解为均值差）在正负50元以内，置信度95%。type为单样本或双样本检验；方法2，在网上搜到一个计算器，可以算出样本量，并附有R代码。(n=(sd*(qnorm(1-alpha/2)+qnorm... 阅8136 转9 评0 公众公开 17-09-08 09:51 |
R语言入门第十讲：使用函数sample进行抽样R语言入门第十讲：使用函数sample进行抽样。接下来介绍第二个函数max/min，以及pmax/pmin。第一行代码我们给原来的向量x添加了一个缺失值“NA”，第二行代码及结果分别显示了在计算x的均值、最大值、最小值时是否带参数na.rm=T的情况。因为R中向量运算时具有“循环原则”，所以在上面的例子中我们也可... 阅26681 转49 评0 公众公开 17-02-09 19:06 |
