require(splines)
?splines
ns(1961:2099)
plot(ns(1961:2099))
plot(ns(1961:2099,df=5))
ns(1961:2099,df=5)
matlines(ns(1961:2099,df=5))
matplot(ns(1961:2099,df=5))
matplot(ns(1961:2099,df=5),type="l")
matplot(ns(1961:2099,knots=1961+(23*(1:5))),type="l")
matplot(1961:2099,ns(1961:2099,knots=1961+(23*(1:5))),type="l")
points(1961+(23*(1:5)),0,pch=15)
points(1961+(23*(1:5)),rep(0,5),pch=15)
matplot(1961:2099,bs(1961:2099,knots=1961+(23*(1:5))),type="l")
require(gam)
?gam
require(mgcv)
?gam
matplot(1961:2099,ns(1961:2099,df=1,type="l")
)
matplot(1961:2099,ns(1961:2099,df=1),type="l")
matplot(1961:2099,ns(1961:2099,df=2),type="l")
cor(ns(1961:2099,df=2))
q()
require(splines)
ns(1961:2099,df=5)
2099-1961
ns(1961:2099,knots=1961+(23*(1:5)))
cor(ns(1961:2099,knots=1961+(23*(1:5))))
matplot(ns(1961:2099,knots=1961+(23*(1:5))),type="l")
matplot(ns(1961:2099,knots=1961+(23*(1:5)),intercept=TRUE),type="l")
matplot(1961:2099,ns(1961:2099,knots=1961+(23*(1:5)),intercept=TRUE),type="l")
ns
require(Hmisc)
require(Hmisc)
?rcs
??rcs
rcspline.eval(1961:2099,knots=1961+(23*(1:5)))
rcspline.eval(1961:2099,knots=1961+(23*(0:6)))
matplot(rcspline.eval(1961:2099,knots=1961+(23*(0:6))),type="l")
matplot(rcspline.eval(1961:2099,knots=1961+(23*(0:6)),pc=TRUE),type="l")
require(rms)
?rcs
rcs
?rms
rcs
rms.trans
rcs(1:10)
demo(rms)
demo(rcs)
demo("rms")
example(rms)
example("rms")
q()
require(rms)
rcspline.eval(1961:2099,knots=1961+(23*(1:5)))
rcspline.eval(1961:2099,knots=1961+(23*(0:6)))
require(rms)
rcs(1:10)
rcs(rcspline.eval(1961:2099,knots=1961+(23*(0:6))))
?rcspline.eval
ns(1961:2099,knots=1961+(23*(0:6)))
ns(1961:2099,knots=1961+(23*(1:5)))
q()
require(Rglimclim)
plot.summary.GLCsim
get.qtab <- plot.summary.GLCsim
edit(get.qtab)
q()
sqrt(0.95*0.05/100)
q()
q()
bibble
q()
x <- 3
q()
5 + 2.6
x <- 5 + 2.6
x
q()
5 + 2.6
x <- 5 + 2.6
x
x <- 8
x
y <- 3.1415
ls()
rm(x)
objects()
rm(y)
rm
x <- c(2.3,1.2,2.4)
x
y <- c(x,9.0,x)
z <- 1:50
z
seq(3,10,2)
z <- c(3,1,4,1,5,9,2,7)
z[3]
z[c(3,2,5)]
x <- c(2,3)
y <- c(1,4,5,6)
2*x
2 + x
y^2
x + y
x*y
4+2*3
2^2*3
4+(2*3)
(4+2)*3
x <- 3
1:x
1:x+2
1:(x+2)
getwd()
setwd("RESEARCH/PROJECTS/VALUE/TRIESTE_2014/PREPARATORY/")
setwd("~/RESEARCH/PROJECTS/VALUE/TRIESTE_2014/PREPARATORY")
for (i in 1:10000) x<- i
for (i in 1:100000) x<- i
for (i in 1:1000000) x<- i
setwd("~/RESEARCH/PROJECTS/VALUE/TRIESTE_2014/DATA")
stop()
source('~/.active-rstudio-document')
?save
save(topo.longs,top.lats,topo.data,file="Topography.rda")
save(topo.longs,topo.lats,topo.data,file="Topography.rda")
source('~/.active-rstudio-document')
source('~/.active-rstudio-document')
5 + 2.6
x <- 5 + 2.6
x
q()
