similar to: failure to access packages

Having given myself carpal tunnel looking for answer to this ... I have a dataset each column of which has 12 rows in it. I created a variable 'z' as follows: z=1:24 Since I have a large number of these plots to make, and they are a bit complex, I want to want to reference the column I want to plot via a variable containing the name of that column. As follows:

I have a set of 52 weekly values, one is the desired high lake stage for a week, and the other is the desired low lake stage for each week. It looks like this: week High Low 1 16 14.5 2 16 14.5 3 15.95 14.45 4 15.84 14.34 5 15.73 14.23 6 15.61 14.11 7 15.5 14 8 15.38 13.88 9 15.25 13.75 10 15.13 13.63 .... 52 15 13.5 I can plot out two nice curves

I issued the following command to obtain the std dev for each month. psd<-numSummary(Sal, groups=month, statistics=c("sd")) which resulted in > psd sd n NA 1 6.930340 9367 2319 2 7.847003 10827 1008 3 5.962308 12988 404 4

I have a file, each column of which is a separate year, and each row of each column is mean precipitation for that month. Looks like this (except it goes back to 1964). month X2000 X2001 X2002 X2003 X2004 X2005 X2006 X2007 X2008 X2009 1 1.600 1.010 4.320 2.110 0.925 3.275 3.460 0.675 1.315 2.920 2 2.960 3.905 3.230 2.380 2.720 1.880 2.430 1.380

My x axis is a series of daily dates (e.g., 01/01/2000, 01/02/2000, etc.) from 2000 to end of 2008. The default only gives me 4 ticks. I want more. Why doesn't this work? sdate<-as.POSIXct(strptime(date,format="%m/%d/%Y")) plot(ppt~sdate,type="l",ylim=c(0,47),col=1,lwd=1,pch=16,ylab="Salinity, psu",xlab="Year",las=1,main="Duck Key

What am I failing to understand here? The script below works fine if the dataset being used is DNase1 <- DNase[ DNase$Run == 1, ] per the example given in help(nlrob). Obviously, I am trying to understand how to use nls and nlrob to fit curves to data using R. #package=DAAG attach(codling) plot(pobs~dose) #next command returns 'step factor reduced below min factor

If you submit these lines, you end up with variable "vare.dis". I want to export vare.dis to csv. Stuck I am. library(vegan,logical.return = TRUE) #return=true verifies package is available library(MASS,logical.return=TRUE) #return=true verifies package is available data(varespec) #varespec is an example data file in the vegan package vare.dis <- vegdist(varespec)

On Mon, Apr 13, 2009 at 12:41 PM, Dan Dube <ddube-at-advisen.com> wrote: > i use tapply and by often, but i always end up banging my head against > the wall with the output. The proposed solution of Dan's problem posted on R-help was: > do.call(rbind,a) When I use this 'solution' I get 'ERROR: second argument must be a list'. So head on wall continues. My

I have 3 columns: flow, month, and monthname, where month is 1-12, and monthname is name of month. I can't get the plot to replace the 1-12 with monthname using ticks.lab. What am I doing wrong? plot(flow~factor(month),xlab="Month",ylab="Total Flow per Month", ylim=c(0,55000), ticks.lab="monthname") Thanks Gregory A. Graves Lead Scientist REstoration

Newbie here~ I've spent a ton of time on this, but have to admit I am stuck. I've tried various combinations of strptime and the package CHRON, but still can't get there. My simple task is that I want to plot salinity (S1 or S2) on y-axis and time on x-axis, but I want to use a lay-understandable representation of time (not Julian). If I use the following command, I can make a

Among those users of Primer, stress values greater than 0.3 are interpreted as "questionable". Using both isoMDS and metaMDS (vegan package), the stress values returned are much higher using my own data and using examples provided in R Help. For example Rstress = 8.3, and the stressplot r2 = 0.99 indicating (to me) that the ordination is OK. I am guessing that the "stress"

Hi !            I am new to R. Can somebody help me in reformatting huge output files ,i.e, rearranging sets of columns in specific order. For example: I have data for three compunds 1, 2 and 3 file1: ID CA1 CA3 CA2 MA2 MA1 MA3 1 14 15 13 7 12 3 2 19 7 12 10 14 5 3 21 12 19 6 8 9   to File 2:   ID CA1 CA2 CA3 MA1 MA2 MA3 1 14 13 15 12 7 3 2 19 12 7 14 10 5 3 21 19 12 8 6 9   or File3: ID

Hi All I have been struggling with this model for some time now and I just can't get it to work correctly. The messages I get when running the code is: DLSODA- Warning..Internal T (=R1) and H (=R2) are such that in the machine, T + H = T on the next step (H = step size). Solver will continue anyway. In above message, R = [1] 0 0 DINTDY- T (=R1) illegal In above message, R = [1]

Hello, I'm developing an ordination using metaMDS (package vegan). The analysis identifies 3 distinct groups that I'd like to define by either adding circles or ellipses to help identify the groups. The data set is a spatial temporal data base depicting change in each of 4 areas over 3 time periods. I can add lines that would link each of the time - x - area trajectories, but

Hello everyone, I'm working with R 2.4.1 on a PC running with XP. Trying to run isoMDS as follows: Gquad.mat <- Gquads[4:10] # extracts only the metric data variables Gquad.dist <- dist(Gquad.mat) Gquad.mds <- isoMDS(Gquad.dist) plot(Gquad.mds$points, type = "n") text(Gquad.mds$points, labels = as.character(Gquads$Quadname), cex = .5) The data.frame

Hello, I'm wondering if there are any packages/functions that can perform conditional sequential gaussian simulation. I'm following an article written by Grunwald, Reddy, Prenger and Fisher 2007. Modeling of the spatial variability of biogeochemical soil properties in a freshwater ecosystem. Ecological Modelling 201: 521 - 535, and would like to explore this methodology.

Currently, I have two networks configured via xml: virsh # net-dumpxml TestNetwork1 <network connections='2'> <name>TestNetwork1</name> <uuid>a76f665a-0196-4edb-81b4-340944a6869c</uuid> <forward dev='p1p1' mode='nat'> <nat> <port start='1024' end='65535'/> </nat> <interface

Hello. I have a question about libvirt remote access. I read carefully how to set up TLS authorization guide and I have done all the steps. Everything works, but for my virtualization infrastructure I need more flexibility. How can I configure a single host to trust more than one certificate authority. Lets say, for example, that I have 10 virtualization hosts and each of them has certificates

I have the following appended to gluster logs at around 100kB of logs per second, on all 10 gluster servers: [2012-06-11 15:08:15.729429] I [dht-layout.c:682:dht_layout_dir_mismatch] 0-sites-dht: subvol: sites-client-41; inode layout - 966367638 - 1002159031; disk layout - 930576244 - 966367637 [2012-06-11 15:08:15.729465] I [dht-common.c:525:dht_revalidate_cbk] 0-sites-dht: mismatching layouts

Hi All Apologies if this is obvious, but I''m a bit flaky around SSL certificates. NB puppet version 0.25.5 We use the brilliant feature of certificates where you can have Alternate DNS names for a certificate which is manifested in the puppet master config file as certdnsnames. All our clients connect to puppet-$ location.example.com, and if $location is down, we can point the CNAME to