similar to: Compact sums in functions definitions

In order to define a function f as: > f <- function(x) (x+1)*(x+2) I want to use the notation: > v = c(1,2) > g <- function(x) prod((v+x)) That apparently works and, for instance, the loop: >for (i in 1:100) { print(f(i)-g(i)) } Produces a sequence of zeros. Nevertheless, if I try to plot the function g by: >t = seq(0,100,1) >plot(t,g(t),type="l") I

Reading a csv file (db.csv) as: yes,full no,full no,empty I can use the command 'levels' to extract the different values appearing for each column as follows: > d<-read.csv("db.csv") > levels(d[,2]) which returns > [1] "empty" "full" while, doing the same with a numerical csv file as: 1,6 0,6 0,7 the same instruction returns an empty

Estimados usarios de R: Tengo una base de datos madre en formato .sav de SPSS y la quiero modificar usando datos de otras base de datos .sav y otra en .csv a las que llamare hijos. No tengo problema en convertirlas en data.frame. Todos los archivos tienen en común una variable que es única. En aquellos casos que les falta un valor a una variable en el archivo madre lo relleno con el valor del

Dear All, I carry out negative binomial estimations using the glm.nb command from the MASS package. Is there a command or a simple procedure for computing marginal effects from a glm.nb fitted object? If these are the same as for a Poisson fitted object (glm), my question remains how to compute them. Thanks in advance for your help. Roberto Patuelli ******************** Roberto Patuelli, Ph.D.

Hi Hadley and everyone, here's a patch for ggplot2 that fixes the behavior of opts(legend.position={left,top,bottom}). If you try the following code in an unmodified ggplot2 options(warn = -1) suppressPackageStartupMessages(library("ggplot2")) data <- data.frame( x = c(1, 2, 3, 4, 5, 6), y = c(2, 3, 4, 3, 4, 5), colour = c(TRUE, TRUE, TRUE, FALSE, FALSE, FALSE))

Hello everybody My collegue and I noticed a strange behaviour of R on different platforms. It's a simple computation, but results are rather different. On Windows XP: > floor(log(8,2)) [1] 3 which is what one should expect. Here's instead the result with Mac OS X (same version, 2.5.0 (2007-04-23)) > floor(log(8,2)) [1] 2 Is it a "bug" in R or in the operating

Dear All, I would like to analyse the residuals from a generalized linear mixed model (GLMM) that I estimated, with random effects, by means of the command glmmPQL, from the MASS package. It is not very clear to me what the actual residuals to analyse are (Y - Yhat): I obtain two columns of residuals, of which the first are population residuals, and the second refer to the grouping used in the

Dear Duncan, Thank you a lot for your feedback. I have a problem though. The branch weights counters overflow in some files and thus I get incorrect numbers. Is there any way to find a workaround for that? Is is supposed to be a known bug or is it something that needs configuration on my part? Again, thank you a lot for your reply. Best Regards, Georgios Zacharopoulos 2016-02-03 18:23

Dear All, I've been using already for a year glm.nb() from the MASS package. But today, R gave me an error message when estimating one of my usual models: > depEsf.nb <- glm.nb(depE ~ manuf00E + corps00E + lngdp00E + lngdp00sqE + > lnpop00E + indshE + scishE + mechshE + elecshE + chemshE + drugshE + > urban_dummyE + aggl_dummyE + + eE1 + eE2 + eE3 + eE4 + eE5 + eE6 + eE7 +

What is the recommended/optimal way to perform aggregates on data frames with chron objects? Here is an example: >raw.data 1 07/09/01 4000 2 07/09/01 2000 3 07/11/01 1000 4 07/13/01 800 5 07/13/01 700 6 07/16/01 600 7 07/17/01 500 I'm trying to construct a function that would first aggregate the data (second column) by day (grouping by the first column) according to a

I am currently trying to use R to construct a regression model to explain output based on temperature. I have combined my output and temp data into a notepad file. there is no problem with loading the data into R. > data.df output temp 1 850 17 2 849 17 3 905 17 4 925 17 5 1043 19 6 1104 20 7 1097 18 8 979 19 9 926 18 10 1133 18 ~~ ~~ 240 1124

Dear Duncan, I am generating branch-weights annotated IR files as described in the documentation of LLVM, using profiling with instrumentation. http://clang.llvm.org/docs/UsersManual.html#profiling-with-instrumentation e.g. llvm-profdata merge -output=$(BENCH).profdata default.profraw > clang -S -emit-llvm -O3 -fprofile-instr-use=$(BENCH).profdata -o > bench.prof.ll bench.c The issue is

Dear ALL, I'm trying to figure out what the percentage effects are in a logistic regression. To be more clear, I'm not interested in the effect on y of a 1-unit increase in x, but on the percentage effect on y of a 1% increase in x (in economics this is also often called an "elasticity"). For example, if my independent variables are in logs, the betas can be directly

Hi all: following the previous discussion, it looks like plot(x) with y=NULL still does not work correctly. If one tries for example plot(1:5) it works, but already for plot(runif(100)) it does not. I posted the proposed correction for plot.POSIXct and plot.POSIXlt before. Please voice your opinions whether the following fix for plot.default could be reasonable? I include the full function and

Dear All, I'm looking to perform an ANOVA between two nested panel fixed effects models. I tried with anova, as well as with waldtest. anova tells me there is no method available for plm objects, while waldtest tells me my models are not nested. I think they are instead. The difference between the two models is that in the second I let the regression coefficients of a given variable variate

Dear Person That Posted This Email, Please be aware that calls for papers are only on-topic for this email list if they are specifically about LLVM. For example, a call for a workshop on LLVM in Scientific Computing would be on topic whereas a general compiler call for papers (like the one below) is not. That said, I hope you receive many strong submissions. Your conference looks

------------------------------------------------------------------------------------------------ CALL FOR PAPERS 15th International Conference on Managed Languages & Runtimes (ManLang'18) September 10-14, 2018, Linz, Austria http://ssw.jku.at/manlang18/ ------------------------------------------------------------------------------------------------ ManLang (formerly PPPJ) is a

Hi! How to find out which plot function is used when i call plot(hist.default(1:10,plot=F)) and all works fine ? The reason why I would like to know it is that after loading some self written R functions > plot(hist.default(1:10,plot=F)) Error in xy.coords(x, y, xlabel, ylabel, log) : x and y lengths differ > traceback() 5: stop("x and y lengths differ") 4: xy.coords(x, y,

Call for Participation http://www.hpdc.org/2012/ The organizing committee is delighted to invite you to *HPDC'12*, the /21st International ACM Symposium on High-Performance Parallel and Distributed Computing/, to be held in *Delft, the Netherlands*, which is a historic, picturesque city that is less than one hour away from Amsterdam-Schiphol airport. HPDC <http://www.hpdc.org> is

Call for Participation http://www.hpdc.org/2012/ The organizing committee is delighted to invite you to *HPDC'12*, the /21st International ACM Symposium on High-Performance Parallel and Distributed Computing/, to be held in *Delft, the Netherlands*, which is a historic, picturesque city that is less than one hour away from Amsterdam-Schiphol airport. HPDC <http://www.hpdc.org> is