Re: [R] Course of R in two days

Who is your audience, please? On Fri, Oct 10, 2008 at 11:27 PM, Lord Tyranus <[EMAIL PROTECTED]> wrote: > Hi > > I have to teach a course of R in two day , I need some an advice about what > must I teach? > > Thanks in advance > >[[alternative HTML version deleted]] > > _

[R] Course of R in two days

Hi I have to teach a course of R in two day , I need some an advice about what must I teach? Thanks in advance [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

Re: [R] Help MLE

Thanks again. Now I want to complicate the model. Supose that Y is a vector with 30 data. Y = c(0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0) Then the r = sum(Y*log(comb)) Is it possible? I want to discover the variables values (alpha12,w_g12,w_u12) f

[R] how to evaluate a cubic Bezier curve (B-spline?) given the four control points

I'm trying to use R to determine the quality of a cubic Bezier curve approximation of an elliptical arc. I know the four control points and I want to compute (x,y) coordinates of many points on the curve. I can't find anything in either the base distribution or CRAN that does this; all the spline-

Re: [R] a question about diparate device quality

Check out: http://www.lanceknadle.com/images/raster_vector.pdf On Fri, Oct 10, 2008 at 12:27 PM, Kyle Matoba <[EMAIL PROTECTED]> wrote: > List, > > I have a question about why plots sent to certain devices can be of so much > worse quality than others. Whenever I plot to x11 or pdf the graphs are

[R] TTR getYahooData?

Would it be possible to get data from uk.finance.yahoo.com instead from finance.yahoo.com. Thanks __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.

Re: [R] Creating GUIs for R

Hi Roger, Bos, Roger wrote: I haven't heard Rpad mentioned yet, so I will mention it. Rpad allows you to run R code inside a browser window. You can many any kind of GUI you want using html forms and then call the R function from the html page. The Apache/Rpad install works on any platform

Re: [R] Altering the cube around a wireframe plot.

On Thu, Oct 9, 2008 at 8:59 AM, Todd Remund <[EMAIL PROTECTED]> wrote: > > I'm trying to create a 3D plot using wireframe with certain parts removed. I > would like to get rid of the part of the outer cube that crosses over the > plot leaving the back two walls and the axes. It would also be us

Re: [R] rgl-snapshot failed (err-msg: "failed")

On 10/10/2008 8:55 AM, Oliver Bandel wrote: Zitat von Duncan Murdoch <[EMAIL PROTECTED]>: On 10/10/2008 8:13 AM, Oliver Bandel wrote: Hello, I tried to use rgl.snapshot and it failed. The error message was not very verbose: == plot3d( motion[[idx+2]], motion[[idx+1]]

[R] bivariate non-parametric smoothing

Hi, I was wondering if there is a function in R which performs bivariate non parametric smoothing which allows for the possibility of including some weights in the smoothing (for each data points in my grid I have some predefined weights that I would like to include in the smoothing). Tha

Re: [R] Creating GUIs for R

Thanks, Roger, It seems simply great to me. I am far from being able to detect the drawbacks of the system, if any. I don't understand why this kind of programs are not of general use... With my best regards, Ricardo Bos, Roger wrote: > I haven't heard Rpad mentioned yet, so I will mention it

Re: [R] plotCI

On 11/10/2008, at 10:15 AM, Caio Azevedo wrote: Hi all, I am using the function "plotCI" with the following command: plotCI(m.residuos.p.2 [1:41],li=m.residuos.p.3 [1:41],ui=m.residuos.p.4 [1:41],lty=1,ylab="") This generates exactly what I want except for the fact that I wanna drawn a

Re: [R] Sampling

On 11/10/2008, at 10:26 AM, Robert W. Baer, Ph.D. wrote: How about: T2=rnorm(1000) temp=list(NULL) for(i in 1:5){temp[[i]]<-sample(T2,40,replace=F)};show(temp) Obviously that *works* and may well be (something like) what ``Alex99'' *should* have done, but what did he actually *do*? T

Re: [R] Sampling

How about: T2=rnorm(1000) temp=list(NULL) for(i in 1:5){temp[[i]]<-sample(T2,40,replace=F)};show(temp) - Original Message - From: "Rolf Turner" <[EMAIL PROTECTED]> To: "Alex99" <[EMAIL PROTECTED]> Cc: Sent: Friday, October 10, 2008 3:36 PM Subject: Re: [R] Sampling On 11/10/2008

[R] plotCI

Hi all, I am using the function "plotCI" with the following command: plotCI(m.residuos.p.2 [1:41],li=m.residuos.p.3 [1:41],ui=m.residuos.p.4 [1:41],lty=1,ylab="") This generates exactly what I want except for the fact that I wanna drawn a line linking the points (m.residuos). How could I do tha

Re: [R] follow-up on Error when reading a SAS transport file (with sasxport.get from Hmisc)

Peter Dalgaard wrote: Jean-Louis Abitbol wrote: I have done what P. Dalgaard has suggested and I don't find a descrepancy between the number of values and the number of labels: there 15 each... Any hint on what might go wrong here ? Actually, I think you got it: > factor(1,c(NA,1:4),c(1:

[R] predicting from a local regression and plotting in lattice

Hi R community, I'm running R 2.7.2 on Windows XP SP2. I'm trying to (1) plot loess lines for each of my groupings using the same color for each group; (2) plot loess predicted values. The first part is easy: data1 <- data.frame(Names=c(rep("Jon",9),rep("Karl",9)),Measurements=c(2,4,16,25,36,4

Re: [R] Sampling

On 11/10/2008, at 7:38 AM, Alex99 wrote: Hi Guys, I am in desperate need with sampling. I am suppose to sample from a dataset. I use the following code: for(i in 1:5){temp[i]<-sample(T2,40,replace=F) + show(temp)} but all the samples are the same, but I want them to be different.any sug

Re: [R] follow-up on Error when reading a SAS transport file (with sasxport.get from Hmisc)

You appear to have a label for an NA value, and that's the problem. I think you need to send an example to the package maintainers (see the posting guide) to help them debug what looks like a bug or undocumented restriction in their code (and it seems this is the same codebase in the two packa

Re: [R] follow-up on Error when reading a SAS transport file (with sasxport.get from Hmisc)

Jean-Louis Abitbol wrote: I have done what P. Dalgaard has suggested and I don't find a descrepancy between the number of values and the number of labels: there 15 each... Any hint on what might go wrong here ? Actually, I think you got it: > factor(1,c(NA,1:4),c(1:5)) Error in factor(1, c

[R] Sampling

Hi Guys, I am in desperate need with sampling. I am suppose to sample from a dataset. I use the following code: > for(i in 1:5){temp[i]<-sample(T2,40,replace=F) + show(temp)} but all the samples are the same, but I want them to be different.any suggestion? Thanks -- View this message in con

Re: [R] follow-up on Error when reading a SAS transport file (with sasxport.get from Hmisc)

I have done what P. Dalgaard has suggested and I don't find a descrepancy between the number of values and the number of labels: there 15 each... Any hint on what might go wrong here ? Here is the output The SAS format from proc contents VISITF

Re: [R] rgl_081.708: rgl.snapshot fails, causing persistent problems

Yihui Xie wrote: Thanks, Duncan, I was encountered with the same problem as Michael, and just now I found the version you provided in your homepage could solve this problem well. Thanks very much, Yihui -- Using Duncan's homepage version also worked for me, and your test helped me slay the f

Re: [R] Help MLE

LFRC yahoo.com.br> writes: > > > Dear Ben Bolker, > > Thanks a lot for your help. > > I have two more questions: > > 1) My goal is maximize the function (> r = Y*log(comb)) but the > parameters found, minimized the function (r = Y*log(comb)). > Oh. Oops. Just change the sign ( r =

Re: [R] glmmPQL

Jean-Baptiste Ferdy univ-montp2.fr> writes: > > Dear all, > > I am experiencing problems with glmmmPQL. I am trying to analyze > binomial data with some spatial autocorrelation. Here is my code and > some of the outputs > > m.1 <- glmmPQL(fixed=cbind(nb_pc_02,I(nb_expr_02-nb_pc_02))~(PSDC99 > +

Re: [R] lag function doesn't work - what am i doing wrong?

Note that there is actually a good reason why this is so. ts can only represent regularly spaced series but the result of a subscripting operation might not be a regularly spaced series so you would not be able to guarantee that it could be a ts series as well. That's the reason for window.ts too

Re: [R] GWR Predictions' standard deviation

Eduardo Marinho gmail.com> writes: > > > Dear all, > > I would like to use a GWR model in order to spatially predict food > insecurity in Africa. I have a georeferenced village data-bases and I've run > a "classic" regression model (taking into account the spatial dependence of > the errors) t

Re: [R] lag function doesn't work - what am i doing wrong?

thanks. i'll try zoo or xts. ts is a p.i.t.a. Gabor Grothendieck wrote: > > lag1resid is a time series, try > > str(residsq) > str(lag1residsq) > > The problem is that when you subscript lag1resid you > don't get a time series out from that. See ?window.ts > or try the zoo or xts packages

Re: [R] number format in plots

Thanks, really it is difficult to find solution Jarek Peter Dalgaard pisze: Jarek Jasiewicz wrote: Hi It is probably simple but how to force in plot command number format: 2, 3, 4 etc. instead of 2e+04, 3e+04 etc. Not all that simple! It's one of those cases where you need to regai

Re: [R] map + some arbitrary locations' heights: some k ind of perspective or contour plot possible?

Werner Wernersen yahoo.de> writes: > > Hi, > > I thought about this but programming it seems rather difficult so I was > wondering if a function exists for > this in R (as most of the times it turns out that it does): > I have a map (shapefile) and for about 50 points on that map (GPS location

Re: [R] Coefficients in a polynomial glm with family poisson/binomial

I don't know what you mean by XCoef x X. But your problem is (as it works if you specify "normal" in a glm) that the functional relationship between your predictors, i.e. Intercept+X+X^2, and Y is not linear for a binomial or a poisson distribution. Generalized linear model implies that the model

Re: [R] lag function doesn't work - what am i doing wrong?

lag1resid is a time series, try str(residsq) str(lag1residsq) The problem is that when you subscript lag1resid you don't get a time series out from that. See ?window.ts or try the zoo or xts packages where subscripting of time series works. On Fri, Oct 10, 2008 at 1:27 PM, Scotty Nelson <[EMAIL

[R] Problems and bugs in vgam()

Hello R-Users, I have recently run into several problems using vgam() in the VGAM package. I am hoping someone might have some solutions... Briefly, I have been trying to fit GAM models for zero-altered negative binomial models. 1. When fitting smoothed parameters (e.g. s(X, df=2)) changin

[R] lag function doesn't work - what am i doing wrong?

I am trying to lag a time series. My data is in a matrix, but I coerce it into a ts object. But when I lag it and then look at the result, nothing has changed. What am I doing wrong? residsq<-resid^2 residsq<-as.ts(residsq) lag1residsq<-lag(residsq,-1) > residsq[1:5] 1

[R] Coefficients in a polynomial glm with family poisson/binomial

Dear R-users When running a glm polynomial model with one explanatory variable (example Y~X+X^2), with a poisson or binomial error distribution, the predicted values obtained from using the predict() function and those obtained from using the coefficients from the summary table "as is" in an eq

[R] help on quantile

Hi, I would like to use the quantile function, but I have some doubts to choose the value of "type" parameter. Does this depend on the shape/structure of my distribution? I cannot unde rstand what mean the expressions Discontinuous sample quantile and Continuous sample quantile, respectively. My

[R] a question about diparate device quality

List, I have a question about why plots sent to certain devices can be of so much worse quality than others. Whenever I plot to x11 or pdf the graphs are of excellent quality. eps is a bit poorer, but of passable quality, and jpeg or tiff are terrible. I tried all manner of parameters, on many

Re: [R] help

Vijaykumar Muley gmail.com> writes: > I like to draw a cluster tree with a hclust and plot command. > But,instead of labels, it is giving number of variables in dendogram. > Could anyone tell me , how to print the labels of variables used for drawing. > Thank you. Please post str(of your data),

Re: [R] Overlay Plots

?matplot On Fri, Oct 10, 2008 at 9:50 AM, mentor_ <[EMAIL PROTECTED]> wrote: > > Hi, > > I have three vectors and I would like to plot them in one window. So R > should overlay them in the > quartz. I achieved this by plotting the first vector: plot(v1) and add the > other two with > points(v2) an

Re: [R] Leap year?

perhaps leapyear() http://finzi.psych.upenn.edu/R/library/pheno/html/leapyear.html --- On Fri, 10/10/08, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > From: [EMAIL PROTECTED] <[EMAIL PROTECTED]> > Subject: [R] Leap year? > To: "r-help" > Received: Friday, October 10, 2008, 10:29 AM > Given

Re: [R] Space between bars in barplot

> I played around with your suggestions to change the appearance of my graph, > but there is still a problem I could not fix. > The vector, which I want to plot, contains 60 entries. After 3 bars I want > to have a large gap > between the next 3 bars. But what always happens is, that R took at th

[R] glmmPQL

Dear all, I am experiencing problems with glmmmPQL. I am trying to analyze binomial data with some spatial autocorrelation. Here is my code and some of the outputs > colnames(d.glmm) [1] "BV" "Longitude" "Latitude""nb_pc_02" "nb_expr_02" [6] "pc_02" "nb_pc_07""nb_expr_

Re: [R] Space between bars in barplot

Hi, thanks again for your help. I played around with your suggestions to change the appearance of my graph, but there is still a problem I could not fix. The vector, which I want to plot, contains 60 entries. After 3 bars I want to have a large gap between the next 3 bars. But what always happens

[R] Overlay Plots

Hi, I have three vectors and I would like to plot them in one window. So R should overlay them in the quartz. I achieved this by plotting the first vector: plot(v1) and add the other two with points(v2) and points(v3). Is there another way to achieve the same result? Cheers -- View this message

[R] Implied GARCH parameters - Error in chol.default(Hessian) : the leading minor of order 1 is not positive definite

Dear R-Users, I'm trying to estimate GARCH parameters implied by options prices (on wednesdays only) by minimizing the average mean squared dollar error between the market and the model price, but I always get the following error message: Error in chol.default(Hessian) : the leading minor of

Re: [R] Help MLE

Dear Ben Bolker, Thanks a lot for your help. I have two more questions: 1) My goal is maximize the function (> r = Y*log(comb)) but the parameters found, minimized the function (r = Y*log(comb)). 2) What this function do? > model2 <- function(p) { > do.call("Model",as.list(p)) > }

Re: [R] Fatal error: unable to restore saved data in .RData

On Fri, 10 Oct 2008, Thompson, David (MNR) wrote: Hello, I may have found the problem, although I don't understand it. When I commented out the following line in my Rprofile.site file the problem _seemed_ to disappear. # ans as shorthand for .Last.value # makeActiveBinding("ans"

Re: [R] Creating a new connection type

Dear Paulo, > Hello. I have an instrument, actually a step motor indexer that I want to > drive from R. It uses an RS-232 connection. I already have a few C > subroutines that I access from R but I would like to be able use R's > connection related functions - readLines, writeLines, cat, read.tabl

Re: [R] Fatal error: unable to restore saved data in .RData

If you are able to create a .Rdata file for which this fails can you record the last few commands before you create the .Rdata file? luke On Fri, 10 Oct 2008, Thompson, David (MNR) wrote: Hello, I am unable to start an R session by double-clicking on an existing .RData file from the Windoze

[R] Creating a new connection type

Hello. I have an instrument, actually a step motor indexer that I want to drive from R. It uses an RS-232 connection. I already have a few C subroutines that I access from R but I would like to be able use R's connection related functions - readLines, writeLines, cat, read.table, etc - directly.

Re: [R] Leap year?

charter.net> writes: > > Given a Date object or simply a year is there an R function to tell me if the it is a leap year or not? I was hoping > for something like 'is.leapyear'. I probably can build my own function (year divisible by 4 etc.) but I > would rather use an existing function if it i

Re: [R] Dump decision trees of randomForest object

If you just want to be able to use the trained RF model in some future R session for prediction on new data, just use save() to save the RF object, and load() it back in the future. If you really want to write your own low-level code for prediction, you can take a look at the predictRegTree() f

Re: [R] ggplot adding points

On Fri, Oct 10, 2008 at 9:18 AM, stephen sefick <[EMAIL PROTECTED]> wrote: > I would like to do the following in ggplot: > what am I missing? > > River.Mile <-c(202, 198, 190, 185, 179, 148, 119, 61) > TSS <- c(1:8) > DOC <- seq(2, by= 0.6, length.out=8) > z <- data.frame(River.Mile, TSS, DOC) > xy

Re: [R] Fatal error: unable to restore saved data in .RData

Hello, I may have found the problem, although I don't understand it. When I commented out the following line in my Rprofile.site file the problem _seemed_ to disappear. # ans as shorthand for .Last.value # makeActiveBinding("ans", function() .Last.value, .GlobalEnv) I decided to d

Re: [R] boxplot and stripchart

On Fri, Oct 10, 2008 at 4:24 PM, Peter Dalgaard <[EMAIL PROTECTED]> wrote: > Göran Broström wrote: >> I want to make illustrations a la boxplot for a data material with >> several groups. Some of the groups are very small, so a boxplot >> doesn't make sense. I would like to use stripchart for these

[R] Fatal error: unable to restore saved data in .RData

Hello, I am unable to start an R session by double-clicking on an existing .RData file from the Windoze file explorer. A dialogue appears with the message "Fatal error: unable to restore saved data in .RData", and I notice in the R Console the last few lines are: "Loading required

Re: [R] ggplot adding points

Dear Stephen, You need the data in long format. Try melt()ing it first. River.Mile <-c(202, 198, 190, 185, 179, 148, 119, 61) TSS <- c(1:8) DOC <- seq(2, by= 0.6, length.out=8) z <- data.frame(River.Mile, TSS, DOC) zMelt <- melt.data.frame(z, id.vars = "River.Mile") library(ggplot2) ggplot(zMelt,

Re: [R] boxplot and stripchart

Göran Broström wrote: > I want to make illustrations a la boxplot for a data material with > several groups. Some of the groups are very small, so a boxplot > doesn't make sense. I would like to use stripchart for these groups. > On the other hand, some of the groups are very large, so for them > s

[R] Leap year?

Given a Date object or simply a year is there an R function to tell me if the it is a leap year or not? I was hoping for something like 'is.leapyear'. I probably can build my own function (year divisible by 4 etc.) but I would rather use an existing function if it is available. Thank you. Kevi

Re: [R] Plot grouped histograms

Another solution would be faclevels<-c("Not","Little","Somewhat","Very", NA) group <- factor(1:6) svdf<-data.frame(value = sample(group, 400, TRUE), category = sample(faclevels,400,TRUE)) library(ggplot2) ggplot(svdf, aes(x = value, fill = category)) + geom_histogram(position = "dodge") HTH, Thi

[R] ggplot adding points

I would like to do the following in ggplot: what am I missing? River.Mile <-c(202, 198, 190, 185, 179, 148, 119, 61) TSS <- c(1:8) DOC <- seq(2, by= 0.6, length.out=8) z <- data.frame(River.Mile, TSS, DOC) xyplot(TSS+DOC~River.Mile, data=z, auto.key=TRUE) thanks -- Stephen Sefick Research Scient

Re: [R] plot-parameter pch without influence when plotting a data-frame

On Thu, Oct 9, 2008 at 8:19 AM, John Kane <[EMAIL PROTECTED]> wrote: > --- On Thu, 10/9/08, Oliver Bandel <[EMAIL PROTECTED]> wrote: > >> From: Oliver Bandel <[EMAIL PROTECTED]> >> Subject: Re: [R] plot-parameter pch without influence when plotting a >> data-frame >> To: "Gerhard Schön, UKE Hambur

Re: [R] Turn factors to numeric

Not an answer but I find that DF <- data.frame(let = letters[1:3], num = 1:3, stringsAsFactors = FALSE) is very handy. Damn it, if I want a factor I'll tell the machine I do :) Or you can also set this option globally with options(stringsAsFactors = TRUE) Gabor Grothendiec

Re: [R] plot-parameter pch without influence when plotting a data-frame

--- On Fri, 10/10/08, Oliver Bandel <[EMAIL PROTECTED]> wrote: > From: Oliver Bandel <[EMAIL PROTECTED]> > Subject: Re: [R] plot-parameter pch without influence when plotting a > data-frame > To: [EMAIL PROTECTED] > > > > > > > what is the result of: > [...] > > > > Try something like >

[R] boxplot and stripchart

I want to make illustrations a la boxplot for a data material with several groups. Some of the groups are very small, so a boxplot doesn't make sense. I would like to use stripchart for these groups. On the other hand, some of the groups are very large, so for them stripchart isn't so good. So my q

[R] Bootstraping problem

__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Credit Scoring Model - SPEC (specificity) and SENS (sensitivity)

Maithili Shiva wrote: Dear R helpers, Hi I am working on credit scoring model using logistic regression. I have main sample of 42500 clentes and based on their status as regards to defaulted / non - defaulted, I have genereted the probability of default. I have a hold out sample of 5000 clien

Re: [R] rgl-snapshot failed (err-msg: "failed")

Zitat von Duncan Murdoch <[EMAIL PROTECTED]>: > On 10/10/2008 8:13 AM, Oliver Bandel wrote: > > Hello, > > > > > > I tried to use rgl.snapshot and it failed. > > > > The error message was not very verbose: > > > > > > > > == > >> > >> plot3d( motion[[idx+2]], motion[[idx+1]], m

Re: [R] Plot grouped histograms

[EMAIL PROTECTED] wrote: r11 -- r16 are variables showing a reason for usage of a product in 6 different situations. Each variable is a factor with 4 levels imported from a SPSS sav file with labels ranging from "not important" to "very important", and NA's for a sample of N = 276. (1) I need

Re: [R] rgl-snapshot failed (err-msg: "failed")

On 10/10/2008 8:13 AM, Oliver Bandel wrote: Hello, I tried to use rgl.snapshot and it failed. The error message was not very verbose: == plot3d( motion[[idx+2]], motion[[idx+1]], motion[[idx]] ) rgl.snapshot(filename="/tmp/shot_01.png", fmt="png") [1] "failed" ==

[R] rgl-snapshot failed (err-msg: "failed")

Hello, I tried to use rgl.snapshot and it failed. The error message was not very verbose: == > > > plot3d( motion[[idx+2]], motion[[idx+1]], motion[[idx]] ) > rgl.snapshot(filename="/tmp/shot_01.png", fmt="png") [1] "failed" > == There was a graphic c

Re: [R] multi-scale singular spectrum analysis (SSA)

wmtsa has some wavelet processing things, but It sounds like you need a very specific implementation. Sorry I would not be of more help- good luck. Stephen On Fri, Oct 10, 2008 at 7:26 AM, Maik Renner <[EMAIL PROTECTED]> wrote: > Hi everybody! > > I am searching for an R implementation of multi-

Re: [R] rgl_081.708: rgl.snapshot fails, causing persistent problems

Thanks, Duncan, I was encountered with the same problem as Michael, and just now I found the version you provided in your homepage could solve this problem well. Regards, Yihui -- Yihui Xie <[EMAIL PROTECTED]> Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http

Re: [R] Creating GUIs for R

I haven't heard Rpad mentioned yet, so I will mention it. Rpad allows you to run R code inside a browser window. You can many any kind of GUI you want using html forms and then call the R function from the html page. The Apache/Rpad install works on any platform and anyone who uses it just ne

[R] multi-scale singular spectrum analysis (SSA)

Hi everybody! I am searching for an R implementation of multi-scale singular spectrum analysis (SSA), which was introduced by Yiou etal (2000) (Data-adaptive wavelets and multi-scale singular-spectrum analysis) and further described in Ghil etal (2002). For SSA alone I found an package recently p

Re: [R] plot-parameter pch without influence when plotting a data-frame

Zitat von John Kane <[EMAIL PROTECTED]>: > --- On Thu, 10/9/08, Oliver Bandel <[EMAIL PROTECTED]> wrote: > > > From: Oliver Bandel <[EMAIL PROTECTED]> > > Subject: Re: [R] plot-parameter pch without influence when plotting > a data-frame > > To: "Gerhard Schön, UKE Hamburg" <[EMAIL PROTECTED]> > >

Re: [R] histogram loses top row with alpha transparency under Windows

Richard.Cotton wrote: > The rectangles being drawn extend higher than the top of the panel. (Your y > axis ranges from 0 to 50, but the bars go up to 100.) > Thankyou - I can also make the bars on the lower panels vanish by tinkering with ylim. > In the top row of plots, depending upon the shap

Re: [R] linear expenditure model

Hi Marie! On Friday 10 October 2008 12:40:23, Marie Vandresse wrote: > I would like to estimate a linear expendire with Systemfit package. > (method: "SUR") If I remember correctly, the "linear expenditure system" (LES) is linear in income but non-linear in the parameters. Hence, you have to est

Re: [R] nls, lattice, and conversion over to ggplot

I just got it and it works wonderfully. thank you both for your help Stephen On Fri, Oct 10, 2008 at 4:01 AM, baptiste auguie <[EMAIL PROTECTED]> wrote: > It worked for me, do you have the latest version of ggplot2 released a few > days ago (ggplot2_0.7) ? > > Baptiste > > On 9 Oct 2008, at 20:55

Re: [R] formal argument "axes" matched by multiple actual arguments

I don't know where you got the idea that 'axes' was a valid argument. From ?plot.Sample ...: further graphical parameters. and 'axes' is *not* a graphical parameter. (It also suggests you call plot, not plot.Sample.) 'xaxt' and 'yaxt' are graphics parameters which control whether axes

[R] help

Hello, R users I like to draw a cluster tree with a hclust and plot command. But,instead of labels, it is giving number of variables in dendogram. Could anyone tell me , how to print the labels of variables used for drawing. Thank you. -- Vijaykumar Yogesh Muley Senior Research Fellow, Dr. A

[R] linear expenditure model

Hi, I would like to estimate a linear expendire with Systemfit package. (method: "SUR") As someone could show me how to define the equations? Thanks. -- Think before you print ! ** Disclaimer: This e-mail may contain

Re: [R] Correlation among correlation matrices cor() - Interpret

On 10-Oct-08 08:07:34, Michael Just wrote: > Hello, > If I have two correlation matrices (e.g. one for each of two > treatments) and then perform cor() on those two correlation > matrices is this third correlation matrix interpreted as the > correlation between the two treatments? > > In my sample

[R] Multi-scale SSA (singular spectrum analysis)

Hi everybody! I am searching for an R implementation of multi-scale SSA (singular spectrum analysis), which was introduced by Yiou etal (2000) (Data-adaptive wavelets and *multi*-*scale *singular-spectrum analysis) and further describe

[R] map + some arbitrary locations' heights: some kind of perspective or contour plot possible?

Hi, I thought about this but programming it seems rather difficult so I was wondering if a function exists for this in R (as most of the times it turns out that it does): I have a map (shapefile) and for about 50 points on that map (GPS locations) I have heights. Is there a function which can m

[R] formal argument "axes" matched by multiple actual arguments

Hi, I'm using the add-on package "sound".I have the following > q<-loadSample("a.wav") > plot.Sample(q,axes=FALSE) this gives me the error. Error in plot.default(sound(s)[1, ], type = "l", col = "red", ylim = c(-1, : formal argument "axes" matched by multiple actual arguments I'm guessing the

Re: [R] R seven years ago

Not to detract from John's very interesting work, but there is a problem with citations as a measure of penetration. When R totally succeeds, there will be no citations. How many citations are there for Excel? However, for the moment, citations are probably a pretty good measure. Using traffic

Re: [R] how to store lme/lmer fit result

> I am building a hierarchical model on a large data set. It can take quite > some time to finish one fit, I was just wondering whether it is possible to > store the fit object (the result) to a file for later (offline) analysis. See ?save. Regards, Richie. Mathematical Sciences Unit HSL --

[R] Antwort: FW: R loops

it did help to create the separate vectors. The problem is, I would like to continue working with each vector on its own and make some further analysis. For that i need to allocate each vector in the filePatterns a unique name, which I just can't do. I'm not so great in R loops and the help

Re: [R] help

Bing Shen wrote: Hello, I'm not sure this is doable but I'm having trouble running my R script with multithreaded capability. With 16x2.93Ghz CPUs available, only one is running with 100%. Any suggestions? Yes: Run 16 R jobs or use some other way to make your job parallel (e.g. the snow fram

[R] Correlation among correlation matrices cor() - Interpretation

Hello, If I have two correlation matrices (e.g. one for each of two treatments) and then perform cor() on those two correlation matrices is this third correlation matrix interpreted as the correlation between the two treatments? In my sample below I would interpret that the treatments are 0.28 co

Re: [R] nls, lattice, and conversion over to ggplot

It worked for me, do you have the latest version of ggplot2 released a few days ago (ggplot2_0.7) ? Baptiste On 9 Oct 2008, at 20:55, stephen sefick wrote: Error in `[.data.frame`(df, , var) : undefined columns selected I got this error in a fresh R session after rerunning all of the comm