search for: gpetris

...C code and calling it from R. (I am also new to shared libraries...) My problem is that my C code uses a C function contained in another library (".a", not ".so" - is that the problem?). This is how I compile the file with the functions I want to call from R: cd /export/home/gpetris/ R SHLIB tryit.c cc -I/export/home/share/R/R-1.4.0/lib/R/include "-I/export/home/gpetris/C" -I/usr/local/include -I/export/home/share/include -KPIC -xO5 -dalign -xlic_lib=sunperf -c tryit.c -o tryit.o cc -G -o tryit.so tryit.o "-L/export/home/gpetris/lib/SUN/32bit -lutil -lrand&q...

...ne a function my.mean that returns 0 for empty vectors) Does anybody have a suggestion for an easy/fast/elegant workaround? Thanks in advance, Giovanni -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________] -.-.-.-.-.-.-.-.-.-...

Hello, Recently updated Ubuntu to 10.04, I have installed the r-base and r-base-dev packages from CRAN with sudo apt-get install r-base But now when I want to start R I get the following error: gpetris at definetti:~$ R /usr/lib/R/bin/exec/R: symbol lookup error: /usr/local/lib/libreadline.so.6: undefined symbol: PC Any clue??? Has anybody seen that before? I have tried to reinstall the readline libraries, both as a ubuntu package and from sources, but the error did not go away.... Thanks in ad...

...year 2003 month 10 day 08 language R Thank you in advance, Giovanni -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________]

I am trying to read an ftable using read.ftable, but I get the following error message: > jobSatTable <- read.ftable("http://definetti.uark.edu/~gpetris/stat5333/jobSatisfaction.dat",skip=2) Error in seek(file, where = 0) : no applicable method for "seek" In addition: Warning messages: 1: no non-missing arguments to max; returning -Inf 2: no non-missing arguments to max; returning -Inf I also tried to play with the argument col.v...

...bstitute(n == k, list(k = k))) > ### this works > plot(0,0, main = substitute(N[eff] == neff, list(neff = n.eff[k]))) > ### this doesn't work > plot(0,0, main = substitute(n == k * N[eff] == neff, list(k = k, neff = n.eff[k]))) Thanks in advance, Giovanni -- Giovanni Petris <GPetris at uark.edu> Associate Professor Department of Mathematical Sciences University of Arkansas - Fayetteville, AR 72701 Ph: (479) 575-6324, 575-8630 (fax) http://definetti.uark.edu/~gpetris/

...similar result (with a tiny size instead of a footnote size, but that is not the point) with knitr I included in the .Rnw file the lines \ifdefined\knitrout \renewenvironment{knitrout}{\begin{tiny}\slshape}{\end{tiny}} \else \fi What I get is something like this: http://definetti.uark.edu/~gpetris/knitr-output.pdf The problem is that the prompt sign '>' did not change either size or shape. So I am probably not using the correct approach. How should I be doing? Also, what do I have customize, and how, in order to get the extra left margin? Final question: how can I get the cont...

...or optimizer + lower = c(1e-7, 0)) \alert<2>{# V must be positive} > mleOut$convergence \alert<3>{# always check this!!!} [1] 0 \end{semiverbatim} How can I get a similar effect using Sweave? Thank you in advance! Best, Giovanni Petris -- Giovanni Petris <GPetris at uark.edu> Associate Professor Department of Mathematical Sciences University of Arkansas - Fayetteville, AR 72701 Ph: (479) 575-6324, 575-8630 (fax) http://definetti.uark.edu/~gpetris/

Is there a way I can have R automatically execute the commands in a source file without ever having to use R interactively? If so, what arguments should I pass to the UNIX call to R? I need to do this to run several R jobs in parallel. An alternative may be to have R and an R package behave as a stand-alone application that can be called from the UNIX command line. Is there any documentation on

...minor 0.1 year 2004 month 11 day 15 language R -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________]

...] 9 7 10 8 5 10 7 7 (55,59] (59,63] 10 5 Thanks in advance, Giovanni -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________]

...ot;, "YYY")) Error in methods::as(data[[i]], colClasses[i]) : no method or default for coercing "character" to "int" Could anybody please tell me what I am doing wrong and how I could read my data into R? Thanks in advance, Giovanni -- Giovanni Petris <GPetris at uark.edu> Associate Professor Department of Mathematical Sciences University of Arkansas - Fayetteville, AR 72701 Ph: (479) 575-6324, 575-8630 (fax) http://definetti.uark.edu/~gpetris/

...8 12.00 1.00 7.71 ... Could anybody tell me what is an efficient way of computing the averages of the simulated values? Thank you in advance, Giovanni -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________] -.-.-.-.-.-.-.-.-.-...

...data(quakes) xyplot> Depth <- equal.count(quakes$depth, number = 8, overlap = 0.1) xyplot> xyplot(lat ~ long | Depth, data = quakes) Bus error -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________]

...data(quakes) xyplot> Depth <- equal.count(quakes$depth, number = 8, overlap = 0.1) xyplot> xyplot(lat ~ long | Depth, data = quakes) Bus error -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________]

On Thu, Jan 21, 2010 at 9:07 AM, Giovanni Petris <GPetris at uark.edu> wrote: > > I have found a good refernce to be "S Programming" by Venables and > Ripley. I'll take a look at this book. Since S and R are not completely the same, there are delicate differences between S and R, which an S book may be confusing if I want to le...

...also on my Mac, and vice versa. I imagine this to be a fairly common problem, so I would like to ask if anybody has suggestions to share about it. Is there a way to make the synchronization automatic? Painless? Thank you in advance for the suggestions. Best, Giovanni -- Giovanni Petris <GPetris at uark.edu> Associate Professor Department of Mathematical Sciences University of Arkansas - Fayetteville, AR 72701 Ph: (479) 575-6324, 575-8630 (fax) http://definetti.uark.edu/~gpetris/

...ots, but I'm not able to find it now... Any suggestions...? Thanks in advance. (Of course, I might have dreamt of such a function...) Best, Giovanni -- __________________________________________________ [ ] [ Giovanni Petris GPetris at uark.edu ] [ Department of Mathematical Sciences ] [ University of Arkansas - Fayetteville, AR 72701 ] [ Ph: (479) 575-6324, 575-8630 (fax) ] [ http://definetti.uark.edu/~gpetris/ ] [__________________________________________________]

...square matrix. Let S be my variance matrix, and S = U D V' be its SVD. In my numerical experiments I always got U = V. Is this necessarily the case? Or I might eventually run into a SVD which has U != V? Thank you in advance for your insights and pointers. Giovanni -- Giovanni Petris <GPetris at uark.edu> Associate Professor Department of Mathematical Sciences University of Arkansas - Fayetteville, AR 72701 Ph: (479) 575-6324, 575-8630 (fax) http://definetti.uark.edu/~gpetris/

...e approach does not work, not either without the label... \begin{Rcode} <<>>= y <- c(53, 57, 66, 67, 72) # count data n <- length(y) s <- sum(y) \label{Rcode:sum} @ \end{Rcode} \end{document} Thank you in advance for your suggestions! Giovanni -- Giovanni Petris <GPetris at uark.edu> Associate Professor Department of Mathematical Sciences University of Arkansas - Fayetteville, AR 72701 Ph: (479) 575-6324, 575-8630 (fax) http://definetti.uark.edu/~gpetris/