similar to: Bug in gap.plot

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

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

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

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

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

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

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

Hi, all, I need to show only values in a barplot from say 4-12. When I do this, values below 4 are not cutoff. I would appreciate any input into how I could cut the bottom part of the bar off. Best wishes, Art Roberts University of Washington Department of Medicinal Chemistry

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

Hi! When running the following code I get the message: Error in eval(expr, envir, enclos) : object "A" not found fm <- function(p,ydata, env = parent.frame()) { #fit model y ~ (b0+b1x1+...+bpxp)*exp(g2plus*z2plus) where bi and g2plus are parameters #create design matrix for linear part fo00 <-

Hi, all, How do I get an inset graph (i.e. graph within a graph)? Your input is greatly appreciated. Best wishes, Art University of Washington Department of Medicinal Chemistry

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 My question is simple - the gap.plot function in the plotrix package allows users to draw graphs that have a broken axis. However, I want to then add a line to the "second" plot, but can't. Eg: twogrp<-c(rnorm(10)+4,rnorm(10)+20) gap.plot(twogrp,rnorm(20),gap.bounds=c(8,16),gap.axis="x",xlab="X values", xtics=c(4,7,17,20),ylab="Y

Hello! Using the package Plotrix I want to do a plot with a broken axis. So far it's working fine but now I want only the x and y axis plotted (the x-axis with a gap in it), but not the axes that are reffered to as axis 3 and 4, and not the lines that additionally mark the gap in the axis. I thought an 'axes=F' in the gap.plot() command would do (as in the usual plot()), but it does

I am trying to make a barplot with a broken axis using gap.barplot (in the indispensable plotrix package). This works well when the data is a vector: > twogrp<-c(rnorm(10)+4,rnorm(10)+20) > gap.barplot(twogrp,gap=c(8,16),xlab="Index",ytics=c(3,6,17,20),ylab="Group values",main="Barplot with gap") But when the data is an array (for a bar plot with multiple

Hi everyone. I'm trying to break the y axis on a plot. For instance, I have 2 series (points and a loess). Since the loess is a "continuous" set of points, it passes in the break section. However, with gap.plot I cant plot the loess because of this (I got the message "some values of y will not be displayed"). Here's my code: library(plotrix); #generate some data x

To Whom It May Concern: I need to use various packages in R and I also want to use Textmate. Problem 1: Textmate doesn't seem to recognize the packages that are already installed on R. Problem 2: When I execute the following command: install.packages("ade4", repos="http://cran.r-project.org", contriburl = contrib.url(

Hello, I try to handle a simple bar-plot, but it turns out to be not as simple as I thought. 1) I have created a .dat-File, e.g. test.dat: DATA DATA-SEM 2.2 0.32 6.2 1.30 12.7 1.61 48.6 3.08 4.1 0.86 4.5 0.32 1.5 1.13 1.2 1.08 The first row is the data represented by bars. The second row deals with the Standard Error of Mean. The lines correspond to time-intervals of experiments. 2) I now