* Problem with R source blocks
@ 2015-06-11 23:12 Thomas S. Dye
2015-06-12 18:50 ` Charles C. Berry
0 siblings, 1 reply; 6+ messages in thread
From: Thomas S. Dye @ 2015-06-11 23:12 UTC (permalink / raw)
To: emacs-orgmode
Aloha all,
Recently, some previously working R code (below) stopped working correctly.
It is not possible to produce a graphic without :session. I get a red
message in the pdf file "arguments imply differing numbers of rows: 86,
79".
With :session, I first get one of two error messages in red in the pdf
file, either "Results must be all atomic, or all data frames" or
"subscript out of bounds".
With the "subscript out of bounds" error, if I run the source code block
again, I get the expected graphic.
With the "Results must be all atomic ..." message, if I re-run the
source code block several times, I see first an error "arguments imply
differing numbers of rows: 67, 60", then on the next run "arguments imply
differing numbers of rows: 86, 79" (like running outside a :session),
then the expected graphic. Subsequent runs all produce the expected
graphic.
If I jump to the session, then print(g) works as expected, without any
error messages.
I haven't the faintest idea what might be happening, or how to debug.
#+begin_src emacs-lisp
(org-version)
#+end_src
#+results:
: 8.3beta
#+begin_src emacs-lisp
(emacs-version)
#+end_src
#+results:
: GNU Emacs 24.5.1 (x86_64-apple-darwin13.4.0, NS apple-appkit-1265.21)
: of 2015-04-18 on tennine-slave.macports.org
All the best,
Tom
,-------------------------------------------------------------------------
| #+name: dated-interval-plot
| #+header: :file ~/temp/test-plot.pdf
| #+header: :var bcal = heiau-dates
| #+header: :var base.font = 11
| #+header: :var theme = "bw"
| #+header: :var xlabel = "Years After Polynesian Discovery"
| #+header: :var ylabel = "Probability"
| #+header: :var rows=4
| #+header: :var cols=4
| #+header: :var xmin=0
| #+header: :var xmax=1000
| #+header: :width 7 :height 5
| #+header: :results graphics
| #+BEGIN_SRC R :session
| library(ggplot2)
| library(ggthemes)
| bcal.df <- data.frame(cal=numeric(0),prob=numeric(0),label=character(0))
| for (i in seq(dim(bcal)[1]))
| {
| t <- read.csv(file=bcal[i,1])
| names(t) <- c("cal","prob")
| t <- cbind(t,label=rep(bcal[i,2], dim(t)[1]))
| bcal.df <- rbind(bcal.df,t)
| }
| switch(theme,
| bw =, BW = theme_set(theme_bw(base_size=base.font)),
| grey =, gray = theme_set(theme_grey(base_size=base.font)),
| tufte = theme_set(theme_tufte(base_size=base.font)),
| calc = theme_set(theme_calc(base_size=base.font)),
| econ =, economist = theme_set(theme_economist(base_size=base.font)),
| excel = theme_set(theme_excel(base_size=base.font)),
| few = theme_set(theme_few(base_size=base.font)),
| five = theme_set(theme_fivethirtyeight(base_size=base.font)),
| google =, gdocs = theme_set(theme_gdocs(base_size=base.font)),
| ## highchart =, hc = theme_set(theme_hc(base_size=base.font)),
| inverted = theme_set(theme_igray(base_size=base.font)),
| solarized = theme_set(theme_solarized(base_size=base.font)),
| solarized2 = theme_set(theme_solarized_2(base_size=base.font)),
| ## solid = theme_set(theme_solid(base_size=base.font)),
| stata = theme_set(theme_stata(base_size=base.font,scheme="s2mono")),
| wsj = theme_set(theme_wsj(base_size=base.font)),
| theme_set(theme_bw(base_size=base.font)))
| g <- ggplot(bcal.df, aes(x=cal, y=prob))
| g <- g + geom_bar(stat='identity') + xlab(xlabel) + xlim(xmin,xmax)
| g <- g + ylab(ylabel) + facet_wrap(~ label, nrow=rows, ncol=cols)
| print(g)
| #+END_SRC
`-------------------------------------------------------------------------
--
T.S. Dye & Colleagues, Archaeologists
735 Bishop St, Suite 315, Honolulu, HI 96813
Tel: 808-529-0866, Fax: 808-529-0884
http://www.tsdye.com
^ permalink raw reply [flat|nested] 6+ messages in thread
* Re: Problem with R source blocks
2015-06-11 23:12 Problem with R source blocks Thomas S. Dye
@ 2015-06-12 18:50 ` Charles C. Berry
2015-06-12 21:55 ` Thomas S. Dye
0 siblings, 1 reply; 6+ messages in thread
From: Charles C. Berry @ 2015-06-12 18:50 UTC (permalink / raw)
To: Thomas S. Dye; +Cc: emacs-orgmode
On Thu, 11 Jun 2015, Thomas S. Dye wrote:
> Aloha all,
>
> Recently, some previously working R code (below) stopped working correctly.
>
> It is not possible to produce a graphic without :session. I get a red
> message in the pdf file "arguments imply differing numbers of rows: 86,
> 79".
>
> With :session, I first get one of two error messages in red in the pdf
> file, either "Results must be all atomic, or all data frames" or
> "subscript out of bounds".
>
> With the "subscript out of bounds" error, if I run the source code block
> again, I get the expected graphic.
>
> With the "Results must be all atomic ..." message, if I re-run the
> source code block several times, I see first an error "arguments imply
> differing numbers of rows: 67, 60", then on the next run "arguments imply
> differing numbers of rows: 86, 79" (like running outside a :session),
> then the expected graphic. Subsequent runs all produce the expected
> graphic.
>
> If I jump to the session, then print(g) works as expected, without any
> error messages.
>
> I haven't the faintest idea what might be happening, or how to debug.
>
The error messages are coming from R. cbind will say:
: Error in data.frame(..., check.names = FALSE) :
: arguments imply differing number of rows: 2, 3
if you try cbind( data.frame(diag(2)), 1:3 )
Also, you are feeding lots of vars to the R src block, so that is a good
place to look for trouble.
You could try to instrument `org-babel-R-evaluate-external-process' and
see if anything jumps out at you.
But if not, ...
To debug, I'd try to replace your block with one that is utterly self
contained - no need for external vars begin passed in from Org mode - and
that invokes ggplot2 to produce a result. Once that works, define :var
args and print them to a temp file from inside that R src block using
capture.output(print(your.var1),...,print(your.varK),
file = "myTemp.Rout")
Once that works, I'd add the first lines of code and print out
intermediate values from inside your for() loop, but before the cbind()
again using capture.output().
I expect by then the issue should be revealed.
HTH,
Chuck
^ permalink raw reply [flat|nested] 6+ messages in thread
* Re: Problem with R source blocks
2015-06-12 18:50 ` Charles C. Berry
@ 2015-06-12 21:55 ` Thomas S. Dye
2015-06-12 22:33 ` Charles C. Berry
0 siblings, 1 reply; 6+ messages in thread
From: Thomas S. Dye @ 2015-06-12 21:55 UTC (permalink / raw)
To: Charles C. Berry; +Cc: emacs-orgmode
Hi Chuck,
"Charles C. Berry" <ccberry@ucsd.edu> writes:
> The error messages are coming from R. cbind will say:
>
> : Error in data.frame(..., check.names = FALSE) :
> : arguments imply differing number of rows: 2, 3
I think the problem is in ggplot2, rather than my for loop:
,---------------------------------------------------------------------------
| [1] "ERROR : arguments imply differing number of rows: 76, 69"
| Warning messages:
| 1: In loop_apply(n, do.ply) :
| Removed 3 rows containing missing values (position_stack).
| 2: In loop_apply(n, do.ply) :
| Removed 1 rows containing missing values (position_stack).
| 3: In loop_apply(n, do.ply) :
| Removed 3 rows containing missing values (position_stack).
| null device
| 1
| > [1] "org_babel_R_eoe"
| > summary(bcal.df$label)
| 50-50-17-01088 Hapaiali`i KAL-1 Makolea KH2-2
| 100 60 82 69 72
| KAL-26 KAL-27 KHO-1 KOL-2 MKI-122
| 79 72 65 63 54
| MKI-123 MKI-124 MKI-125 KH1-3 KAL-24
| 62 69 65 65 57
| Cook
| 41
`---------------------------------------------------------------------------
Note that the warnings removed 7 rows, which is the difference pointed
out by the error.
I don't know what position_stack might be, but it doesn't map directly
to the data frame I've constructed, which lacks an entry with 76 rows.
I'm guessing this brouhaha has to do with facet_wrap.
In any case, Org Babel seems to be behaving correctly.
I've had this code around for years and have run it without problem many
hundreds of times.
Is ess a potential source of the problem? I upgraded that recently.
Thanks for your help,
Tom
--
Thomas S. Dye
http://www.tsdye.com
^ permalink raw reply [flat|nested] 6+ messages in thread
* Re: Problem with R source blocks
2015-06-12 21:55 ` Thomas S. Dye
@ 2015-06-12 22:33 ` Charles C. Berry
2015-06-12 23:54 ` Thomas S. Dye
0 siblings, 1 reply; 6+ messages in thread
From: Charles C. Berry @ 2015-06-12 22:33 UTC (permalink / raw)
To: Thomas S. Dye; +Cc: emacs-orgmode
On Fri, 12 Jun 2015, Thomas S. Dye wrote:
> Hi Chuck,
>
> "Charles C. Berry" <ccberry@ucsd.edu> writes:
>
>> The error messages are coming from R. cbind will say:
>>
>> : Error in data.frame(..., check.names = FALSE) :
>> : arguments imply differing number of rows: 2, 3
>
> I think the problem is in ggplot2, rather than my for loop:
>
> ,---------------------------------------------------------------------------
> | [1] "ERROR : arguments imply differing number of rows: 76, 69"
> | Warning messages:
> | 1: In loop_apply(n, do.ply) :
> | Removed 3 rows containing missing values (position_stack).
> | 2: In loop_apply(n, do.ply) :
> | Removed 1 rows containing missing values (position_stack).
> | 3: In loop_apply(n, do.ply) :
> | Removed 3 rows containing missing values (position_stack).
> | null device
Those messages come from plyr, which ggplot2 uses.
[snip]
>
> In any case, Org Babel seems to be behaving correctly.
>
> I've had this code around for years and have run it without problem many
> hundreds of times.
>
> Is ess a potential source of the problem? I upgraded that recently.
>
I don't think so.
I'd rewrite the src block with
save.image(file="setup-for-ggplot.RData")
before the graphics commands.
Restart emacs and execute the src block.
Start a fresh R session, then type
load("setup-for-ggplot.RData")
and type/paste in all the graphics commands and see what happens. If the
error repeats and you can verify that each object is as you expected it to
be, then this is not an Org mode problem.
HTH,
Chuck
^ permalink raw reply [flat|nested] 6+ messages in thread
* Re: Problem with R source blocks
2015-06-12 22:33 ` Charles C. Berry
@ 2015-06-12 23:54 ` Thomas S. Dye
2015-06-13 0:26 ` Thomas S. Dye
0 siblings, 1 reply; 6+ messages in thread
From: Thomas S. Dye @ 2015-06-12 23:54 UTC (permalink / raw)
To: Charles C. Berry; +Cc: emacs-orgmode
"Charles C. Berry" <ccberry@ucsd.edu> writes:
> I'd rewrite the src block with
>
> save.image(file="setup-for-ggplot.RData")
>
> before the graphics commands.
>
> Restart emacs and execute the src block.
>
> Start a fresh R session, then type
>
> load("setup-for-ggplot.RData")
>
> and type/paste in all the graphics commands and see what happens. If
> the error repeats and you can verify that each object is as you
> expected it to be, then this is not an Org mode problem.
Done. I can verify this is not an Org mode problem.
I'll report it to the ggplot2 folks.
Thanks for your help.
All the best,
Tom
--
Thomas S. Dye
http://www.tsdye.com
^ permalink raw reply [flat|nested] 6+ messages in thread
* Re: Problem with R source blocks
2015-06-12 23:54 ` Thomas S. Dye
@ 2015-06-13 0:26 ` Thomas S. Dye
0 siblings, 0 replies; 6+ messages in thread
From: Thomas S. Dye @ 2015-06-13 0:26 UTC (permalink / raw)
To: Charles C. Berry; +Cc: emacs-orgmode
tsd@tsdye.com (Thomas S. Dye) writes:
> I'll report it to the ggplot2 folks.
Too late. The bug in plyr was corrected a week ago.
> library(devtools)
> devtools::install_github("hadley/plyr")
> library(plyr)
> library(ggplot2)
> ...
Thanks again for your help.
All the best,
Tom
--
Thomas S. Dye
http://www.tsdye.com
^ permalink raw reply [flat|nested] 6+ messages in thread
end of thread, other threads:[~2015-06-13 0:26 UTC | newest]
Thread overview: 6+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2015-06-11 23:12 Problem with R source blocks Thomas S. Dye
2015-06-12 18:50 ` Charles C. Berry
2015-06-12 21:55 ` Thomas S. Dye
2015-06-12 22:33 ` Charles C. Berry
2015-06-12 23:54 ` Thomas S. Dye
2015-06-13 0:26 ` Thomas S. Dye
Code repositories for project(s) associated with this public inbox
https://git.savannah.gnu.org/cgit/emacs/org-mode.git
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).