similar to: Probability of data values form empirical distribution

Hello,   Could someone help me please and to tell how to get the probability from empirical DENSITY (not parametric) for each data value (R function). For example, for normal distribution there is such a function like:   “dnorm(q, mean = 0, sd = 1, lower.tail = TRUE, log.p = FALSE)”   I need the same function only for the empirical DENSITY function (which does not correspond to any typical

Hello,   What function can I use for matrices addition? I couldn’t find any information about it in the manual or in the internet. (A+B suits, when the number of matrixes is small, function sum() doesn’t suit for matrices addition, because it sums all variables in the matrices and produces as an answer single number, not a matrix).   Best regards, Lina [[alternative HTML version

Dear Sirs, I would like to ask you, what function can I use for matrices addition? I couldn't find any information about it in the manual or in the internet. (A+B suits, when the number of matrixes is small, function sum() doesn't suit for matrices addition, because it sums all variables in the matrices and produces as an answer single number, not a matrix). I would be very thankful for

I ran into a problem with titles on graphs. I wanted a graph with multiple subplots, with each having a title that involved both a Greek letter and an identifier for each graph. Below is a simplified version of code to do this. The graph appears fine, with the first graph having "i=1" in the title, and the second graph having "i=2" in the title. However, when I resize the

I am trying to call a precompiled C function that uses a struct as one of it's arguments. I could write a wrapper function in C, but I was hoping there is some way to pack fields into an array of type raw that could be passed directly to the function. Here is some more detail. The C struct is simple, but has mixed types: struct STRUCT1 { long type; long nx; double

There is a contradiction between what the help page says and what constrOptim actually does with the constraints. The issue is what happens on the boundary. The help page says The feasible region is defined by ?ui %*% theta - ci >= 0?, but the R code for constrOptim reads if (any(ui %*% theta - ci <= 0)) stop("initial value not feasible") The following example

Ben Bolker gives some reasons why Numerical Recipes may be problematic as a starting point for R codes. CUP did a masterful job of marketing, but the license is restrictive as the links he gives points out. In some tests, I've also noted that some of the algorithms are less than stellar e.g, convergence tests in one or two optimization routines. Should we have a wiki item to help people find

I have developed an R package that works under Win32, but when I attempt to build it on Win64, I get ERROR: loading failed for 'x64' More precisely, I developed and tested the package under Win32 and it works. But when I move to a 64 bit Windows 7 (Home Premium) system, and attempt to build both 32 bit and 64 bit packages, the 32 bit package seems to build, but the 64 bit build

Hello: The example "integrate(dnorm,0,20000)" says it "fails on many systems". I just got 0 from it, when I should have gotten either an error or something close to 0.5. I got this with R 2.12.0 under both Windows Vista_x64 and Linux (Fedora 13); see the results from Windows below. I thought you might want to know. Thanks for all your work in creating

Dear all, I am trying to use the rstable(n, alpha, beta, gamma = 1, delta = 0, pm = c(0, 1, 2))) function to generate positive stable random numbers. For positive stable distribution, beta==1 and alpha is in (0,1), which defines random variables with support (0, infinity). So, I used rstable(100, 0.5, 1) for an example. I found that this gives me some negative numbers. For example, >

Hi I am new to writing C code and am trying to write an R extension in C. I have hit a wall with F77_CALL(dgemm) in that it produces wrong results. The code below is a simplified example that multiplies the matrices Ab and Bm to give Cm. The results below show clearly that Cm is wrong. Am= 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 Bm= 1 1 1 1 1

Hi, The `quietly` argument of `require` is documented as follows: quietly: a logical. If ?TRUE?, no message confirming package attaching is printed, and most often, no errors/warnings are printed if package attaching fails. However: > require(foo, quietly=TRUE) Warning message: In library(package, lib.loc = lib.loc, character.only = TRUE, logical.return = TRUE, :

Hello, How do I import data from Excel onto R? I am using S-Plus Envstat module. I need to use some of the data that already exist in the Envstat and put it into R to make graphs, find basic informations like mean, median, standard deviation, etc. I have been reading the help links, but I don't see anything in reference to this. Please help me! -Maher Lina

Dear list! I have run into a problem that seems very simple but I can't find any solution to it (have searched the internet, help-files and "An introduction to R" etc without any luck). The problem is the following: I would like to create a data.frame with two components (columns), the first component being a matrix and the second component a vector. Whatever I have tried so far, I

Where does the value 2.2e-16 come from in p-values for chisq tests such as those reported below? > Anova(cm.mod2) Analysis of Deviance Table (Type II tests) Response: Freq LR Chisq Df Pr(>Chisq) B 11026.2 1 < 2.2e-16 *** W 7037.5 1 < 2.2e-16 *** Age 886.6 8 < 2.2e-16 *** B:W 3025.2 1 < 2.2e-16 *** B:Age 1130.4 8 < 2.2e-16 *** W:Age 332.9 8 < 2.2e-16 *** --- Signif.

If we do add an argument to get(), then it should be named consistently with the ifnotfound argument of mget(). As mentioned, the possibility of a NULL value is problematic. One solution is a sentinel value that indicates an unbound value (like R_UnboundValue). But another idea (and one pretty similar to John's) is to follow the SYMSXP design at the C level, where there is a structure that

Hi! Does anyone know how to do the test for goodness of fit of a logistic model (in rms package) after running fit.mult.impute? I am using the rms and Hmisc packages to do a multiple imputation followed by a logistic regression model using lrm. Everything works fine until I try to run the test for goodness of fit: residuals(type=c("gof")) One needs to specify y=T and x=T in the fit. But

Well, I'm getting a warning (not an error) when the package doesn't exist. I interpreted "most often" to mean that suppressing warnings/errors is why you'd most often use this argument, as most packages don't emit startup messages. And technically there isn't a problem with attaching the package, since we don't even try to attach packages that don't exist.

Hello all, this is my first post to the mailing list, so I hope this is the right way to ask a question. I''ve installed xen on my laptop (Asus N51Vf, specs here: http://www.asus.com/product.aspx?P_ID=0Tc42s23bOdbbvY5), as I''m currently at university and don''t have my proper desktop hardware to play with. The laptop supports IOMMU/Vt-d and has one PCI-e Geforce GT130M

Hi, I'm trying to compile a simple program that only printf("hello") using llvm-gcc with -emit-llvm flag, but I'm getting the following error stream. Any ideas what could be the problem? Thanks, Tal In file included from /usr/include/stdio.h:28, from simple.c:1: /usr/include/features.h:324:26: error: bits/predefs.h: No such file or directory