similar to: R killed on Solaris

Hi, are there any function to plot the empirical probability distribution function? I just don't want to reinvent the wheel... Best wishes, Ott -- Ott Toomet PhD Student Dept. of Economics ?rhus University Building 322 Universitetsparken 8000 ?rhus C Denmark otoomet (a) econ au dk ph: (+45) 89 42 20 27 ------------------------------------------- (o_ (*_ (O_

Hi folks, I am planning to write some more time-consuming matrix manipulations in c++. What is the experience with the existing c++ matrix libraries? Do you have some recommendations? Are some libraries more compatible with R than the others? All suggestions welcome! Best, Ott

Can anyone explain me the following behavior: > 1:2/1 [1] 1 2 -- makes sense > 1:2/matrix(1,1,1) [1] 1 2 -- makes sense > 1:2/data.frame(a=1) a 1 1 -- why is this different? Best, Ott -- Ott Toomet Visiting Researcher School of Information Mary Gates Hall, Suite 095 University of Washington Seattle, WA 98195 [[alternative HTML version deleted]]

Hi R lovers! I would like to know how you can write comments inside the code of a function is it latex style % ? or any other language style? thanks a lot ************************************************************************* Ce message et toutes les pieces jointes (ci-apres le "message") sont confidentiels et etablis a l'intention exclusive de ses destinataires. Toute

Do I mess up something or is this a bug? If I define an S4 object that contains "factor", all the tests indicate that it is a factor but as.character.factor() complains of it being a non-factor... > setClass("Foo", contains="factor") > a <- new("Foo", factor(1:3)) > a Object of class "Foo" [1] 1 2 3 Levels: 1 2 3 > class(a) [1]

Hi, I want to calculate expected likelihood over a 2D discrete distribution, something like \sum_k \sum_l p_{kl} L(v_k, v_l) It is trivial to write a cycle like for(k in 1:K) for(l in 1:L) sum <- sum + p[k,l]*L(v[k], v[l]) But is there a more clever way for R? Best wishes, Ott

Dear R-people, are there any way to change the size of legend in barplot? I have tried various versions of cex, both as par(cex.*= ) and barplot(..., cex.*= ). So long without success. Sincerely, Ott > version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1

hist() is designed so that the total area sums to 1. You should build you desired behavior using a barchart. ? David Sent from my iPhone > On Aug 5, 2023, at 11:50 PM, Ott Toomet <otoomet at gmail.com> wrote: > > ?Sorry if this topic has been discussed earlier. > > Currently, hist(..., log="y") fails with > >> hist(rexp(1000, 1), log="y")

Hi, I have a following vector: > smallch [1] 0.0652840 0.1181300 0.0319370 0.0155700 0.0464110 0.0107850 [7] 0.0158970 0.0375900 0.0603090 0.0310300 0.0105920 0.0540580 [13] -0.0177740 0.0039393 Pretty (R 1.5.1) has problems with zero: > pretty(smallch) [1] -2.000000e-02 -3.469447e-18 2.000000e-02 4.000000e-02 6.000000e-02 [6] 8.000000e-02 1.000000e-01 1.200000e-01

Hi, I have a following vector: > smallch [1] 0.0652840 0.1181300 0.0319370 0.0155700 0.0464110 0.0107850 [7] 0.0158970 0.0375900 0.0603090 0.0310300 0.0105920 0.0540580 [13] -0.0177740 0.0039393 Pretty (R 1.5.1) has problems with zero: > pretty(smallch) [1] -2.000000e-02 -3.469447e-18 2.000000e-02 4.000000e-02 6.000000e-02 [6] 8.000000e-02 1.000000e-01 1.200000e-01

Hi, I have a problem with spacing in XML files when reading them with xmlTreeParse. I don't know the exact specification of xml but according what I have red before it should work. consider a tiny test.xml file: <?xml version="1.0"?> <fields> <v1>1 </v1> <v2> 2 </v2> <v3> 3</v3> </fields> i.e. I have three fields v1, v2 and

Full_Name: Ott Toomet Version: 1.5.1 OS: Linux (RH 7.1) Submission from: (NULL) (192.38.43.138) Hi, Try to give an impossible file name (command to pipe through) for dev.print: > plot(1:4) > dev.print(file="|llllll") I got the response sh: llllll: command not found Process R broken pipe at Wed Aug 28 14:42:27 2002 After which R was killed (or at least vanished from the

Hi everybody, I am using gzipped files quite extensively and cannot figure out how I can close (or reuse) connections. Example: > for(i in 1:300) {cat(i, "\n"); a <- readLines(gzfile("file.gz"), n=1) } 1 2 ... 46 47 Error in gzfile("~/tyyq/andmebaasiq/ETU/ETU1997.tsv.gz") : all connections are in use To use just gzfile without any open/close should

Dear R-people, I have to compute C - -(pnorm(B)*dnorm(B)*B + dnorm(B)^2)/pnorm(B)^2 This expression seems to be converging to -1 if B approaches to -Inf (although I am unable to prove it). R has no problems until B equals around -28 or less, where both numerator and denominator go to 0 and you get NaN. A simple workaround I did was C <- ifelse(B > -25, -(pnorm(B)*dnorm(B)*B

Sorry if this topic has been discussed earlier. Currently, hist(..., log="y") fails with > hist(rexp(1000, 1), log="y") Warning messages: 1: In plot.window(xlim, ylim, "", ...) : nonfinite axis=2 limits [GScale(-inf,2.59218,..); log=TRUE] -- corrected now 2: In title(main = main, sub = sub, xlab = xlab, ylab = ylab, ...) : "log" is not a graphical

I have installed R 1.6 for Linux from source, but unfortunately, I didn't check if the Xwindow headers were available, so I have installed it without Xwindows capabilities. Thus, I want to uninstall it and install it again properly with the required features, but I don't know how to uninstall R when built from source. I would appreciate any help on this. Thanks very much in advance!

Hi R-devl, I'm still unable to force opening an *interactive* non-Rstudio platform-specific plot device on *unix* systems. dev.new() add a new argument 'noRStudioGD' in R 3.1.1. Thank you. It works for me when using RStudio on Windows, but on the unix system it opens a pdf device instead of an interactive device when using an interactive RStudio session (with R_DEFAULT_DEVICE

I would like to create a barchart that has vertical instead of horizontal bars. My goal is something similar to Figure 6.10 (p. 166) in "The Basics of S and S-plus," except for the orientation of the bars. I don't think bwplot() or xyplot() will work. Any suggestions? Scott Waichler Hydrology Group/Environmental Technology Division Battelle Pacific Northwest National Laboratory

Hi, I am trying to study R internal behaviour. So long, I have not succeeded to access the value of R strings from C. I use: void salvesta_tabel( SEXP data_frame, SEXP file ) { printf( "nimi %d\n", (R_CHAR)( file)); } and from the R side: salvesta.tabel <- function (x, file = "") { .Call( "salvesta_tabel", x, file) } When calling

Full_Name: Ott Toomet Version: 2.6.0, 2.5.x OS: debian etch, lenny Submission from: (NULL) (80.235.63.243) When building a package with 'R CMD build name_of_directory" using "et_EE.UTF-8" locale, I get the following: siim at tancredi:~/tyyq/econ/micEcon$ R CMD build trunk * checking for file 'trunk/DESCRIPTION' ... OK * preparing 'trunk': * checking