search for: calcualt

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®.

...t_id'', :include => [:order, :product]) I want to get back a collection of order_items, GROUPed BY product_id. It is important for me to get back order_items, as then I can get to the product using the belongs_to association to get the name of the product, plus the business rules to calcualte the price of the product. Thanks Joerg -- Posted via http://www.ruby-forum.com/.

...0% 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%.   My query is (1) if there is any R package which will calcualte ytm or (2) is there any method in 'R' which can solve the above equation.   Thanking you all in advance   Regards   Madhavi       Your Mail works best with the New Yahoo Optimized IE8. Get it NOW! http://downloads.yahoo.com/in/internetexplorer/ [[alternative HTML version deleted]]

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$LE...

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]]

...s 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 XMEAN=XMEAN+X(J) END DO XMEAN=XMEAN/N RETURN end subroutine my_xmean When I call this DLL from R, it gets loaded properly but the values of XMEAN calcualted are way off: x <- 1:6 > .Fortran("my_xmean",as.double(X),as.integer(length(X)),double(1)) [[1]] [1] 1 2 3 4 5 6 [[2]] [1] 6 $xmean [1] 5.336073e-315 Can someone please let me what is causing this huge difference??? Thank you. Ravi -- View this message in context: http:/...

#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, ncol=20) ) w=12 f...

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 known observation. I guess the idea would be to start from the begining of my series and, for example, if there is no timestamp for t+5min add that time with value NA. Than I could use locf() function to...

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, S...

...ot;two", and "three". This is the same result as second <- list(one=list(),two=list(),three=list()) Is there a way to exploit lapply to convert a list such as first to a list of empty lists such as second? Question 2: In a parallel routine using foreach and doMPI have to put calcualted values into a list of lists such as list(one=one, two=two, three=three) It uses the same (large) list of variable names as first in question 1. Is there a simpler way to accomplish this with a lapply statement? I've tried several permutations but no luck! Thanks a bunch. Jim...

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 jesper [[alternative HTML version deleted]]

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. Thanks all [[alternative HTML version deleted]]

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

...st" ... ..$ Docs : chr [1:5677] "1" "2" "3" "4" ... - attr(*, "class")= chr [1:2] "TermDocumentMatrix" "simple_triplet_matrix" - attr(*, "Weighting")= chr [1:2] "term frequency" "tf" SVD is calcualted using: > tdm_matrix <- as.matrix(tdm2) > svd_out<-svd(tdm_matrix) Error: cannot allocate vector of size 767.7 Mb In addition: Warning messages: 1: In matrix(0, n, np) : Reached total allocation of 3583Mb: see help(memory.size) 2: In matrix(0, n, np) : Reached total allocation of...

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...

...uot; Thanks for you attention. Bernardo. > options(STERM='iESS', editor='emacsclient') > rm(list = ls()) > > source("fgenIGLD.R") #RNG of IGLD(mu,b,sigma) > > source("fsYUTH.R") #the fucntions Y,U,Tbbarran and H necessary to calcualte the hat{b} > > source("fMLEPFP.R") #the maximun likelihood profile estimation of IGLD(mu,sigma) > > source("f04P.R") #it to compute the mle, jacknife and bootstrap > > source("fprofileloglik06.R") #the log-likelihood profile of IGLD(...

...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).....

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! --