search for: cade

Just noticed that I get a similar error about object 'kronecker' in "Matrix" package when trying to load "lme4". So this is a more pervasive problem. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade@usgs.gov tel: 970 226-9326 ----- Forwarded by Brian S Cade/BRD/USGS/DOI on 11/30/2012 02:07 PM ----- From: Brian S Cade/BRD/USGS/DOI To: R-help@r-project.org Date:...

...riable (dbc), i.e, shift them down by 1 row replacing the first row with NA, within groups defined by locid nested within population. I thought I could do something using by(csss3,list(locid, population), function) but don't seem to be having any success. Any suggestions?? Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade@usgs.gov tel: 970 226-9326 [[alternative HTML version deleted]]

I am trying to implement an ANOVA on a pair of quantile regression models in R. The anova.rq() function performs a basic check to see whether the models are nested, but I think this check is failing in my case. I think my models are nested despite the anova.rqlist() function saying otherwise. Here is an example where the GLM ANOVA regards the models as nested, but the quantile regression ANOVA

...ind a problem that I know how to fix, I'ld like to fix it - once. I suspect I'm missing something pretty obvious about how to edit these but my search of the online documentation, etc. didn't suggest any simple approach like I use to do. Any help would be welcome. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade@usgs.gov tel: 970 226-9326 [[alternative HTML version deleted]]

...(Windows, R 1.91)? I've several times started functions that for various reasons are not executing properly and it would be nice to stop them without killing the R session. I've been using taskmgr to end the R session but I then lose all objects created during the session. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade at usgs.gov tel: 970 226-9326

...on 70F221 1730 ESE 20 Hardwood ACMA Class1 0 4 Modern Cameron 70F221 1730 ESE 20 Hardwood AECA Class1 0 5 Modern Cameron 70F221 1730 ESE 20 Hardwood ARME Class1 0 6 Modern Cameron 70F221 1730 ESE 20 Conifer CADE Class1 15 7 Modern Cameron 70F221 1730 ESE 20 Hardwood CELE Class1 0 8 Modern Cameron 70F221 1730 ESE 20 Hardwood CONU Class1 0 9 Modern Cameron 70F221 1730 ESE 20 Conifer JUCA Class1 0 10 Modern Cameron 70F221...

...s them on her version of R (not sure what number hers is) on her Mac, she gets these warnings constantly. I've checked some records manually by doing the algebra and the predict.glm() function is working correctly incorporating the factor levels on my machine. Any thoughts??? Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: cadeb at usgs.gov <brian_cade at usgs.gov> tel: 970 226-9326 [[alternative HTML version deleted]]

...the same on the 2 machines. If so, ignore the above. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Tue, Dec 5, 2017 at 12:17 PM, Cade, Brian <cadeb at usgs.gov> wrote: > I am helping a student with some logistic regression analyses and we are > getting some strange inconsistencies regarding a warning about factor > levels being dropped when running predict.glm(, newdata = ournewdata) on > the logistic regressio...

2012/3/23 Vladi Belperchinov-Shabanski <cade at datamax.bg> > > hello, > Hi, I'm adding -upsuser list for info > from the source documentation I figured you are dev coordinator. > I'm sorry if I got it wrong and please forward to the correct person > if possible, thanks! > I'm indeed leading NUT, but...

...? I have not been able to relate them to the $coefficients, $fitted.values, $resid, or $y completely. For example, if I estimate the simplest model with just an intercept term, the $effects for the n observations differ from $y by a constant quantity except for one of the values. Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: cadeb at usgs.gov <brian_cade at usgs.gov> tel: 970 226-9326 [[alternative HTML version deleted]]

...71 -0.1455 -0.0784 42.81 <.0001 n 1 0.0001 0.0163 -0.0318 0.0320 0.00 0.9941 Scale 0 1.0000 0.0000 1.0000 1.0000 NOTE: The scale parameter was held fixed. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade@usgs.gov tel: 970 226-9326 [[alternative HTML version deleted]]

Your understanding isn't similar to mine. Mine says robust/resistant methods are for data with heavy tails, not heteroscedasticity. The common ways to approach heteroscedasticity are transformation and weighting. The first is easy and usually quite effective for dose-response data. The second is not much harder. Both can be done in R with nls(). Andy From: Quin Wills > > I am

...oading required package: SparseM Error : object ?kronecker? is not exported by 'namespace:methods' Error: package ?SparseM? could not be loaded or Error: object ?kronecker? is not exported by 'namespace:methods' During startup - Warning message: unable to restore saved data in C:\CADESTUFF\DATA\BarryNoon\.RData > Brian Brian S. Cade, PhD U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade@usgs.gov tel: 970 226-9326 [[alternative HTML version deleted]]

Dear moderator, If possible I would like to send in the following question for R-help: I am analyzing a small data set using PGLS with phylogenetic uncertainty taken into account and thereby including 100 potential phylogenetic tree scenarios. I've managed to run models on all of the different trees and performed model averaging to get parameter estimates for the intercept and most of the

Maybe it''s all in my head, but I''m having the strangest little issue ever. My pricetags in my rails app are always whole numbers. Rails seems to completely ignore my decimal input. So as a consequence, 10.41 becomes 10.00, always. Here''s some worthwhile stuff: within my 004_create_products.rb file: :price, :decimal, :percision => 2, :default => 0 And calling

...quot; out and look nice, but then the symbols have been decomposed into their constituent elements and are screwed up. I'm guessing somewhere in R's graphing parameters/controls there might be a solution but I've yet to find it. Any suggestions would be welcome. Brian Brian S. Cade U. S. Geological Survey Fort Collins Science Center 2150 Centre Ave., Bldg. C Fort Collins, CO 80526-8818 email: brian_cade@usgs.gov tel: 970 226-9326 [[alternative HTML version deleted]]

...ription #---------??????????????????????? ----------------------------- ?DEPT.M?????????????????????????? :DEPTH ?GRDE.GAPI??????????????????????? :GAMMA FROM DENSITY TOOL ?CODE.G/C3??????????????????????? :COMPENSATED DENSITY??? ? # # # ~A? Depth?????? GRDE?????? CODE?????? LSDU?????? BRDU?????? CADE?????? DENL?????? DENB?????? ADEN?????? VL2F?????? VL4F?????? VL6F?????? VL2A ??? 51.000??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25 ??? 50.990??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25??? -999.25???...

Hi I am using the quantreg package for median regression for a large series of subsets of data. It works fabulously for all but one subset. When it reaches this subset, R takes the command and never responds. I end up having to kill R and restart it. It appears to be something with the particular data subset, but I can't pinpoint the problem. Here are some details Operating system:

Hi, I am currently trying to do a GLMM on a dataset with percent cover of seagrass (dep. var) and a suite of explanatory variables including algal (AC) and epiphyte cover (EC), rainfall, temperature and sunshine hours. M2=glmer(SG~AC+EC+TP+SS+RF+(1|Location/fSi/fTr), family=binomial,data=data,nAGQ=1) As the dependent variable is percent cover, I used a binomial error structure. I also have a

I am perfectly aware that this question is not an R question, at least not yet, but I have not succeeded in finding what I want in other ways, so ... What I am looking for are two algorithms, preferabley in Pascal, but other languages may do. For (a) systematic (complete) permutations for grouped data with unequal group sizes, and (b) random permutations for the same kind of data. I know