This repository has been archived on 2023-08-20. You can view files and clone it, but cannot push or open issues or pull requests.
yap-6.3/library/r_session/r_demo.pl
2010-10-16 14:54:48 +01:00

185 lines
5.9 KiB
Prolog

:- use_module(library(lists)). % member/2.
:- use_module(library(readutil)). % read_line_to_codes/2.
:- nl, nl.
:- ( r_bin(Rbin) ->
write( 'Will be using the R found at: ' ), nl,
write( Rbin ), nl, nl
;
write( 'This session cannot locate an R executable. Please register the location' ), nl,
write( 'of the R executalbe using r_bin/1 before you can run the demos.' ), nl, nl
).
:- write( 'Demo predicates for R (r_session) package.' ), nl.
:- write( 'See r_demo_1/0,...,r_demo_10/0.' ), nl.
:- write( 'The goal r_demo/0 is a shorthand for r_demo_1/0,...,r_demo_7/0 which are the main demos.' ), nl, nl.
:- write( 'r_demo_all/0 and r_demo_clauses/0 for r_demo_1,...,r_demo_10.' ), nl.
:- write( 'which include demos for some non-basic features.' ), nl.
:- write( 'You need to look at the sources before running r_demo_8,9 and 10.' ).
:- nl, nl.
r_demo :-
nl, nl,
Rdemos = [r_demo_1,r_demo_2,r_demo_3,r_demo_4,r_demo_5,r_demo_6,r_demo_7],
r_demo( Rdemos, false ).
r_demo_all:-
nl, nl,
Rdemos = [r_demo_1,r_demo_2,r_demo_3,r_demo_4,r_demo_5,r_demo_6,r_demo_7,r_demo_8,r_demo_9,r_demo_10],
r_demo( Rdemos, false ).
r_demo_clauses :-
nl, nl,
Rdemos = [r_demo_1,r_demo_2,r_demo_3,r_demo_4,r_demo_5,r_demo_6,r_demo_7,r_demo_8,r_demo_9,r_demo_10],
r_demo( Rdemos, true ).
r_demo( Rdemos, Clauses ) :-
member(Wh, Rdemos ),
write( doing-Wh ), nl,
( Clauses == true ->
write( 'Clauses: ' ), nl,
findall( Wh-Body, (clause(Wh,Body),
portray_clause((Wh:-Body)), nl), _ )
;
true
),
( call(Wh) ->
true
;
write( 'Demo ended with failure.' ), nl
),
nl, nl,
fail.
r_demo( _Rdemos, _ ) :-
write( done ), nl.
r_demo_1 :-
write( 'Demo: basic vector interactions.' ), nl, nl,
r_open,
r_in( x <- c(10.4, 5.6, 3.1, 6.4, 21.7) ),
( r_out( print(x), Lines ), r_lines_print( Lines ), fail; true ),
r_print( x ),
r_in( x ),
r_in( (y <- c(6,5,4,3,2,1); y) ), % The extra paranthesis are only
% needed for Yap.
r_in( Z <- c(10.4, 5.6, 3.1, 6.4, 21.7) ),
write( z(Z) ), nl,
r_close.
r_demo_2 :-
write( 'Demo: plots (screen and postscript).' ), nl, nl,
r_open,
r_in( y <- rnorm(50) ),
r_print( y ),
r_in( Y <- y ),
write( y(Y) ), nl,
r_in( x <- rnorm(y) ),
r_print( x ),
r_in( X <- x ),
write( x(X) ), nl,
r_in( x11(width=5,height=3.5) ),
r_in( plot(x,y)),
write( 'Press Return to continue...' ), nl,
current_prolog_flag( version, V ),
( integer(V) -> User = current_input % SWI Prolog
; User = user ),
read_line_to_codes( User, _ ),
r_in( 'dev.off()' ),
r_in( 'postscript(file="x_vs_y.eps")' ),
r_in( plot(x,y)),
r_in( 'dev.off()' ),
r_close,
!, % Swi leaves a backtracking point at read_line_to_codes/2
write( 'Check that file x_vs_y.eps has been created.' ), nl.
r_demo_3 :-
write( 'Demo: aliases.' ), nl, nl,
r_open( [alias(mamonaku)] ),
( current_r_session(Alias),write(session(Alias)),nl, fail; true ),
r_in( mamonaku, x <- c(10.4, 5.6, 3.1, 6.4, 21.7) ),
r_print( mamonaku, x ),
r_close( mamonaku ).
r_demo_4 :-
write( 'Demo: history.' ), nl, nl,
r_open,
( r_history(A,B), write(history(A,B)), nl, fail; true ),
r_in( x <- c(10.4, 5.6, 3.1, 6.4, 21.7) ),
( r_out( print(x), Lines ), r_lines_print( Lines ), fail; true ),
r_print( x ),
( r_history(C,D), write(h(C,D)), nl, fail; true ),
r_close.
r_demo_5 :-
write( 'Demo: calls to R functions.' ), nl, nl,
r_open,
r_in( i <- 0:14 ),
r_print( i ),
r_in( I <- i ),
write( 'I'(I) ), nl,
r_in( x <- i/10 ),
r_in( y <- c(176.547,194.2552,218.5462,441.3706,795.786,1190.8606,1321.0128,1326.4694,1437.3068,1364.6906,1343.768,1513.7298,1553.8264,1554.1748,1549.399) ),
r_print( (integrate(splinefun(x,y), 0.2, 0.4)) ),
r_close.
r_demo_6 :-
write( 'Demo: copying output and error on to file.' ), nl, nl,
r_open( [copy('rec_both.txt',both)] ),
r_in( x <- c(10.4, 5.6, 3.1, 6.4, 21.7) ),
r_print( x ),
write( 'Check that file rec_both.txt has been created.' ), nl,
r_close.
r_demo_7 :-
write( 'Demo: error on R.' ), nl, nl,
r_open( [at_r_halt(restart)] ),
r_in( x <- c(10.4, 5.6, 3.1, 6.4, 21.7) ),
r_print( x ),
r_print( y ),
r_print( x ),
r_close.
%%% Cut-off
r_demo_8 :-
write( 'Demo: reinstate on halt.' ), nl,
write( 'This is no longer valid.' ), nl, nl,
r_open( [at_r_halt(reinstate)] ),
r_in( x <- c(10.4, 5.6, 3.1, 6.4, 21.7) ),
r_print( y ),
% here slave dies
% and r_session tries to restar_demo it and replay all commands.
r_print( x ),
% succeeds
r_close.
/* change 192.168.0.* to a host in your domain before running the following. */
r_demo_9 :-
write( 'Demo: ssh on a machine with R on a different location.' ), nl, nl,
r_open( [ssh('192.168.0.3')] ),
r_in( I <- 0:14 ),
write( 'I'(I) ), nl,
r_close.
r_demo_10 :-
write( 'Demo: ssh on a machine with explicit set of the remote R location.' ),
nl, nl,
r_bin( '/usr/local/users/nicos/local/bin/R' ),
r_open( [ssh('192.168.0.3')] ),
r_in( I <- 0:14 ),
write( 'I'(I) ), nl,
r_close.
/*
% You can replace any of the above r_open/0,1, with one of the following
r_open( [with(restore)] ).
% do not use --no-restore on the R flags
r_open( [copy(copied,both)] ).
% copy both input and output to file copied
r_open( [at_r_halt(restar_demo),alias(mamunaku),copy(copied_in,in)] ).
% copy input to file copied_in
r_open( [at_r_halt(restar_demo),alias(mamunaku),copy(copy_out,out)] ).
% copy output to file copied_out
*/