Re: [R] Multiple CPU HowTo in Linux?

2010-09-14 Thread Johnson, Cedrick W.

 ?multicore perhaps

On 09/14/2010 10:01 AM, Edwin Groot wrote:

Hello all,
I upgraded my R workstation, and to my dismay, only one core appears to
be used during intensive computation of a bioconductor function.
What I have now is two dual-core Xeon 5160 CPUs and 10 GB RAM. When I
fully load it, top reports about 25% user, 75% idle and 0.98 short-term
load.
The archives gave nothing helpful besides mention of snow. I thought of
posting to HPC, but this system is fairly modest WRT processing power.
Any pointers of where to start?
---
#Not running anything at the moment

sessionInfo()

R version 2.11.1 (2010-05-31)
x86_64-pc-linux-gnu

locale:
  [1] LC_CTYPE=en_GB.UTF-8   LC_NUMERIC=C
  [3] LC_TIME=en_GB.UTF-8LC_COLLATE=en_GB.UTF-8
  [5] LC_MONETARY=C  LC_MESSAGES=en_GB.UTF-8
  [7] LC_PAPER=en_GB.UTF-8   LC_NAME=C
  [9] LC_ADDRESS=C   LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_GB.UTF-8 LC_IDENTIFICATION=C

attached base packages:
[1] stats graphics  grDevices utils datasets  methods   base


loaded via a namespace (and not attached):
[1] tools_2.11.1
---
$ uname -a
Linux laux29 2.6.26-2-amd64 #1 SMP Sun Jun 20 20:16:30 UTC 2010 x86_64
GNU/Linux
---
Thanks for your help,
Edwin


__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] R time series analysis

2010-09-05 Thread Johnson, Cedrick W.
 You also may want to look at auto.arima in the 'forecast' package, 
which will return the best ARIMA model based on AIC/AICc/BIC values


?auto.arima

hth
c

On 09/05/2010 06:02 PM, Stephan Kolassa wrote:

Hi,

basically, you know 5 periods later. If you use a good error measure, 
that is.


I am a big believer in AIC for model selection. I believe that arima() 
also gives you the AIC of a fitted model, or try AIC(arima1).


Other ideas include keeping a holdout sample or some such.

I'd recommend looking at a time series textbook.

HTH,
Stephan


Am 05.09.2010 22:37, schrieb lord12:


How do you evaluate the predictive models? For example if I have:

arima1 = arima(training, order = c(1,1,1))
arima2 = arima(training, order = c(0,0,0))
x.fore = predict(arima1, n.ahead=5)
x.fore1 = predict(arima2, n.ahead = 5)

How do I know which arima model is better for prediction?


__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide 
http://www.R-project.org/posting-guide.html

and provide commented, minimal, self-contained, reproducible code.


__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] online documentation only?

2010-08-30 Thread Johnson, Cedrick W.

try help.start()

that starts a local help process (within R) and open your browser to 
that local location.


-c

On 08/30/2010 03:53 PM, Bob McCall wrote:


Greetings:

I recently installed R 2.11.1 for windows. It seems that there is only
online help now. Is there any way to get the local docs? I don't have
always-on high-speed internet.

thanks,
Bob


__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Please sigh me out from the R-help. Thank you!

2010-08-29 Thread Johnson, Cedrick W.

https://stat.ethz.ch/mailman/listinfo/r-help

At the bottom of the webpage.

-c

On 08/29/2010 03:58 AM, Kaigang Li wrote:






[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Axes=F and plotting dual y axes

2010-07-28 Thread Johnson, Cedrick W.
Howdy. Been running into a bit of trouble with plotting. Seems that 
axes=F is not working. Whenever I plot (either a dataframe or xts/zoo 
series) and I set axes=F along with xlab/ylab= I still get the default 
axes printed in my chart. Consider this:



#Create some sample data, both 50 units of blah
series2 = c(1:50)
series1 = rep(25:74)

testdf1 = as.data.frame(series1)
testdf1$series2 = series2

As a note, I converted my original xts/zoo dataset into a dataframe 
thinking it could be weirdness on the part of that. I just did this here 
to have something reproducible since it's not feasible to put my entire 
original dataset in the email.


plot(testdf1[,1], main=Woo, col=rich12equal, xlab=, ylab=, axes=F, 
type='l')


It still displays the x and y axes. What I'm trying to do:

par(new=T)

plot(testdf1[,2], ...)

and from there I can manually label the axes, title, etc.

Any idea what's going on? I'm using R 2.10.1 on ubuntu (I probably need 
to upgrade via debian packages soon)



Regards,
CJ

ps: I am trying to accomplish this, as seen here: 
http://blog.earlh.com/index.php/2009/07/multiple-y-axes-in-r-plots-part-9-in-a-series/

 -- but axes=F is giving me trouble.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Axes=F and plotting dual y axes

2010-07-28 Thread Johnson, Cedrick W.
That worked. Stupid me forgot that I had the stock ticker 'F' assigned 
in my workspace.


Well.. guess I'll hit myself with a 2x4 now.. Thanks for your help guys..


-c

On 07/28/2010 12:37 PM, Ben Bolker wrote:

Johnson, Cedrick W.cedrickat  cedrickjohnson.com  writes:



Howdy. Been running into a bit of trouble with plotting. Seems that
axes=F is not working. Whenever I plot (either a dataframe or xts/zoo
series) and I set axes=F along with xlab/ylab= I still get the default
axes printed in my chart.


   A quick guess: what happens if you try axes=FALSE instead? (i.e., have
you assigned 'F' a value in your workspace?)

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.