[R] spatial analysis using quickPCNM

ompute PCNMs = 0.82 sec Error in if (temp2.test[1, 5] <= alpha) { : argument is of length zeroTiming stopped at: 1.06 0.05 1.19 I do not understand the error message coming up and would appreciate some advice. Andy -- Andrew Halford Ph.D Research Scientist (Kimberley Marine Parks) Dep

[R] multiple graphs on one plot

d(mf[[1]])),type="l", col="dark gray") Even trying to add an abline to the plot doesn't work. What am I missing? cheers Andy -- Andrew Halford Ph.D Senior Coastal Fisheries Scientist Pacific Community | Communauté du Pacifique CPS – B.P. D5 | 98848 Noumea, New Caledoni

Re: [R] multiple graphs on one plot

uot;All Measured > Crabs", col=c("dark gray", "light gray"), > breaks=seq(90,210, by=10),beside=TRUE,space=c(0,0.5)) > legend("topright", c("Females", "Males"), fill=c("dark gray", "light > gray")) > l

[R] defining group colours in a call to rda

r 2 groups and leaves the rest uncoloured. The call to MI_fish_all.mrt$where gives... [1] 3 3 8 6 6 9 5 5 9 3 8 6 9 6 5 9 5 3 8 6 9 6 5 9 5 3 3 8 6 6 9 5 5 9 6 9 5 9. These are the split groupings for all 39 sites in the analysis and there are 5 numbers corresponding to 5 final leaves in th

[R] Fwd: defining group colours in a call to rda

-- Forwarded message - From: Abby Spurdle Date: Wed, Aug 5, 2020 at 3:07 PM Subject: Re: [R] defining group colours in a call to rda To: Andrew Halford Hi Andrew, Perhaps you want this: cols <- rep_len (c ("red", "green", "blue", "aq

[R] incompatible dimensions error

al=nrow(fish05),xvmult=100,which=4,pca=TRUE) The tree is produced no problem but it wont produce a PCA. I am just keen to understand what this error means as I dont see anything unusual about the dataset used, notwithstanding the data is rather messy. Andy -- Andrew Halford Ph.D Senior C

[R] extracting x and y coordinates from rnorm

nother package I am using. Andy -- Andrew Halford Ph.D Adjunct Research Scientist University of Guam & Curtin University Ph: +61 (0) 468 419 473 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://st

Re: [R] grImport/ghostscript problems

\Users\ahalford\AppData\Local\Temp\Rtmp6BOVDe\fileffc30613d6 -sstdout=fish.ps.xml capturefish.ps' Any advice appreciated. Andy -- Andrew Halford Ph.D Adjunct Research Scientist University of Guam & Curtin University Ph: +61 (0) 468 419 473 [[alternat

Re: [R] grImport/ghostscript problems

27; What am I doing wrong here? Andy On 14 August 2013 12:10, Andrew Halford wrote: > Hi Listers > > I have been trying to import a .ps graphic file into R using the grImport > package but I keep getting the following error message > > Error in PostScriptTrace("fish.ps&quo

[R] loop for calculating 1-se in rpart

0));x=printcp(fit1);x[which.min(x[,"xerror"]),"nsplit"]}) > hist(cp50,main="optimal splits for tree",xlab= "no. of optimal tree splits", ylab= "frequency") Any help appreciated. Andy -- Andrew Halford Ph.D Adjunct Researc

[R] repeating an analysis

ts after 50 runs. Any help appreciated. Andy -- Andrew Halford Associate Researcher Marine Laboratory University of Guam Ph: +1 671 734 2948 [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/lis

Re: [R] repeating an analysis

il Spector > Statistical Computing Facility > Department of Statistics > UC Berkeley > spec...@stat.berkeley.edu > > > > &

Re: [R] repeating an analysis

tatistical Computing Facility > Department of Statistics > UC Berkeley > spec...@stat.berkeley.edu > > > > > > On Wed, 13 Oct 2010, Andrew Halford wrote: > > Hi All, >>

[R] penalized regression analysis

Loglikelihood = -154.1055 L1 penalty = 4944.889 at lambda1 = 356.0856 L2 penalty = 234.7781 at lambda2 = 3.458605 > coefficients (fitfinal) (Intercept)BC POC EXPFI 4.685739e+01 2.074521e-01 1.079459e-01 -1.373058e-05 -2.295339e

[R] goodness-of-fit test

0 Warning message: In summary.goodfit(gf) : Chi-squared approximation may be incorrect Am I getting caught out because of zero counts or frequencies in my data? Andy -- Andrew Halford Ph.D Associate Research Scientist Marine Laboratory University of Guam Ph: +1 671 734 2948 [[altern

[R] ordiellipse

x27;m assuming that 2 points is not enough to perform the relevant calculations here, however I would like to plot one if I could, if only for the sake of pictorial consistency. Andy -- Andrew Halford Ph.D Associate Research Scientist Marine Laboratory University of Guam Ph: +1 671 734

[R] which multivariate regression?

ccount for variation in one of the EVs affecting others. I appreciate any help I can get, -- Andrew Halford Ph.D Associate Research Scientist Marine Laboratory University of Guam Ph: +1 671 734 2948 [[alternative HTML version deleted]] __ R-help@r-p

[R] reflecting a PCA biplot

ecies"),type="n",scaling=3) # how do I get the plot to draw up the blank display based on the reversed site and species scores? Any help appreciated. cheers Andy -- Andrew Halford Ph.D Associate Research Scientist Marine Laboratory University of Guam Ph: +1 671 734 2948

[R] plot 3d info in 2d

axes would be larger while sites that are actually further away in 3 d space would be proportionally smaller. any help/advice appreciated Andy -- Andrew Halford Ph.D Associate Research Scientist Marine Laboratory University of Guam Ph: +1 671 734 2948 [[alternative HTML version deleted

[R] plotting issues with PCA

isplay="species",scaling=3)[omanfish.mrt.indval$pval<=0.05,] > graph <- plot(fish.pca,display=c("sites","species"),type="n",scaling=3) > plotcolor <- c("red","green","blue&quo