search for: calcualtions

Dear list members: I searched the R-help for packages to calculate power for an ANCOVA problem I have. I have found power.t.test, power.anova.test. But it seems that I can not found one with ANCOVA. I have two datasetsets with variables: univariate response(one data with continous response and one with 0 and 1), treatment(two levels), covariates(x1,x2,x3). I would appreciate your help. Tony

Does anyone have a good method for calculating Jaccard coefficients now that the dissimilarity() function is no longer an option? Wen Gu John Jay College of Criminal Justice445 West 59 StreetNew York, NY 10029 wgu@gc.cuny.edu _________________________________________________________________ Express your personality in color! Preview and select themes for Hotmail®.

I have orders, order_items, and products. I want to collate several orders so that I can get a SUM of quantities ordered for each product etc. Can I say something like (the below gives an error on :sum, and ignores :group) OrderItem.find(:all, :sum => ''quantity'', :group => ''product_id'', :include => [:order, :product]) I want to get back a

Dear R helpers   I am calculating the 'Yield to Maturity' for the Bond with following characteristics.   Its a $1000 face value, 3 year bond with 10% annual coupon and is priced at 101. The yield to maturity can be calculated after solving the equation -   1010 = [100 / (1+ytm)]  + [100 / (1+ytm)^2] + [ 1100 / (1 + ytm)^3]   This can be solved by trial and error method s.t. ytm = 9.601%.

Dear All, I have a matrix of size 10000 x 50. I would like to calculate all possible pair-wise correlation coefficient (5x10^7 combinations) using cor(). How can I efficiently calcualte and save the result in a matrix? Thanks in advance. Kind regards, Ezhil

...ector and paste together the object name and pass it on to a function. Here an (hopefully) instructive example #Data Example gnuff<-list() gnuff$IHD$LE<-66 gnuff$LUNG$LE <-55 #This is the list, where I collect data for different diseases at the second level of the list #Now I want to do calcualtions just for these two diseases and the sub-list "LE" within these diseases nam <- c("LUNG","IHD") for(i in 1:2) x[i] <- paste("gnuff",nam[i],"LE",sep="$") /2 x #So I try to paste the name of the object which I mean (gnuff$IHD$LEand...

Dear All, In the N-dimensional space, give a data point A and a curve f, how to write the explicit expression for calculating the minimal distance between A and f? Or have to use some nonlinear optimization method to calcualte it? Thanks for your point. Fred [[alternative HTML version deleted]]

I have the following FORTRAN code converted to a DLL: ! my_xmean.f90 ! ! FUNCTIONS/SUBROUTINES exported from my_function.dll: ! my_function - subroutine ! subroutine my_xmean(X,N,XMEAN) ! Expose subroutine my_function to users of this DLL ! !DEC$ ATTRIBUTES DLLEXPORT,C,REFERENCE,ALIAS:'my_xmean_'::my_xmean ! Body of my_function DOUBLE PRECISION X(N) XMEAN=0D0 DO J=1,N

#using "rollapply" to calculate a moving sum or running sum? #I am tryign to use rollapply to calcualte a moving sum? #I tried rollapply and get the error message #"Error in seq.default(start.at, NROW(data), by = by) : # wrong sign in 'by' argument" #example: mymatrix <- ( matrix(data=1:100, nrow=5, ncol=20) ) mymatrix_cumsum <- ( matrix(data=NA, nrow=5,

Hello, I have a series of intraday (high-frequency) price data in the form of POSIX timestamp followed by the value. I sucesfuly loaded that into "its" package object. I would like to create from it a regularly spaced time series of prices (for example 1min, 5min, etc apart) so i could calcualte returns. There is an interpolation function locf() that for timestamp with value NA uses last

Hi,I try to calcualte AIC or Loglik to GARCH model,But the Packege Tseries do not deal with them.How can I calculate AIC or Loglike to GARCH Model By Packege Tseries? Thanks. ____________________________________________________ Free Internet Access NOW! In Alexandria, Ismaileya, Suez, Portsaid, Hurgadha, Sharm Banha, Shebin El-Kom, Damietta, Tanta, Zagazig, Mansoura, Damanhour, Assyout, Qena

Question 1: I have a long list of variable names such as first <- c("one","two","three") and what I want to do is create a list of lists ... where the names of each of overall lists components are "one","two", and "three". This is the same result as second <- list(one=list(),two=list(),three=list()) Is there a way to exploit

Hi I would like to use the newton raphson method to find the root if the equation x^3-0.165*x+0.0003993 without using any readliy available program in r but instead by writing my own code and loop. the problem is that i really cant understand how to write the loop so that it keeps using the last calcualted values. if anyone could help me or give me some tips i would deeply appriciate it thanks

Hi, i'm beginner in Bayesian methods, I'm reading the documentation about dlm package and kalman filters, I'm looking for a example of transformation of ARIMA in a state space equivalent to use the dlm package and calcualte the hyperparameters. Someone can help me about it?. If it's possible with a arima(1,0,1) example, or more complex model. While I have more examples best for me.

Hello, I am a complete newbie with Ruby and Rails. I am working on learning it by writing a project I was going to write in PHP using Rails. I am having a difficult time finding clear information on how I would do this. I need to take a date out of a table in my database and given another day calculate the next 28 day anniversery of the original date. This is used to calculate a delivery

I am trying to perform Singular Value Decomposition (SVD) on a Term Document Matrix I created using the 'tm' package. Eventually I want to do a Latent Semantic Analysis (LSA). There are 5677 documents with 771 terms (the DTM is 771 x 5677). When I try to do the SVD, it runs out of memory. I am using a 12GB Dual core Machine with Windows XP and don't think I can increase the memory

Any suggestions on applications that involve alot of calcualtions on a fairly large data set? My app uses a set of raw data ~5k, applies some default/override rules on the raw data and does some calculations on the data in combination with a list of assumptions. A ranked list along with detailed metrics is generated. The end user can manipulate some of the...

Hello R Community, I've been run the following codes. However, I've been getting an unusual segfault that I'm unable to trace its origin. Please give me a light to decipher the "caught segfault" Thanks for you attention. Bernardo. > options(STERM='iESS', editor='emacsclient') > rm(list = ls()) > > source("fgenIGLD.R") #RNG

...ation a number of times for a set of different initial conditions. If I run them in a batch, say 100 initial conditions, the whole procedure tends to run slower with each new additional condition. I must tell you that not much of information is stored, just pooled results over 1000 cycles and final calcualtions for each initial set. I avoided FOR loops as musch as I could... I have just two: one for the cycles and one for the initial conditions, everything else is vectorized. What could I do to improve the procedure? I am thinking to additionaly vectorize the last for loop (the initial conditions)... wo...

Hello everyone, does anybody know how to increase the decimal places that R uses to calculate something. I think that in default R uses 6 decimal places but I need 12. Thank you very much! Best regards, Maja! --