[R] Data manipulation

2007-10-14 Thread Klaus Friis Østergaard
Hi, I have a data set which is like this I write as the CSV I import: Sample;Hole;Feature;Value 1;5;x;4,2334 1;5;y;3,3434 1;5;r;0,1080 1;10;x;5,2526 1;10;y;4,3434 1;10;r;0,1080 with 98 sample and 10 different holes. These are measured values. Now I also have a list of nominel values: Hole

Re: [R] Need some help

2007-10-14 Thread Daniel Nordlund
> -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf > Of azzza > Sent: Sunday, October 14, 2007 10:21 PM > To: r-help@r-project.org > Subject: [R] Need some help > > > Hi! > I'm taking a course that requires some programming background, but I'm a > complete

[R] Need some help

2007-10-14 Thread azzza
Hi! I'm taking a course that requires some programming background, but I'm a complete novice in the field. when asked to generate a list of 20 uniform random numbers, is it alright if I put in >randu, and just copy-paste the first 20 numbers?? Or is there, as I suspect, a better way of calling ou

Re: [R] R API - optim

2007-10-14 Thread Wang Chengbin
The attached code is compiled properly. The error happens at the linking stage. It seems that R shared library is not built by configure/make. 2007/10/14, Dirk Eddelbuettel <[EMAIL PROTECTED]>: > > > On 14 October 2007 at 05:58, Wang Chengbin wrote: > | I am trying to use the R API to call optim

Re: [R] how to use normalmixEM to get correct result?

2007-10-14 Thread Tatiana Benaglia
Hi, you're having problems with the initial values for the EM. If you don't give initial values for lambda, mu, and sigma, the function tries to generate them, which is causing the errors. I'd try to input initial values based on your data. Do the histogram of log.x (you don't need to sort

[R] some question about partial prediction in survival

2007-10-14 Thread coldeyes.Rhelp
Hi there: i got a problem to get the prediction from a model recently. for example if i use a survival analysis to predict the risk. i use the code like below: i found the the prediction is not equal to (coef * x + coef * sex) , could someone help me with why this happened? and can someone ex

[R] how to use normalmixEM to get correct result?

2007-10-14 Thread bgchen
Dear R-Users, I have a large number of data(54000) and the field of data is 50 to 2.0e9. I want to use normalmixEM (package:mixtools) to fit them in finite mixture narmal distributions,but get some mistakes.I don't know which steps make the error. I have used the following functions

[R] specifying models ":" and %in% and soforth

2007-10-14 Thread Daniel Lakeland
If a and b are factor variables are there any differences between ~ a:b ~ b:a ~ a %in% b ~ b %in% a playing with model.matrix doesn't show any differences in the resultant matrices, which is what I expect based on my understanding if what's being specified. I'm hoping that someone can clarify

[R] significance for a random effect in Mixed Model ANOVA

2007-10-14 Thread Nathaniel Street
In a number of cases I want to use mixed-model ANOVA tests where I am interested in whether both the fixed and random effects (and their interactions) are significant. If I use this example >library(nlme) >data(Orthodont) >anova(lme(distance ~ age + Sex, data = Orthodont, random = ~ 1)) I get th

Re: [R] cross tabulation

2007-10-14 Thread James Reilly
On 14/10/07 5:39 PM, Louis Martin wrote: > I am looking for an efficient and fast way for doing this: > > for (v in 1:dw) {ok <-list1[v,] > > for (z in 1:5) {kk = ok[z] ; > for (o in 1:t) {if (kk== list[o]) n<- n+1}} > show(n) } > > list1 : a data.fra

Re: [R] ggplot2 geom_hline error

2007-10-14 Thread Jim Nikelski
Hi Hadley (and List), Your suggested geom_hline work-around did the trick. Also, thanks for the info regarding axis labels. My world is a happy place again. Enjoy your weekend ;) -Jim On 10/14/07, hadley wickham <[EMAIL PROTECTED]> wrote: > On 10/14/07, Jim Nikelski <[EMAIL PROTECTED]> wrot

Re: [R] Dendrogram

2007-10-14 Thread jim holtman
assign then names to the rows in the dataframe: row.names(d) <- c("Italy", "France", "Spain", ...) On 10/14/07, Alessandra Trimarchi <[EMAIL PROTECTED]> wrote: > Dear r-help mailing list, > here there is a data-frame of 24 observations and 7 variables. >W3 D3P3 T3 U3 V3 R3 > 1 31.3

Re: [R] ggplot2 geom_hline error

2007-10-14 Thread hadley wickham
On 10/14/07, Jim Nikelski <[EMAIL PROTECTED]> wrote: > Hello R-List > > I've run into a minor problem with ggplot2. In particular, I get > an error message when trying to plot a horizonal line through zero. > Here is what I do: > > plt.df <- subset(xfm.df, select=c(stage, subj, tX)) > plt <- g

[R] Dendrogram

2007-10-14 Thread Alessandra Trimarchi
Dear r-help mailing list, here there is a data-frame of 24 observations and 7 variables. W3 D3P3 T3 U3 V3 R3 1 31.3 52 50200 1.67 59.4 55 57 2 33.0 61 45684 0.83 62.0 64 85 3 29.6 55 30455 0.60 56.6 45 52 4 28.0 59 8858 0.96 22.1 40 45 5 29.1 53 8334 0.63 42.9 46 51 6 23.1 79

[R] ggplot2 geom_hline error

2007-10-14 Thread Jim Nikelski
Hello R-List I've run into a minor problem with ggplot2. In particular, I get an error message when trying to plot a horizonal line through zero. Here is what I do: plt.df <- subset(xfm.df, select=c(stage, subj, tX)) plt <- ggplot(pp.df, aes(x=stage, y=tX, group=subj)) + geom_line() plt plt

Re: [R] reading data and Variance Ratio Tests for Weak-form Market Efficiency

2007-10-14 Thread John Kane
--- Jasim Al-AJmi <[EMAIL PROTECTED]> wrote: > Dear Rhelp > > I am first time user of R and I have the following > questions: > > 1. How I can load my data (it is a matrix 10x1618) > with the titles of the variables in the first row? You probably should give us some more infromation about th

[R] repeated measures - aov, lme, lmer - help

2007-10-14 Thread juli pausas
Dear all, I'm not very sure on the use of repeated measures in R, so some advice would be very appreciate. Here is a simple example similar to my real problem (R 2.6.0 for windows): Lets supose I have annual tree production measured in 9 trees during 3 years; the 9 trees are located in 3 different

Re: [R] Make playwith a default graphic device

2007-10-14 Thread Ronaldo Reis Junior
Em Dom 14 Out 2007, Prof Brian Ripley escreveu: > You seem to be under the misapprehension that playwith() is a graphics > device: it is not. > > To set a graphics device as the default, use options(device=). I am not > sure what 'more functionalities' you are looking for: the playwith() uses > th

Re: [R] Question re matplot

2007-10-14 Thread jim holtman
You can always add the ylab with 'mtext' and reduce the font size if it is too big. On 10/13/07, David Kaplan <[EMAIL PROTECTED]> wrote: > Yea, sorry, that was a typo when I copied into my emal. > > Here it is again > matplot(battingagg$X, battingagg[, c("HR","RBI","X2B", "BB", > "R", "SB")], type

Re: [R] Adjusting for heaping in data

2007-10-14 Thread Ted Harding
On 14-Oct-07 08:33:41, Thomas Frööjd wrote: > Hi R users. I am new to the community and have got myself into > a little problem. It does not look as though it was yourself who got you into this problem! You have been given the bathwater along with the baby. > I have a dataset of birth weights rec

[R] ggplot2: ordering categorial data

2007-10-14 Thread Christoph Krammer
Hello again, everytime I think I got something to work, the next issue comes up... I have the following data.frame, I want to visualize: > data_rb tld spam1 spam2 share 1 ca 826436 73452 0.0889 2 org 470550 25740 0.0547 3 de 156042 15531 0.0995 4 com 140753 7527 0.

Re: [R] Make playwith a default graphic device

2007-10-14 Thread Felix Andrews
Hi Ronaldo, As Professor Ripley explained, playwith is a function, not a device. It uses a Cairo-based device from the cairoDevice package. I too would like to have `playwith` as the default behavior for some plots, but there is no standard way to do that. It can not be implemented as a device bec

Re: [R] Make playwith a default graphic device

2007-10-14 Thread Prof Brian Ripley
You seem to be under the misapprehension that playwith() is a graphics device: it is not. To set a graphics device as the default, use options(device=). I am not sure what 'more functionalities' you are looking for: the playwith() uses the cairoDevice package for its graphics device, and ther

[R] Adjusting for heaping in data

2007-10-14 Thread Thomas Frööjd
Hi R users. I am new to the community and have got myself into a little problem. I have a dataset of birth weights recorded by nurses at a delivery clinic in an developing country. The weights are entered in KiloGrams with one decimal. However there is substantial heaping at each 500g when lookin