search for: carrillo

 Hi, What's the equivalent of "Save to File" from the R console File menu on an R routine? Just trying to capture the whole R console into a text file when my code fails. Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA http://www.fws.gov/redbluff/rbdd_jsmp.aspx [[alternative HTML version deleted]]

...n(name = 'Número\nde\nespecies',  colours = rainbow(6), breaks = c(2, 4, 6, 8, 10, 12, 14, 16, 18, 20))+  xlab("Longitud") + ylab("Latitud") + opts(axis.text.x = theme_text(size = 8, vjust = 1)) +  opts(axis.text.y = theme_text(size = 8, hjust = 1))        Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA http://www.fws.gov/redbluff/rbdd_jsmp.aspx > >From: Manuel Spínola <mspinola10@gmail.com> >To: Felipe Carrillo <mazatlanmexico@yahoo.com> >Sent: Sat, April 2, 2011 11:2...

...R-Forge.R-project.org) ?Warning: unable to access index for repository http://R-Forge.R-project.org/bin/windows/contrib/2.10 Warning message: In getDependencies(pkgs, dependencies, available, lib) : ? package ?gridextra? is not available I would like to download the binary for windows Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

...s trying i<- nrows(weekly$a) i but returns 7 when it should be 4. Thanks a b c d 27.000 27.000 1.569 0.013 160.000 27.000 1.632 0.013 146.000 27.000 1.830 0.015 70.000 27.000 2.475 0.019 156.000 27.000 27.000 Felipe D. Carrillo Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA ____________________________________________________________________________________ Be a better friend, newshound, and

...but I want seq to start at 27 by=2, but when it reaches 51 start with with number 1 to 25. is this possible. I can do the basics of seq() but I can't figure how to do this one. This is how I want my sequence to look like: 27 29 31 33 35 37 ............51 1 3 5 7 9 11 13 ...........25 Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

...(week,value/1000)) + geom_line(aes(colour=variable))+ opts(legend.position="none") + facet_wrap(~variable,ncol=2) + opts(title="Winter") + labs(y="Number of fish X 1,000",x="WEEK")) @ \caption{Weekly estimates.} \label{figure:ggplot1} \end{figure} ? Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA ? Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

...te dates with the month and day? I tried the code below: seq(as.Date("7-1","%B%d"),as.Date("6-30","%B%d"),by="week") and got this?error message: Error in seq.int(0, to - from, by) : 'to' must be finite Thanks for any pointers ? Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

...t;- read.table(textConnection(x), header = TRUE) p <- ggplot(DF,aes(x=Location)) + geom_bar() p + geom_text(aes(y=Location),label=sum(count)) # Error because count doesn't exist in dataset What should I use instead of 'count' to be able to sum the number of Locations? Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

...f my data is skewed then use wilcox.test. Something like the pseudo code below: fishlength <- c(35,32,37,39,42,45,37,36,35,34,40,42,41,50) if fishlength= "normally distributed" then t.test(fishlength) else wilcox.text(fishlength) I hope this isn't very confussing Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

...0.0000 \\ male & $-$0.4094 & 0.0482 & $-$8.49 & 0.0000 \\ langarts & $-$0.0147 & 0.0013 & $-$11.34 & 0.0000 \\ \hline \end{tabular} \end{center} \end{table} I would appreciate it if someone can point me to some examples about how to do this. Thanks Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

...### predict for another 2 weeks newD <- data.frame(week = 1:52);newD newD$pred_wt <- predict(M_model, newD) newD plot(pred_wt ~ week, newD, pch = 4, col = "red", ylab = "Weight", xlab = "Week") with(weightData, points(week, weight,col='blue')) Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA http://www.fws.gov/redbluff/rbdd_jsmp.aspx [[alternative HTML version deleted]]

...13716163737 0.000145848 8.59635E-07 11742149934 0.000155914 . . . 8.9031E-07 10039878212 0.0001513 9.1802E-07 10665936164 0.000146999 8.52681E-07 9130907448 0.000157145 . . . 9.10279E-07 8543438113 0.000148139 8.4587E-07 7313876766 0.000158363 . . . 8.71193E-07 7769949102 0.000153862 Felipe D. Carrillo Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

Hi all: How can the R graphics window be customized programmatically? Either minimized,maximized or change the size of the default that ships with R. Thanks Felipe D. Carrillo Fishery Biologist US Fish & Wildlife Service Red Bluff, California 96080

....0122 2.392 0.0146 2.244 0.0175"), header=TRUE) anova(lm(Efficiency~percentQ,data=dfr)) Analysis of Variance Table Response: Efficiency Df Sum Sq Mean Sq F value Pr(>F) percentQ 1 0.00014432 0.00014432 4.3774 0.04561 * Residuals 28 0.00092312 0.00003297 Felipe D. Carrillo Fishery Biologist US Fish & Wildlife Service California, USA ____________________________________________________________________________________ Looking for last minute shopping deals?

...usedLevels() doesn't seem to work here ? dropUnusedLevels() ggplot(firstyear,aes(year,fatPerc)) + geom_boxplot() + geom_jitter() # code below also should drop levels but it doesn't #data.frame(lapply(firstyear, function(x) if (is.factor(x)){ factor(x)} else{x})) str(firstyear) ? Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

...010194 0.0000449 417 133067 29 0.012567 0.0000449 82 10184 30 0.014695 0.0000449 413 68893 30 0.010634 0.0000449 373 103688 30 0.013934 0.0000449 914 272719 30 0.014107 0.0000449 213 28473 30 0.014339 0.0000449 651 147016 30 0.014375 0.0000449 521 102429 30 0.016373 0.0000449 979 229104 Felipe D. Carrillo Fishery Biologist US Fish & Wildlife Service California, USA ____________________________________________________________________________________ Be a better friend, newshound, and

...e max Total for 2012   x <- subset(winter,Year==2012 & Total==max(Total));b  # How come one line doesn't work?     # It works if I subset the year first and then get the Total max from it   xx <- subset(winter,Year==2012)  xxx <- subset(xx,Total==max(Total));xxx   xxx Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA http://www.fws.gov/redbluff/rbdd_jsmp.aspx [[alternative HTML version deleted]]

...r{webgreen}{rgb}{0, 0.5, 0} % less intense green \definecolor{webblue}{rgb}{0, 0, 0.5} % less intense blue \definecolor{webred}{rgb}{0.5, 0, 0} % less intense red \usepackage[all]{hypcap} \pdfbookmark[1]{Contents}{table} \title{BROOD-YEAR 2007 WINTER CHINOOK JUVENILE PRODUCTION} \author{Felipe D. Carrillo} \date{January 5, 2009} \pagestyle{headings} \parskip 7.2pt \begin{document} \setkeys{Gin}{width=1.1\textwidth} \maketitle \tableofcontents The hyperlink of table 1 doesn't work!!!! but the hyperlink of table 2 works OK \newpage \listoftables \listoffigures \section{Introduction} \label{sec...

...at the weight will be by week 10 based on my weight values and make a line plot of all the weights(including the predicted values). I have two questions: 1- Should the predicted values be linear or exponential? 2- Is the predict.arima function appropriate to do this? Thanks in advance. Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife Service California, USA

Dear r-users,   I have daily rainfall data from 1971 to 2000.  I would like to extract november and december data only.  I would also like to do column bind for november and december, therefore I would like to delete 31 December from december data so that the length of november and december are the same.  Hope somebody can help me.  I tried this below:   > kuantan.dt.1 <-