similar to: Using as.polynomial() over a matrix

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 I've written a function that uses Rmpi to perform a calculation in parallel. It works fine, but I'm trying to improve efficiency in terms of memory usage and the amount of data being passed back and forth between mater and slaves. Calculations are performed on a symmetrical matrix in order to zero-out some of the cells. In the parallel

Hi All I am getting the above mentioned error when I run the code below. I don't know why because I have implemented the function and I'm calling the function with a parameter. I'm obviously missing the plot ... Can someone perhaps point out the error of my ways? Error: > out<-ode(y=init, times=times, func=G1999, parms=parms, method="lsoda") Error in m_Na(v) : attempt

void inverse_mdct_slow(float *buffer, int n) { =A0=A0 int i,j; =A0=A0 int n2 =3D n >> 1; =A0=A0 float *x =3D (float *) malloc(sizeof(*x) * n2); =A0=A0 memcpy(x, buffer, sizeof(*x) * n2); =A0=A0 for (i=3D0; i < n; ++i) { =A0=A0=A0=A0=A0 float acc =3D 0; =A0=A0=A0=A0=A0 for (j=3D0; j < n2; ++j) =A0=A0=A0=A0=A0=A0=A0=A0 // formula from paper: =A0=A0=A0=A0=A0=A0=A0=A0 //acc +=3D n/4.0f *

Hi All I downloaded the binary package of Samba 2.2.2 on Friday from one of the Samba mirror sites, I think Belgium, and installed it on a Sun E450 running Solaris 2.6. Today everything has fallen flat on its face. Unfortunately I am on holiday this week so not able to deal with the problem directly, but via telephone calls. Apparently what is happening is that smbd is gobbling up all the

Dear Wiza[R]ds, I have the following data in a data.frame. I need to do an anova with multiple comparison but I don't know how to factor the groups for analysis. There are 3 groups, 1,2 and 3 labelled in column 1. Help appreciated with thanks in advance. Group SI Sif SG Io I2 lol2 1 9.08e-05 9.08e+00 0.060842287 1.798556446 32.574500

A prebuilt binary. Needs non-free Watcom compiler to build. Signed-off-by: Lubomir Rintel <lkundrak at v3.sk> --- dosutil/mdiskchk.com | Bin 8106 -> 0 bytes 1 file changed, 0 insertions(+), 0 deletions(-) delete mode 100644 dosutil/mdiskchk.com diff --git a/dosutil/mdiskchk.com b/dosutil/mdiskchk.com deleted file mode 100644 index

Hi, I have a casting/formatting question on hist.POSIXt: The histogram plot from POSIXct works perfect (with help of Prof. Ripley -thanks!). When processing the hist(plot=FALSE) output and then plotting the results over the x-axis (bins) coming from hist(), I lose the date/time labels, getting instead integers displayed. Trying to cast the $breaks with as.POSIXct gives silly results with

Hi all, I am trying to run Weibull PH model in R. Assume in the data set I have x1 a continuous variable and x2 a categorical variable with two classes (0= sick and 1= healthy). I fit the model in the following way. Test=survreg(Surv(time,cens)~ x1+x2,dist="weibull") My questions are 1. Is it Weibull PH model or Weibull AFT model? Call: survreg(formula = Surv(time, delta) ~ x1

I am having trouble understanding results I'm getting back from coxph doing a recurrent event analysis. I've included the model below and the summary. In some cases, with minor variations, the Robust variance and Wald tests are significant, but the individual covariates may or may not be significant. My main question is: If Wald and robust tests both take into account the

Hi All , Can anybody tell me if there's any way to get the summarized anova values.Now i will explain what i mean , when i say "*summarized*". Below you can see the anova table of recmeanC1 with rest* all* i.e from recmeanC2 to i15(predictors),as shown in table. Df Sum Squares Mean Square F value Significance [Pr(>F)] recmeanC2 1 89.272 89.272

Hello, I have some trouble using step() and stepAIC() functions. I'm predicting recruitment against several factors for different plant species using a negative binomial glm. Sometimes, summary(step(model)) or summary(stepAIC(model) does not select the best model (lowest AIC) but just stops before. For some species, step() works and stepAIC don't and in others, it's the opposite.

Hello I am an R newbie, and am coming against a couple of problems when I try and apply the scale function across all the rows of a matrix. ------------------------- #I have a matrix "dt_l". > str(dt_l) num [1:40962, 1:885] 3.04 4.1 3.4 3.58 3.77 ... #I want to convert the values in each row of this matrix into standardised values (mean 0, sd 1). I believe the command for doing

Hi R community! I've fitted three mixed-effects regression models to a thousand bootstrap samples (case-resampling regression) using the lme4 package in a custom-built for-loop. The only output I saved were the inferential statistics for my fixed and random effects. I did not save any output related to the performance to the machine learning algorithm used to fit the models (REML=FALSE).

Hallo, I have this dataset with repeated measures. There are two within-subject factors, "formant" (2 levels: 1 and 2) and "f2 Ref" (25 levels: 670, 729, 788, 846, 905, 1080, 1100, 1120, 1140, 1170, 1480, 1470, 1450, 1440, 1430, 1890, 1840, 1790, 1740, 1690, 2290, 2210, 2120, 2040, 1950), and one between-subject factor, lang (2 levels:1 and 2). The response variable

Hi, I can't apply the summary function when I use de armaFit (fArma). Can you help? fit<-armaFit(~arma(1,0),data=age55) fit Title: ARIMA Modelling Call: armaFit(formula = ~arma(1, 0), data = age55) Model: ARIMA(1,0,0) with method: CSS-ML Coefficient(s): ar1 intercept -0.564184 0.001190 Description: Mon Feb 16 15:00:54 2009 by user: f003553

Hi, I am working on problem 2 of Chapter 8 in Data Analysis and Graphics Using R and don't know how to approach the second half of the question: In the data set (an artificial one of 3121 patients, that is similar to a subset of the data analyzed in Stiell et al., 2001) head.injury, obtain a logistic regression model relating clinically.important.brain.injury to other variables. Patients

Hi, I am running the following model: fit1.full <- coxme(Surv(age_sym1, sym1) ~ sex + lifedxm*sex + (1|famid), data=bip.surv) I would like to extract the AIC from that object to calculate the AICC. However, when I look at str(fit1.full) and summary(fit1.full) (pasted below) I don't see anything that would allow me to get pull the AIC out from that object. Is there a way to retrieve the

Dear R help-list, I am trying to calculate Monmoniers algorithm. I have 3 data sets: coordinates (coord), distance matrix (comdist) and geographical distance matrix (geodist). I consistently get the error " dimnames does not match array extent". However, when I use the function "structure" I get the following results: > str(coord) 'data.frame': 114 obs. of 3

Hi! I'm new to the list. I have a problem with mounting my home directory since my PC crashed. I hope that I can get some help on this list as I don't know much of ext3 myself. The mount command for my /home gives me the following output: # mount /dev/sda6 /mnt/tmp/ mount: wrong fs type, bad option, bad superblock on /dev/sda6, missing codepage or other error In some cases useful info

I have a samba version 3.024 server joined to a ADS as a "Workstation or server" acting as a file/print server. It is a Debian 4.0 amd64 up to date with stable. ftp works just fine so not network, must be samba config. vary slow file transfer when copying file from one computer but copy from two or more at the same time speed is 10 times faster. You will be copying a 600MG file with one