Re: [R] Calculating the trading days

Hi Ravi, See if the following helps: = Lines - Date1 Date2 29-Dec-06 25-Jan-07 29-Dec-06 25-Jan-07 29-Dec-06 25-Jan-07 2-Jan-07 25-Jan-07 2-Jan-07 25-Jan-07 2-Jan-07 25-Jan-07 DF - read.table(con- textConnection(Lines), skip = 1) close(con); names(DF) -

Re: [R] get a vector with filenames with a certain extension in a folder

milton ruser wrote: Hi there, try this: setwd(/home/mcr) dir(path=., pattern=.txt$) OR dir(path=/home/mcr/., pattern=.txt$) $ means finished with .txt A word of warning. The pattern is a regular expression, so the above means that the name ends with an arbitrary character followed by

Re: [R] Comprehensive power analysis/sample size package in R?

Many thanks for your helpful suggestions and the detailed feedback! We will have a look at your suggestions before buying the (quite expensive) PASS software. -Udo Quoting Tobias Verbeke tobias.verb...@openanalytics.be: Frank E Harrell Jr wrote: Greg Snow wrote: I don't know of a single

Re: [R] ggplot2: geom_errorbarh()

when I use size, it sets the global size of the whiskers (I'd rather say the thickness of the line) but for sure it does not set the width of the whisker. Try to use the option 'width' with errorbar() not errorbarh(), you'll see that it works perfectly, it sets the width of the whisker. It just

[R] Remembering a value in multiple calls of a function

Hello, I tried this pseudo-generator style experiment (see below). The - operator assigns to in the calling environment which would be the environment of getN. Yet when the function incr is returned, isn't this environment lost? Also the print displays GlobalEnv, yet the globalenv does not have

Re: [R] Renviron.in error when compiling R 2.9.1 under Debian.

Milton, On 16 July 2009 at 18:08, milton ruser wrote: | I'm trying to compile R 2.9.1 on my home directory under debian, | and as I need to play with RPy, I tryed compile R using: To use RPy on Debian (or Ubuntu) you can also just try sudo apt-get install python-rpy which has been

Re: [R] Problems generating image from tiff file

Hi Mehdi, PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. I don't mean to rude, but you seem to lack basic knowledge regarding the use of spatial data in R. Starting to use R is hard, but please ask

Re: [R] Remembering a value in multiple calls of a function

On Fri, Jul 17, 2009 at 9:15 AM, Saptarshi Guhasaptarshi.g...@gmail.com wrote: Hello, I tried this pseudo-generator style experiment (see below). The - operator assigns to in the calling environment which would be the environment of  getN. Yes. Yet when the function incr is returned, isn't

Re: [R] Getting the C-index for a dataset that was not used to generate the logistic model

Professor Harrell, Thanks for your lightning-fast reply. It was extremely helpful, and pointed me exactly to where I needed to go to solve my problem. For others reading, my problem was that I was incorrectly dealing with the validation data. I tried to do this to generate predictions from a

Re: [R] Remembering a value in multiple calls of a function

'startgiven' is in the environment of your function 'incr' (which is what your function 'a' becomes). It might be more transparent to define your function with an enclosing local environment explicitly, for example: b - local({ + startgiven - 0 + function(n) { + if(!missing(n)) { +

[R] how to evaluate character vector within pnorm()

Hi, I'm trying to evaluate a character vector within pnorm. I have a vector with values and names x = c(2,3) names(x) = c(mean, sd) so that i tried the following temp = paste(names(x), x, sep = =) #gives # temp #[1] mean=2 sd=3 #Problem is that both values 2 and 3 are taken as values for

Re: [R] Remembering a value in multiple calls of a function

Thanks to all for the helpful explanation. Regards Saptarshi On Fri, Jul 17, 2009 at 12:42 AM, bill.venab...@csiro.au wrote: 'startgiven' is in the environment of your function 'incr' (which is what your function 'a' becomes).  It might be more transparent to define your function with an

Re: [R] how to evaluate character vector within pnorm()

Dear Thomas, You don't need parse() X - c(mean = 2, sd = 3) pnorm(0, mean = X[mean], sd = X[sd]) Have a look at library(fortunes) fortune(106) HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en

Re: [R] how to evaluate character vector within pnorm()

Here is one way: x - c(2,3) names(x) - c(mean, sd) do.call(pnorm, c(list(0), as.list(x))) Bill Venables http://www.cmis.csiro.au/bill.venables/ -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Thomas Roth (geb. Kaliwe)

[R] Solving two nonlinear equations and two knowns

Dear R users, I have two nonlinear equations, f1(x1,x2)=0 and f2(x1,x2)=0. I try to use optim command by minimize f1^2+f2^2 to find x1 and x2. I found the optimal solution changes when I change initial values. How to solve this? BTW, I also try to use grid searching. But I have no information

Re: [R] Solving two nonlinear equations with two knowns

yhsu6 wrote: Dear R users, I have two nonlinear equations, f1(x1,x2)=0 and f2(x1,x2)=0. I try to use optim command by minimize f1^2+f2^2 to find x1 and x2. I found the optimal solution changes when I change initial values. How to solve this? BTW, I also try to use grid searching.

Re: [R] curve from a formula with ggplot2

Hi, It works perfectly! Though I'm still stuck on the height of the whiskers with errorbarh() :) Regards/Cordialement Benoit Boulinguiez -Message d'origine- De : hadley wickham [mailto:h.wick...@gmail.com] Envoyé : dimanche 12 juillet 2009 11:09 À : Benoit Boulinguiez Cc :

[R] Package norm has been removed. What to use for Maximum Likelihood Missing Data Imputation?

Hello, I apologize if an answer to my questions is available, or if I submitted this question incorrectly. I have read the mailing lists, as well as the R Project and CRAN homepages. However, I may have missed something. I noticed the package 'norm' has been removed. Its page

Re: [R] Package norm has been removed. What to use for Maximum Likelihood Missing Data Imputation?

The package is ORPHANED (i.e. no active maintainer) and does not pass the checks for recent versions of R anymore. If you think the package is very useful for you, you may want to become its new active maintainer. Just take the sources from the archive, fix and update the package as you like

Re: [R] DataFrame help

Don MacQueen wrote: As others noted, you can use the built in function colSums, but you said you're writing your own. Given what you've got so far, that makes the issue one of structuring the output. Try csum - function(m) { a = data.frame(m) s = lapply(a,sum) unlist(s) }

Re: [R] Package norm has been removed. What to use for Maximum Likelihood Missing Data Imputation?

I appreciate the info! I am mostly interested with ML missing data imputation. Is there another R module which has such a function that is well regarded? Thank you. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Package norm has been removed. What to use for Maximum Likelihood Missing Data Imputation?

On Fri, 17 Jul 2009, Daniel Abbott wrote: Hello, I apologize if an answer to my questions is available, or if I submitted this question incorrectly. I have read the mailing lists, as well as the R Project and CRAN homepages. However, I may have missed something. I noticed the package 'norm'

Re: [R] Transformation of data!

On Fri, Jul 17, 2009 at 10:29 AM, Andriy Fetsunfet...@googlemail.com wrote: I want to perform some sort of transformation on all the elements in the matrix I have posted and that I have only presented those 3 elements as an example of how the transformation will affect those 3 elements. Do

Re: [R] ggplot2: geom_errorbarh()

I am still failing to understand exactly what you want to set. Maybe you want to set the height of the whiskers? Hadley On Fri, Jul 17, 2009 at 9:07 AM, Benoit Boulinguiezbenoit.boulingu...@ensc-rennes.fr wrote: when I use size, it sets the global size of the whiskers (I'd rather say the

Re: [R] Package norm has been removed. What to use for Maximum Likelihood Missing Data Imputation?

Achim Zeileis wrote: On Fri, 17 Jul 2009, Daniel Abbott wrote: Hello, I apologize if an answer to my questions is available, or if I submitted this question incorrectly. I have read the mailing lists, as well as the R Project and CRAN homepages. However, I may have missed something. I

Re: [R] Weibull Prediction?

On Jul 17, 2009, at 12:52 AM, Sean Brummel wrote: This conversation is digressing... Here is my question, what does the value from the predict function mean? AS STATED IN MY CODE. Try running the example that I gave. predict(model,type=c(response))[1]=? Look at the help page

Re: [R] SOLVED: how to evaluate character vector within pnorm()

That what i was looking for :-) bill.venab...@csiro.au schrieb: Here is one way: x - c(2,3) names(x) - c(mean, sd) do.call(pnorm, c(list(0), as.list(x))) Bill Venables http://www.cmis.csiro.au/bill.venables/ -Original Message- From: r-help-boun...@r-project.org

[R] pass an array from a perl script to a R script

Hi, All, Many thanks for letting me know there is no references and dereferencing in R, but I saw a library called (ref) on cran, but when I use library(ref), I was told there is no such library in R. So if I can not passing an array from a perl script to a R script by its reference, how to

[R] Count data categories from table

Hi there, I have a relatively simple question, though, I couldn't find a solution for it so far. I have a table with 1000 entries and columns containing information about different parameters for each entry. What I want to do is group all parameters from one of the columns [e.g. if all 1000

Re: [R] Package norm has been removed. What to use for Maximum L

On 17-Jul-09 12:12:01, Uwe Ligges wrote: Achim Zeileis wrote: On Fri, 17 Jul 2009, Daniel Abbott wrote: Hello, I apologize if an answer to my questions is available, or if I submitted this question incorrectly. I have read the mailing lists, as well as the R Project and CRAN homepages.

Re: [R] Count data categories from table

On Jul 17, 2009, at 5:47 AM, Miroslav Nikolov wrote: Hi there, I have a relatively simple question, though, I couldn't find a solution for it so far. I have a table with 1000 entries and columns containing information about different parameters for each entry. What I want to do is group

Re: [R] Package norm has been removed. What to use for Maximum L

Follow-up: On 17-Jul-09 12:38:27, Ted Harding wrote: On a point of information: The licence in question: License: The software may be distributed free of charge and used by anyone if credit is given. It has been tested fairly well, but it comes with no guarantees and the

Re: [R] pass an array from a perl script to a R script

On Jul 17, 2009, at 4:44 AM, xin liu wrote: Hi, All, Many thanks for letting me know there is no references and dereferencing in R, but I saw a library called (ref) on cran, but when I use library(ref), I was told there is no such library in R. An error from that function may indicate

Re: [R] Package norm has been removed. What to use for Maximum L

On Fri, Jul 17, 2009 at 8:54 AM, Ted Hardingted.hard...@manchester.ac.uk wrote: Follow-up: On 17-Jul-09 12:38:27, Ted Harding wrote: On a point of information: The licence in question: License: The software may be distributed free of charge and used          by anyone if credit is given. It

Re: [R] Package norm has been removed. What to use for Maximum L

On 17-Jul-09 13:01:46, Gabor Grothendieck wrote: On Fri, Jul 17, 2009 at 8:54 AM, Ted Hardingted.hard...@manchester.ac.uk wrote: Follow-up: On 17-Jul-09 12:38:27, Ted Harding wrote: On a point of information: The licence in question: License: The software may be distributed free of charge

Re: [R] Count data categories from table

aggregate? Something like this should work although it is not very elegant. mydata - data.frame(aa=rep(letters[1:10],2), bb=rnorm(20, 5,1)) ss - aggregate(mydata[,2],list(aa=mydata$aa), sum) pie(ss[,2])

Re: [R] dse model setup help

I think the problem here is that you seem to be trying to specify a non-innovations form model, for which both Q and R need to be specified, but you have only specified Q. The code guesses whether you are specifying an innovations model or an non-innovations model based on whether you specify

Re: [R] Entire Organization Switching from SAS to R - Any experience?

Kel Lam lamk...@gmail.com wrote: My institute has been heavily dependent on SAS for the past while, and SAS is starting to charge us a very deep amount for license renewal. Since we are a non-profit organization that is definitely not sustainable. The team is brainstorming possibility of

[R] Solving two nonlinear equations with two knowns

---BeginMessage--- My question is as follows Y~N(mu2,sigma2^2), i.e. Y has cdf F2 X generates from: X=F2^{-1}(u)+a*|u-tau|^b*I(u0.75), where u~U(0,1) Given tau=0.75, I want to find a and b such that E(X)-E(Y)=2 and Var(X)/Var(Y)=3 I try optim and grid seaching and get different results,

Re: [R] Solving two nonlinear equations with two knowns

My question is as follows Y~N(mu2,sigma2^2), i.e. Y has cdf F2 X generates from: X=F2^{-1}(u)+a*|u-tau|^b*I(u0.75), where u~U(0,1) Given tau=0.75, I want to find a and b such that E(X)-E(Y)=2 and Var(X)/Var(Y)=3 I try optim and grid seaching and get different results, any solution?

Re: [R] Passing additional arguments through '...'

Thanks everyone for your suggestions, they've been very helpful. I think I've just about solved my problem, but not quite yet. As suggested, I've used sys.call and match.call. Here is what I have so far: fun.tester - function(abc,def,ghi,jkl,...){ mf - match.call(expand.dots = FALSE) sc -

[R] R: extract data.frames from a list

Dear useRs and developeRs, I am struggling with a simple but not obviously solvable issue. Suppose I have the following list of data.frames called 'tmp': a - data.frame(a=rnorm(10),b=letters[1:10]) (tmp - list(a,a[1:4,],a[1:7,])) It is known that all data.frames in this list have the same

[R] Arules questions. I need some help please

Question 2a) I am also working with arules package and I have the following problem let suppose the matrix b like: b-matrix(c(1,1,1,1,1,1,0,0,1,1,1,1,0,0,1,1,0,1,1,1,1,1,1,1),nrow=6) rownames(b)=c(T1, T2, T3, T4, T5, T6) colnames(b)=c(It1, It2, It3, It4) bt-as(b, transactions) rules-apriori(bt,

Re: [R] Plotting hourly time-series data loaded from file using plot.ts

Yep, it's my mistake while manipulating the data. Now it works fine except there is a warning message coming out: Warning message: In zoo(rval, ix) : some methods for “zoo” objects do not work if the index entries in ‘order.by’ are not unique I checked the faq and some other documents, it

[R] How to do adjust for sex, age, genotype for a data

Hello R gurus, I am biologist doing biomarker research and I have a data set where I have 6 proteins and close to 3000 samples, i have to look for differences between disease(Y) and controls(N) along with genetic risk, genotypes, sex and other demographic info available. however i do not know

Re: [R] R: extract data.frames from a list

try this, do.call(rbind, tmp) baptiste 2009/7/17 Angel Spassov anspas...@googlemail.com Dear useRs and developeRs, I am struggling with a simple but not obviously solvable issue. Suppose I have the following list of data.frames called 'tmp': a - data.frame(a=rnorm(10),b=letters[1:10])

Re: [R] R: extract data.frames from a list

do.call(rbind, tmp) On Fri, Jul 17, 2009 at 10:55 AM, Angel Spassov anspas...@googlemail.comwrote: Dear useRs and developeRs, I am struggling with a simple but not obviously solvable issue. Suppose I have the following list of data.frames called 'tmp': a -

Re: [R] How to do adjust for sex, age, genotype for a data

1Rnwb wrote: Hello R gurus, I am biologist doing biomarker research and I have a data set where I have 6 proteins and close to 3000 samples, i have to look for differences between disease(Y) and controls(N) along with genetic risk, genotypes, sex and other demographic info available. however i

Re: [R] Entire Organization Switching from SAS to R - Any experience?

I should elaborate the situation a bit more. We store our data in UNIX and have been using UNIX SAS for our work. My Biostat dept has 40 SAS users from which at most 10 also use R. The Epi/Grad Students/Investigators combine for another 30-40 not-so-frequent SAS users let alone R. So we are

Re: [R] R: extract data.frames from a list

Right on target! Very appreciated. Thanks. Cheers, AS On Fri, Jul 17, 2009 at 4:39 PM, Henrique Dallazuanna www...@gmail.comwrote: do.call(rbind, tmp) On Fri, Jul 17, 2009 at 10:55 AM, Angel Spassov anspas...@googlemail.comwrote: Dear useRs and developeRs, I am struggling with a simple

Re: [R] Plotting hourly time-series data loaded from file using plot.ts

Need reproducible code and data to respond. As per the last line to every message in r-help, reduce the file to the smallest subset that exhibits the error and post it with self contained code that will read it in and reproduce the error message. On Fri, Jul 17, 2009 at 10:29 AM,

Re: [R] Help with Conditional Logit

On Thu, 16 Jul 2009, Noah Silverman wrote: Thanks for the quick reply. I found that page and do have the survival package installed. My problem is that I don't understand a few things: 1) What do I input for the formula field? (In RapidMiner, the functions don't have this input.) 2)

Re: [R] Passing additional arguments through '...'

Thanks for your suggestions, they were quite helpful. I've just about solved the problem, but now I just need to know where the additional arguments passed through '...' are stored (in order to make my own matching scheme I need to be able to change the value of those to the correct values upon

Re: [R] Question regarding package submission to CRAN

It usually gets uploaded within a day if it passes all the checks run by Kurt Hornik. Then the Windows version gets built within the next 24 hours as well. Vik Jim Lemon wrote: Gates, Michael wrote: ... Members, I recently submitted my first package to the submissions ftp site on

Re: [R] Passing additional arguments through '...'

On Fri, 17 Jul 2009, escher2079 wrote: Thanks for your suggestions, they were quite helpful. I've just about solved the problem, but now I just need to know where the additional arguments passed through '...' are stored (in order to make my own matching scheme I need to be able to change the

[R] Fisher's exact test

Hi, I'm trying to run Fisher's Exact test on the data below, but I'm not sure how interpret the data shown. Can someone tell me what this is saying? Looking at the numbers it should be that there's no significant difference between the HDL and LDL, but a p-value of 1 seems high. Is the low value

Re: [R] Passing additional arguments through '...'

Fair enough. I'm not going to post all the code because it's distracting, but the following example function is demonstrative of the problem: exfun-function(...){ print(x) } I'm aware that leaving making the only inputs additional arguments is bad form, but this is merely an example. So here

[R] log error for use on axis

Hi, I am new to R plot. I am trying to scale my y axis in log. When I do this I receive the following error Error in axis(side = side, at = at, labels = labels, ...) : CreateAtVector [log-axis()]: axp[0] = 0 0! In addition: Warning messages: 1: In plot.window(...) : nonfinite axis

Re: [R] Fisher's exact test

James Allsopp-2 wrote: Hi, I'm trying to run Fisher's Exact test on the data below, but I'm not Bound Unbound HDL27 9 LDL 8 2 fisher.test(data) At odds of 3 and 4 with one count =2, even the best test won't help. Dieter -- View this message in context:

Re: [R] Solving two nonlinear equations with two knowns

Kate, This is a difficult problem, mainly because your function is not deterministic. Run the function a couple of times with the same parameter values and you will get a different value each time. Obviously, you cannot optimize such functions in the ususal sense. However, you can rewrite the

Re: [R] Solving two nonlinear equations and two knowns

yhsu6 wrote: I have two nonlinear equations, f1(x1,x2)=0 and f2(x1,x2)=0. I try to use optim command by minimize f1^2+f2^2 to find x1 and x2. I found the optimal solution changes when I change initial values. How to solve this? BTW, I also try to use grid searching. But I have no

[R] c-index validation from Design library

Hi Group, I have a question about obtaining the bias-corrected c-index using validate from the Design library. As an example, consider the example from help page: library(Design) ?validate.lrm n - 1000 age- rnorm(n, 50, 10) blood.pressure - rnorm(n, 120, 15) cholesterol-

Re: [R] Solving two nonlinear equations with two knowns

yhsu6 wrote: ## #R code: mu2=0.4 sigma2=4 tau=0.75 mean.diff=2 var.ratio=3 #Use optim: parameter-function(c) { a-c[1] b-c[2] u-runif(1) Y-qnorm(u,mean=mu2,sd=sigma2) u-runif(1) X-qnorm(u,mean=mu2,sd=sigma2)+a*abs(u-tau)^b*(utau)

Re: [R] Passing additional arguments through '...'

On Fri, 17 Jul 2009, escher2079 wrote: Fair enough. I'm not going to post all the code because it's distracting, but the following example function is demonstrative of the problem: exfun-function(...){ print(x) } I'm aware that leaving making the only inputs additional arguments is bad form,

[R] poly objects as data frame columns

Dear UseRs, I just learnt that the number of columns of a data frame is not always what I thought it to be, and I wonder where I should have learnt about this. Consider the following example: dat - data.frame(X1=1:10, X2=LETTERS[1:10]) ncol(dat) ## evaluates to 2 (of course) dat$X1poly

Re: [R] Fisher's exact test

On 7/17/2009 12:13 PM, James Allsopp wrote: Hi, I'm trying to run Fisher's Exact test on the data below, but I'm not sure how interpret the data shown. Can someone tell me what this is saying? Looking at the numbers it should be that there's no significant difference between the HDL and LDL, but

[R] Data mining, R and MySQL ...

The question is: you have your data in MySQL, build your model in R, but now want to use the model to score your MySQL data on an on-going basis, what to do? MySQL users frequently use R for data mining and to build statistical models. They benefit from the RMySQL package which builds an

Re: [R] How to do adjust for sex, age, genotype for a data

then what will be the other factors needed to be adjusted and whether I should adjust or use them as covariates. Finally how these analysis will be done in R Harrell, Frank E wrote: 1Rnwb wrote: Hello R gurus, I am biologist doing biomarker research and I have a data set where I have 6

[R] dot plot with several points for 2 categories

I'm trying to wean myself off the very limited capabilities of Excel and Oo. Currently, I am trying to make a plot showing several values for 2 categories in a dot blot (see http://www.nabble.com/file/p24538360/Picture%2B1.png Picture+1.png except that the x axis should contain the category not

Re: [R] Entire Organization Switching from SAS to R - Any experience?

On Jul 17, 2009, at 9:57 AM, Kelvin Lam wrote: I should elaborate the situation a bit more. We store our data in UNIX and have been using UNIX SAS for our work. My Biostat dept has 40 SAS users from which at most 10 also use R. The Epi/Grad Students/Investigators combine for another

Re: [R] poly objects as data frame columns

Dataframes are lists. Look at dat with str and you will see that the third column (actually the third list element) is a matrix. It's not hard to find the documentation. If you read the documentation on the help page for data.frame you should see this: If a list or data frame or matrix is

Re: [R] dse model setup help

I thought that I wanted the non-innovations form but wanted R and w(t) to be zero. I thought leaving out R would give that. What I'm trying to do is estimate a transfer function noise model. I noticed that you suggested to another person to use ARMA() and one can estimate the forecast function of

Re: [R] dot plot with several points for 2 categories

On Fri, Jul 17, 2009 at 7:17 PM, jaregisuck...@mpi-cbg.de wrote: I'm trying to wean myself off the very limited capabilities of Excel and Oo. Currently, I am trying to make a plot showing several values for 2 categories in a dot blot (see http://www.nabble.com/file/p24538360/Picture%2B1.png

Re: [R] dot plot with several points for 2 categories

Hi, On Jul 17, 2009, at 1:17 PM, jaregi wrote: I'm trying to wean myself off the very limited capabilities of Excel and Oo. Currently, I am trying to make a plot showing several values for 2 categories in a dot blot (see http://www.nabble.com/file/p24538360/Picture%2B1.png Picture+1.png

Re: [R] Entire Organization Switching from SAS to R - Any experience?

We use SAS and R here (a biostat department and consulting unit), in part because there are some things SAS does that R doesn't. In particular, we use SAS proc nlmixed with custom likelihood functions. R has similar capability but does not allow custom likelihood; the authors say adding it would

Re: [R] log error for use on axis

NDC/jshipman wrote: Hi, I am new to R plot. I am trying to scale my y axis in log. When I do this I receive the following error Error in axis(side = side, at = at, labels = labels, ...) : CreateAtVector [log-axis()]: axp[0] = 0 0! In addition: Warning messages: 1: In plot.window(...)

Re: [R] log error for use on axis

On Jul 17, 2009, at 11:27 AM, NDC/jshipman wrote: Hi, I am new to R plot. I am trying to scale my y axis in log. When I do this I receive the following error Error in axis(side = side, at = at, labels = labels, ...) : CreateAtVector [log-axis()]: axp[0] = 0 0! In addition: Warning

Re: [R] dse model setup help

I reworked this using the innovations form and it seems to work. But none of the parameters were changed after estimation. summary lists: $counts function gradient 1 1 It guess it calculates the function but doesn't optimize the parameters. I just used: tfn.est -

Re: [R] poly objects as data frame columns

David, thanks. Your explanation does not quite fit, though, as it refers to using function data.frame, while I assigned the new column with $-. poly() does return an object of classes poly and matrix, not model.matrix, and handing a poly object to function data.frame does behave like I would

Re: [R] Solving two nonlinear equations with two knowns

Thanks a lot Berend and Ravi!!! I am thinking to use integrate command to get theoretical E(X), E(Y), Var(X), and Var(Y). Then the result will not depend on u anymore. After that, I will use dfsane command in BB package. Kate Original message Date: Fri, 17 Jul 2009 09:51:34 -0700

Re: [R] poly objects as data frame columns

On Jul 17, 2009, at 3:24 PM, Ulrike Grömping wrote: David, thanks. Your explanation does not quite fit, though, as it refers to using function data.frame, while I assigned the new column with $-. poly() does return an object of classes poly and matrix, not model.matrix, But

Re: [R] Entire Organization Switching from SAS to R - Any experience?

I differ with Marc in one way. It is amazing what people can learn when you create an emergency for them to do so. Frank Marc Schwartz wrote: On Jul 17, 2009, at 9:57 AM, Kelvin Lam wrote: I should elaborate the situation a bit more. We store our data in UNIX and have been using UNIX

Re: [R] c-index validation from Design library

Juliet Hannah wrote: Hi Group, I have a question about obtaining the bias-corrected c-index using validate from the Design library. As an example, consider the example from help page: library(Design) ?validate.lrm n - 1000 age- rnorm(n, 50, 10) blood.pressure - rnorm(n, 120, 15)

[R] How do I run the by function on an object of class by???

I wrote a function with some loops and avoided a third loop by using the by function. The output object is an array (or list) of class by that looks like this: StationCover.df Empirical.df$Quadrat: VFFF9-21 Station Shrub Tree Woody Invasive Herb Litter Bare [1,] 0.0 11

Re: [R] poly objects as data frame columns

David Winsemius wrote: On Jul 17, 2009, at 3:24 PM, Ulrike Grömping wrote: David, thanks. Your explanation does not quite fit, though, as it refers to using function data.frame, while I assigned the new column with $-. poly() does return an object of classes poly and matrix,

[R] package to do inverse probability weighting in longitudinal data

Hi there, I have a dataset from a longitudinal study with a lot of drop-out. I want to implement the inverse probability weighting method by Robins 1995 JASA paper Analysis of semiparametric regression models for repeated outcomes in the presence of missing data. Does anyone know if there is

[R] Help with R.oo

Hi! I'm trying to learn about object oriented R, as it seems like it would be very useful. I'm going over an example from the documentation, and I'm very confused: http://www1.maths.lth.se/help/R/R.oo/ [assume you've called library (R.oo)] setConstructorS3(SavingsAccount,

Re: [R] Problems generating image from tiff file

Paul, Thank you very much for your explanation. I have only been using R for GIS purposes for about a week, so my expertise in the field is shallow. Your explanation helped me out and I was able to plot the tif file simply by using the spplot command. The reason I was projecting it was that I

Re: [R] poly objects as data frame columns

On Jul 17, 2009, at 5:25 PM, Ulrike Grömping wrote: David Winsemius wrote: On Jul 17, 2009, at 3:24 PM, Ulrike Grömping wrote: David, thanks. Your explanation does not quite fit, though, as it refers to using function data.frame, while I assigned the new column with $-. poly() does

Re: [R] Help with R.oo

On Fri, Jul 17, 2009 at 11:41 PM, Jim Nemeshnem...@broad.mit.edu wrote: Hi!  I'm trying to learn about object oriented R, as it seems like it would be very useful. I'm going over an example from the documentation, and I'm very confused: http://www1.maths.lth.se/help/R/R.oo/ [assume you've

Re: [R] Help with R.oo

On Sat, Jul 18, 2009 at 3:21 AM, Henrik Bengtssonh...@stat.berkeley.edu wrote: On Fri, Jul 17, 2009 at 11:41 PM, Jim Nemeshnem...@broad.mit.edu wrote: Hi!  I'm trying to learn about object oriented R, as it seems like it would be very useful. I'm going over an example from the documentation,

[R] how to run the R script in background in Windows XP?

Hello ,everyone. I am a fresh user of R. I wrote several several R script in R with the version 2.9.1 in Windows XP. But I only can run these scrpit via R-GUI. Can anyone tell me how to run the scripts in background without openning R-gui.In another word, run the R-script and get the result

[R] apply/return and reuse

Hi, Is it possible to make something like the following code actually work? My goal in this example would be that I'd see results like 1 1 10100 2 10100 10200 3 10200 10300 4 10300 10400 In real usage the function would obviously do a lot more work, but the question I

Re: [R] apply/return and reuse

See ?Reduce On Fri, Jul 17, 2009 at 11:10 PM, Mark Knechtmarkkne...@gmail.com wrote: Hi,   Is it possible to make something like the following code actually work? My goal in this example would be that I'd see results like 1   1   10100 2   10100   10200 3   10200   10300 4   10300  

Re: [R] how to run the R script in background in Windows XP?

Check out: ?Rscript On Fri, Jul 17, 2009 at 10:14 PM, Jie TANGtotang...@gmail.com wrote: Hello ,everyone.  I am a fresh user of R. I wrote several several R script in   R with the version 2.9.1 in Windows XP. But I only can run these scrpit via R-GUI.  Can anyone tell me how to run the

Re: [R] How to do adjust for sex, age, genotype for a data

On Fri, 17 Jul 2009, 1Rnwb wrote: then what will be the other factors needed to be adjusted It is NOT an exaggeration to say that hundreds of research papers, dozens of books, and many dissertations have been written on how to go about answering that question in one context or another.