similar to: control data symbols in pairs function

Is there a "canned" function in R for finding the standard error of the mean? I have tried > sem <- function(x) c(mean =mean(x), + SEM = stdev(x)/sqrt(length(x))) > sem(pnet.lai) Error in sem(pnet.lai) : couldn't find function "stdev" It looks like there is no stdev function in R Thanks, Kirk Kirk R. Wythers email: kwythers at umn.edu University of

I just compiled R on my freebsd box. I don't see anything when firing it up with -g or --gui or --gui=X11, etc... only the standard > prompt (see below). Is that what I should expect? I'm running the kde desktop if that matters. Any suggestions from those who use R on freebsd? Thanks Kirk bash-2.05a$ R R : Copyright 2001, The R Development Core Team Version 1.3.1 (2001-08-31)

I am struggling to get JMP installed. When I try launching the installer with "wine /mn/cdrom/Launch.exe", I get the JMP splash screen. When I use the mouse to click "install" on the splash screen, I get a dialog box with the error "Setup failed to launch installation engine: (0x80070005). The console output reads as follows: [kwythers@lorax kwythers]$ wine

I am trying to load the GRASS library and am getting the following error: > library(GRASS) Loading required package: akima Error: package 'akima' could not be loaded In addition: Warning message: There is no package called 'akima' in: library(pkg, character.only = TRUE, logical = TRUE, lib.loc = lib.loc) This is on R 2.0.0, library(GRASS) is installed, and update.packages()

I have the following problem. I have some multidimensional data points "x" and a curve "fit" fitted to these points. How can I combine R> pairs(x) and R> pairs(fit,panel=lines) in one plot, so that I can see how good the curve fits the data? A command like R> pairs(x, panel=function(x,y) {points(x,y); lines(lowess(x,y))}) does not work, since I fit the curve in all

I have a set of modeled climate data recorded at irregular intervals. The format of the data is such that there are monthly measurements for the years 2000, 2020, 2050, 2080, etc. Therefore I have 12 regular records, a skip of some number of years, then 12 more monthly records, another skip, and so on.... I created a dataframe from .txt with the read.table() command. For starters I need

Dear R users, subject: How to filter XY pairs of inacurate gps position along track, taking into account the time index to not mix track from different days in one average track or subject as: How to filter XY pairs of inacurate gps position: (latitude, longtitude) to get estimated correct position before calculating track i.e. How to estimate (time series? )? How to filter XY pairs of

R-helpers: Is there a simple way to extract every nth element from a very long vector? [The vector I want to sample from is an object returned by lowess(). It is so long (11,628 pairs of elements) that it is causing non-R-related memory problems elsewhere. I don't see a better way other than sampling the output returned from lowess.] Thanks in advance. Barry Cooke

Hello, can anyome tell me how to access the full script of the panel.smooth function so that I can change the thickness of the smoothing line or its colour ? All I could access is : > panel.smooth function (x, y, col = par("col"), bg = NA, pch = par("pch"), cex = 1, col.smooth = "red", span = 2/3, iter = 3, ...) { points(x, y, pch = pch, col = col, bg =

Hello Masters, I run the loess() function to obtain local weighted regressions, given lowess() can't handle NAs, but I don't improve significantly my situation......, actually loess() performance leave me much puzzled.... I attach my easy experiment below #------SCRIPT---------------------------------------------- #I explore the functionalities of lowess() & loess() #because I have

I am trying to change the attributes of the lowess lines fit to an xyplot command, but have been unsuccessful in my search of the online help. Right now, both the points and lowess line come out in the same color (blue). I am unsure how I can change the properties of the lowess line separately. xyplot ( SnowLineElevation ~ Year | Model, data = data, ylim = c(0,1800), type =

I need some help understanding how on of the example data sets is formatted in the basic R installation. If I load the Mona Loa CO2 data, with the command: > data(co2) I can view the data with: > co2 And the data are in the form of 11 rows labeled as years (1994-2004) and 12 columns labeled (Jan - Dec). This structure appears to be a dataframe, however, if I type the command

Hi all, I am new in R tool. I would like to plot two dataset in in fig. Here is what I did for both a and b data sets jpeg(file="a.jpeg") dat<-read.table('a', header=F, sep=',') dim(dat) y<-dat[,1] y<-y[!is.na(y)] plot(y);lines(lowess(y, f=0.05), col = ("red"), lwd=5) dev.off Two questions: 1. How I can save this lowess smooth data? 2. Once I

As I've mentioned a number of times. I find it very useful to have lowess() become a generic function so that a lowess.formula() can be defined. Below is a patch that makes both changes, as well as updating the corresponding help documentation. Gregory R. Warnes Manager, Non-Clinical Statistics Pfizer Global Research and Development Tel: 860-715-3536 ? DESCRIPTION ? Makefile ?

Hi, I have a scatter plot of the variables GNI and Lifeexp (Gross National Income and Life Expectancy, both metric). So I plotted them and I want to add a regression line and a lowess line. I use lowess and not loess because I have missing values. My code: plot(GNI,Lifeexp) abline(lm(Lifeexp~GNI), col="red") y.loess<-loess(Lifeexp~GNI,na. action = na.exclude)

I'm using this code to plot a smoothed line. These two columns of data really represent 4 groups and I'd like to plot a separate line for each group but have them all in the same plot. The R-Docs for lowess do not seem to indicate some type of "GROUPS=var_name" option. What would be the syntax for this? plot(AWGT ~ lipid ) lines(lowess(lipid , AWGT, f=.8)) -- Dean

The lowess() help page refers to documentation in "src/appl/lowess.doc". This was moved to "src/library/stats/src/lowess.doc" in 2007. This patch fixes it: Index: src/library/stats/man/lowess.Rd =================================================================== --- src/library/stats/man/lowess.Rd (revision 79137) +++ src/library/stats/man/lowess.Rd (working copy) @@

Hi, Are there any functions in R that could be used to estimate the phase-shift between two semi-sinusoidal vectors? Here is what I have tried so far, using the spectrum() function -- possibly incorrectly: # generate some fake data, normalized to unit circle x <- jitter(seq(-2*pi, 2*pi, by=0.1), amount=pi/8) # functions defining two out-of-phase phenomena f1 <- function(x)

Hi, A lowess method is defined in our package for one of our S4 classes. To explicitely define the generic this is being used if (!isGeneric("lowess")) setGeneric("lowess", useAsDefault = lowess) This works fine for many other methods inherited from various R packages. In this case a warning is issued in R CMD check: * checking R code for possible problems ... NOTE

Dear all, To draw a lowess line on a plot was a piece of cake; to draw a loess line, however, seems not that easy. Is the loess plotting implemented at all in relation to the loess function, or do I have to look in add-on packages? Thanks, Marcin