similar to: LDA fuction

I would like to integrate the area under a curve without any smoothing or the like- just on the raw numbers. I looked at integrate() but it requires a function which I assume means something like x+x^2+x^3 is there a built in function in R for this? #let's say x <- seq(1:50) y <- seq(1:50) plot(y~x) # the are would be 1250 # I would like to be able to do this but on more complicated

Hello all, I am struggling with sending voicemail as an attachement in Email. When i have given the email like someone at gmail.com it is delivering to my gamil account perfectly(of course to spam folder). But when i given the email like someone at mycompanymail.com it is not delivering to my company email account.. What should i do ? Actually my company is using a third party email server..

#is there a way to get NA in the table of descriptive statistics instead of the function stopping Thank you in advance #data x.f <- structure(list(Site = structure(c(9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L), .Label = c("BC", "HC", "RM119", "RM148", "RM179", "RM185",

library(zoo) library(chron) t1 <- chron("1/1/2006", "00:00:00") t2 <- chron("1/31/2006", "23:45:00") deltat <- times("00:15:00") tt <- seq(t1, t2, by = times("00:15:00")) d <- sample(33:700, 2976, replace=TRUE) sin.zoo <- zoo(d,tt) #there are ninety six reading in a day d.max <- rollapply(sin.zoo, width=96, FUN=max)

Dear users I am fitting a Generalized Additive Mixed Models (gamm) model to establish possible relationship between explanatory variables (water temperature, dissolved oxygen and chlorophyll) and zooplankton data collected in the inner and outer estuarine waters. I am using monthly time-series which are auto-correlated. In the case of the inner waters, I have applied satisfactoryly (by

The following C code : #include <stdio.h> #include <stdlib.h> int TESTE2( int parami , int paraml ,double paramd ) { int varx=0,vary; int nI =0; //varx= parami; if( parami > 0 ) { varx = parami; vary = varx + 1; } else { varx = vary + 1; vary = paraml; } varx = varx + parami + paraml; for( nI = 1 ; nI <= paraml; nI++) { varx =

hi all why does the following not work??? this was someone elses code and i couldnt explain why it doesn't work. m=matrix(c(0,0),2,1) v=matrix(c(1,0,0,1),2,2) Y=function(X1,X2,mu=m,V=v) { X=matrix(c(X1,X2),2,1) a=(1/((2*pi)*sqrt(det(V))))*exp((-0.5)*(t(X-mu)%*%solve(V)%*%(X-mu))) a[1] } x1=seq(-1,1) x2=x1 Z=outer(x1,x2,FUN="Y",mu=m,V=v) persp(x1,x2,Z) my code:

I'm not able to create the proper syntax to specify a lattice bwplot() for only one of two conditioning factors. The syntax that produces a box plot of each of the two conditioning factors is: bwplot(quant ~ param | era, data=mg.d, main='Dissolved Magnesium', ylab='Concentration (mg/L)') What I've tried unsuccessfully are: bwplot(quant ~ param |

The following C test program was compiled using LLVM with -O3 option and MSVC with /O2. The MSVC one is about 600 times faster than the one compiled with the LLVM. We can see that the for loop in MSVC assembler is solved in the optimization pass more efficiently than that in LLVM. Is there an way to get a optimization result in LLVM like that of the MSVC? Manoel Teixeira #include

Dear all, Being a newbie to R, I've trawled through many old posts on this list looking for a solution to my problem, but unfortunately couldn't quite figure it out myself. I'd be very grateful if someone here on this list could perhaps help me out. I have a lattice plot with several panels and would like to add some text next to the y-axis on the right hand side of each row of

Hola a todos, Partiendo de: library(ggplot2); library(data.table) datos <- data.table(Date = seq.Date(as.Date("2017-01-01"), as.Date("2017-01-10"), by = "day"), V1 = rnorm(10), V2 = rnorm(10)) Quiero crear una función para pintar un gráfico cuyos argumentos sean

I don't find how to do what I need to do in Dalgaard or 'R Cookbook', so I'm asking here. I have a data frame with water chemistry data and I want to start exploring these data. There are three factors (site, date, chemical) associated with each measurement. The data frame looks like this: > summary(chemdata) site_id.sample_date.param.quant

R users, I have df like this a <- data.frame( indx = 1:20, var1 = rep(c("I20", "I40", "A50", "B60"), each=5), var1_lab= rep(c("cat", "dog", "mouse", "horse"), each=5), var2 = rep(c("B20", "X40", "D50",

Hi Everyone, I want to generate some random numbers according to some empirical distribution. Therefore I am looking for the inverse of an empirical cumulative distribution function. I haven't found any in R. Can anyone give a pointer? Thanks, Geoffrey _______________________________________________________=0A= =0A= =0A= The information in this email or in any file attached

I'm having a function of the form 1> f<-function(x){ 1+ 1+ return(x^p) 1+ 1+ } ,and I would like to integrate it with respect to x, where p should be any constant. One way would be to set a value for p globally and then call integrate function: p=2 integrate(f, lower = -1, upper = 1) However, I would like to use 'integrate' inside a function, so I could call it passing

Hi, I am having an issue when creating a postscript file from RGL window. It seems to cut off some of the axis labels. Here is the code I am using. I created a 3D plot using RGL_0.71 with R 2.5 on Windows XP. z1<-c(5,4,1,4.5,2,3,2,1,1) z2<-c(6,8,7,7.5,5,3.5,4,1,1) z3<-c(3,2,4,7,3,4.5,6,2,3) x1<-seq(1,9) x2<-seq(1,9) x3<-seq(10,18) y1<-seq(8,0) y2<--1*y1

I add a function abs64 to support abs for 64 bit integer -- 此致 礼 罗勇刚 Yours sincerely, Yonggang Luo -------------- next part -------------- A non-text attachment was scrubbed... Name: llvmabsfault.patch Type: application/octet-stream Size: 800 bytes Desc: not available URL: <http://lists.llvm.org/pipermail/llvm-dev/attachments/20090809/0d1aa27d/attachment.obj>

hi, i just ported the patch of David Hinkle http://www.mail-archive.com/asterisk-users@lists.digium.com/msg49781.html to the current cvs version of asterisk. the theory is that the leds of supervisioned extensions are blinking until a call is established whereafter the leds should be constantly lit. however it's not working. the asterisk server is sending the following xml notify to the

On 06/08/18 13:48, m.roth at 5-cent.us wrote: > Frank Cox wrote: >>>> so if it would work, replace shortname with short and short1? >> >> With all of this hokey-pokey surrounding licensing and mac addresses, I >> wonder if this outfit is actually still in compliance with the terms of >> their license for this software, whatever it may be? >> >> If

Hi everyone, I'm working on a script trying to use foreach %dopar% but without success, so I manage to run the code with foreach %do% and looks like this: The code is part of a MCMC model for projects valuation, returning the most important results (VPN, TIR, EVA, etc.) of the simulation. foreach (simx = NsimT, .combine=cbind, .inorder=FALSE, .verbose=TRUE) %do% { MCPVMPA = MCVAMPA[simx]