Re: [R] plot problem

On Mon, 21 Aug 2006, Daniil Ivanov wrote: > Ok, what is wrong with a following code: > > # remove all the present objects > rm(list = ls()) > > # load the libraries we need > library(gstat) > > data(meuse) > vgm1 <- variogram(log(zinc)~1, ~x+y, meuse) > plot(vgm1) > dev.copy2eps(file="fig2.eps"

[R] Clique technique-Package

Dear R Users, I am looking clique a graph technique, to identify the values from the dataset. I tried with help.search("graph") it show's the graphics related stuff Is there any package that i can use to find Clique in a dataset? Thanks in Advance JJ -- Lecturer J. Joshua Thomas

Re: [R] Clique-Method-Package-Help

Dear R Users, I am looking clique a graph technique, to identify the values from the dataset. I tried with help.search("graph") it show's the graphics related stuff Is there any package that i can use to find Clique in a dataset? Thanks in Advance JJ [[alternative HTML

Re: [R] Grid Points

That's the default. See the relation subargument to scales if you want them different. e.g. library(lattice) y <- c(601:700, seq(6510,7000, by=10)) x <- c(601:700, 601:650) g <- rep(1:2, c(100, 50)) xyplot(y ~ x | g) On 8/20/06, Anupam Tyagi <[EMAIL PROTECTED]> wrote: > Thanks. How do I retain

Re: [R] Permutations with replacement

Thanks, David. That worked fabulously! Here is the R code for the hypercube test example: ## begin R code library(combinat) x <- rep(3,3) # for partitions of 3 units into the three classes {1,2,3} hcube(x, scale=1, transl=0) ### end R code For th

Re: [R] plot problem

Ok, what is wrong with a following code: # remove all the present objects rm(list = ls()) # load the libraries we need library(gstat) data(meuse) vgm1 <- variogram(log(zinc)~1, ~x+y, meuse) plot(vgm1) dev.copy2eps(file="fig2.eps",horizontal=T) dev.off() it plots nothing but from the R console

Re: [R] Variance Components in R

