[R] Fit a 3-Dimensional Line to Data Points

Hi Experts, I am new to R, and was wondering how to do 3D linear regression in R. In other words, I need to Fit a 3-Dimensional Line to Data Points (input). I googled before posting this, and found that it is possible in Matlab and other commercial packages. For example, see the Matlab link: http:

Re: [R] Create an R package

Vista works fine for developing packages. You might look at http://batchfiles.googlecode.com which has some batch files for use with Vista and R. On Thu, Jul 24, 2008 at 9:13 PM, FScottDahlgren <[EMAIL PROTECTED]> wrote: > > Hello all! > > I'm new to the forum, and relatively new to R. I've tr

[R] Installation error for RCurl in Redhat enterrpise 5

I am getting the following error while trying to install the RCurl library. I have checked that the curl and the libcurl.so.3 is already installed in the /usr/bin > install.packages("RCurl") --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done trying URL '

Re: [R] simple random number generation

Duncan Murdoch recently (24 July 2008) posted code providing a much more elegant and efficient approach. (Subject line: truncated normal). Also there is apparently a function in the msm package which will do this for you. cheers, Rolf Turner On 25/07/2008, at 2:43

[R] Bug in gap.plot

Hi, all I am trying to make a plot with a axis break and I want the whole plot to be line, not points. However, when I execute the following command half of the graph is points and the other lines. gap.plot(Xdata, Ydata,gap=c(5,6),gap.axis="x",type="l") I think it might be a bug in plotri

Re: [R] simple random number generation

Since the standard normal distribution goes to infinity in both directions, you can't have random normal constrained to +-1.5. You can have _truncated_ standard normal, though, if that's really what you want. +-1.5 is/are the normal quantiles at pnorm(c(-1.5,1.5)). So if we generate 500 uniform

[R] nlminb--lower bound for parameters are dependent on each others

Hello I'm trying to solve two sets of equations (each set has four equations and all of them share common parameters) with nlminb procedure. I minimize one set and use their parameters as initial values of other set, repeating this until their parameters become very close to each other. I have

Re: [R] simple random number generation

Or, as suggested by Duncan Murdoch, qnorm(runif(500,pnorm(-1.5),pnorm(1.5))) --- On Fri, 25/7/08, jim holtman <[EMAIL PROTECTED]> wrote: > From: jim holtman <[EMAIL PROTECTED]> > Subject: Re: [R] simple random number generation > To: "dxc13" <[EMAIL PROTECTED]> > Cc: r-help@r-project.org > Rec

[R] Create an R package

Hello all! I'm new to the forum, and relatively new to R. I've tried searching a couple of times for my question and came up empty handed. This may be due to my lack of knowledge of keywords. I'm planning on putting a package together, and I'm only about half way through. I'm running Windows

Re: [R] simple random number generation

Generate more than you need and then just keep 500 of them: > x <- rnorm(3000) # make sure we have enough > y <- x[(x > -1.5) & (x < 1.5)][1:500] # keep 500 0f them > hist(y) > On Thu, Jul 24, 2008 at 7:35 PM, dxc13 <[EMAIL PROTECTED]> wrote: > > useR's, > > I want to randomly generate 500 numb

[R] Rolling range and regression calculations

Hi everyone, I want to calculate a min and max (i.e. range) on a rolling time frame of 50 periods for a matrix with the number of periods extending along the row direction. So for 300 periods, there will be 6 time frame windows per column, and 6 min max pairs. I then want to o/p these pairs to a

[R] simple random number generation

useR's, I want to randomly generate 500 numbers from the standard normal distribution, i.e. N(0,1), but I only want them to be generated in the range -1.5 to 1.5. Does anyone have a simple way to do this? Thanks, dxc13 -- View this message in context: http://www.nabble.com/simple-random-numb

[R] Relying on R to generate html tables with hyperlinks and pictures ?

Any packages that can do that ? _ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the postin

Re: [R] How to preserve the numeric format and digits ?

