similar to: books on Time serie

Hi, I have managed to format my data into a single datframe consisting of two AsIs response and predictor dataframes in order to supply the plsr command of the pls package for principal components analysis. When I execute the command, however, I get this error: > fiber1 <- plsr(respmat ~ predmat, ncomp=1, data=inputmat,validation="LOO") Error in model.frame.default(formula =

Hello, Is possible to place a label on a multipanel figure outside the boundaries of any single one of the existing graphs in a panel? Specifically, I have a function that creates a panel with quartz() and then fills it with 6 graphs consisting of exactly the same plot for 6 different groups. I would like to add one title for the entire figure at the top of the figure, but I cannot figure out

Hi all! Maybe someone could help me with the following. I know this hasn't directly to do with ecology but I'm also using glm. I have a list of 16 genes and 10 samples. The samples are of two types, 4 Ctrl and 6 Diseased. If, labelInd<-as.factor(c(rep("0",4),rep("1",6))) genes.glm<-glm(labelInd ~ ., family=binomial, data=mat) beeing "mat" the 10x16

Hi, How do I use the string content of a string variable to reference a data frame of the same name? I want to do the typical tasks of 1) building a name with a string variable and using the string variable to create a data frame (or any object) whose name is the string value of the variable and 2) pass on a string to a function as a parameter, and then use that string to refer to an existing

Hello everyone, I have a couple of fairly simple questions (I hope) the answers to which I cannot find through the documentation at the moment. 1. I would like to delete the a row from a matrix if a certain elimination criterion is met. I am familiar with x <- x[-7,] (to remove row 7, for example). Are there any other means of removing an entire row? 2. Is there a single command that

Hello everyone, I am trying to prune a data frame for partial least squares analysis. I need to delete an entire row if one cell in the row contains a NA. Presently, I am running a loop that is supposed to extract the rows that are full of numbers into a second data frame and skips the rows that contain a single NA value. I want to know if there is a simple way to determine if a row

Dear R users, I have two simple questions here, and hope someone can help me on this. Thanks in advance. 1. I have a list object lst=list(a1=matrix(rnorm(4),2,2), a2=matrix(rnorm(4),2,2),a3=matrix(rnorm(4),2,2)). Here I only use three elements for illustration, and in fact the length of lst, n, is unknown in advance. I want to define an object for each element of this lst, and the objects have

Dear all, I am struggling with calculation of Schoenfeld residuals of my Cox Ph models. Based on the formula as attached, I calculated the Schoenfeld residuals for both non tied and tied data, respectively. And then I validated my results with R using the same data sets. However, I found that my results for non-tied data was ok but the results for tied data were different from R's. How

hi sorry if this has been discussed before, but I'm wondering why the scaled Schoenfeld residuals do not follow the defining formula for obtaining them from the ordinary Schoenfeld residuals, but are instead offset by the estimated parameter values. e.g. library(survival) attach(ovarian) sv<-Surv(futime,fustat) f1<-coxph(sv~age+ecog.ps) f1

Hi all, Although I sent the mail to Piotr, the author of FSelector, it should be better to ask here to let others know. Yanwei Begin forwarded message: From: Yanwei Song <yanwei.song@gmail.com> Date: September 7, 2011 4:41:58 PM EDT To: p.romanski@stud.elka.pw.edu.pl Subject: FSelector and RWeka problem Dear Piotr, Thanks for developing the FSelector package for us. I'm a new

Hi, I have a coxph model like coxph(Surv(start, stop, censor) ~ x + y, mydata) I would like to calculate the Schoenfeld residuals for the null, i.e the same model where the beta hat vector (in practical terms, the coeff vector spat out by summary()) is constrained to be all 0s --all lese stays the same. I could calculate it by hand, but I was wondering if there is a way of doing it with

Well.. plans are there from a long time to rewrite SciViews completely and make it platform independent (to work on Linux/Unix and MacOS X, as well as Windows). I have done some work in this direction when time permitted, but I am pretty busy with other work. During the holidays, I will continue to work in this direction. I will try to package a first running version of SciViews compatible with

I'm useing Samba (2.0.7) on an Compaq TRUE64 (Digital OSF) Server. Everything is working fine, but the dfree-script does not work. I think it's not even called. The corosponding entry in my smb.conf is: [global] server string = Samba %v Server, xxxx workgroup = xxxx keep alive = 30 dead time = 15 dfree command = /usr/local/samba/bin/dfree The script it self works fine

Thanks Philippe for your answer! It doesn't work for me, but I guess I have something wrong with Tcl/Tk. Ivan -- Ivan CALANDRA Universit? de Bourgogne UMR CNRS/uB 6282 Biog?osciences 6 Boulevard Gabriel 21000 Dijon, FRANCE +33(0)3.80.39.63.06 ivan.calandra at u-bourgogne.fr http://biogeosciences.u-bourgogne.fr/calandra -------- Message original -------- Sujet: Re: [R] choose folder

Dear R-Sig-Jobs members, A PhD position is available at the COMPLEXYS Research Institute of the University of Mons (UMONS), Belgium (http://www.complexsystems.be/). The PhD position is part of a joint research project on "Ecological Studies of Open Source Software Ecosystems (ECOS)". This interdisciplinary project between biology and computer science aims to use software repository

Dear all, I'm trying to use the "optim" function to replicate the results from the "glm" using an example from the help page of "glm", but I could not get the "optim" function to work. Would you please point out where I did wrong? Thanks a lot. The following is the code: # Step 1: fit the glm clotting <- data.frame( u =

Hello, Does anybody know if there is somewhere in R a function to calculate the Lyapunov exponent in a time series? Thanks, Philippe Grosjean .......................................................<??}))><.... ) ) ) ) ) ( ( ( ( ( Prof. Philippe Grosjean \ ___ ) \/ECO\ ( Numerical Ecology of Aquatic Systems /\___/ ) Mons-Hainaut University, Pentagone / ___ /( 8, Av. du

Hi I agree with your comments and of course you identifie the problems regarding the implementation of your suggestions. Of course we'd need bidireccional communication to implement it. I think it's impossible to do it statically. What I want with NEdit is a very ligth editor (which jedit is not) that highlights the most important language syntax, so that It becomes easier to write

> From: Philippe Grosjean > > Duncan, > > I agree totally with you on all points, now that we clarified our > respective ideas. I am afraid I probably agree also with your last > point, from a theoretical point-of-view ("I still think we need more > glue and am working on that while we continue to experiment with the > design of GUIs for the next 5 years

Hello, I am looking for more information about lazy loading introduced in R 2.0.0. Doing ?lazyLoad I got some and there is a 'see also' section that points to 'makeLazyLoading'... But I cannot reach this page. My problem is: I recompiled a library that uses a lot of functions from other libraries (of course I can give details if needed). I load it in my computer: library(svGUI),