Dear Harold and others, I have changed the syntax for lmer() and used this one: require(lme4) gt <- read.table("gt5.txt") sink("GT output.txt") attach(gt) system.time(fm <- lmer(RATING ~ 1 +(1|CHAIN) +(1|SECTOR) +(1|RESP) +(1|ASPECT) +(1|ITEM) +(1|SECTOR*RESP) +(1|SECTOR*ASPECT) +(1|SECTOR*ITEM) +

Re: [R] plot problem

Hi, Ok, thanks to all. Problem was with class of variogram > class(vgm1) [1] "gstatVariogram" "data.frame" If I fix it manually to class(vgm1) <- "gstatVariogram" everything runs as it should. Thanks, Daniil. On 8/21/06, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > On Sun, 20 Aug 2006, Dan

Re: [R] unquoting

Thank you, Brian, Peter and Gabor Brian has what want. My heading was a bit misleading. I was looking for a function that would, in logicians' terms, convert 'mention' into 'use'. (This usually goes along with the story about the importance of knowing the difference between a lion and "lion".)

Re: [R] plot problem

On Sun, 20 Aug 2006, Daniil Ivanov wrote: > Hello. > > I'm pretty much new to R and I'm trying to produce some figures. What have you been reading to get the ideas below? People new to R do not tend to use dev.next ... indeed experienced users very rarely use it. > It seems to me, that R has

[R] plot problem

Hello. I'm pretty much new to R and I'm trying to produce some figures. It seems to me, that R has some asynchronous way of plotting figures. When I run this code: #constructs the semivariogram of SC1929 vgm1 <- variogram(SC1929~1,~U+V,puerto.map$att.data) # trying to make new plot dev.set(which

Re: [R] C compile problem on Ubuntu linux

Manuel Castejón Limas wrote: > Hello, > I've just compiled Hmisc ok under dapper. > I think you need to further install some packages. > Have you installed libc6-dev? > I would start installing the build-essential package. > Best wishes > Manuel Thanks Manuel, apt-get install build-essential solve

[R] C compile problem on Ubuntu linux

Under Ubuntu dapper, after installing packages gcc and g77, under platform i486-pc-linux-gnu arch i486 os linux-gnu system i486, linux-gnu status major2 minor2.1 year 2005 month12 day 20 svn rev 36812 language R I get an error when trying to update.packages('Hmis

Re: [R] Simulate p-value in lme4

You've raised a very interesting question about testing a fixed-effect factor with more than 2 levels using Monte Carlo. Like you, I don't know how to use 'mcmcsamp' to refine the naive approximation. If we are lucky, someone else might comment on this for us. Beyond this, you are

Re: [R] Grid Points

Thanks. How do I retain the same scale of grid.points from one panel to next even if the scale of the data changes? For example: c(seq(601:700),seq(6510,7000, by=10)) ~ seq(601:700) | gl(2,50). --- Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > Try this. gl(2,50) is such that the first 50 poi

Re: [R] fit the series data

try this: > x<-c(1,2,5,6) > y<-c(5,3,6,8) > xy<-data.frame(x,y) > xy x y 1 1 5 2 2 3 3 5 6 4 6 8 > new.df <- data.frame(x=seq(max(xy$x)), y=rep(0, max(xy$x))) > new.df x y 1 1 0 2 2 0 3 3 0 4 4 0 5 5 0 6 6 0 > new.df$y[xy$x] <- xy$y > new.df x y 1 1 5 2 2 3 3 3 0 4 4 0 5 5 6 6 6 8 > On 8/

Re: [R] string-to-number

On Sat, 2006-08-19 at 10:25 -0600, Mike Nielsen wrote: > Wow. New respect for parse/eval. > > Do you think this is a special case of a more general principle? I > suppose the cost is memory, but from time to time a speedup like this > would be very beneficial. > > Any hints about how R programm

[R] fit the series data

Dear friends, suppose my dataset *xy* : xy 1 5 2 3 5 6 6 8 -generated the data-- x<-c(1,2,5,6) y<-c(5,3,6,8) xy<-data.frame(x,y) --- I want to fit the gap in x with the corresponding y=0, I use the following programs to generate

Re: [R] split a y-axis to show data on different scales

Look at oma= and mar= parameters to par for controlling the space when using mfrow=. e.g. opar <- par(oma = c(6, 0, 5, 0), mar = c(0, 5.1, 0, 2.1), mfrow = c(2,2)) for(i in 1:4) plot(1:10) par(opar) On 8/20/06, Anupam Tyagi <[EMAIL PROTECTED]> wrote: > I think information can be enhanced by usin

Re: [R] Grid Points

Try this. gl(2,50) is such that the first 50 points are series 1 and the second 50 points are series 2. The scales= argument defines the positions of the tick marks and the xlim= argument defines the x axis limits. The layout puts the panels on top of each other rather than side by side. strip

Re: [R] Variance Components in R

Harold, I have tried the following syntax: > fm <- lmer(RATING ~ CHAIN*SECTOR*RESP +(1|CHAIN*SECTOR*RESP), gt) > summary(fm) Linear mixed-effects model fit by REML Formula: RATING ~ CHAIN * SECTOR * RESP + (1 | CHAIN * SECTOR * RESP) Data: gt AIC BIClogLik MLdeviance REMLdevianc

Re: [R] split a y-axis to show data on different scales

I think information can be enhanced by using different scaled graphs next to each other. mfrow() created too much space, there may be no need to again draw the x-axis. It can be very useful to have different scales of the same data presented next to each other, in addition to the main graph. So I t

Re: [R] Variance Components in R

Harold, I have tried to adapt your syntax and got some problems. Some responses from lmer: On this one, I have tried to use "1" as a grouping variable. As I understood from Bates (2005), grouping variables are like nested design, which is not the case. > fm <- lmer(RATING ~ CHAIN*SECTOR*RESP +(CHA

[R] Grid Points

How do I put grid points (not grid lines) as the base layer of an xyplot? Is there a way to vary the interval at which x and y grid points are placed? Is it possible to start a graph so that Y axis begins at 500 and ends at 800? I am only interested in focusing on the relative distance between t

Re: [R] Variance Components in R

Reading Bates' article on R News, I see that random effects require a grouping variable. As, by convention, all variables in G-studies are supposed random, what could be a grouping variable in that case? I see that the model I wrote before (if ever ran...) would take all effects as fixed. Is it po

Re: [R] unquoting

?get I really think this has nothing to do with `quoting', rather to do with evaluating variables from their names. At first I though you were looking for noquote(), which does unquote in the conventional sense. > noquote(names(AF)[2]) [1] Second > get(names(AF)[2]) [1] 3 4 On Sun, 20 Aug 2006

Re: [R] unquoting

Murray Jorgensen <[EMAIL PROTECTED]> writes: > I would like a function to strip quotes off character strings. I should > work like this: > > > A <- matrix(1:6, nrow = 2, ncol=3) > > AF <- as.data.frame(A) > > names(AF) <- c("First","Second","Third") > > AF >First Second Third > 1 1

Re: [R] unquoting

Try these get(names(AF)[2]) AF["Second"] # this one different than the rest AF[["Second"]] AF[, "Second"] AF$Second On 8/20/06, Murray Jorgensen <[EMAIL PROTECTED]> wrote: > I would like a function to strip quotes off character strings. I should > work like this: > > > A <- matri

[R] unquoting

I would like a function to strip quotes off character strings. I should work like this: > A <- matrix(1:6, nrow = 2, ncol=3) > AF <- as.data.frame(A) > names(AF) <- c("First","Second","Third") > AF First Second Third 1 1 3 5 2 2 4 6 > names(AF)[2] [1] "Second" >

Re: [R] how to the p-values or t-values from the lm's results

try the following: data <- data.frame(matrix(rnorm(900), ncol = 9)) names(data) <- c("y", paste("x", 1:8, sep = "")) logr <- lm(y ~ . - 1, data) a <- summary(logr) coef(a) coef(a)[, 3:4] coef(a)[, "t value"] coef(a)[, "Pr(>|t|)"] I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D.

Re: [R] Query: how to modify the plot of acf

Try constructing the acf plot using the traditional plot tools. Then you can do what you like with it. Eg if your model is called model.lme, then something like this should work: acf.resid <- ACF(model.lme, resType = "n") my.lags <- acf.resid$lag > 0.5 plot(acf.resid$lag[my.lags], acf.resid$ACF[

Re: [R] how to the p-values or t-values from the lm's results

Hi there Zhang, While there might be a better way... an ugly but generic way of accessing this type of information is to use str() and a little experimentation... here is a little history() of what I did to find it... a str(a) str(logr) a[[1]] a[[2]] a[[3]] a[[4]] a[[4]][[1]] a[[4]][1,] a[[4]][,4]