小编Adr*_*ian的帖子

数据来自这里

library(nlme)
dat0 <- read.table("aids.dat2",head=T)
dat1 <- dat0[dat0$day<=90, ]   # use only first 90-day data
dat2 <- dat1[!apply(is.na(dat1),1,any),]  # remove missing data 

# Next, let's treat the data as longitudinal (or grouped) data 
aids.dat <- groupedData(lgcopy ~ day | patid, data=dat2)

# A NLME model fit, with random effects on all 4 parameters 
start <- c(10,0.5,6,0.005)  # starting value 

aids.dat$log10copy = log10(aids.dat$lgcopy)

nlme.fit <- nlme(log10copy ~ exp(p1-b1*day) + exp(p2-b2*day + 1),
                 fixed = list(p1 ~ 1, b1 ~ 1, …

当我有data.frame对象时,我可以简单地做View(df),然后我会看到data.frame一个漂亮的表(即使我看不到所有的行,我仍然知道我的数据包含哪些变量).但是当我有一个list对象时,同样的命令不起作用.当列表很大时,我不知道列表是什么样的.我试过head(mylist)但我的控制台根本无法一次显示所有信息.查看R中的大型列表有什么有效的方法？

>dput(dummy)
c(NA, "MGM", NA, "M", NA)

>dummy
[1] NA    "MGM" NA    "M"   NA 

>ifelse(dummy == "NA", 0, 1)
NA  1 NA  1 NA

我有一个字符向量,我想将NA值替换为0,其他所有值都替换为1.但是,我不知道这些NA字符到底是什么？正如你在上面看到的,他们不是"NA".当我尝试时,dput(dummy[1])我得到了NA_character_.我该怎么做才能让ifelse语句识别出那些NA？

    string = "ABC3JFD456"

假设我有上面的字符串,我希望找到字符串中的第一个数字并存储其值.在这种情况下,我想存储值3(因为它是字符串中的第一个出现的数字).grepl("\\d", string)只返回一个逻辑值,但不会告诉我任何关于第一个数字的位置或内容的信息.我应该使用哪个正则表达式来查找第一个数字的值？

> df = data.frame(A = c(1, 2, 3), B = c(3, 2, 2), C = c(3, 2, 1)); df
  A B C
1 1 3 3
2 2 2 2
3 3 2 1
> df2 = data.frame(A = c(1, 2, 3), B = c(1, 2, 3), C = c(1, 2, 3)); df2
  A B C
1 1 1 1
2 2 2 2
3 3 3 3

我想知道 data.frame 中的所有列是否都相同。对于df，它应该是 FALSE，而对于 ，df2它应该是 TRUE。

v1 = c(2, 2.01, 2.02, 2.03, 2.04, 2.05, 2.06, 2.07, 2.08, 2.09, 2.1, 
  2.11, 2.12, 2.13, 2.14, 2.15, 2.16, 2.17, 2.18, 2.19, 2.2, 2.21, 
  2.22, 2.23, 2.24, 2.25, 2.26, 2.27, 2.28, 2.29, 2.3, 2.31, 2.32, 
  2.33, 2.34, 2.35, 2.36, 2.37, 2.38, 2.39, 2.4, 2.41, 2.42, 2.43, 
  2.44, 2.45, 2.46, 2.47, 2.48, 2.49, 2.5, 2.51, 2.52, 2.53, 2.54, 
  2.55, 2.56, 2.57, 2.58, 2.59, 2.6, 2.61, 2.62, 2.63, 2.64, 2.65, 
  2.66, 2.67, 2.68, 2.69, 2.7, 2.71, 2.72, 2.73, 2.74, 2.75, 2.76, 
  2.77, …

我有一个简单的x, ydata.frame.

mydata <- data.frame(days = 1:96, risk = c(5e-09, 5e-09, 5e-09, 1e-08, 4e-08, 6e-08, 9e-08, 1.5e-07, 4.2e-07, 
                                           7.2e-07, 1.02e-06, 1.32e-06, 1.66e-06, 2.19e-06, 2.76e-06, 3.32e-06, 
                                           3.89e-06, 4.55e-06, 5.8e-06, 7.16e-06, 8.51e-06, 9.85e-06, 1.138e-05, 
                                           1.396e-05, 1.672e-05, 1.947e-05, 2.222e-05, 2.521e-05, 2.968e-05, 
                                           3.439e-05, 3.909e-05, 4.378e-05, 4.894e-05, 5.697e-05, 6.546e-05, 
                                           7.392e-05, 8.236e-05, 9.16e-05, 0.00010573, 0.00012063, 0.00013547, 
                                           0.00015025, 0.00016642, 0.00019127, 0.00021743, 0.00024343, 0.00026924, 
                                           0.00029818, 0.00034681, 0.00039832, 0.00044932, 0.00049976, 0.0005451, 
                                           0.00056293, 0.00057586, 0.00058838, 0.0006005, 0.00061562, 0.00065079, 
                                           0.00068845, 0.00072508, 0.00076062, 0.00079763, 0.00084886, 0.00090081, 
                                           0.0009507, 0.00099844, 0.00104427, 0.00108948, …

Ancestry <- data.frame(Race = c("European", "African American", "Asian", "Hispanic", "Other"), Proportion = c(40, 30, 10, 15, 5))
Ancestry %>% 
ggplot(aes(y = Proportion, fill = Race)) + 
  geom_bar(stat="identity", colour="white")

运行上面的代码给我以下错误：

Warning in min(x, na.rm = na.rm) :
  no non-missing arguments to min; returning Inf
Warning in max(x, na.rm = na.rm) :
  no non-missing arguments to max; returning -Inf
Warning in min(diff(sort(x))) :
  no non-missing arguments to min; returning Inf
Warning in x - width/2 :
  longer object length is not a …

我在 Windows 中使用 Rstudio。没有红色八角形供我点击。我试过在控制台中按 ESC 和 Ctrl + C 和 Ctrl + Z 但这些都没有用。

library(ggplot2)
df2 <- data.frame(supp=rep(c("VC", "OJ"), each=3),
                dose=rep(c("D0.5", "D1", "D2"),2),
                len=c(6.8, 15, 33, 4.2, 10, 29.5))
head(df2)
ggplot(data=df2, aes(x=dose, y=len, fill=supp)) +
  geom_bar(stat="identity")

我有一个简单的堆叠条形图，我想手动更改颜色。更具体地说，我想翻转用于的颜色fill = supp（即 OJ 的青色）。我已经尝试向其中添加一个color = ...参数，geom_bar但这只是勾勒出条形图而不是将它们着色。