similar to: Inverse t-distribution

Manuel: I changed your variable names from x to 'long' and y to 'lat' on the riqueza_out.csv file. The code below should do what you want. Also, since the legend title is kind of long, I broke it down into three lines so you can see more plot area. I am cc'ing the other groups so more people use it if needed. library(rgdal) library(ggplot2) library(sp) library(maptools)

Hi all: Can someone help me count the number of rows with values in colum "a" only. assume the name of my dataframe is "weekly" I was 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

 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]]

? Hi: I get?a couple of warnings? when trying to download gridExtra: install.packages("gridExtra",repos=http://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

Hi: I am using Sweave and texi2dvi to generate a LaTeX document but can't find the way to hide the graphics while the R chunks are being executed. I thought results=hide would do it but that't not the case. If I do: \begin{figure}[h] <<figA=true,echo=F,fig=T,results=hide>>= a < rnorm(1000) plot(a) @ \caption{Weekly estimates.} \label{figure:ggplot1} \end{figure} The

HI: Could someone help me with the seq function? I have a range of values starting from 1 to 52 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.

Hi: I have a dataframe named 'spring' and I am trying to add a new variable named 'IdDate' This line of code works fine: spring$idDate <- seq(as.Date("2008-07-01"),as.Date("2009-06-30"),by="week") But I don't want to hardcode the year because it will be used again the following year Is it possible to just generate dates with the month and

Hello all: Is there a way to detect in R if a dataset is normally distributed or skewed without graphically seeing it? The reason I want to be able to do this is because I have developed and application with Visual Basic where Word,Access and Excel "talk" to each other and I want to integrate R to this application to estimate confidence intervals on fish sizes (mm). I basically want to

Hi: Can anyone advice me on how to loop and perform a calculation through all the columns. here's my data xd<- c(2.2024,2.4216,1.4672,1.4817,1.4957,1.4431,1.5676) pd<- c(0.017046,0.018504,0.012157,0.012253,0.012348,0.011997,0.012825) td<- c(160524,163565,143973,111956,89677,95269,81558) mydf<-data.frame(xd,pd,td) trans<-t(mydf) trans I have these values that I need to

Hi: I need help with geom_text(). I would like to count the number of Locations and put the sum of it right above each bar. x <- "Location Lake_dens Fish Pred Lake1 1.132 1 0.115 Lake1 0.627 1 0.148 Lake1 1.324 1 0.104 Lake1 1.265 1 0.107 Lake2 1.074 0 0.096 Lake2 0.851 0 0.108 Lake2 1.098 0 0.095 Lake2 0.418 0 0.135 Lake2 1.256 1 0.088 Lake2 0.554 1 0.126 Lake2 1.247 1 0.088

 Hi:  Using nls how can I increase the numbers of iterations to go beyond 50.  I just want to be able to predict for the last two weeks of the year.  This is what I have:  weight_random <- runif(50,1,24)  weight <- sort(weight_random);weight weightData <- data.frame(weight,week=1:50)                          weightData plot(weight ~ week, weightData) M_model <- nls(weight ~ alpha +

HI: Is it possible to use the RODBC package within MS Access. I have been using from R but was just wondering if it could be used along with R(D)COM. Something like this: Dim dbs As DAO.Database Dim rst As DAO.Recordset Dim myApp As StatConnector Set myApp = New StatConnector myApp.GetErrorText myApp.Init "R" myApp.EvaluateNoReturn "library(RODBC)" myApp.EvaluateNoReturn

Hi all: I am having trouble dropping levels, got a few hints online?without success. Please consider the dataset below: ?I was under the inpression that subset(......drop=TRUE) would work but it doesn't library(ggplot2) ??? library(hmisc) x <- structure(list(first = c(38.2086, 43.1768, 43.146, 41.8044, 42.4232, 46.3646, 38.0813, 40.0745, 40.4889, 38.6246, 40.2826, 41.6056, 34.5353,

Hi: After a couple of days trying to synchonize or connect MiKTex and R I'm getting a little frustrated. My knowledge about creating and running batch files is kind of limited so, I was wondering if someone could explain in plain english how to do this task. I'm able to run MiKtex by itself and accomplish what I am trying to do, but I want to be able to use the Sweave function from within

Hi everyone: I did this before with R and I can't remember how. I got some fish forklength values FL number 34 4 35 3 36 7 I need to breakdown the FL by the number of fish with the same length like this: 34 1 34 1 34 1 34 1 35 1 35 1 35 1 36 1 36 1 36 1 36 1 36 1 36 1 36 1 Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish

Hi: Is my first attempt to try to download from github. Nothing happens by clicking on the 'download' button. Could anyone give me a hint on how to get all the files from the link below? Thanks http://github.com/hadley/ggplot2-bayarea/tree/0a8bf71dea38cfbf2d928eb713d24dfd928359fc Felipe D. Carrillo Supervisory Fishery Biologist Department of the Interior US Fish & Wildlife

Hi: Can't find?a way to convert?from shortDate to LongDate format. I?got: 3/10/10 that I want to convert to March 10, 2010. I am using: \documentclass[11pt]{article} \usepackage{longtable,verbatim} \usepackage{ctable} \usepackage{datetime} \title{my title} \begin{document} ? % Convert date \mmmmddyyyydate\3/10/10 end{document}? My report is changing every two weeks so I?will eventually use

Hi Hadley: I'm practicing writing a document using MikTex but can't make qplot to execute. It works when using 'plot' though. Is this a normal behaviour? \documentclass[9pt]{article} \title{ggplot2 example} \begin{document} \maketitle \section*{Examples of using ggplot2} The goal is to be able to import ggplot2 graphics into the annual report. Hadley Wickham has done a great

Hi: Using this dataframe with quite long column headers, how can I wrap the text so that the columns are narrower. I was trying to use strwrap without success. Thanks reportDF <- structure(list(IDDate = c("3/12/2010", "3/13/2010", "3/14/2010", "3/15/2010"), FirstRunoftheYear = c("33 (119 ? 119)", "n (0 ? 0)", "893 (110 ?

Hello: I have 5 years of weekly passage data and want to predict fish passage for the following year. I don't have a model to use to predict data for the sixth year. Can I somehow still predict based on these five years? I just want to see on the graph what the predicted year would look like and how those?new values are generated. Gracias? ? ?fall <- structure(list(week = c(48L, 49L, 50L,