search for: grave

Displaying 20 results from an estimated 933 matches for "grave".

Did you mean: gave
2008 Feb 02
ARCH LM test for univariant time series
Hi, Does anyone know if R has a Lagrange multiplier (LM) test for ARCH effects for univariant time series? Thanks! -- Tom [[alternative HTML version deleted]]
2005 Jul 06
Plotting confidence intervals for lme
Hello and sorry to disturb. I'm trying to plot the confidence intervals for the fixed effects of a lme. I want to obtain graphically, if it is possible, a bar with Estimate, upper and lower CI for each level of the factors. I know how to do for a lm model but for a lme one, I tried with plot(intervals(...)) and plot(ci(...)) from the gmodels package but it doesn't work well. Thanks
2005 Nov 11
no package 'Matrix' at the repositories
...t the repositories I installed lme4, maps, mapproj, CircStats, scatterplot3d, gregmisc, Hmisc without problems. To confirm, 'library(lme4)' produced the following error: Error: package 'Matrix' required by 'lme4' could not be found What do you suggest? Spencer Graves p.s. I get the same result using several different (US) mirrors. -- Spencer Graves, PhD Senior Development Engineer PDF Solutions, Inc. 333 West San Carlos Street Suite 700 San Jose, CA 95110, USA spencer.graves at <> Tel: 408-938-4420 Fax: 4...
2005 Sep 28
R-code for binormla distribution
Dear users, does any one have a code (S or R) to compute the binormal distribution (or the upper its quadrant area) other than the pmvnorm. Thanks -- Nabil Channouf etudiant en Ph.D. Bureau 3733 Departement d'Informatique et de Recherche Operationnelle (D.I.R.O.) Universite de Montreal, C.P. 6128, succ. Centre-Ville, Montreal, H3C 3J7 Tel.: (514) 343-6111, poste 1796 Fax.: (514) 343-5834
2003 Jun 01
compositional data: percent values sum up to 1
again, under another subject: sorry, maybe an all too trivial question. But we have power data from J frequency spectra and to have the same range for the data of all our subjects, we just transformed them into % values, pseudo-code: power[i,j]=power[i,j]/sum(power[i,1:J]) of course, now we have a perfect linear relationship in our x design-matrix, since all power-values for each subject sum
2005 Aug 12
General expression of a unitary matrix
Hi, all, Does anybody got the most general expression of a unitary matrix? I found one in the book, four entries of the matrix are: (cos\theta) exp(j\alpha); -(sin\theta)exp(j(\alpha-\Omega)); (sin\theta)exp(j(\beta+\Omega)); (cos\theta) exp(j\beta); where "j" is for complex. However, since for any two unitary matrices, their product should also be a unitary matrix. When I
2005 Jul 08
Finite Mixture Models with logistic regression
Do we have any R package that can do analysis on finite mixture model with logistic regression? Thanks Faith Feng Gao Dept. of Statistics Virginia Tech. Email: fgao at
2005 Aug 17
power of a matrix
Dear all, I have a population with three age-classes, at time t=0 the population is: <- c(1,0,0) I have a transition matrix A which denotes "fertility" and "survival": A <- matrix(c(0,1,5, 0.3,0,0, 0,0.5,0), ncol=3, byrow=TRUE) To obtain the population at t=1, I calculate: A %*% To obtain the population t=2, I calculate: A %*% (A %*% ... and so
2009 Feb 09
Michael Graves post
Michael Grave just posted a question about surround conferences. 3908&id=564633430&index=0 I didn't see it posted on the ast-list, what do you think? Does something like this have potential? I'...
2005 Sep 19
FDR analyses: minimum number of features
Dear List, We are planning a genotyping study to be analyzed using false discovery rates (FDRs) (See Storey and Tibshirani PNAS 2003; 100:9440-5). I am interested in learning if there is any consensus as to how many features (ie. how many P values) need to be studied before reasonably reliable FDRs can be derived. Does anyone know of a citation where this is discussed? Bill Dupont William D.
2004 Feb 26
RE: system.time(), sys.time() etc
> From: Spencer Graves [mailto:spencer.graves at] > > Martin says, "This is another instance of S-Plus following R > behind and doing it incompatibly [with a reason?] ... ." > > This is one example of a major issue in "how to wage and win a > standards war",...
2005 Jun 15
Kalman Filtering?
...2) Does anyone have any comments on this? The data for 1892, 1894, and 1897 seem slightly outside the 90% confidence interval in Durbin and Koopman (2001, Fig. 2.1) but just inside the figure created by this code. Am I doing this correctly? Is there an easier way? Thanks, spencer graves
2018 Feb 20
deparseDots to get names of all arguments?
...s(...)? ????? I'm able to get the name of the first argument but not the second: deparseDots <- function(...){ ? deparse(substitute(...)) } a <- 1 b <- 2 deparseDots(a, b) [1] "a" ????? I'd like to get c('a', 'b'). ????? Thanks, ????? Spencer Graves > sessionInfo() R version 3.4.3 (2017-11-30) Platform: x86_64-apple-darwin15.6.0 (64-bit) Running under: macOS High Sierra 10.13.3 Matrix products: default BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versi...
2005 Jul 11
exact values for p-values - more information.
Hi there, If I do an lm, I get p-vlues as p-value: < 2.2e-16 This is obtained from F =39540 with df1 = 1, df2 = 7025. Suppose am interested in exact value such as p-value = 1.6e-16 (note = and not <) How do I go about it? stephen
2004 Jul 24
help(package)->sink()", while displaying mclustInfo on my screen. After a few more failed attempts, I gave up trying to be clever and copied the text from the screen into Word. Is this a bug, that sink does not capture the output of objects of class packageInfo? Best Wishes, spencer graves
2004 Aug 28
removing invariant columns from a matrix
I'm looking for an efficient way of removing zero-variance columns from a large matrix. Any suggestions? Thanks, - Moises [[alternative HTML version deleted]]
2013 Feb 04
rJava works with 32-bit but not 64
...States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] rJava_0.9-3 -- Spencer Graves, PE, PhD President and Chief Technology Officer Structure Inspection and Monitoring, Inc. 751 Emerson Ct. San Jos?, CA 95126 ph: 408-655-4567 web:
2013 May 20
Objects created by more than one data call?
...t argument for a function. With LazyData, that NOTE disappears. However, then I get, "Warning: objects 'Hstarts', 'Hstarts', 'MedExp' are created by more than one data call". What do you suggest I do to fix this problem? Thanks, Spencer Graves > sessionInfo() R version 3.0.0 (2013-04-03) Platform: i386-w64-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=English_United States.1252 [2] LC_CTYPE=English_United States.1252 [3] LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C [5] LC_TIME=English_United States.1252 attached ba...
2018 Feb 21
deparseDots to get names of all arguments?
...sin(3:1), z=stop("Oops"), "untagged arg") $y 1/(1:10) $z stop("Oops") [[3]] [1] "untagged arg" > names(.Last.value) [1] "y" "z" "" Bill Dunlap TIBCO Software wdunlap On Tue, Feb 20, 2018 at 2:36 PM, Spencer Graves < spencer.graves at> wrote: > Hi, All: > > > How can I get the names of all the arguments in dots(...)? > > > I'm able to get the name of the first argument but not the second: > > > > deparseDots <- function(...){ &gt...
2004 Jun 04
Error() term in glm model formula
Hello, My data are numbers of trees in plots sampled in a number of forest stands. Some stands were subjected to a treatment, others not. Several plots were sampled per stand to get a better idea of what the stand means were, but replication is really at the stand level. Therefore I think this is a split-plot design. I would like to know whether the treatment affected the number of trees, so: