如何在xyplot / panel函数中使用predict？
在xyplot / panel函数中使用predict的方法如下：
- 首先，确保你已经加载了适当的R包，例如lattice或ggplot2。
- 创建一个数据集，包含你要绘制的x和y变量。
- 使用xyplot函数创建一个基本的散点图或线图。例如，使用xyplot(y ~ x, data = dataset)创建一个简单的散点图。
- 在panel函数中使用predict函数来添加预测线或曲线。panel函数是xyplot函数的一个参数，用于自定义每个面板的绘图。
- 在panel函数中，使用predict函数来生成预测值。你需要提供一个适当的模型对象和要预测的新数据。
- 使用lines函数将预测线或曲线添加到图形中。例如，使用lines(x, predicted_values)将预测值绘制为线条。
下面是一个示例代码：
# 加载所需的包
library(lattice)
# 创建一个数据集
x <- 1:10
y <- x^2
dataset <- data.frame(x, y)
# 创建散点图
xyplot(y ~ x, data = dataset,
main = "Scatter Plot with Predicted Line",
xlab = "X", ylab = "Y",
type = c("p", "g"))
# 在panel函数中使用predict
panel <- function(x, y, ...) {
panel.xyplot(x, y, ...)
# 创建线性模型
lm_model <- lm(y ~ x, data = dataset)
# 生成预测值
predicted_values <- predict(lm_model, newdata = data.frame(x = x))
# 添加预测线
panel.lines(x, predicted_values, col = "red")
}
# 添加panel函数到xyplot中
xyplot(y ~ x, data = dataset, panel = panel)这个例子中，我们创建了一个简单的散点图，并使用线性回归模型生成了预测线。在panel函数中，我们首先创建了一个线性模型，然后使用predict函数生成了预测值。最后，我们使用lines函数将预测线添加到图形中。
请注意，这只是一个示例，你可以根据你的需求和数据类型进行适当的修改和扩展。
相关·内容
1回答
我已经定义了这个功能：
fresults <- function(svrFit,testSet,var) {
library(caret)
if(substr(deparse(substitute(svrFit)), 1, 3)!="svr") {
options(digits=3)
tiff(filename = paste("Predicted vs. Residuals (Train)", "_", deparse(substitute(svrFit)), ".tiff", sep=""
浏览 2提问于2014-09-10得票数 1
回答已采纳
我通过绘制和拟合数据
> xyplot(Amplification~Voltage,
data = before_database[before_database$Serial_number==912009913,],
grid = TRUE,
main = "SN = 912009913",
panel = function(x, y) {
panel.xyplot(x, y)
fit<-lm(y ~ poly(x,2))
浏览 7提问于2017-07-05得票数 0
回答已采纳
1回答
我想在ellipse colors和width中更改lattice xyplot。下面是一个示例：
library(lattice)
xyplot(Sepal.Length ~ Petal.Length, groups=Species,
data = iris, scales = "free",
par.settings = list(superpose.symbol = list(pch = 18, cex = 0.9,col = c("green", "orange","brown"),superp
浏览 0提问于2016-10-04得票数 1
回答已采纳
3回答
我这里有点麻烦，请帮帮我。我有这些数据
set.seed(4)
mydata <- data.frame(var = rnorm(100),
temp = rnorm(100),
subj = as.factor(rep(c(1:10),5)),
trt = rep(c("A","B"), 50))
这个适合他们的模型
lm <- lm(var ~ temp * subj, data = mydata)
我想用格子来绘制结果
浏览 1提问于2012-01-13得票数 3
回答已采纳
1回答
我是一个格包中的木作图的初学者，我正在尝试可视化我的数据框架。我想给出x轴上的连续变量和y轴上的一对因子变量(换句话说，y轴上的嵌套变量)的值。
问:我如何向y轴添加自定义条形标签？
我阅读了手动文档，在没有找到解决方案的情况下搜索了这个问题。我希望得到的是这样的东西
这有可能得到这样的图解用木偶吗？
更新
事实上，ggplot2可能是一个不错的选择。然而，通过y~x1\x2的公式，xyplot似乎更方便，而我不知道如何制作这个图。
这里是一个可重复的例子
method<-rep(c("XGBM", "LGBM", "WGBM",
浏览 0提问于2019-04-02得票数 0
回答已采纳
1回答
我知道我可以用panel.xyarea从latticeExtra填补该地区的地块任何颜色。在没有在type中定义xyplot参数的情况下，这种填充将遵循默认的type="p"路径。
library(lattice)
library(latticeExtra)
data <- data.frame(time=1:24,value=rnorm(24))
xyplot(value~time, data,
panel=function(x,y,...){
panel.xyarea(x,y,...)
panel.x
浏览 2提问于2014-12-10得票数 9
回答已采纳
1回答
我已经开始转换从基本R图形到格子图形。现在我有了一些数据和一个适合于这些数据的nls()，我正在试图将两者都绘制到一个图上。因此，作为一个简单的R基图形代码：
require(lattice)
# Data
V <- c(45.5513,57.9145,66.3092,72.3338,79.7007,92.5738,107.9735,153.7357,210.8857)
E <- c(-0.05403,-0.20890,-0.24562,-0.25644,-0.26019,-0.25213,-0.23181,-0.17106,-0.10860)
# Guesses
V0g <
浏览 1提问于2014-04-03得票数 0
回答已采纳
我试图在lattice图中自动添加多个模型评估参数。我使用以下代码
library(lattice)
library(tidyverse)
library(hydroGOF)
#Calculation of model evaluation parameters
summ <- iris %>%
group_by(Species) %>%
summarise(Rsq = cor(Sepal.Length, Petal.Length)^2,
RMSE = RMSE(Sepal.Length, Petal.Length),
浏览 2提问于2021-08-23得票数 2
回答已采纳
嗨，我试着用下面的代码向格木素图中添加一条中线.这不管用..。
library(lattice)
temp <- data.frame(c(1,1,1,1,1,1,2,2),
c(1,1,2,2,1,1,1,1),
c(1,1,1,2,1,2,1,2),
c(0,0,0,150,150,0,0,0),
c(3,10,4,15,14,2,4,5))
names(temp) <- c("instNo","NoOfAC","NoReleaseLevels","trExt",&#
浏览 3提问于2014-01-24得票数 1
回答已采纳
我正在尝试用lattice绘制一个图，其中点的颜色和形状都会根据数据中的变量发生变化。当我自己设置col和pch值，并且不将panel.xyplot放在ifelse语句中时，这样做效果很好： #Setup
library(cluster)
library(lattice)
xyplot(V8 ~ V7,
data = flower,
V8_var = flower$V8,
panel = function(x, y, ..., V8_var, subscripts) {
col <- ifelse(V8_
浏览 20提问于2021-02-03得票数 0
回答已采纳
我想添加不同的文字到每个面板在格中的木素图。
res<- xyplot(CumSpec ~ CumTotal | Site, data=data1, index.cond=list(c(1,2,3)),layout = c(3,1,1), aspect = 1,
axis=axis.overlap, origin=0, xlab="Total number of individuals", ylab="Total number of species",
between = list(x = 0),
s
浏览 6提问于2014-05-21得票数 6
回答已采纳
我用xyplot创建了一个图，它必须显示两组点，它们是回归线。
xyplot(log(Vegetati)~log(Reprod),
type=c("p", "r"),
group=Espece,
data=plantes,
panel = "panel.superpose",
auto.key =list(
points = FALSE,
columns=2),
xlab="log(modules végétati
浏览 0提问于2013-09-16得票数 1
1回答
这里有一个简单的问题，你肯定已经遇到了，但却让我很头疼……
我有一个这样的数据帧：
set.seed(3)
mydata <- data.frame(var = rnorm(100,20,1),
temp = sin(sort(rep(c(1:10),10))),
subj = as.factor(rep(c(1:10),5)))
我需要为每个subj做一个散点图，这不是问题，但是...我想要的是替换网格中的条带，并为每个图添加一个标签。我用下面的代码设法做到了这一点，但我仍然不满意……
xyplot(v
浏览 0提问于2011-12-15得票数 9
回答已采纳
我有一个数据集，在两个日期x的值为y，对应于两个组gp。另外，对于这些值中的每一个，我都有一个图形n，我想在x轴上添加它。下面是生成数据和图形的代码： set.seed(1234)
data0 <- expand.grid(gp = c("A","B")
, x = 1:2
)
data0$y <- round(runif(4),2)
data0$n <- round(100*runif(4),0)
data0
# gp x y n
# 1 A 1 0.1
浏览 22提问于2019-10-07得票数 0
回答已采纳
1回答
我已经加载了晶格包。然后我运行：
> xyplot(yyy ~ xxx | zzz, panel = function(x,y) { panel.lmline(x,y)}
这将生成绘图面板，显示回归线，而不显示xyplot。我在做panel.lmline的时候并没有完全理解它是如何做到的。我知道有一个数据参数，数据是什么，知道我有3个变量xxx，yyy，zzz
浏览 0提问于2012-10-19得票数 13
回答已采纳
1回答
我在数据框上使用xyplot，并希望向面板函数提供不是(x，y，...)的数据。参数，但数据框的一些附加列(例如下面示例中的k)：
library(lattice)
mydata <- data.frame(xxx = 1:20,
yyy = rbinom(20,20,.5),
uuu = gl(2,10),
k = rnorm(20))
xyplot( formula = yyy ~ xxx | uuu, data = mydata,
pan
浏览 2提问于2013-04-16得票数 8
回答已采纳
1回答
是否有人用CombinedXYPlot类型的地块来实现JSplitPane中的每个地块，而不是在一个JPanel中实现所有的情节？这样，可以在不影响其他绘图的情况下修改每个地块的一个大小维度(例如，如果每个地块位于分裂窗格中，则可以通过拖动JSplitPane分隔器手动调整每个地块的垂直高度)。
示例代码如下：
import java.awt.Color;
import java.awt.EventQueue;
import java.awt.FlowLayout;
import java.awt.Panel;
import java.util.Random;
import javax.sw
浏览 4提问于2014-05-18得票数 1
回答已采纳
我正在尝试构建一个xyplot，它根据一个额外的分类变量的不同值包含不同的背景颜色。使用panel.xblocks (package: latticeExtra)函数获得重复的背景着色是没有问题的，但是到目前为止，我还没有找到一种方法来为xyplot中的不同子图实现不同的着色。 JD <- c(seq(0,19, 1), seq(0,19, 1))
VAR <- c(rnorm(20, mean=10, sd=1), rnorm(20, mean=10, sd=1))
CATEG <- c(rep("A", 5), rep("B", 15),
浏览 14提问于2019-06-05得票数 0
1回答
我正在尝试为我的情节创建一个定制的黄土面板功能。基本上，它应该做的与简单的panel.loess或type = "smooth"相同。原因是以后我想让事情变得更复杂一些，这些事情很难用上面的方法实现。然而，它失败了。
这里有一个MWE (其中一些是基于Sarkar的格子书第235页的例子)：
library(lattice)
set.seed(9)
foo <- rexp(100)
bar <- rexp(100)
thing <- factor(rep(c("this", "that"), times = 50))
d.f
浏览 4提问于2015-03-12得票数 1
回答已采纳
1回答
这是我在书中看到的代码：
print(xyplot(lower + upper ~ x, data=NPL.bands(dataset$variable),
panel = function(x,y,...) {
panel.xyplot(x,y,...)
panel.curve(pnorm(x,7.22,1.14), lwd=1.5)
}))
因此，我开始阅读以下内容:为该数据集绘制置信区间的上端和下端。
但在那之后我很困惑。我知道面板功能可以对面板中的任何内容进行更专门的更改。
代码似乎是这样写的：“我们在这个面板中放置了一个xy函数，
浏览 4提问于2016-11-25得票数 1
回答已采纳
我试图将广义加性模型gam()从mgcv包合并到xyplot()函数或R中的lattice包中的coplot()函数。
通过选择臭氧数据，可以在中找到这些数据。
这是我的内核平滑代码。
ozonedata=ozone
Temperature=equal.count(ozonedata$temperature,4,1/2)
Wind=equal.count(ozonedata$wind,4,1/2)
xyplot(ozone^(1/3) ~ radiation | Temperature * Wind, data = ozonedata, as.table =
浏览 2提问于2014-06-01得票数 3
回答已采纳
2回答
在此图中，分组由颜色表示。
xyplot(mpg ~ disp, data = mtcars,
type = c('p', 'g'), auto.key = TRUE,
groups = cyl)
在下面的图中，我也想通过颜色来表示分组，就像上面的例子一样。
xyplot(mpg ~ disp, data = mtcars, auto.key = TRUE,
groups = cyl,
panel = function(x, y, ...) {
panel.xyplot(x, y, type = 'g', ..
浏览 2提问于2014-08-03得票数 2
1回答
我想在网格中自定义标题...特别是对于splom，我想去掉标题“散点图矩阵...”在图形的底部。
splom(
iris[,1:4],
upper.panel = function(x, y, ...) {
panel.xyplot( x, y )
l <- lm(y~x)
panel.xyplot( x,fitted(l), type="l" )
l <- lowess(x,y)
panel.xyplot( l$x, l$y, col="red", t
浏览 5提问于2011-02-24得票数 2
回答已采纳
我正在尝试使用下面的代码创建具有多个面板的格子xyplot。问题是，它在所有面板上添加了所有中线(即，不按主题分组)。有没有办法对每个受试者的中位线进行子集？
下面是一个使用PKPDdatasets包中的Theoph数据集的可重现的示例。
# devtools::install_github("dpastoor/PKPDdatasets")
library(PKPDdatasets)
library(lattice)
library(latticeExtra)
library(dplyr)
Theoph_Data <- Theoph
Theoph_Data1 <
浏览 0提问于2016-08-07得票数 1
好的，我一直在寻找一个函数来帮助我改变xyplot中回归线的大小。我尝试过panel.abline(reg = coef, size = 2)，但它不起作用。
听到就是我的代码，
Plot.col <- brewer.pal(8,"Set1")[8][cut(Std.Change, c(-2.5,2.5), label = FALSE)]
Plot.ord <- rev(order(Std.Change))
coef <- coef(lm(Std.ModernC ~ Std.Change, data = Std.DataReg))
xyplot(Std.Mo
浏览 0提问于2012-03-03得票数 4
回答已采纳
1回答
我在理解xyplot()函数在R中创建绘图时遇到了困难，下面有一个R代码示例，它确实创建了一个很好的绘图
install.packages("mice")
library("mice")
data <- airquality[, c("Ozone", "Solar.R")]
# Applies regression imputation to Ozone w.r.t Solar.R
imp <- mice(data, method = "norm.predict", seed = 1,
浏览 2提问于2022-05-09得票数 0
回答已采纳
2回答
我不知道如何在特定的点阵图中处理存储在data.frame中的对象。在第二个图中，我得到下面的错误消息。有没有可能让它正常工作？
require(lattice)
require(latticeExtra)
data<-data.frame(a=I(list(1,2,3)),b=factor(1:3))
ecdfplot(~a|b,data=data
,layout=c(1,3)
,panel=function(x,...){
print(x[[1]])
浏览 0提问于2012-12-15得票数 0
回答已采纳
1回答
我使用wght列中列出的权重对不同的组(保存在grp列中)执行线性回归。这工作得很好。 by_time_frame <- group_by(df, grp)
weighted_model <- do(by_time_frame, tidy(lm(y~x, weights=wght, data=.))) 但是，我无法在散点图上绘制每组的(加权)回归线： xyplot(y ~ x | grp, data=df, main='Original Data'),
panel = function(x, y) {
panel.xyplot(
浏览 24提问于2019-10-18得票数 0
2回答
我有一个“长”数据帧定义为：
q <- data.frame(Indicator.Code=factor(),Year=numeric(),Value=numeric())
并且我尝试在单个xyplot中为每个不同的Indicator.Code绘制作为年份的函数的值，如下所示
xyplot( Value~Year,data=q,group=Indicator.Code)
到现在为止还好。现在，我尝试添加与线性回归相对应的行
rlm(q$Value[q$Indicator.Code==a]~q$Year[q$Indicator.Code==a])
用于Indicator.Code的所有值。
浏览 1提问于2013-05-03得票数 4
我想在晶格的一个面板上叠加多个组，并想要独立的回归线。
很容易得到多个面板，每个面板都有一条回归线，通过使用一个条件因子：
xyplot(
Petal.Width ~ Petal.Length | Species,
data = iris,
panel = function(x, y, ...) {
panel.xyplot(x, y, ...)
panel.abline(lm(y~x), col='#0080ff')
},
grid = TRUE
)
打印叠加xyplot中所有点的单个回归也是相当容易的：
xyplot(
P
浏览 4提问于2013-01-07得票数 6
回答已采纳
我正在尝试添加一个线性趋势线到散点图中，它已经被panel.smoothScatter平滑了。
基于其他帖子，我尝试添加type=c("p","r")，但这没有任何作用。我还尝试将面板设置为调用两者的函数，但结果却是一个空白的图：
#plots smoothed scatter without trend line
p2 <- xyplot(y~x, grid=T, type=c("p","r"), panel=panel.smoothScatter)
#plots nothing
p2 <- xyplot(y~x,
浏览 12提问于2014-10-01得票数 1
回答已采纳
1回答
我想用线条来画一个小贴画。举个例子，我使用的是虹膜数据，我知道实际的情节并没有真正的意义:)
library(lattice)
panel.iris <- function(x,y){
panel.xyplot(x,y)
panel.lines(x,y)
}
xyplot(Sepal.Length ~ Sepal.Width| Petal.Width, data = iris, groups= Species,
auto.key=list(space="right"), panel=panel.iris)
这给了我以下数字：
但实际上我希望
浏览 5提问于2017-09-14得票数 3
回答已采纳
1回答
样本数据/示例：
x <- log(seq(1,11,1))
y <- seq(0,1,0.1)
plot <- xyplot(y~x, type="l")
update(plot, panel = function(...) {
panel.xyplot(...)})+as.layer(xyplot(y[3]~x[3], type=c("p"), lwd=3, cex=3, pch=3))
现在，假设我对函数x输出值高于2.5的点有一个新的度量，类似于数字10，我必须将初始的地轴保持在y=0:1和x=0:2,5上。
如何以图形方式显示该值超出
浏览 1提问于2013-11-21得票数 0
回答已采纳
1回答
我想创建一个包含两个面板的图形，每个面板包含两条ROC曲线。为了介绍我的(失败)方法，我生成了包含四种方法的真正标签的数据帧(每个方法对应于)
N <- 20
TF <- rep(c(0,1),each=N/2)
pred <- method <- true <- NULL
for (imethod in 1 : 4){
pred <- c(pred,seq(-1,1,length.out=N) + rnorm(N) )
method <- c(method,rep(imethod,N))
true <- c(true,TF)
}
浏览 2提问于2014-07-09得票数 0
回答已采纳
我想创建一系列x-y散点图，其中y始终是相同的变量，x是我要检查它们是否相关的变量。作为示例，让我们使用mtcars数据集。
我对R相对较新，但正在变得更好。
下面的代码可以工作，列表图表包含所有图表，除了X轴显示为"x"，并且我希望它是变量的名称。我尝试了许多xlab=的组合，但我似乎不能理解
如果我使用names(data)，我会看到我想要使用的名称。我想我想在apply的第一次迭代中引用names(data)的第一次，在第二次迭代中引用第二次，等等，我该怎么做呢？
下一步是将它们一起打印在一个网格中，我假设lapply或sapply可以使用打印功能-我也很欣赏这个想法，只
浏览 0提问于2013-01-15得票数 3
回答已采纳
2回答
我想使用网格来创建多个面板散点图，并使用分组因子通过每个面板添加平滑线，这是一个虚拟数据帧
x <- rep(1:10, 4)
a <- as.factor(rep(1:4, 10))
b <- as.factor(sort(rep(1:2, 20)))
y <- rep(NA, 80)
df <- data.frame(x, y, a, b)
df$y[df$b=="1"] <- df$x[df$b=="1"]+df$x[df$b=="1"]^0.5
df$y[df$b=="2"] &l
浏览 0提问于2013-04-18得票数 4
回答已采纳
我有类似数据的100+文件遵循了几乎相同的趋势。我已经找到了他们的最佳人选，但现在我想把这与理论上的论点相比较。换句话说，我想用实验数据为我所生成的最佳拟合曲线生成一个方程；该方程适用于任意随机值，在一个特定的范围内，并产生与以前相似的曲线，当然也有一些误差。
代码：
set.seed(42)
x <-sort(round(runif(10,0,53))) ## random x values
y <- runif(10,0,400) ## random y values
data1 <- data.frame(y=y,x=x) ##
浏览 1提问于2015-11-18得票数 3
回答已采纳
5回答
基于这个伟大的问题：
在晶格中如何做到这一点？
plot(rnorm(120), rnorm(120), col="darkblue", pch=16, xlim=c(-3,3), ylim=c(-4,4))
points(rnorm(120,-1,1), rnorm(120,2,1), col="darkred", pch=16)
points(c(-1,-1.5,-3), c(4,2,0), pch=3, cex=3)
xspline(c(-1,-1.5,-3), c(4,2,0), shape = -1)
下面是类似的数据，它们的格式更适合lat
浏览 0提问于2013-01-13得票数 7
回答已采纳
1回答
我使用网格包来绘制图片使用下面的csv文件。 csv file 我使用了下面的代码，得到了预期的图片。 library(lattice)
data_1 <- read.table("./stackoverflow.csv", header=T, sep=",")
data <- na.omit(data_1)
xyplot(a ~ b | c, data = data, panel = function(x, y){
panel.xyplot(x, y)
}
) ? 此外，我使用了下面的代码，并得到了下一个预期的图片。 librar
浏览 24提问于2021-07-02得票数 0
回答已采纳
我正在使用Growthrate包为我的数据生成增长率曲线的参数估计。我已经进行了回归，并查看了生成的曲线图，我对数据很满意，但我想在ggplot2中重现以下曲线图。 图1：Multiplot of a regression for each group:treatment combo 我想要每个组的回归线的多线图:处理组合，但包括我在(即logistic，gompertz，gompertz2等)中对其执行的所有回归。到目前为止，我有： library(growthrates)
####Using logistic regression to fit the data across mutl
浏览 11提问于2020-01-27得票数 1
回答已采纳
1回答
当我尝试在日志尺度散点图上使用panel.text (或ltext)添加文本时，我遇到了这个问题。我知道正常的情节，你会用：
attach(mtcars) # using the mtcars as an example:
mytext<- paste("text here")
xyplot(mpg~wt, scales=list(cex=.8, col="red"),
xlab="Weight", ylab="Miles per Gallon",
main="MGP vs Wei
浏览 0提问于2014-08-16得票数 1
回答已采纳
我有一个实验数据的点阵图，它有4个数据集(即“Config”中的4个不同的子集。下面)，每个数据集根据配置细分，然后还根据风向"Dir“在每个测试类型中分组。我可以根据分组"Dir“在所有图上设置点的形状，但不能根据"Config”设置颜色，例如第一个图全是红点，第二个图全是蓝点等等。我在网格中使用的是xyplot，该图的代码是：
xyplot(ach[,"F"]~ach[,"Ar_sqrt"] | Config., type=c("p"),ach, groups=Dir,
auto.key=list(TRUE,
浏览 0提问于2014-06-18得票数 1
我有一个点阵图，显示三个系列，标有A，B，C。在图表中，我使用panel.text添加了最大值。但是，我怎样才能把这个信息添加到条子中，就在条名旁边呢？(如在顶层:最大value=2.61)
#Libraries used:
library(lattice)
#Data: Create three random walks:
A <- c(rnorm(99), cumsum(rnorm(99)))
B <- c(rnorm(99), cumsum(rnorm(99)))
C <- c(rnorm(99), cumsum(rnorm(99)))
#combi
浏览 2提问于2013-09-25得票数 1
回答已采纳
我有一个简单的格子散点图。两个小组-男性/女性。ID是每个参与者的唯一编号。var1是一个总测试时间。Mean.values是两个数字的向量(性别的手段)。
没有点，包括一个最佳拟合线，所以我想要的是绘制一个趋势线的平均值在每个面板。这两个小组有不同的方式，比如说男性=1分钟，女性=2分钟。
xyplot(var1 ~ ID|Gender, data=DF,
group = Gender,
panel=function(...) {
浏览 3提问于2015-08-01得票数 1
回答已采纳
1回答
我的笔记本电脑的基本信息：
SO: Windows 7 x64
R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
问题：
当我运行以下代码时：
x <- rnorm (100)
z <- z + rnorm(100)
f <- gl(2,50,labels =c("Groups 1" , "Groups 2"))
xyplot (z ~ x | f,
panel = function (x, z, ...) {
panel.xyplot(x,z, ...)
浏览 3提问于2014-01-22得票数 3
回答已采纳
我在创建一个动物园对象的xyplot时遇到了一些麻烦。
我已经创建了一个简化的示例：
z <- zoo(cbind(a=1:4,b=11:14,c=10:13,d=5:8,e=10:7,f=1:4), 1991:1994)
我使用以下代码来创建一个多面板xyplot：
numLbl <- 4
xx <- seq(from = min(time(z)), to = max(time(z)), length.out = numLbl)
xyplot(z[,c(1,2,4,6)],scales = list(x = list(at = time (z), rot = 90)), l
浏览 2提问于2012-01-23得票数 4
1回答
我想绘制几条分位数回归线(quantreg包)，如下所示：
library(quantreg)
data(engel)
attach(engel)
plot(income, foodexp, cex = 0.25, type = "n",
xlab = "Household Income", ylab = "Food Expenditure")
points(income, foodexp, cex = 0.5, col = "blue")
abline(rq(foodexp ~ income, tau =
浏览 0提问于2013-03-23得票数 2
回答已采纳
例如，我应该在这个表达式中添加什么：
xyplot(a~b, groups=abc, data=super,
scales=list(x=list(log=TRUE),y=list(log=TRUE)),
panel = panel.grid)
为了有100而不是10^2？
浏览 0提问于2012-02-15得票数 4
回答已采纳
我正在尝试在lattice图形包中生成一个自定义的xyplot。它是一个由两条线组成的图，并绘制了范围限制。
# Raw dataframe in long format
df <- data.frame(
Response = c(runif(100), rnorm(100)),
Trial = c(rep("A", 100), rep("C", 100)),
Year = 1:10,
Rep = rep(1:10, each = 10))
# Aggregate the data (take mean/min/max
浏览 0提问于2018-07-06得票数 2
当在包含panel.rug()的格子中生成宽图时，凹凸不平的边距中的线条长度在y轴中比x轴长：
library(lattice)
png(width=800, height=400)
xyplot(Fertility ~ Education, swiss, panel = function(x, y,...) {
panel.xyplot(x, y, col=1, pch=16)
panel.rug(x, y, col=1, end= ...)})
dev.off()
我希望那些x-和y-轴的地毯线是相同的长度，不管一个情节的形状(注意:现在地毯线将只有相同的长度，当情节是正
浏览 3提问于2015-10-19得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
实时音视频
对象存储
云直播活动推荐
腾讯云开发者