paste("A", format(m, scientific=FALSE, trim=TRUE), "B", sep="") -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Daren Tan Sent: Thursday, July 24, 2008 5:47 PM To: [EMAIL PROTECTED] Subject: [R] How to preserve the numeric format and digits ? Instead of

Re: [R] How to preserve the numeric format and digits ?

Like this? options (scipen=100) m <- c(4, 500) paste("A", m, "B", sep="") Regards, Miltinho astronauta brazil On 7/24/08, Daren Tan <[EMAIL PROTECTED]> wrote: > > > Instead of > > > m <- c(4, 500) > > paste("A", m, "B", sep="") > [1] "A4e+08B" "A5e+10B" > > I wa

[R] How to preserve the numeric format and digits ?

Instead of > m <- c(4, 500) > paste("A", m, "B", sep="") [1] "A4e+08B" "A5e+10B" I want "A4" and "A500" __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] ggplot2 help

You might want to take a look at http://had.co.nz/ggplot2/ There are many examples of how to use this function. And, there's a book coming out soon. hope this helps david freedman Edna Bell wrote: > > Hi yet again! > > Thank you for being patient with me. > > Is there a "how to" for ggplot

Re: [R] how to export ".xls" file with colorful cells?

on 07/23/2008 08:29 PM XinMeng wrote: Hi all: I wanna know how to expoort ".xls" file. And in the exported ".xls" file, how to set different cells with different color? Thanks a lot! My best. You don't indicate what operating system you are using and that is important as most of the relevant

Re: [R] Installing R packages in Textmate

Perhaps this has something to do with it... getOption("pkgType") [1] "mac.binary" whereas you have "MacBinary". But binary is the default in OS X, so why specify it at all? Is there some reason why you need to specify both repos and contriburl in order to get ade4? A simple install.p

[R] ANNOUNCE: rapache 1.1.0 release

Hello, First presented at DSC2005, rapache is a project supporting web application development using the R statistical language and environment and the Apache web server. The current release runs on UNIX/Linux and Mac OS X operating systems. http://biostat.mc.vanderbilt.edu/rapache/ rapache

[R] Odp: sum each row and output results

Hi [EMAIL PROTECTED] napsal dne 23.07.2008 19:43:33: > > Hello, > > I have the following data frame (DF): > > V5V5.1V5.2 V5.3 V5.4 V5.5 > 2 -5890.18905 -6019.84665 -6211.06545 -6198.9353 -6616.8677 -6498.7183 > 3 -5890.18905 -6019.84665 -6211.06545

[R] Installing R packages in Textmate

To Whom It May Concern: I need to use various packages in R and I also want to use Textmate. Problem 1: Textmate doesn't seem to recognize the packages that are already installed on R. Problem 2: When I execute the following command: install.packages("ade4", repos="http://cran.r-project.o

[R] Integrating R and Textmate

To Whom It May Concern: I need to use various packages in R and I also want to use Textmate. Problem 1: Textmate doesn't seem to recognize the packages that are already installed on R. Problem 2: When I execute the following command: install.packages("ade4", repos="http://cran.r-project.o

Re: [R] still more dumb questions

