[R] Minimization of the distance

2010-12-09 Thread bluesky
I just contect R,and still learn how to write the code. I have a problem with argmin sum d(pi,p)/n for example I have 3 points (a1,b1)(a2,b2)(a3,b3) ,then I want to find p(x,y) make sure that (sqrt((a1-x)^2+(b1-y)^2)+sqrt((a2-x)^2+(b2-y)^2)+sqrt((a3-x)^2+(b3-y)^2))/3 is the minimum. -- View

Re: [R] ReadWrite.xls problem

2010-12-09 Thread haruo0409
Hello tkdweber. I'm also annoyed at same problem. I installed xlsReadWriter today and entered x <- read.xls("data.xls",sheet=1) But I got Error Message: 以下にエラー .Call("ReadXls", file, colNames, sheet, type, from, rowNames, : 引数の個数(11)が不正です。10 個が ReadXls に対しては必要です (It's Japanese.Its Engli

Re: [R] Question on ARIMA Prediction

2010-12-09 Thread Dennis Murphy
Hi: Usually, if you get no response from the list, it's because either (a) your question is vague or ill-defined; (b) you supplied no reproducible example to illustrate what you tried, what was the error and what you expected to get (per the instructions in the Postin

Re: [R] RES: Barplot with "Independent" Lines Y axis

2010-12-09 Thread Dennis Murphy
Hi: Like Peter Ehlers, I'm not a big fan of multiple response variables on a page with different y-axis scaling, but if you have to do it, try not to let one graphical metaphor interfere with/obscure/dominate the other (e.g., bars and lines). This is my attempt, both with connecting lines and poin

Re: [R] Question on ARIMA Prediction

2010-12-09 Thread Jiakang Lu
Does anyone have similar experience before? I'd appreciate any suggestion/advice from the community. On Wed, Dec 8, 2010 at 4:05 PM, Jiakang Lu wrote: > Dear all, > > I'm new to R and time series analysis. I'd appreciate if you could shed > light on my problem. > > Here is what I have been tryin

Re: [R] Projecting data on a world map using long/lat

2010-12-09 Thread Greg Snow
The symbols function can add circles to an existing plot with the diameter based on another variable. You may need to project points if you draw a map with a projection. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -O

Re: [R] Delete observations with a frequency < x

