我正在尝试下载动物园包以使用时间序列:
install.packages("zoo")
但我收到以下消息:
--- Please select a CRAN mirror for use in this session ---
Warning: unable to access index for repository http://cran.cnr.Berkeley.edu/bin/windows/contrib/2.12
Warning: unable to access index for repository http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.12
Warning messages:
1: In open.connection(con, "r") :
unable to connect to 'cran.r-project.org' on port 80.
2: In getDependencies(pkgs, dependencies, available, lib) :
package ‘zoo’ is not available
我尝试了几个不同的镜子,并且我一直得到同样的错误......动物园在大多数镜子中都没有,或者出了什么问题？
如何for并行运行循环(所以我可以使用我的Windows机器上的所有处理器),结果是3维数组？我现在的代码需要大约一个小时才能运行,如下所示:
guad = array(NA,c(1680,170,15))
for (r in 1:15)
{
name = paste("P:/......",r,".csv",sep="")
pp = read.table(name,sep=",",header=T)
#lots of stuff to calculate x (which is a matrix)
guad[,,r]= x #
}
我一直在寻找相关问题并认为我可以使用,foreach但我找不到将矩阵组合成数组的方法.
我是并行编程的新手,所以任何帮助都将非常感谢!
我有一个每日动物园(xts),有几十年的数据,格式如下:
head(almorol)
1973-10-02 1973-10-03 1973-10-04 1973-10-05 1973-10-06 1973-10-07
183.9 208.2 153.7 84.8 52.5 35.5
我想绘制冬季数据(12月,1月和2月的整月).我找到了xts的子集,所以我想我可以使用以下方法提取所有的Decembers:
x<-apply.yearly(almorol, FUN=last(almorol, "1 month"))
然后在Jan和Feb做类似的事情,但是我收到以下错误:
Error in get(as.character(FUN), mode = "function", envir = envir) :
object 'FUN' of mode 'function' was not found
我可以单独使用apply.yearly和,last(almorol, "1 month")但是当我将它们组合起来时它不起作用.有没有人知道如何对这3个月的时间序列进行子集化？谢谢你的帮助!
我有2个对称矩阵(矩阵的数学意义),一个位置之间的距离(位置用4位数字编码:2030,2059,2095 ......),如下所示:
2030 2059 2095 ...
2030 NA 59328 68464
2059 59328 NA 37196
2095 68464 37196 NA
...
和另一个与位置之间的相关性:
2030 2059 2095...
2030 1.0000000 0.4651804 0.6185849
2059 0.4651804 1.0000000 0.4428746
2095 0.6185849 0.4428746 1.0000000
...
我需要在相关性与距离的关系图中组合这两个矩阵,但不知道如何在R中进行,并且考虑到我有80多个位置我不想手动执行!有谁知道这样做的方法？
谢谢!
我试图从一系列文件中绘制每月每天的长期(几年)平均降水量(pp)的直方图.每个文件都包含从不同位置收集的数据(并且具有不同的代码).我的每个文件都是这样的:
X code year month day pp
1 2867 1945 1 1 0.0
2 2867 1945 1 2 0.0
...
我使用以下代码:
files <- list.files(pattern=".csv")
par(mfrow=c(4,6))
for (i in 1:24) {
obs <- read.table(files[i],sep=",", header=TRUE)
media.dia <- ddply(obs, .(day), summarise, daily.mean<-mean(pp))
codigo <- unique(obs$code)
hist(daily.mean, main=c("hist per day of month", codigo))
}
我在标题中得到了24个直方图,其中包含24个不同的代码,但不是24个不同位置的24个不同直方图,我得到相同的直方图24次(有24个不同的标题).谁能告诉我为什么？谢谢!
我正在用R制作一些图表,我正在将它们应用于Word.我正在将它们作为元文件处理,但Word似乎无法应对它们.R中复制图形的另一个选项是位图,但是当我使用它时,单词中图形的质量很糟糕.
我看到了一些关于改变本网站分辨率的答案,但前提是我保存了我想避免的图表.有没有办法改变复制图形的分辨率？
谢谢,
SBG