similar to: apply'ing cut to a data.set column

Hi all, I'm new with R (and S), and relatively new to statistics (I'm a computer scientist), so I ask sorry in advance if my question is silly. My problem is this: I have a (sample of a) discrete time stochastic process {X_t} and I want to estimate Pr{ X_t | X_{t-l_1}, X_{t-l_2}, ..., X_{t-l_k} } where l_1, l_2, ..., l_k are some fixed time lags. It will be enough for me to compute

Full_Name: Massimo santini Version: 1.0.1 OS: Linux Submission from: (NULL) (159.149.147.89) I've just dl the .rpm of R 1.0.1, installed it and then launched insatll.packages("akima") after re-launching R [--vanilla], the command library(akima) has the only effect to coredump the program. Here is the output of gdb and strace... --- GDB (commands given can be seen from the

Hello, I have a big table with 3 columns and 103918 rows. This is the example,            time             species                 dbh 5 1 4.9377297 575 1 11.64127213 575 1 109.8182438 575 1 8.029809521 5 1 24.32501874 575 1 4.895992119 575 1 11.40567637 575 1 2.795090562 575 1 21.79281837 575 1 52.57476174 575 1 27.7290919 575 1 3.23262083 575 2 19.30612651 575 1

Dear R users, I have used the following function (in blue) aiming to find the linear regression between MOE and XLA and nesting my data by Species. I have obtained the following results (in green). model4<-lme(MOE~XLA, random = ~ XLA|Species, method="ML")summary(model4) Linear mixed-effects model fit by maximum likelihood Data: NULL         AIC     BIC   logLik  -1.040187 8.78533

How can locate in my program an error of that kind (Error: subscript out of bounds)? -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch

Hello! I have several variables. Each of them has a different distribution. I was thinking to use a Generalized Linear Model, glm(), but I need to introduce the family. Do you know if R has any tests for matching data to any distribution ( I am aware of shapiro.test). All the best, -- BAS [[alternative HTML version deleted]]

Hello list, I'm having problems with newest BIOS updates from Gigabyte to my Mother Boards, they don't want to boot syslinux, extlinux y usb-hdd mode. Then I need to know if someone else are having same problem, and can send me the Mother Board versions and bios versions to report as a bios bug. this is from the support, I'm awaiting for answer.

Dear all, Can anyone help me, my R software can not run a nested linear regression by using the lme funcion. The message that appears is  Error: could not find function "lme" I already downloaded and loaded the package, please see below. Thank you in advance for any help! Nadia. > data<-read.csv("/Users/nadiasan1/Desktop/MOE and MOR.csv")> attach(data)>

Dear R users, I have to call fortran program from within R (R 2.8.1 on ubuntu 8.10 machine). Suppose I have a fortran code like this (this is only a toy model, my working model is far more complex, but input/output is similar) DOUBLE PRECISION FUNCTION model(times, alfa, beta) DOUBLE PRECISION alfa, beta, times model=beta*sin(times)+alfa*cos(times) END FUNCTION which

Hello. I don't understand a particular output of portfolio.optim (tseries). I have 4 assets and the portfolio.optim returns an asset with weight equals to zero. If I do a portfolio.optim with 3 assets, without the asset with weight equals to zero, it returns a completely different result. That's I would expected the same weights as the run with 4 assets. Below the code. Thanks in

Hi, I am looking how to setup speex codec in codecs.conf for 8 Kbps and 6 Kbps. In config file are many parameters for setting. I don't know what is need to change for narrowbad like 8 Kbps and 6 kbps. Any suggestion? [speex] ;0-10 quality => 4 ;0-10 complexity => 4 ; true / false enhancement => true ; true / false vad => false ; true / false vbr => false ; 0 = off, otherwise,

Hi there, > > I am using the function pgls in the CAPER package, everything seems to run > fine except for one of my variables. > When using this variable as continuous the anova works without problem, > but if using the same variable as factor (in fact, this is what it is), the > program returns this answer: > > Error in terms.formula(formula, data = data) : >

Hi all I am analyzing a data set containing information about the behaviour of marine molluscs on a vertical wall. Since I have replicate observations on the same individuals I was thinking to use the geepack library. The data are organised in a dataframe with the following variables Date = date of sampling, Size = dimensions (mm) Activity duration of activity (min) Water = duration of

Hi everyone, I can't figure out how to format numbers to have a certain number of significant figures (as opposed to decimal places) without using scientific notation and including trailing zeros if necessary. e.g. I would like to achieve the following: 0.800001 ---> 0.8000 123.4567 ---> 123.4 0.1234567 ---> 0.1234 7.654321 ---> 7.654 7654321 --->

Suppose I have some data >x <- rnorm(1000); >y <- x*x; then try to cut it into 2 chunks, >c <- cut(y, breaks=2); >summary(y) Min. 1st Qu. Median Mean 3rd Qu. Max. 6.879e-06 9.911e-02 3.823e-01 9.499e-01 1.297e+00 8.342e+00 >summary(c) (-0.00833,4.17] (4.17,8.35] 958 42 Is that the correct behavior? Why is the left hand

tabulate has .C("R_tabulate", as.integer(bin), as.integer(length(bin)), as.integer(nbins), ans = integer(nbins), PACKAGE="base")$ans The implementation of R_tabulate has if(x[i] != R_NaInt && x[i] > 0 && x[i] <= *nbin) and so copes with (silently drops) NA. Perhaps the .C could have NAOK=TRUE? This is useful in apply'ing tabulate to

I came to report this same bug and found it already in the trash, but I slightly disagree with that assessment. If it's not a bug, then perhaps it's a feature request. Comments at the end. On Mon, May 14, 2007, Duncan Murdoch wrote: >>On 13/05/2007 8:46 PM, scott.wilkinson at csiro.au wrote: >> >> In the example below round() does not report to the specified number of

Hi all I have this problem: In my database .dta, called "data" I have five rows data<-read.dta("C:\\2_CO_mmobile_ALL_Rid.dta") # From this database I wuold like to create another coppie<-c(data[4:length(data)]) but I find this # Length of original data length(data[,4]) 5 RIGHT!! # Length of new data

Hi, the problem I am reporting relates to RH Linux 6.0 I have successfully mounted a disk from a windows NT server machine, using smbmount as follows: $ smbmount //machine/service -c 'mount /mnt/NTgroups -u local_owner' -U remote_username and everything looked fine. Now I have just noticed that when "touching" or creating files everything behaves strangely, that is the date is

Hello everybody: I have just started using R, specifically the clustering packages. Despite the vastness of the functionality of the system, I have found it very flexible and intuitive. It has been great fun to experiment with it. Nevertheless, I am under tight time constraints to complete a cluster analysis of some data and need to come up with a "quick" solution for the following