On 25/07/2008, at 9:58 AM, Edna Bell wrote: Hello. When I type in "median", I get the following: median function (x, na.rm = FALSE) UseMethod("median") How do I look at the code, please? It seems that median() is ``generic'' (although with only one --- the default --- ``method''.

[R] still more dumb questions

Hello. When I type in "median", I get the following: > median function (x, na.rm = FALSE) UseMethod("median") > How do I look at the code, please? Thanks again, Edna __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help P

Re: [R] incrementing for loop by 2

you can define the 'seq' argument of for() accordingly, e.g., for(i in seq(0, 20, 2)) print(i) I hope it helps. Best, Dimitris Dimitris Rizopoulos Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/33689

Re: [R] Interactive plot

Googling MDI gives a wikipedia entry in the first place which should explain things. David Scott On Thu, 24 Jul 2008, Duncan Murdoch wrote: On 7/24/2008 8:24 AM, Megh Dal wrote: Thanks for this. Can you please explain me what are those MDI and SDI and the diff in laymans language? Not re

Re: [R] NAs - NAs are not allowed in subscripted assignments

This may not be exactly what the OP wanted, but it may be easier to plow the original vectors a and b thru something like > alogic <-rle(!is.na(as.logical(a))) which returns a bunch of lengths and strings of TRUEs and FALSEs . Then play with alogic and blogic (created similarly to alogic), co

Re: [R] R-package install

Hi, Mark, Thank you for your response. No, I don't know how to compile it. Do I need to use unix to compile? After compiling, can I use it in Windows R? Attached is the package. Cindy On 7/24/08, Mark Difford <[EMAIL PROTECTED]> wrote: > > > Hi Cindy, > > >> Hi, I have a R package, which is a .ta

[R] OT: course in research administration

I'm currently enrolled in the distance education masters degree program in statistics at Texas A&M University. I'm hoping the program will help me develop skills useful in advancing my primary care research endeavors. Having a great time so far, learning a lot. I have one semester in which to ta

[R] R2WinBUGS problem

Dear friends - I'm on winXP, R 2.71 - I have with some help dveloped this multivariate normal model, which gives very plausible results in WinBUGS even without any initial values specified. However, when I then try to run the same model via the bugs function in R2WinBUGS with inits specified as

Re: [R] sem & testing multiple hypotheses with BIC

Dear Donald, My apologies for the slow response: I was out of town and taking a break from r-help when you posted your message. The sem() function does not automatically fit alternative models. If you're considering removing all possible combinations of a subset of paths, then you should be able

Re: [R] Should this PDF render correctly without font embedding?

On 25/07/2008, at 8:29 AM, Jeffrey Horner wrote: A professor here at Vanderbilt sent me the following code. Each of the text strings should right justify against the center vertical bar, but because of font issues it doesn't. I understand that there are workarounds, but I was just curious

Re: [R] error with .dll file and RGtk2

On Thu, Jul 24, 2008 at 7:39 AM, paulandpen <[EMAIL PROTECTED]> wrote: > Hi all, > > I am getting the following error message. > > Does somebody know what needs to happen here? > > I have tried re-installing the RGtk2 package and also downloading a .dll > file and installing it in the RGtk2 file f

Re: [R] Should this PDF render correctly without font embedding?

On 24-Jul-08 20:29:09, Jeffrey Horner wrote: > A professor here at Vanderbilt sent me the following code. Each of the > text strings should right justify against the center vertical bar, but > because of font issues it doesn't. I understand that there are > workarounds, but I was just curious if

Re: [R] Should this PDF render correctly without font embedding?

Mac OS X 10.3.5 R 2.7.1 works fine in the quartz window. Screwed up in Preview with the pdf device On Thu, Jul 24, 2008 at 4:29 PM, Jeffrey Horner <[EMAIL PROTECTED]> wrote: > A professor here at Vanderbilt sent me the following code. Each of the text > strings should right justify against th

Re: [R] Question: how to build a subset to do separate calculations

sorry disreguard On Thu, Jul 24, 2008 at 4:45 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > or sapply as I have recently found to be useful. > > > On Thu, Jul 24, 2008 at 4:33 PM, Erik Iverson <[EMAIL PROTECTED]> > wrote: > >> Spencer - >> >> Spencer wrote: >> >>> >>> Dear R Experts, >>> >>> I

Re: [R] Question: how to build a subset to do separate calculations

Erik Iverson wrote: Spencer - Spencer wrote: Dear R Experts, I am trying to create several subsets that I want to use as separate datasets. After separate subsets are created for each country, I want to do some calculations for each. I did the following: data <- data.frame(na.omit(read.

Re: [R] Question: how to build a subset to do separate calculations

or sapply as I have recently found to be useful. On Thu, Jul 24, 2008 at 4:33 PM, Erik Iverson <[EMAIL PROTECTED]> wrote: > Spencer - > > Spencer wrote: > >> >> Dear R Experts, >> >> I am trying to create several subsets that I want to use as separate >> datasets. After separate subsets are creat

Re: [R] R-package install

Hi Cindy, >> Hi, I have a R package, which is a .tar.tar file. And it is said to be >> source code for all >> platforms. ... I am wondering if I can use this package in Windows R. If it is source code you would first need to compile it to binary format before you can use it. Can you do that? W

[R] Should this PDF render correctly without font embedding?

A professor here at Vanderbilt sent me the following code. Each of the text strings should right justify against the center vertical bar, but because of font issues it doesn't. I understand that there are workarounds, but I was just curious if this was consistent across all platforms. On linux

Re: [R] Cairo graphics

> What are Cairo graphics and how to find out about them, please? http://www.google.de/search?hl=de&q=cairo+graphics cu Philipp -- Dr. Philipp Pagel Lehrstuhl für Genomorientierte Bioinformatik Technische Universität München Wissenschaftszentrum Weihenstephan 85350 Freising, Germany ht

Re: [R] Question: how to build a subset to do separate calculations

Spencer - Spencer wrote: Dear R Experts, I am trying to create several subsets that I want to use as separate datasets. After separate subsets are created for each country, I want to do some calculations for each. I did the following: data <- data.frame(na.omit(read.spss("trial.sav"))) at

Re: [R] Subsetting an array by a vector of dimensions

Returning to the old question ... > Is it possible to subset an n-dimensional array by a > vector of n dimensions? On Sat, Jul 12, 2008 at 1:34 AM, Wolfgang Huber <[EMAIL PROTECTED]> wrote: > > do.call("[", list(x,1,2,TRUE)) > [1] 3 9 15 21 > > See also http

[R] SPR experiment: using lmer, transforming data, collinearity, and using a covariable

Dear R users, I try to analyse a self paced reading experiment where I have two fixed variables (Relativiser, Attachment), two random variables (Participant, ItemNbr), and one covariable (BiasValue). The dependant variable is RT, reading time region by region. My aim is to remove the variance i

[R] Question: how to build a subset to do separate calculations

Dear R Experts, I am trying to create several subsets that I want to use as separate datasets. After separate subsets are created for each country, I want to do some calculations for each. I did the following: data <- data.frame(na.omit(read.spss("trial.sav"))) attach (data) Netherlands <-

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

--- On Thu, 7/24/08, Michal Figurski <[EMAIL PROTECTED]> wrote: > From: Michal Figurski <[EMAIL PROTECTED]> > Subject: Re: [R] Coefficients of Logistic Regression from bootstrap - how to > get them? > To: > Cc: "r-help@r-project.org" > Received: Thursday, July 24, 2008, 11:02 AM > Greg and a

[R] Cairo graphics

Hi again. Sorry to be a pest. What are Cairo graphics and how to find out about them, please? Thanks yet again, Edna __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-projec

Re: [R] Dividing by 0

On 25/07/2008, at 5:24 AM, Robert Baer wrote: I'm trying to calculate the percent change for a time-series variable. Basically the first several observations often look like this, x <- c(100, 0, 0, 150, 130, 0, 0, 200, 0) and then later in the life of the variable they're are generally

[R] ggplot2 help

Hi yet again! Thank you for being patient with me. Is there a "how to" for ggplot2, please? I would like to look at it, but have no idea where to start, please. Thanks, Edna Bell __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] factor question

Hi Edna, Because I am "always" subsetting, I keep the following function handy mydata[] <- lapply(mydata, function(x) if(is.factor(x)) x[,drop=T] else x) This will strip out all factor levels that have been dropped by a previous subsetting operation. For novice users of R (though I am not sugge

Re: [R] factor question

Try: fac1[3:9][drop=T] On Thu, Jul 24, 2008 at 4:19 PM, Edna Bell <[EMAIL PROTECTED]> wrote: > Hi! > > Suppose I have a factor: > >> fac1 <- factor(rep(c("dog","cat","tree"),c(2,3,4))) >> fac1 > [1] dog dog cat cat cat tree tree tree tree > Levels: cat dog tree >> length(fac1) > [1] 9 > > No

[R] factor question

Hi! Suppose I have a factor: > fac1 <- factor(rep(c("dog","cat","tree"),c(2,3,4))) > fac1 [1] dog dog cat cat cat tree tree tree tree Levels: cat dog tree > length(fac1) [1] 9 Now I want to get rid of the dogs: > fac1 <- fac1[3:9] > fac1 [1] cat cat cat tree tree tree tree Levels: cat do

Re: [R] NAs - NAs are not allowed in subscripted assignments

On Thu, Jul 24, 2008 at 10:39:34AM -0700, Nordlund, Dan (DSHS/RDA) wrote: [...] > Yes, it does help. I was misunderstanding how logical values are > used for indexing. I assumed incorrectly that a value would be > returned only if the index expression evaluated as TRUE. It would > seem that the

[R] Multinomial logistic regression with non-numerical data.

Hi all, I have a set of data which looks similar to this (only way bigger!) linktyperegister F D a E D b R D c A T d D T a F D b . . . . . . and would like to

Re: [R] Coconut benchmark for R?

I wonder what such a benchmark test would be good for. Relevant and speaking tests in Global Optimization are formulated as models in AMPL, GAMS, GMPK, MPS, etc. Rewriting a model in matrix form can be a tedious manual task and I have not yet seen it. Very few optimization packages in R accept s

Re: [R] [Fwd: Re: Coefficients of Logistic Regression from bootstrap - how to get them?]

Michal Figurski wrote: Thank you Frank and all for your advices. Here I attach the raw data from the Pawinski's paper. I have obtained permission from the corresponding Author to post it here for everyone. The only condition of use is that the Authors retain ownership of the data, and any public

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Michal Figurski wrote: Thank you all for your words of wisdom. I start getting into what you mean by bootstrap. Not surprisingly, it seems to be something else than I do. The bootstrap is a tool, and I would rather compare it to a hammer than to a gun. People say that hammer is for driving na

[R] Problem with GLS dwtest function

Hi, I want to test for independence in my GLS model fitp2, but when I try to use the dwtest function in the lmtest library, I get the error message "Error in terms.default(formula) : no terms component". The model and data set are below. Any suggestions would be really helpful! Thanks a lot in ad

Re: [R] Calculating Betweenness - Efficiency problem

Dear Gabor, I am really sorry about the file attachment. As you might have figured out I am quite new to the forum interaction techniques. I will keep your suggestion in mind. Thanks for taking the time to test the data I sent you. I found out where exactly the problem is. The surprising p

Re: [R] incrementing for loop by 2

on 07/24/2008 07:59 AM Clark, Jennifer H wrote: Dear List-serv members: How can I structure a for loop so that it will increment the counter by two using R? I am just looking for a simple example that shows where the incrementing part should go. Thanks! It would be helpful to have an idea of w

Re: [R] How to get rule number in arules

As always, 'str(rules)' will provide you with a view on the internal structure of 'rules'. By the way, 'rules' is of class 'rules', see the "rules-class" entry in the ARULES help pages, with more explanations on the meaning of attributes. For example, '[EMAIL PROTECTED]' will show support, confid

Re: [R] Help with which()

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 [EMAIL PROTECTED] wrote: | Hi guys, | | Thank you all for looking at this, and I'm sorry - I was making a silly | mistake. | I was using a previous version of the data table in R, while looking at | a new | file in excel (its easier for me to view the

Re: [R] NAs - NAs are not allowed in subscripted assignments

> -Original Message- > From: Gabor Csardi [mailto:[EMAIL PROTECTED] > Sent: Thursday, July 24, 2008 9:59 AM > To: Nordlund, Dan (DSHS/RDA) > Cc: r-help@r-project.org > Subject: Re: [R] NAs - NAs are not allowed in subscripted assignments > > On Thu, Jul 24, 2008 at 09:30:54AM -0700, Nordl

[R] R-package install

Hi, I have a R package, which is a .tar.tar file. And it is said to be source code for all platforms. And the author said it should install on any system (but he didn't know about Windows). I am wondering if I can use this package in Windows R. Thank you, Cindy [[alternative HTML version d

[R] Parallel Processing and Linear Regression

Does anybody have any suggestions regarding applying standard regression packages lm(), hccm(), and others within a parallel environment? Most of the packages I've found only deal with iterative processes (bootstrap) or simple linear algebra. While the latter might help, I'd rather not program th

Re: [R] Dividing by 0

I'm trying to calculate the percent change for a time-series variable. Basically the first several observations often look like this, x <- c(100, 0, 0, 150, 130, 0, 0, 200, 0) and then later in the life of the variable they're are generally no more 0's. So when I try to calculate the percent c

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

What are the arguments against fidelity of this concept to scientific validity? The concept of predictive performance was devised by one of you, biostatisticians - not me! I accept the authoritative view of the person that did it, especially because I do understand it. When I think of it, ex

[R] Dividing by 0

I'm trying to calculate the percent change for a time-series variable. Basically the first several observations often look like this, x <- c(100, 0, 0, 150, 130, 0, 0, 200, 0) and then later in the life of the variable they're are generally no more 0's. So when I try to calculate the percent c

Re: [R] NAs - NAs are not allowed in subscripted assignments

On Thu, Jul 24, 2008 at 09:30:54AM -0700, Nordlund, Dan (DSHS/RDA) wrote: [...] > > > a <- c(rep(seq(1,4),4),NA,NA) > > > b <- c(rep(seq(1,2),7),NA,NA,1,2) > > > > Andreas, > > > > what is wrong with > > > > a[ (a < 2 | a > 3) & b==1 ] <- NA > > > > ? Isn't this what you want? > > [...] > >

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

To quote (or as nearly so as I can) Einstein's famous remark: "Make everything as simple as possible ... but no simpler" Moreover, "as possible" here means "maintaining fidelity to scientific validity," not "simple enough for me to understand." So I don't think a physicist can explain relativist

Re: [R] Help with which()

Ben Bolker wrote: Erik Iverson biostat.wisc.edu> writes: See http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f naw3 duke.edu wrote: I'm using which to find the position of a value in my data table [snip] For example, when I do: where

Re: [R] Help with which()

On 24-Jul-08 15:35:45, Erik Iverson wrote: > See > > http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-the > se-numbers-are-equal_003f That's not an applicable response: The difference between 3573.1 and 3573.15 is far too large for them to be confused by "==". Thr confusion on

Re: [R] NAs - NAs are not allowed in subscripted assignments

> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Gabor Csardi > Sent: Thursday, July 24, 2008 8:39 AM > To: Kunzler, Andreas > Cc: r-help@r-project.org; [EMAIL PROTECTED]; > [EMAIL PROTECTED] > Subject: Re: [R] NAs - NAs are not allowed in subscripte

Re: [R] Help with which()

Erik Iverson biostat.wisc.edu> writes: > > See > > http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f > > naw3 duke.edu wrote: > > I'm using which to find the position of a value in my data table [snip] > For example, when I do: > > > > where1<

Re: [R] What is wrong with this contrast matrix?

Dear Ted, Thanks for your help; is there any straightforward way to construct an orthogonal contrast matrix by hand? Of course contrasts(as.factor(letters[1:6])) would also do the job, but user-defined contrasts would make more sense. Ideally, I would like to have a function that "tests" for

Re: [R] incrementing for loop by 2

Clark, Jennifer H wrote: > Dear List-serv members: > How can I structure a for loop so that it will increment the counter by two > using R? I am just looking for a simple example that shows where the > incrementing part should go. Thanks! > > for (i in seq(, , by=2)) { ... } vQ ___

Re: [R] What is wrong with this contrast matrix?

On 24-Jul-08 15:30:57, Christoph Scherber wrote: > Dear all, > I am fitting a multivariate linear model with 7 response variables and > 1 explanatory variable. > > The following matrix P: > > P <- cbind( > c(1,-1,0,0,0,0,0), > c(2,2,2,2,2,-5,-5), > c(1,0,0,-1,0,0,0), > c(-2,-2,0,-2,2,2,2), > c(-2

Re: [R] NAs - NAs are not allowed in subscripted assignments

> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Kunzler, Andreas > Sent: Thursday, July 24, 2008 7:45 AM > To: [EMAIL PROTECTED] > Cc: r-help@r-project.org; [EMAIL PROTECTED] > Subject: Re: [R] NAs - NAs are not allowed in subscripted assignments >

Re: [R] NAs - NAs are not allowed in subscripted assignments

On Thu, Jul 24, 2008 at 04:45:14PM +0200, Kunzler, Andreas wrote: > Hy Richie, thank you for the quick response. > > Unfortunately my problems hold on. > > Once again: 2 vectors (numeric) including NAs > My intention: I want to replace the values of vector a that are smaller > than 2 and larger

[R] What is wrong with this contrast matrix?

Dear all, I am fitting a multivariate linear model with 7 response variables and 1 explanatory variable. The following matrix P: P <- cbind( c(1,-1,0,0,0,0,0), c(2,2,2,2,2,-5,-5), c(1,0,0,-1,0,0,0), c(-2,-2,0,-2,2,2,2), c(-2,1,0,1,0,0,0), c(0,-1,0,1,0,0,0)) should consist of orthogonal elemen

Re: [R] ggplot question

On Thu, Jul 24, 2008 at 2:23 AM, Williams Scott <[EMAIL PROTECTED]> wrote: > I am trying to do something simple with ggplot. I wish to draw a density > plot split by group, and fill each group with a different colour (and > each with an alpha =0.25). I have tried a number of variations of the > fol

Re: [R] error with .dll file and RGtk2

On Fri, 25 Jul 2008, paulandpen wrote: Hi all, I am getting the following error message. Does somebody know what needs to happen here? I have tried re-installing the RGtk2 package and also downloading a .dll file and installing it in the RGtk2 file folder Error in dyn.load(x, as.logical(lo

Re: [R] Problem with scatterplot3d example

Duncan Murdoch wrote: On 7/24/2008 8:02 AM, Tom La Bone wrote: I tried to run the following example from section 4.1.4 of the "Scatterplot3d - an R package for Visualizing Multivariate Data" vignette and got an error on the part that plots the regression plane: library(scatterplot3d) data

Re: [R] Help with which()

See http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f [EMAIL PROTECTED] wrote: Hi, I'm using which to find the position of a value in my data table, and it is returning the correct position and the position of another number that differs by an ex

[R] Help with which()

Hi, I'm using which to find the position of a value in my data table, and it is returning the correct position and the position of another number that differs by an extra decimal place. For example, when I do: where1<-which(Operons==3573.1,arr.ind=TRUE) it returns the position of that number and

Re: [R] Coefficients of Logistic Regression from bootstrap - how to get them?

Greg and all, Just another thought on bias and variability. As I tried to explain, I perceive this problem as a very practical problem. The equation, that is the goal of this work, is supposed to serve the clinicians to estimate a pharmacokinetic parameter. It therefore must be simple and al

Re: [R] NAs - NAs are not allowed in subscripted assignments

Hy Richie, thank you for the quick response. Unfortunately my problems hold on. Once again: 2 vectors (numeric) including NAs My intention: I want to replace the values of vector a that are smaller than 2 and larger than 3 into NAs only in case vector b equals 1 a <- c(rep(seq(1,4),4),NA,NA) b

Re: [R] How to delete duplicate cases?

this works cno = c(rep(1342,times=3),rep(2568,times=2)) rank= c(.23,.14,.56,.15,.89) df1 = data.frame(cno,rank)[order(cno,rank),] cnou= unique(cno) ind = match(cno,cnou) where = tapply(rank,ind,length) where = cumsum(as.numeric(where)) df1[where,] regards, PF 2008/7

[R] [Fwd: Re: Coefficients of Logistic Regression from bootstrap - how to get them?]

Thank you Frank and all for your advices. Here I attach the raw data from the Pawinski's paper. I have obtained permission from the corresponding Author to post it here for everyone. The only condition of use is that the Authors retain ownership of the data, and any publication resulting from the

Re: [R] How to delete duplicate cases?

on 07/24/2008 09:00 AM Daniel Wagner wrote: Dear R users, I have a dataframe with lot of duplicate cases and I want to delete duplicate ones which have low rank and keep that case which has highest rank. e.g df1 cno rank 1 13420.23 2 13420.14 3 13420.56 4 25680.1

[R] error with .dll file and RGtk2

Hi all, I am getting the following error message. Does somebody know what needs to happen here? I have tried re-installing the RGtk2 package and also downloading a .dll file and installing it in the RGtk2 file folder Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to

[R] Error - unable to load shared library

Hi, I'm trying to use fastICA package but I only get an error message like following. > library(fastICA) Error in dyn.load(file, DLLpath = DLLpath, ...) : unable to load shared library '/usr/lib/R/library/fastICA/libs/fastICA.so': /usr/lib/R/library/fastICA/libs/fastICA.so: undefined symbol:

Re: [R] How to delete duplicate cases?

Daniel - First, use order() to arrange the data.frame into an appropriate format. Then, use duplicated() with the negation operator to get rid of the duplicated values. Daniel Wagner wrote: Dear R users, Â I have a dataframe with lot of duplicate cases and I want to delete duplicate ones

Re: [R] Weighted variance function?

On Thu, 2008-07-24 at 02:25 +0530, Arun Kumar Saha wrote: > There is a R function to calculate weighted mean : weighted.mean() under > stats package. Is there any direct R function for calculating weighted > variance as well? Here are two ways; weighted.var() is via the usual formula and weighted.

Re: [R] How to delete duplicate cases?

Try this: aggregate(x$rank, list(cno=x$cno), max) On 7/24/08, Daniel Wagner <[EMAIL PROTECTED]> wrote: > Dear R users, > > I have a dataframe with lot of duplicate cases and I want to delete duplicate > ones which have low rank and keep that case which has highest rank. > e.g > > > df1 > cno

Re: [R] How to delete duplicate cases?

Dear Daniel, Try this: x=read.table(textConnection("cno rank 1 13420.23 2 13420.14 3 13420.56 4 25680.15 5 25680.89"),header=TRUE,sep="") x[cumsum(tapply(x$rank,x$cno,which.max)),] cno rank 3 1342 0.56 5 2568 0.89 HTH, Jorge On Thu, Jul 24, 2008 at 10:00 AM, D

Re: [R] how know the status of particular process

Kurapati, Ravichandra (Ravichandra) wrote: Hi I have the process id, How can I know whether it is in active state or not This is not an R question. Find out how to do it in your operating system (whatever that is), and do the same in R using system(). Duncan Murd

[R] How to delete duplicate cases?

Dear R users,   I have a dataframe with lot of duplicate cases and I want to delete duplicate ones which have low rank and keep that case which has highest rank. e.g   > df1   cno      rank 1  1342    0.23 2  1342    0.14 3  1342    0.56 4  2568    0.15 5  2568    0.89

[R] incrementing for loop by 2

Dear List-serv members: How can I structure a for loop so that it will increment the counter by two using R? I am just looking for a simple example that shows where the incrementing part should go. Thanks! Jennifer -- Jennifer Hayes Clark Assistant Professor Department of Political Science Un

  1   2   >