similar to: object not found inside function

hi all i have a function that uses two inputs, say xdata and ydata. An example is the following, simple1<-function(xdata,ydata) { ofit<-lm(ydata~xdata) list(ofit) } say i use arbitray number for xdata and ydata such that D = x1 x2 y 1 1 10 2 6 6 3 10 7 x<-D[,1:2] and y<-D[,3] if one uses these inputs and rund the program we get the following: >simple(xdata=x,ydata=y)

Hi there, I don't have much experience with fitting at all and I'd like to get some advice how to use the "weights"-argument with nls correctly. I have created some data with a sigmoidal curve shape. Each y-Value was generated by the mean of three values. A standard deviation was calculated too. Now, I'd like to weight the data points respective to its standard

Dear R-users and experts I want to create to analyse my data which looks like follows:( I have show only 8 variables but original variables much more number >1000) *sub* *ydata* *X1* *X2* *X3* *X4* *X5* *X6* *X7* *X8* 1 12 1 1 1 2 1 1 1 1 2 13 2 2 1 2 2 1 1 1 3 11 1 1 1 2 1 2 1 2 4 12 1 1 2 1 1 2 2 2 5 14 1 2 2 1 1 2 2 2 6 12 2

Hello all, I have been using R for about 3 weeks and I am frustrated by a problem. I have read R in a nutshell, scoured the internet for help but I either am not understanding examples or am missing something completely basic. Here is the problem: I want to plot data that contains dates on the x axis. Then I want to fit a line to the data. I have been unable to do it. This is an example of

I have a set of data in the form x1, y1, z1 x1, y2, z2 ... x1, yN, zN x2, y1, z(N+1) x2, y2, z(N+2) ... x2, yN, z(2N) ...and so on... xM, yN, val(M*N) I have been trying to figure out how to get R to use this data in a persp plot. So far the only thing that I can figure out to do is to break the data file into three different files. The first file contains the x-coordinate data: x1 x2

Airon, I used R2.4.0 on a Windows XP (SP2) (not Chinese :-)) and it still works: > data = read.csv("u:/airon.csv") > xdata = data$Adj..Close > modwt.la8 = modwt(xdata, "la8", n.level=6) > summary(modwt.la8) Length Class Mode d1 1467 -none- numeric d2 1467 -none- numeric d3 1467 -none- numeric d4 1467 -none- numeric d5 1467 -none- numeric d6 1467

Good afternoon My name is Joanne Lenehan, I am a post grad at UNE using R version 2.9.0 I came across the Tweedie package in old R help posts and was interested in giving it a go for some data. The data is below and also attached as BaregroundLitterLogs Site Treatment Graze Dam Plot Time Bare Litter Logs 1 C remote yes 1 A 0 2 0 1 C remote yes 2 A 0 15 0 1 HE

Dear all: I have the following codes: Xdata<-c(2,3,8,9,10) Ydata<-1:5 gap.plot(Xdata, Ydata,gap=c(5,6),gap.axis="x",type="o") However, the type='o' seems only work on the first part of gap plot, the second half of the plot always just points, you can not add lines on that part, any help will be highly appreciated. I would like to have these two parts of

Hi: I have encountered problems with imodwt and universal.thresh.modwt and cannot find any reference in R Search. Hope someone can give me some ideas: Starting with modwt.la8 <- modwt(xdata, "la8", n.level=6) <-- this seems to work fine (1) ydata <- imodwt(modwt.la8) will always give ydata as numeric(0) (no values) instead of being a time series data with

Dear list members. I have problems with the usage of cv.glm from the boot package. Here are some parts of the script I wanted to use: data <- read.table("selected_2D.csv", header=TRUE, sep=",") ? glm.fitted <- glm("ydata$ y ~ 1 + density + vsurf_ID6 + vsurf_S ", data=data) error <- cv.glm(data=data, glm.fitted, K=6) ydata$y is a separate data set, where

Hi, all I am trying to make a plot with a axis break and I want the whole plot to be line, not points. However, when I execute the following command half of the graph is points and the other lines. gap.plot(Xdata, Ydata,gap=c(5,6),gap.axis="x",type="l") I think it might be a bug in plotrix. I would greatly appreciate your input. If there is another way to do it, I

Airon, I don't think you have to find an English computer 'cause the following must work in your Chinese one :-) Let me explain. First of all, change your lines to xdata <- ckhdat$Adj..Close[1:1447] #names(ckhdwt.la8) <- c("w1", "w2", "w3", "w4", "w5","w6", "v6") note the # sign, i.e., DO NOT change the names

Greg, Two thoughts: 1. It might be possible that 'vars' is a reserved word of sorts and if you change the name of your vector RF might be happier 2. A way that works for me is to call importance as follows: sel.imp <- importance(sel.rf, class=NULL, scale=TRUE, type=NULL) and then use the 'names' of the imp data frame to be absolutely clear to RF you are talking about the

Hello to everyone, I am trying to compute a mean and show it to the screen and then save to later be used as a boxplot.I have used the following code: dat<-read.table(file="C:\\Documents and Settings\\Owner\\My Documents\\Yeast\\Yeast.txt",header=T,row.names=1) file.show(file="C:\\Documents and Settings\\Owner\\My Documents\\Yeast\\Yeast.txt") x<-dat[2,23:46]

I've spent the weekend or so trying to get my head around Theora and I'm just not having any luck. First of all, I really don't have any interest in immersing myself in video encoding--I just have one simple thing I'd like to do--so I haven't spent any time reading about the theory or specs or anything like that. I'm attempting to write a "render-to-theora"

I've been playing around with some ways to tighted up our AVX specification and have hit upon a nice way to reduce a bunch of code. Unfortunately, right now TableGen can't handle it. Here's a simple example of what I want to do: class Data<string n, int v> { string Name = n; int Value = v; } // Define some objects usable as arguments. def X : Data<"X", 1>;

Hi all, I am trying to print some generated text on plot which is comprised of text and the value of a variable. Thanks to the response from Joshua Wiley <jwiley.psych@gmail.com> on the thread "How to format R superscript 2 followed by "=" value, I can now get R^2 = value. However with the following code: ggplot(blData, aes(x = xData, y = yData)) + geom_point(aes(colour

----------------------- *** Disclaimer *** ----------------------- This e-mail and its contents are subject to the SA Reserve Bank's Disclaimer and Confidentiality Clause, which can be viewed at: http://www.reservebank.co.za/disclaimer Should you be unable to access the link provided, please send a blank e-mail to Disclaimer@resbank.co.za ----------------------- *** Disclaimer ***

Hello, I am building a package to read in a pretty common atmospheric data format and create a SKEW-T, log p diagram (also very common for the atmospheric community). The problem is this: How should one address the (ASCII) data files so that the examples find them? I have read through the "Writing R Extensions" document and still have a problem. There seems to be some sort of

Hi, P-values in Durbin-Watson test obtained through the use of functions available in packages "lmtest" and "car" are different. The difference is quite significant. function "dwtest" in "lmtest" is much faster than "burbinwatson" in "car". Actually, you can take a nap while the latter trying to calculated Durbin-Watson test. My question