2010-12-09 Thread Seeliger . Curt
mathijsdevaan wrote on 12/09/2010 04:21:54 PM: > I have two columns with data (both identifiers - it's an affiliation list) > and I would like to delete the rows in which the observations in the second > column have a frequency < 5 in the entire second column. Example: > > 1 a > 1 b > 1

[R] Projecting data on a world map using long/lat

2010-12-09 Thread mathijsdevaan
Hi, I have a dataset (CSV) with some counts of firms located around the globe. Each count is assigned to the longitude and latitude of the specific location. Now I want to plot these counts on a world map using dots (size of dots represent the count). I have been unable to find any info on whethe

Re: [R] Delete observations with a frequency < x

2010-12-09 Thread mathijsdevaan
Hi Phil, That worked perfectly! Thanks Mathijs -- View this message in context: http://r.789695.n4.nabble.com/Delete-observations-with-a-frequency-x-tp3081226p3081264.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.

Re: [R] Delete observations with a frequency < x

2010-12-09 Thread Phil Spector
Suppose this is your data frame: df = data.frame(x=c(1,1,1,2,2,2),y=c('a','b','c','a','b','d')) df x y 1 1 a 2 1 b 3 1 c 4 2 a 5 2 b 6 2 d df[!table(df$y)[df$y] < 2,] x y 1 1 a 2 1 b 4 2 a 5 2 b Note that this will only work properly if y is a factor or character variable. If y was num

Re: [R] help on timeseries

2010-12-09 Thread Dennis Murphy
Hi: This is more or less the same idea as the other respondent, but uses function rle() instead. The function: findDates <- function(df) { u <- rle(df$Signals) n <- length(u$values) ends <- with(u, cumsum(lengths)[values == -1]) begins <- with(u, cumsum(lengths)[values == 1]

[R] Delete observations with a frequency < x

2010-12-09 Thread mathijsdevaan
Hi, I have two columns with data (both identifiers - it's an affiliation list) and I would like to delete the rows in which the observations in the second column have a frequency < 5 in the entire second column. Example: 1 a 1 b 1 c 2 a 2 b 2 d Let's say, I would like to

Re: [R] RES: Barplot with "Independent" Lines Y axis

2010-12-09 Thread Peter Ehlers
On 2010-12-09 14:22, Rodrigo Aluizio wrote: Oh sorry. An example say lots more than words. The data below, when submitted to twoord.plot return the mentioned error. Rain are bars and Salt lines, the bars appear and the error occurs with the salt data. If you get an error complaining about need

Re: [R] Number of dimension in Multidimensional Scaling

2010-12-09 Thread Michael Bedward
Just to add to Michael F's comments: I've looked for that elbow many a time but never found it :) Admittedly, I typically deal with fairly noisy, ecological data, but I think it's a mistake to try to identify the "optimal" number of dimensions. Better instead to concentrate on a "useful" number, i

Re: [R] Reading SQL Server Tables using RODBC

2010-12-09 Thread pfruin
Paul Miller wrote: > > Hello Everyone, >   > I'm still quite new to R and am having trouble reading SQL Server Tables > using RODBC. I've looked though the RODBC documentation as well as > material I found online and in the book "Data Manipulation with R," but I > just can't quite seem to get th

Re: [R] order matrix by column position

2010-12-09 Thread Folkes, Michael
Thanks Gabor, I should have paid closer to the examples within the function help. I discovered the same approach minutes after submission. :\ problem solved! ___ Michael Folkes Salmon Stock Assessment Canadian Dept. of Fisheries & Oceans

Re: [R] order matrix by column position

2010-12-09 Thread Gabor Grothendieck
On Thu, Dec 9, 2010 at 5:24 PM, Folkes, Michael wrote: > Hello all, > I'm trying to grasp a way to order a matrix by giving order() only a vector > of the columns it can look to for sorting. > The approach has to be generic such that I can feed order() a vector of > variable length (ie varying #

[R] order matrix by column position

2010-12-09 Thread Folkes, Michael
Hello all, I'm trying to grasp a way to order a matrix by giving order() only a vector of the columns it can look to for sorting. The approach has to be generic such that I can feed order() a vector of variable length (ie varying # of columns to sort by). x<-rep(1,9) y<-c(1, 1, 2, 2, 2, 1, 2, 2

[R] RES: Barplot with "Independent" Lines Y axis

2010-12-09 Thread Rodrigo Aluizio
Oh sorry. An example say lots more than words. The data below, when submitted to twoord.plot return the mentioned error. Rain are bars and Salt lines, the bars appear and the error occurs with the salt data. Point RainSalt Fev/03 365.6 13 Mar/03 235 18 Abr/03 115.1 18 Mai/03 47.4

Re: [R] Encoding problem - I fails to read Hebrew text from online

2010-12-09 Thread Matt Shotwell
Tal, OK, let me clarify my understanding. The original and decoded file are text, encoded by UTF-8. In the original file, there are HTML `entities' that represent UTF-8 Hebrew characters. In the decoded file, the entities are converted to UTF-8 characters. The question is how to convert these ent

Re: [R] Bivariate kernel density bandwidth selection

2010-12-09 Thread Gaj Vidmar
I'm faar from knowledgeable on R or kernel density estimation (or many other statistical things for that matter), but allow me a simplistic suggestion: I remember that in a package I co-authored quite some time ago (chplot -- maybe you might have a look at what it does) we use bkde2D from the

Re: [R] How does Sweave write to a file?

2010-12-09 Thread Duncan Murdoch
On 09/12/2010 2:07 PM, Arnau Mir Torres wrote: Hello. I need that Sweave writes R output to a file. To do it, I put the following into the foo.Rnw file: \documentclass{article} \usepackage{Sweave} \begin{document} ... <>= x<- rnorm(100) y<- rnorm(100) sink("foo.txt") summary(lm(y~x)) sink()

Re: [R] Plotting 3d surfaces

2010-12-09 Thread Duncan Murdoch
On 09/12/2010 11:20 AM, Uwe Wolfram wrote: Dear Subscribers, I am using R for quite a while nowadays on Ubuntu 10.04 LTS. I a using R for doing my statistics. Furthermore I am using it as a tool to generate the graphics for my publications. I am currently working on a project which involves nls

[R] attributable cost estimation using aggregate data

2010-12-09 Thread n . mitsakakis
Hello, I am facing with an unusual problem of using aggregate data in order to estimate the attributable cost of a disease, for different stages. My data set consist of mean and std estimates of the cost outcome corresponding to strata coming from cross-classification of a set of factors

Re: [R] Barplot with "Independent" Lines Y axis

2010-12-09 Thread Greg Snow
Without seeing a reproducible example we cannot be sure, but my guess is that you are letting twoord.plot set the limits and the function does not remove missing values, if you specify rylim and or lylim arguments specifically, then it should not run into the problem you are seeing (I hope). If

Re: [R] Plotting 3d surfaces

2010-12-09 Thread John Fox
Dear Ben, Though still accessible through the Rcmdr, the scatter3d() function has been moved to the car package. Best, John On Thu, 9 Dec 2010 19:41:12 + (UTC) Ben Bolker wrote: > Barry Rowlingson lancaster.ac.uk> writes: > > > > > On Thu, Dec 9, 2010 at 4:20 PM, Uwe Wolfram in-chemn

Re: [R] lapply getting names of the list

2010-12-09 Thread David Winsemius
On Dec 9, 2010, at 2:21 PM, David Winsemius wrote: On Dec 9, 2010, at 12:44 PM, Sashi Challa wrote: Hello All, I have a toy dataframe like this. It has 8 columns separated by tab. NameSampleIDAl1 Al2 X Y R Th rs191191A1 A B 0.9

Re: [R] Plotting 3d surfaces

2010-12-09 Thread Ben Bolker
Barry Rowlingson lancaster.ac.uk> writes: > > On Thu, Dec 9, 2010 at 4:20 PM, Uwe Wolfram in-chemnitz.de> wrote: > > I am currently working on a project which involves nls-fits of three > > dimensional surfaces such as ellipsoids or even more complex. I have > > been searching R help and the m

Re: [R] Number of dimension in Multidimensional Scaling

2010-12-09 Thread Michael Friendly
On 12/9/2010 7:26 AM, Petar Milin wrote: Hello! Very often one can hear that MDS usually ends with two-dimensional solution. Of course, there are methods, like Scree-test (proposed by Kruskal and Wish, 1981), to determine optimal number of dimensions. However, I am trying to find references to th

Re: [R] Encoding problem - I fails to read Hebrew text from online

2010-12-09 Thread Tal Galili
Hi Matt, Thanks for having a look at this. I just spent some time looking around and couldn't find any R function to decode decimal HTML code. Do you (or someone else on the list) knows how to program this sort of thing? (is there a formula for the translation? p.s: For it to work on my end I a

[R] RES: Barplot with "Independent" Lines Y axis

2010-12-09 Thread Rodrigo Aluizio
Thank you for the function suggestion, works nicely for complete data vectors. Just another question. When using the twoord.plot I reached another issue, it seems that this function can't handle NAs in one of the variables (the data of one month for another variable is missing). The lines and barpl

[R] lapply getting names of the list

2010-12-09 Thread Sashi Challa
Thanks a lot Joshua, that works perfectly fine. I could not think to lapply on the names instead of data itself. I don't now notice SampleID names in the column names. Thanks for your time, -Sashi -Original Message- From: Joshua Wiley [mailto:jwiley.ps...@gmail.com] Sent: Thursday, Dece

Re: [R] lapply getting names of the list

2010-12-09 Thread David Winsemius
On Dec 9, 2010, at 12:44 PM, Sashi Challa wrote: Hello All, I have a toy dataframe like this. It has 8 columns separated by tab. NameSampleIDAl1 Al2 X Y R Th rs191191A1 A B 0.999 0.090.780.090 abc928291 A1 B

Re: [R] Condional Density Plot from different data

2010-12-09 Thread Deepayan Sarkar
On Thu, Dec 9, 2010 at 10:11 AM, Doran, Harold wrote: > Perhaps I found a solution as: > > datA$type <- gl(1, nrow(datA), label = 'PopA') > datB$type <- gl(1, nrow(datB), label = 'PopB') > dat <- rbind(datA, datB) Or, slightly shorter dat <- make.groups(datA, datB) densityplot(~scores|condition,

Re: [R] [R-sig-Debian] X11 --> Wayland under Ubuntu WAS: Re: RGL crashes

2010-12-09 Thread Alex Mandel
On 12/09/2010 04:02 AM, Rainer M Krug wrote: > On 12/09/2010 12:53 PM, Duncan Murdoch wrote: >> On 09/12/2010 5:52 AM, Matthew Dowle wrote: >>> Weyland is the project to remove X11 from Linux. >>> >>> http://en.wikipedia.org/wiki/Wayland_(display_server) >>> >>> Ubuntu chiefs have said they support

[R] How does Sweave write to a file?

2010-12-09 Thread Arnau Mir Torres
Hello. I need that Sweave writes R output to a file. To do it, I put the following into the foo.Rnw file: \documentclass{article} \usepackage{Sweave} \begin{document} ... <>= x <- rnorm(100) y <- rnorm(100) sink("foo.txt") summary(lm(y~x)) sink() @ ... \end{document} When I run: R CMD Sw

[R] Bivariate kernel density bandwidth selection

2010-12-09 Thread Glen Sargeant
I've been trying to implement bivariate kernel density estimation. For data like mine, function "kde" from package "ks" with bandwidth matrix derived by function "Hscv" seems like a very good choice. Unfortunately, Hscv seems unmanageably slow except for very small sample sizes (up to a few hund

Re: [R] Encoding problem - I fails to read Hebrew text from online

2010-12-09 Thread Matt Shotwell
Tal, It looks like the data you received has HTML special hex characters. That is, 'ש' is just an ASCII HTML representation of a hex character. It's not encoded in a special manner. The trick is to substitute the HTML encoded hex character for its binary representation, or "decode" the character

[R] RODBC sqlSave question

2010-12-09 Thread Felipe Carrillo
Hi: I am getting the following error message when using the sqlSave function: "Error in sqlSave(myDB,myset)" table "myset" already exists" I want this table to be populated with new data everytime I execute sqlSave. The documentation says that: "If the table exists and has the appropriate structur

Re: [R] Trouble Loading doBy and coin Packages

2010-12-09 Thread Adam Carr
Hi Tal: No I have not tried this. I will do it this evening and we'll see what happens. Thanks for the suggestion. Adam From: Tal Galili Cc: r-help@r-project.org Sent: Thu, December 9, 2010 12:29:20 PM Subject: Re: [R] Trouble Loading doBy and coin Packages

Re: [R] lapply getting names of the list

2010-12-09 Thread Joshua Wiley
Hi Sashi, On Thu, Dec 9, 2010 at 9:44 AM, Sashi Challa wrote: > Hello All, > > I have a toy dataframe like this. It has 8 columns separated by tab. > > Name    SampleID        Al1     Al2     X       Y       R       Th > rs191191        A1      A       B       0.999   0.09    0.78    0.090 > abc9

Re: [R] Condional Density Plot from different data

2010-12-09 Thread Doran, Harold
Perhaps I found a solution as: datA$type <- gl(1, nrow(datA), label = 'PopA') datB$type <- gl(1, nrow(datB), label = 'PopB') dat <- rbind(datA, datB) densityplot(~scores|condition, group = type, dat) > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-projec

Re: [R] Reshape Columns

2010-12-09 Thread Henrique Dallazuanna
Try this: c(rbind(ColA, ColB)) On Thu, Dec 9, 2010 at 3:30 PM, Ross, Stephanie wrote: > Hello, > > I have a general formatting question. I have two columns of data: > > ColA <- c("m", "m", "m", "m") > ColB<- c("d","d","d","d") > > And I would like to reorder them into a new column that looks lik

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread Jon Erik Ween
Just for clarification: This all has to do with the Wechsler Memory Scale revised (WMS-R). The task is digit span forward and backwards (digit span total). A raw score is recorded based on subject performance and translated into an age-adjusted standard score in one table and the standard score

Re: [R] Reshape Columns

2010-12-09 Thread Erik Iverson
One way: as.vector(t(cbind(ColA, ColB))) Ross, Stephanie wrote: Hello, I have a general formatting question. I have two columns of data: ColA <- c("m", "m", "m", "m") ColB<- c("d","d","d","d") And I would like to reorder them into a new column that looks like this: ColC<- c("m","d","m","d",

Re: [R] prevent line breaks with sink()

2010-12-09 Thread Greg Snow
You can set options(width= ) with a sufficiently large number and you should not see the wrapping. If you don't want wrapping or the initial [1], then consider using cat instead of print (or implied print), cat can even bypass the need for sink if you just want to send vector contents to a file

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread Jon Erik Ween
Thanks to David and William for helpful comments. I'm not sure if the list will accept attachments, but am trying with this. Z_example.txt is a pared-down sample of the target table in text format. What is needed is a z-score based on the DSF and DSB fields relative to the age field. The second

Re: [R] Barplot with "Independent" Lines Y axis

2010-12-09 Thread Greg Snow
Look at the twoord.plot function in the plotrix package, but be sure to read the note on the help page, then reread it and take its advice if you decide to stick with this type of plot. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.

Re: [R] how to use diff() with different variables?

2010-12-09 Thread Seeliger . Curt
Luca wrote on 12/09/2010 09:38:07 AM: > ... > What I am trying to do is to build another variable fine1 that > should contain the lagged value for "fine", that is: > > xfine fine1 > 1 A 2010-12-09 07:57:33 NA > 2 B 2010-12-09 08:05:00 2010-12-09 07:57:33 > 3 C 2010

[R] prevent line breaks with sink()

2010-12-09 Thread Duncan.Ramsey
Hello. I have a script that writes vectors to a txt file using sink(). But depending on the console window's size, sometimes you will get output like: [1] 38 84 272 287 305 348 411 413 425 426 427 457 513 517 569 583 784 [18] 866 933 996 Instead of: [1] 38 84 272 287 305 348 411 413 425 426

Re: [R] Plotting 3d surfaces

2010-12-09 Thread Barry Rowlingson
On Thu, Dec 9, 2010 at 4:20 PM, Uwe Wolfram wrote: > Dear Subscribers, > > I am using R for quite a while nowadays on Ubuntu 10.04 LTS. I a using R > for doing my statistics. Furthermore I am using it as a tool to generate > the graphics for my publications. > > I am currently working on a project

[R] Reshape Columns

2010-12-09 Thread Ross, Stephanie
Hello, I have a general formatting question. I have two columns of data: ColA <- c("m", "m", "m", "m") ColB<- c("d","d","d","d") And I would like to reorder them into a new column that looks like this: ColC<- c("m","d","m","d","m","d","m","d") Thank you! Stephanie ___

Re: [R] VARMA

2010-12-09 Thread Garten Stuhl
I used estVARXls from dse, but it doesn’t run. Can anybody give me a simple example? model <- estVARXls(c) R gives me: Fehler in x$input : $ operator is invalid for atomic vectors c is a matrix with 2 columns with measured values. Whats my mistake? Best, Thomas 2010/12/9 Giovanni P

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread William Dunlap
> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Jon Erik Ween > Sent: Thursday, December 09, 2010 8:27 AM > To: David Winsemius > Cc: r-help@r-project.org > Subject: Re: [R] set dataframe field value from lookup table > > Sorr

Re: [R] Dynamically build functional expression and taking derivative

2010-12-09 Thread Henrique Dallazuanna
Try this: D(parse(text = foo), "x") On Thu, Dec 9, 2010 at 2:59 PM, Oleksandr Dyklevych wrote: > Hi! > > I'm faced the next problem: > suppose I have a vector of some coefficient and for some purpose from these > coefficients I need dynamically build a functional expression (because, the > numb

[R] lapply getting names of the list

2010-12-09 Thread Sashi Challa
Hello All, I have a toy dataframe like this. It has 8 columns separated by tab. NameSampleIDAl1 Al2 X Y R Th rs191191A1 A B 0.999 0.090.780.090 abc928291 A1 B J 0.3838 0.3839 0.028 0.888 abcnab A1

[R] Condional Density Plot from different data

2010-12-09 Thread Doran, Harold
I'm not certain I am using the lattice plot correctly here. Below is reproducible code. Suppose I have two data frames, such as: set.seed(1234) datA <- data.frame(condition = gl(3, 100), scores = c(rnorm(100), rnorm(100, 1,1), rnorm(100, 2,1))) datB <- data.frame(condition = gl(3, 1000), scores

[R] how to use diff() with different variables?

2010-12-09 Thread Luca Meyer
Hi, I first should say I am new to R. I have searched without success the R-archives to see if I could find an answer to what I am about to ask you. My dataset is like: xfine 1 A 2010-12-09 07:57:33 2 B 2010-12-09 08:05:00 3 C 2010-12-08 20:42:00 ... that is: 'data.frame':

[R] survival: ridge log-likelihood workaround

2010-12-09 Thread Damjan Krstajic
Dear all, I need to calculate likelihood ratio test for ridge regression. In February I have reported a bug where coxph returns unpenalized log-likelihood for final beta estimates for ridge coxph regression. In high-dimensional settings ridge regression models usually fail for lower values of

Re: [R] Odp: 2 plots with different scales in the same graphical window

2010-12-09 Thread Greg Snow
Just be sure to read the note on the help page for twoord.plot, then reread it, then consider alternatives (including the one mentioned in the note), then if you still want to do this, follow the advice in the note. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcar

Re: [R] Trouble Loading doBy and coin Packages

2010-12-09 Thread Tal Galili
I Adam, Have you tried deleting the package files and then reinstalling them from a different CRAN mirror? Tal Contact Details:--- Contact me: tal.gal...@gmail.com | 972-52-7275845 Read me: www.talgalili.com (Hebrew) | www.bios

[R] Plotting 3d surfaces

2010-12-09 Thread Uwe Wolfram
Dear Subscribers, I am using R for quite a while nowadays on Ubuntu 10.04 LTS. I a using R for doing my statistics. Furthermore I am using it as a tool to generate the graphics for my publications. I am currently working on a project which involves nls-fits of three dimensional surfaces such as e

Re: [R] Encoding problem - I fails to read Hebrew text from online

2010-12-09 Thread Tal Galili
I am bumping this question in the hopes that someone might be able to advise. This Hebrew and R business is not as smooth as I had hoped... Thanks, Tal Older massage: On Tue, Dec 7, 2010 at 2:30 PM, Tal Galili wrote: > Hello all, > > # I am trying to read the text in this URL: > u <- > http://

[R] Mixed Model analysis instead of repeated measures ANOVA: how to call lmer or lme instead of aov?

2010-12-09 Thread peter_petersen
Dear R Community! I try to detect differences in the the reaction times of different experimental conditions. The treatment has three levels and was manipulated within each subject. Each subject is tested e.g. 100 times for each treatment. Usually, in psychophysics, data like this is analyzed by

[R] Dynamically build functional expression and taking derivative

2010-12-09 Thread Oleksandr Dyklevych
Hi! I'm faced the next problem: suppose I have a vector of some coefficient and for some purpose from these coefficients I need dynamically build a functional expression (because, the number of elements in the vector of coefficients might differ each time). So, I'm trying to do like this (l

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread David Winsemius
Offlist comments No reply needed .. This is just for emphasis and clarification. On Dec 9, 2010, at 11:19 AM, Jon Erik Ween wrote: David I see how findInterval is a more elegant way of doing 1). I'd need to change the indices in the lookup table, as findInterval(36, c(0, 17, 19, 24, 29, 3

Re: [R] Statistics of count data for differential expression

2010-12-09 Thread Ben Bolker
Vishal Thapar gmail.com> writes: > This is more of a general statistics question that I am not so sure about so > I need some thoughts from you experts out there. I have Chip-Seq read count > data for 2 different treatments (T1 and T2) and my final aim is to calculate > differential expression be

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread David Winsemius
On Dec 9, 2010, at 11:27 AM, Jon Erik Ween wrote: Sorry, I should have included the error I get when using the initial vesion of step 2): Error in `$<-.data.frame`(`*tmp*`, "DSTz", value = list(Age7 = c(-1.55, : replacement has 20 rows, data has 955 In addition: Warning message: In DSTz

[R] Barplot with "Independent" Lines Y axis

2010-12-09 Thread Rodrigo Aluizio
Hi list. I'm plotting pluviometric (Rain) data as a barplot, and then adding the salinity variable to this plot as lines. Obviously as these Y scales are completely different the salinity appears at the lower part of the graph extremely compacted. I need to plot the line at the exactly same area of

[R] Statistics of count data for differential expression

2010-12-09 Thread Vishal Thapar
Hi All, This is more of a general statistics question that I am not so sure about so I need some thoughts from you experts out there. I have Chip-Seq read count data for 2 different treatments (T1 and T2) and my final aim is to calculate differential expression between them. I also have Chip-Seq r

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread Jon Erik Ween
Sorry, I should have included the error I get when using the initial vesion of step 2): Error in `$<-.data.frame`(`*tmp*`, "DSTz", value = list(Age7 = c(-1.55, : replacement has 20 rows, data has 955 In addition: Warning message: In DSTzlook[, 1] == df$DSF + df$DSB : longer object length is

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread Jon Erik Ween
David I see how findInterval is a more elegant way of doing 1). I'd need to change the indices in the lookup table, as > findInterval(36, c(0, 17, 19, 24, 29, 34, 44, 54, 64, 69, 74, 79, 84, 89) ) [1] 6 should be 7, not 6. The age range for the 7th column 35-44. But that's easy. I can't see ho

Re: [R] VARMA

2010-12-09 Thread Giovanni Petris
Package dse does. HTH, Giovanni On Wed, 2010-12-08 at 17:45 +0100, Garten Stuhl wrote: > Hi all, > > > > I want to estimate parameters from a VARMA(p,q)-Modell. > > > > The equations of the model or the model structures is given by: > > > > Xt=beta1+beta2*Xt-1+beta3*Yt-1+epsilon1 > > Y

Re: [R] problem on Matrix package~~

2010-12-09 Thread Douglas Bates
On Thu, Dec 9, 2010 at 5:37 AM, 朱曼 wrote: > Hi,all! > I encountered a problem "not-yet-implemented method for %*% > " when I tried "diagr %*% design_mat", where > "diagr<-diag(vr,100)", vr is a vector with 100 elements, and design_mat is a > sparse matrix with dimension 100*677000 constructed by

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread David Winsemius
On Dec 9, 2010, at 10:51 AM, Jon Erik Ween wrote: Thanks David What I am trying to do is set up a script that assigns z-scores to a large dataframe (2500x300, but has Age in years and test scores as columns.) from a published table of age-corrected standard scores on this cognitive test.

Re: [R] String to array

2010-12-09 Thread Henrique Dallazuanna
Try this: library(gsubfn) strapply("11 - 23", "\\d{1,3}", simplify = as.numeric) On Thu, Dec 9, 2010 at 12:24 PM, Romildo Martins wrote: > Hello, > > how convert x in xarray (numbers)? > > > x > [1] "0 - 13" > > y > [1] "11 - 23" > > z > [1] "220 - 9" > > xarray > [1] 0 13 > > yarray > [1] 11 2

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread David Winsemius
On Dec 9, 2010, at 10:33 AM, David Winsemius wrote: On Dec 9, 2010, at 9:34 AM, Jon Erik Ween wrote: Hi This is (hopefully) a bit more cogent phrasing of a previous post. I'm trying to compute a z-score to rows in a large dataframe based on values in another dataframe. Here's the scrip

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread Jon Erik Ween
Thanks David What I am trying to do is set up a script that assigns z-scores to a large dataframe (2500x300, but has Age in years and test scores as columns.) from a published table of age-corrected standard scores on this cognitive test. 1) The age intervals in the lookup table are given and n

[R] Odp: 2 plots with different scales in the same graphical window

2010-12-09 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 09.12.2010 16:25:06: > Dear useRs, > I have two sets of data that I would like to plot in the same window, but their > ranges are really different, e.g. > > a <- c(0.2, 0.5, 0.8, 0.3, 0.4, 0.5, 0.2, 0.2, 0.3) > b <- c(100, 2000, 800, 200, 100, 50, 4,

Re: [R] set dataframe field value from lookup table

2010-12-09 Thread David Winsemius
On Dec 9, 2010, at 9:34 AM, Jon Erik Ween wrote: Hi This is (hopefully) a bit more cogent phrasing of a previous post. I'm trying to compute a z-score to rows in a large dataframe based on values in another dataframe. Here's the script (that does not work). 2 questons, 1) Anyone know of a

[R] Using Lagsarlm

2010-12-09 Thread Saswati Neogi
I'm trying to use the spdep package to calculate this: y = rho W y + e I don't want to use explanatory variables, just the lag from the dependent variable. How would I code this? __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

[R] 2 plots with different scales in the same graphical window

2010-12-09 Thread Rafael Moral
Dear useRs, I have two sets of data that I would like to plot in the same window, but their ranges are really different, e.g. a <- c(0.2, 0.5, 0.8, 0.3, 0.4, 0.5, 0.2, 0.2, 0.3) b <- c(100, 2000, 800, 200, 100, 50, 4, 2, 0) If I do plot(a, ty='l'); points(b, ty='l') I won't be able to see sequen

Re: [R] error in lrm( )

2010-12-09 Thread David Winsemius
On Dec 9, 2010, at 8:06 AM, 笑啸 wrote: Dear Sir or Madam: I am a doctor of urology,and I am engaged in developing a nomogram of bladder cancer. May I ask for your help on below issue? I set up a dataset which include 317 cases. I got the Binary Logistic Regression model by SPSS.And then

[R] problem using Matrix package

2010-12-09 Thread 朱曼
Hi,all! after executing : library(Matrix) spr_mat<-sparseMatrix(mydataframe$r,mydataframe$c,x=mydataframe$v) diagr<- diag(vr,nrow=ins_num) A<-diagr%*%spr_mat I encountered a error "not-yet-implemented method for %*% " could someone show me some tips? many hanks! best regards! June

Re: [R] String to array

2010-12-09 Thread Jorge Ivan Velez
Try f <- function(string) as.numeric(strsplit(string, "- ")[[1]]) f(x) f(y) f(z) HTH, Jorge On Thu, Dec 9, 2010 at 9:24 AM, Romildo Martins <> wrote: > Hello, > > how convert x in xarray (numbers)? > > > x > [1] "0 - 13" > > y > [1] "11 - 23" > > z > [1] "220 - 9" > > xarray > [1] 0 13 > > ya

Re: [R] String to array

2010-12-09 Thread Keith Jewell
"Romildo Martins" wrote in message news:aanlktinbiaexcobzyqdbtr62xr9q=kjvwaazaqi-k...@mail.gmail.com... > Hello, > > how convert x in xarray (numbers)? > >> x > [1] "0 - 13" >> y > [1] "11 - 23" >> z > [1] "220 - 9" >> xarray > [1] 0 13 >> yarray > [1] 11 23 >> zarray > [1] 220 9 > > > > Than

Re: [R] String to array

2010-12-09 Thread Ben Bolker
Romildo Martins gmail.com> writes: > how convert x in xarray (numbers)? > > > x > [1] "0 - 13" > > y > [1] "11 - 23" > > z > [1] "220 - 9" > > xarray > [1] 0 13 > > yarray > [1] 11 23 > > zarray > [1] 220 9 Is as.numeric(unlist(strsplit("0 - 13","-"))) what you want?

[R] set dataframe field value from lookup table

2010-12-09 Thread Jon Erik Ween
Hi This is (hopefully) a bit more cogent phrasing of a previous post. I'm trying to compute a z-score to rows in a large dataframe based on values in another dataframe. Here's the script (that does not work). 2 questons, 1) Anyone know of a more elegant way to calculate the "rounded" age value t

Re: [R] Complete newbie - create microarray image from data file?

2010-12-09 Thread Ben Bolker
wocket yahoo.co.uk> writes: > > > Hey > > So sorry to be a total newbie, but i'm just finding my feet with R. > > I heard on the grapevine i could recreate a scanned microarray image, or at > least get a good graphic of it from a just a data file. > > I have .txt files for illumina beadarray

[R] String to array

2010-12-09 Thread Romildo Martins
Hello, how convert x in xarray (numbers)? > x [1] "0 - 13" > y [1] "11 - 23" > z [1] "220 - 9" > xarray [1] 0 13 > yarray [1] 11 23 > zarray [1] 220 9 Thanks, RMB __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Calculating odds ratios from logistic GAM model

2010-12-09 Thread David Winsemius
On Dec 9, 2010, at 7:14 AM, denis.ay...@unibas.ch wrote: Dear R-helpers I have a question related to logistic GAM models. Consider the following example: # Load package library(mgcv) # Simulation of dataset n <- 1000 set.seed(0) age<- rnorm(n, 50, 10) blood.pressure <- rnorm(n,

Re: [R] Sequence generation in a table

2010-12-09 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 09.12.2010 13:24:17: > Dear Sirs, > > I understand these already are numeric values. Sir, Basically I am working on > Value at Risk for the Bond portfolio using the historical simulation and for > this I need to find out Marked to Market (MTM) value

[R] error in lrm( )

2010-12-09 Thread 笑啸
Dear Sir or Madam: I am a doctor of urology,and I am engaged in developing a nomogram of bladder cancer. May I ask for your help on below issue? I set up a dataset which include 317 cases. I got the Binary Logistic Regression model by SPSS.And then I try to reconstruct the model (lrm(REC

Re: [R] Sequence generation in a table

2010-12-09 Thread Henrique Dallazuanna
Try this: sapply(sapply(c(257, 520, 110), seq, by = -100), '/', 360) On Thu, Dec 9, 2010 at 10:24 AM, Vincy Pyne wrote: > Dear Sirs, > > I understand these already are numeric values. Sir, Basically I am working > on Value at Risk for the Bond portfolio using the historical simulation and > f

Re: [R] Sequence generation in a table

2010-12-09 Thread Liviu Andronic
On Thu, Dec 9, 2010 at 1:24 PM, Vincy Pyne wrote: > yy <- lapply(c(257, 520, 110), seq, to=0, by=-100) > > yy/360, I get following error. > > Error in yy/360 : non-numeric argument to binary operator > > On the other hand, > > yy[[1]]/365  fetches me > > [1] 0.7138889 0.436 0.158 >

[R] Number of dimension in Multidimensional Scaling

2010-12-09 Thread Petar Milin
Hello! Very often one can hear that MDS usually ends with two-dimensional solution. Of course, there are methods, like Scree-test (proposed by Kruskal and Wish, 1981), to determine optimal number of dimensions. However, I am trying to find references to this two-dimensional gold-standard. Can

Re: [R] Sequence generation in a table

2010-12-09 Thread Vincy Pyne
Dear Sirs, I understand these already are numeric values. Sir, Basically I am working on Value at Risk for the Bond portfolio using the historical simulation and for this I need to find out Marked to Market (MTM) value suing the Present Value of the coupon payments for each Bonds (here as an ex

[R] Complete newbie - create microarray image from data file?

2010-12-09 Thread wocket
Hey So sorry to be a total newbie, but i'm just finding my feet with R. I heard on the grapevine i could recreate a scanned microarray image, or at least get a good graphic of it from a just a data file. I have .txt files for illumina beadarrays but no images cos the service we used didn't send

[R] Calculating odds ratios from logistic GAM model

2010-12-09 Thread Denis . Aydin
Dear R-helpers I have a question related to logistic GAM models. Consider the following example: # Load package library(mgcv) # Simulation of dataset n <- 1000 set.seed(0) age<- rnorm(n, 50, 10) blood.pressure <- rnorm(n, 120, 15) cholesterol<- rnorm(n, 200, 25) sex<-

Re: [R] Sequence generation in a table

2010-12-09 Thread Petr PIKAL
Hi r-help-boun...@r-project.org napsal dne 09.12.2010 12:41:47: > Dear Sir, > > Sorry to bother you again. Sir, the R code provided by you gives me following output. > > > yy <- lapply(c(257, 520, 110), seq, to=0, by=-100) > > yy > [[1]] > [1] 257 157 57 > > [[2]] > [1] 520 420 320 220 120

  1   2   >