我试图弄清楚为什么我在使用ddply时收到错误消息.
示例数据:
data<-data.frame(area=rep(c("VA","OC","ES"),each=4),
sex=rep(c("Male","Female"),each=2,times=3),
year=rep(c(2009,2010),times=6),
bin=c(110,120,125,125,110,130,125,80,90,90,80,140),
shell_length=c(.4,4,1,2,.2,5,.4,4,.8,4,.3,4))
bin7<-ddply(data, .(area,year,sex,bin), summarize,n_bin=length(shell_length))
错误消息:.fun(piece,...)出错:缺少参数"by",没有默认值
我昨天收到此错误消息.我重新启动R并重新编写代码,一切都很好.今天早上我再次收到错误消息并重新启动R并没有解决问题.
我还尝试运行一些示例代码并得到相同的错误消息.
# Summarize a dataset by two variables
require(plyr)
dfx <- data.frame(
group = c(rep('A', 8), rep('B', 15), rep('C', 6)),
sex = sample(c("M", "F"), size = 29, replace = TRUE),
age = runif(n = 29, min = 18, max = 54)
)
# Note the use of the '.' function to allow
# group and sex to be used without quoting
ddply(dfx, .(group, sex), summarize,
mean …我使用聚合函数来获得因子级别的范围.我正在尝试重命名列,但聚合函数的输出没有min和max作为单独的列.
# example data
size_cor <- data.frame(SpCode = rep(c(200, 400, 401), 3),
Length = c(45, 23, 56, 89, 52, 85, 56, 45, 78))
# aggregate function
spcode_range <- with(size_cor, aggregate(Length, list(SpCode), FUN = range))
输出:
spcode_range
Group.1 x.1 x.2
1 200 45 89
2 400 23 52
3 401 56 85
数据结构:
str(spcode_range)
'data.frame': 3 obs. of 2 variables:
$ Group.1: num 200 400 401
$ x : num [1:3, 1:2] 45 23 56 89 52 85
dim(spcode_range)
[1] 3 …我使用格子和格子额外绘制五个不同区域的观察值和预测值.我使用xyplot绘制观察值,然后使用格子额外的as.layer函数来添加预测线.我想在图表中添加一个图例,但没有任何运气.
以下是两个示例数据集以及用于绘图的代码.
示例数据集一.我只包括两个方面.
example1 <-
structure(list(model_predict = c(10, 25, 95, 23, 56, 70, 56,
45, 25, 50), Shell_Height = c(27, 33, 115, 25, 46, 50, 35, 35,
23, 45), SAMS_region_2015 = structure(c(1L, 1L, 1L, 1L, 1L, 2L,
2L, 2L, 2L, 2L), .Label = c("DMV", "LI"), class = "factor")), .Names = c("model_predict",
"Shell_Height", "SAMS_region_2015"), row.names = c(NA, -10L), class = "data.frame")
示例数据集二.也只包括两个方面.
example2 <-
structure(list(Meat_Weight = c(15, 27, 100, 15, 60, 75, 50, 37,
28, 60), Shell_Height = c(25, 30, 110, …我正在尝试将 A 和 B 标签添加到我使用 grid.arrange 函数组合的两个点阵图中。我无法让 B 标签出现在我想要的位置。我希望字母 A 位于顶部图的左侧和上方，而字母 B 位于底部图的左侧和上方。我尝试将行代码更改为不同的数字，但没有成功。
以下是一些示例数据：
data1<-data.frame(x=c(1:60),y=rep(c("A","B"),each=30),w=rep(c("V1","V2","V3"),times=20))
Plotting Code
require(lattice);require(grid.Extra)
p1<-bwplot(x~y, data = data1, groups = w,
pch = "|", box.width = 1/6, par.settings =mytheme ,scales=list(x=list(cex=1.2),y=list(cex=1.2)),
ylab=list("Number of Scallops per Station",fontsize=15),xlab=list("Strata",fontsize=15),
panel = panel.superpose,
panel.groups = function(x, y, ..., cex,group.number,col) {
panel.bwplot(x + (group.number-1.5)/6,y,...)
})
p2<-barchart(x~y|w,type="count",layout=c(3,1),data=data1,
col="light gray", xlab="Binned Number of Scallops per Station",ylab="Count",horizontal=F,
strip = strip.custom(bg="white",strip.levels = T),as.table=T,origin=0,family="A",
scales=list(relation="same",alternating=1,x=list(cex=1.2),y=list(cex=1.2)),
par.settings = list( grid.pars =
list(fontfamily = 'serif',cex=1.2)),
panel=function(x,y,...){
panel.barchart(x,y,...)})
grid.arrange(p1,p2)
mtext("A",3,line=2,cex=1.2,at=0,family="serif") …我有一个带小数点的数字向量.我想在数字前面添加一个0,在小数位前只有一个整数,同时在有一个的值中保留尾随0.例如,5.466应阅读05.466并43.770应阅读43.770.
这是数据
longM<-c(18.846, 18.906, 5.466, 11.19, 6.894, 7.578, 5.13, 53.868, 57.216,
0.228, 41.892, 57.576, 44.19, 29.292, 16.452, 48.306, 16.224,
28.566, 32.328, 25.908, 9.318, 7.332, 38.838, 4.812, 27.588,
20.922, 30.804, 43.872, 7.734, 41.856, 3.222, 47.364, 5.754,
21.792, 43.494, 50.82, 48.618, 4.524, 5.454, 48.744, 36.504,
44.328, 38.61, 39.084, 38.52, 39.504, 4.206, 47.928, 15.996,
42.99, 46.116, 26.496, 44.148, 40.284, 43.296, 34.248, 22.242,
29.838, 31.38, 27.12, 22.278, 21.36, 20.142, 15.93, 54.726, 12.072,
13.05, 11.886, 1.242, 54.42, …