similar to: S speedup in R

Greetings, I am in the process of moving several large sites to a Tyan+BQ+CentOS box (built from the Nuonce ISO) and I see that special characters are not showing up in the browser properly but do show up in the source "view source" properly. -- very strange. Example: on the web page --> Heizungsbau M?ller in the source --> Heizungsbau M?ller When I view this on the existing

Hi, I am running simulations that does multiple comparisons to control. For each simulation, I need to model 7 nls functions. I loop over 7 to do the nls using try if try fails, I break out of that loop, and go to next simulation. I get warnings on nls failures, but the simulation continues to run, except when the internal call (internal to nls) of the chol2inv fails.

ggplot2 =================================== ggplot2 is a plotting system for R, based on the grammar of graphics, which tries to take the good parts of base and lattice graphics and avoid bad parts. It takes care of many of the fiddly details that make plotting a hassle (like drawing legends) as well as providing a powerful model of graphics that makes it easy to produce complex multi-layered

ggplot2 =================================== ggplot2 is a plotting system for R, based on the grammar of graphics, which tries to take the good parts of base and lattice graphics and avoid bad parts. It takes care of many of the fiddly details that make plotting a hassle (like drawing legends) as well as providing a powerful model of graphics that makes it easy to produce complex multi-layered

Hi everybody, I am reading the Lomb paper (Lomb, 1976) and I found an interesting equation, and I wish to resolve it using R. I am wondering if anybody has a hint. The equation is: Yis[i] = a*cos((2*pi/T)*(Times[i] - Tau)) + b*sin((2*pi/T)*(Times[i] - Tau)) ... (1) Where T and Tau are constants. I know the "Times" and "Tis" values (in fact these values come from a Time

Hello, Please copy all replies directly to me (my account is having difficulty with receiving lists these days). This is primarily a programming question, but the specifics regard start values for a nonlinear regression (if you have suggestions on alternative ways to obtain start values they are welcome as well). I'm using nls to estimate a nonlinear time series equation of the form:

Full_Name: Edward Huang Version: 2.6.1 OS: Mac OS 10.5.1 Leopard Submission from: (NULL) (71.198.106.232) I'm trying to run quantile regression on my data. I just couldn't make it work. The same dataset ran okay on STATA 10, tho. Would you please take a look at it? Here is the error message: *** caught segfault *** address 0x3ff00008, cause 'memory not mapped' Traceback:

Dear all, I am trying to compare the estimated coefficients of a quantile regression model between two different samples. It is a Wald test, but I cannot find one way to do that in R.The samples are collected conditional on a specific characteristic and I would like to test whether such characteristic indeed affect the estimators. The problem in the test anova.rq is that the response variable

Hello. I am a beginer in R and I can not implement the KhmaladzeTest in the following command. Please help me!!!!!!!!!!! PD: I attach thw results and the messages of the R program R : Copyright 2006, The R Foundation for Statistical Computing Version 2.3.1 (2006-06-01) ISBN 3-900051-07-0 R es un software libre y viene sin GARANTIA ALGUNA. Usted puede redistribuirlo bajo ciertas

Hi, I'd like to make a legend with a mix of mathematical symbol (tau), numeric variable and character variables.I have tried : types<-c("Type 1","Type 2","Type 2") tau<-c(1,3,2) legend(x="topright",legend=paste(types,"tau=",expression(tau))) but it doesn't work: the 'tau' symbol is not written in its 'symbol

Hi all, I would like to ask a question related to Kendall package. I ran Kendall (x,y) and saw the results. But I am not sure which tau values R reported. I have ties in my data set, so I want tau-b. Can anybody tell how Kendall package is calculating tau values? I have looked at the package PDF, but I could not find any useful information. As long as I see from the following link, there

Hello there Can anybody please tell me if there is any package in R to solve the following 4 nonlinear equations with 4 unknowns: alpha*exp(20/sigma)+ beta*exp(21/tau) = 2 alpha*exp(22/sigma)+ beta*exp(9/tau) = 4 alpha*exp(10/sigma)+ beta*exp(30/tau) = 6 alpha*exp(40/sigma)+ beta*exp(39/tau) = 5 where alpha = exp(lambda/sigma) beta= exp(delta/tau) I need to estimate lambda, sigma, delta, tau

Dear Rxperts, Running the following code: ======================================================= twlo=10; twhi=20; wt=154; vd=0.5; cl=0.046; tau=6; t=3; F=1; wtkg <- wt/2.2 # convert lbs to kg vd.pt <- wtkg * vd # compute weight-based vd (L) cl.pt <- wtkg * cl # compute CL (L/hr) k <- cl.pt/vd.pt # compute k (hr^-1) cave <-

Hello, I would like to use mapply to avoid using a loop but for some reason, I can't seem to get it to work. I've included copies of my code below. The first set of code uses a loop (and it works fine), and the second set of code attempts to use mapply but I get a "subscript out of bounds" error. Any guidance would be greatly appreciated. Xj, Yj, and Wj are also lists, and s2,

Hi all, I am using quantile regression with svy design. I want to retrieve summary regression statistics (std error, p-value), since I don't have any in my output: Commands: clus1_d<- svydesign(id=~cd002_co, weights=~wtper, strata=~str, data=data) bclus1<-as.svrepdesign(clus1_d,type="bootstrap",replicates=100) fit1<-

Hi I want to write a C function for the R Code below and call it with .Call: SimPoisson <- function(lambda,tgrid,T2M) #Simulation eines Poissonprozesses { NT <- 0 Ni <- rep(0,length(tgrid)) tau <- 0 sign <- 0 if(lambda != 0) { i=1 j=1 while (1) { EVar <- rexp(1,lambda) sign <- sign + EVar if (sign > T2M) { break } tau[i] <- sign

I would like to use greek "tau" as a symbol of variable to integrate over in plotmath expression(integral(f(tau)*dtau, 0,t)) but nothing seems to work. I tried d{\tau}, d\tau, etc., without any success Is it possible? How can I accomplish this? Best regards, Ryszard -------------------------------------------------------------------------- Confidentiality Notice: This message is

(BlueQuartz+CentOS4 server) I installed an SSL cert on my server this week. SSL seems to work fine in the site except for scripts. :( this works... http://secure3.pdcweb.net https://secure3.pdcweb.net And this script works... http://secure3.pdcweb.net/test.cgi But this does not work... https://secure3.pdcweb.net/test.cgi ...It just displays the code :( Any ideas? - Bill

Dear expeRts, I would like to create a scatter plot matrix with splom(). The lower panel should contain some additional information about the samples shown in the upper panel plot, see the splom() call below. Now two questions came up: (1) The lower panels show "tau" and "alpha" on top of each other. How can I plot *three* expressions on top of each other? I tried several

Dear R family I have a following question. Suppose I have a matrix as follows, for instance: tau= 0 1 0 0 0 0 0 1 0 0 0 0 0 1 0 0 0 0 0 1 1 0 0 0 0 I want to calculate (-m) power of tau, for example, m=893. When I run tau^2, the outcome is just tau. Any help would be appreciated. thanks in advance. best moohwan