Re: [R] Non-linear Regression best-fit line

Hi: Perhaps the self-starting functions may be helpful. See ?selfStart. There are self-starting functions for both the logistic and Gompertz models (SSlogis and SSgompertz, respectively). Go through the examples to see how they work. HTH, Dennis On Fri, Jun 17, 2011 at 2:14 PM, Sean Bignami wro

Re: [R] Bartlett's Test of Sphericity

cortest.bartlett() in the psych package. I've never seen a non-significant Bartlett's test. Jeremy On 17 June 2011 12:43, thibault grava wrote: > Hello Dear R user, > > I want to conduct a Principal components analysis and I need to run two > tests to check whether I can do it or not. I found

[R] double integral calculation

a=[0.1,0.2,0.1,0.3,0.4] b=[0.2,0.3,0.1,0.2,0.5] c=[1,1,1,1,1] log(c+a-x*b) where x=unknown scale variable. int=$$log(c+a-x*b)dadb, where $ denotes integral sign. Actually, how could I calculate the integral's approximation? double summation? best, moohwan ___

Re: [R] graphsheet and export.graph equivalents in R?

On Fri, 17 Jun 2011, Michael Karol wrote: R Experts I'm currently using an S+ script of the following format and would like to convert it to R. The script opens a graphsheet with an associated name, plots something (in this case a boxplot) and then exports the contents of the graphsheet of

Re: [R] can this sequence be generated easier?

Here ia an idea that might be useful to adapt fixedSumCombinations <- function(N, terms) if(terms == 1) return(N) else if(terms == 2) return(cbind(0:N, N:0)) else { X <- NULL for(i in 0:N) X <- rbind(X, cbind(i, Recall(N-i, terms-1))) X }

Re: [R] How to join matrices of different row length from a list

No. --- Jeff Newmiller The . . Go Live... DCN: Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Research Engineer (Solar/Batteries O.O#. #.O#. with /Software/Embedded Controllers) .OO#. .OO#. rocks...1k --

[R] can this sequence be generated easier?

I have 'x' variables that I need to find the optimum combination of, with the constraint that the sum of all x variables needs to be exactly 100. I need to test all combinations to get the optimal mix. This is easy if I know how many variables I have - I can hard code as below. But what if I don't

Re: [R] Bartlett's Test of Sphericity

The formula for the chi-square value is: -( (n-1) - (2*p-5)/6 )* log(det(R)) where n is the number of observations, p is the number of variables, and R is the correlation matrix. The chi square test is then performed on (p^2-p)/2 degrees of freedom. So you can compute it by hand. Or you can use t

Re: [R] How to join matrices of different row length from a list

Can this be used in matplot(x,y,)? where x and y have matching rows. However, each column may have different rows or length. Thanks. Sheng -- View this message in context: http://r.789695.n4.nabble.com/How-to-join-matrices-of-different-row-length-from-a-list-tp3177212p3607158.html Sent from t

Re: [R] Conditional Correlation

This question was just answered yesterday in this post: http://r.789695.n4.nabble.com/Correlations-by-subgroups-td3599548.html#a3600553 One solution: x<-c(1,1,1,1,1,2,2,2,2,2) y<-rnorm(10) z<-y+rnorm(10) by(data.frame(y,z),factor(x),cor) HTH, Daniel Mateus Rabello wrote: > > Hi, > How c

[R] Conditional Correlation

Hi, How can I accomplish this in R. Example: I have the following data.frame: data <- data.frame(x=c(1,2,3,4,5,6,5,3,7,1,0,4,8),y=c(1,2,1,2,2,2,1,1,1,2,2,2,2),z=c(5,8,4,3,4,1,6,3,3,6,3,5,7)) Supposing that data$y is a factor, I would like to find the Spearman correlation between data$x and

[R] Vim-R-Plugin issue : Python interface must be enabled to run Vim-R-Plugin

I am trying to get the Vim-R-Plugin to work with gvim and R on Windows 7. When I open a .R file in VIM, it complains and says ""Python interface must be enabled to run Vim-R-Plugin." I have installed pywin32 for python 2.7, and added the follo

Re: [R] sink() and library() calls cause script to stop when called from php

chmod 777 on the destination folder fixed the sink() issue. Still don't know why library() was causing the script to stop though. - Original Message - From: "Yoni Teitelbaum" To: r-help@r-project.org Sent: Friday, June 17, 2011 5:41:33 AM Subject: [R] sink() and library() calls caus

Re: [R] rle on large data . . . without a for loop!

Justin gmail.com> writes: > > I think need to do something like this: > > > > dat<-data.frame(state=sample(id=rep(1:5,each=200),1:3, 1000, > > replace=T,prob=c(0.7,0.05,0.25)),V1=runif(1,10,1000),V2=rnorm(1000)) > > brown bag... ... its friday and im sleepy!... dat<-data.frame(id=rep(1:5,each

[R] Inconsistent results from var.get.nc in RNetCDF

Hello - I am having trouble extracting data from NetCDF data files using RNetCDF. The data files each have 3 dimensions (longitude, latitude, and a date) and 3 variables (latitude, longitude, and a climate variable). Here is some of the output from print.nc for clarity: - dimensions: m

[R] Using MCMC sampling to estimate p values with a mixed model

Hi everyone, Apologies if this is a silly question but I am a student and this is my first time using R so I am still trying to educate myself on commands, models e.t.c I have a mixed model with four dichotomous fixed factors and subject as a random factor (as each person completed four vignettes

[R] graphsheet and export.graph equivalents in R?

R Experts I'm currently using an S+ script of the following format and would like to convert it to R. The script opens a graphsheet with an associated name, plots something (in this case a boxplot) and then exports the contents of the graphsheet of the assigned name to an EMF file. I've bee

[R] Non-linear Regression best-fit line

I am trying to fit a curve to a cumulative mortality curve (logistic) where y is the cumulative proportion of mortalities, and t is the time in hours (see below). Asym. at 0 and 1 > y [1] 0. 0.04853859 0.08303777 0.15201970 0.40995074 0.46444992 0.62862069 0.95885057 1. [10] 1.0

Re: [R] Size argument in sample function

you can do something like this x<-rnorm(20,0,1) time=c(1,2,3) sapply(1:length(time),function(t) sample(x,time[t])) Weidong Gu On Fri, Jun 17, 2011 at 9:59 AM, alfredo wrote: > Hi All, > > I'd like to randomly sample a vector N times, where each successive random > sample increases in size. I ha

[R] Bartlett's Test of Sphericity

Hello Dear R user, I want to conduct a Principal components analysis and I need to run two tests to check whether I can do it or not. I found how to run the KMO test, however i cannot find an R fonction for the Bartlett's test of sphericity. Does somebody know if it exists? Thanks for your hel

Re: [R] rle on large data . . . without a for loop!

Justin Haynes gmail.com> writes: > > I think need to do something like this: > > dat<-data.frame(state=sample(id=rep(1:5,each=200),1:3, 1000, > replace=T,prob=c(0.7,0.05,0.25)),V1=runif(1,10,1000),V2=rnorm(1000)) brown bag... dat<-data.frame(id=rep(1:5,each=200),state=sample(1:3, 1000, replac

[R] rle on large data . . . without a for loop!

I think need to do something like this: dat<-data.frame(state=sample(id=rep(1:5,each=200),1:3, 1000, replace=T,prob=c(0.7,0.05,0.25)),V1=runif(1,10,1000),V2=rnorm(1000)) rle.dat<-rle(dat$state) temp<-1 out<-data.frame(id=1:length(rle.dat$length)) for(i in 1:length(rle.dat$length)){ temp2<-

Re: [R] Size argument in sample function

Hi, There is nothing sacrosanct about taking a random sample of size 1 then 2 then 3 versus just taking one random sample of size 6 (with replacement, if you like) and then just use the first element, elements 2:3 and elements 4:6. HTH, Josh On Jun 17, 2011, at 6:59, alfredo wrote: > Hi All

Re: [R] is this a bug?

On 6/17/2011 2:24 PM, (Ted Harding) wrote: And the extra twist in the tale is exemplified by this mini-version of Albert-Jan's first example: DF<- data.frame(A=c(1,2,3)) DF$B<- c(4,5,6) DF$C<- c(7,8,9) DF # A B C # 1 1 4 7 # 2 2 5 8 # 3 3 6 9 DF$D<- DF["A"]/DF["B"]

Re: [R] is this a bug?

And the extra twist in the tale is exemplified by this mini-version of Albert-Jan's first example: DF <- data.frame(A=c(1,2,3)) DF$B <- c(4,5,6) DF$C <- c(7,8,9) DF # A B C # 1 1 4 7 # 2 2 5 8 # 3 3 6 9 DF$D <- DF["A"]/DF["B"] DF # A B CA # 1 1 4 7 0.25 # 2 2 5 8

Re: [R] is this a bug?

df$varname is a column of df. df["varname"] is a one-column df containing that column. df[["varname"]] is a column of df (same as df$varname). df[,"varname"] is a column of df (same as df$varname). df[,"varname",drop=FALSE] is a one-column df (same as df$varname). df$newVarname <- df["varname"

Re: [R] combining strings

What do you want to happen when both are NA? what do you want to happen if both have values? -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-bounce

[R] is this a bug?

Hello, Is the following a bug? I always thought that df$varname <- does the same as df["varname"] <- > df <- data.frame(weight=round(runif(10, 10, 100)), sex=round(runif(100, 0, 1))) > df$pct <- df["weight"] / ave(df["weight"], df["sex"], FUN=sum)*100 > names(df) [1] "weight" "sex""pct"

Re: [R] Polygon question

Does this do what you want? x <- abs(rnorm(100)) tt <- 1:100 m <- mean(x) par(mfrow=c(2,1)) yy <- c(0,3)# y-limit plot(tt,x,type="l",ylim=yy) abline(h=m) clip(0,100,0,m) polygon( c(1,tt,100), c(m,x,m), col='red' ) -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare

Re: [R] prediction intervals

I am not an expert in time series (that is why I referred you to the task view rather than give my own inexpert opinion). I do remember from a textbook that covered the basics of time series that prediction 2 time points ahead was different from plugging in the next time estimate and predicting

Re: [R] plot the y-axis upside down

That's fantastic! Thanks. > -Original Message- > From: Marc Schwartz [mailto:marc_schwa...@me.com] > Sent: Friday, June 17, 2011 3:47 PM > To: Owen, Jason > Cc: 'r-help@r-project.org' > Subject: Re: [R] plot the y-axis upside down > > On Jun 17, 2011, at 2:31 PM, Owen, Jason wrote: > > >

Re: [R] plot the y-axis upside down

On Jun 17, 2011, at 2:31 PM, Owen, Jason wrote: > Hello, > > I need to create a scatterplot where the y-axis is upside down. If I have > non-negative bivariate data in objects "x" and "y," then the operation > > plot(x, -y) > > gives me the figure I want -- a mirror image of plot(x, y) -- exc

[R] profile plot in R

Hi friends, I have a matrix with following format. group var1 var2 ...varN c1 group1 1.2399 1.4990-1.4829 c2 group4 0.8989 0.7849.1.8933 ... ... c100 group10 . I want to draw a profile plot of each condition c1 to c100, which rows in above matrix and each line represen

Re: [R] Factor labels

?plotmath On Fri, Jun 17, 2011 at 6:32 AM, Hugo Loyola wrote: > Please, can you help me with the following? I want to include a Greek letter > as part of a factor label that I need to use as factor on a xyplot. Probably > this was already solved so can I get the place where I can look for? > > Th

[R] plot the y-axis upside down

Hello, I need to create a scatterplot where the y-axis is upside down. If I have non-negative bivariate data in objects "x" and "y," then the operation plot(x, -y) gives me the figure I want -- a mirror image of plot(x, y) -- except that the y-values (coordinates) are negative, which I don't

Re: [R] Specify ID variable in daisy{cluster}

You need to use hhid as the rownames for housing.cluster rather than including it as a variable in the data.frame: housing.cluster <-data.frame(htypec1, afforcr1, resyrc1, crowdcc1, chprbos1) rownames(housing.cluster) <- hhid Then it will not be included in the cluster analysis but will be used t

[R] observation limits on correlog and mantel?

Hello, I am working with two datasets with ~5,200 observations. Apparently, neither the correlog function (ncf package) nor the mantel function (ecodist package) can handle these datasets. Does anyone have any experience about the maximum number of observations these functions

[R] Server question

Hi A client of mine has asked me to investigate the installation of R-software. Could anyone tell me whether the software works only on a client machine or whether it sits on a server with clients attaching to it? Not immediately clear from the docs. Best Oliver -- Oliver Jone

Re: [R] Porting "unmaintained" packages to post R 2.10.0 era

> -Original Message- > > > You can install to any library, even a temporary one that > does not need to be in .libPaths() of your regular R > installation. So you won't do any harm to your sytsem. > The principal point I was making is that "the manuals" are not very informative on

[R] Errorbars

Hi, Does anyone know how i get rid of the marker point in my error bars? The bars default function (it seems) is to have a dot in the middle of the bar, however i don't want it there. I am using Hmis to draw my errorbars and code is a<-as.vector(tapply(Sporangia,list(Host,Isolate),mean)) b<-as.v

[R] Size argument in sample function

Hi All, I'd like to randomly sample a vector N times, where each successive random sample increases in size. I have realised that the function sample does not take vectors for arguments. For example, x<-rnorm(20,0,1) sample(x,c(1,2,3)) ## will only return one random sample of size 1. The trick

[R] Will R CMD build delete the .Rnw file?

Hi, I used R CMD build to build a package but realized that the .Rnw file was deleted after the .tar.gz file was generated. Is it supposed to be the case? If not, any idea as to why it happened? Thanks, YL __ R-help@r-project.org mailing list https:

Re: [R] R in MAC add many extra <´s

If the problem is that in the history R adds a new newline avery time, you should go to Options and activate "clear input lines" (or something similar). marco t. On Thu, Jun 16, 2011 at 7:38 PM, Jonathan Daily wrote: > Do you mean prompts (the > that indicates R is waiting for input)? > Your iss

Re: [R] Porting "unmaintained" packages to post R 2.10.0 era

> -Original Message- > From: r-help-boun...@r-project.org > [mailto:r-help-boun...@r-project.org] On Behalf Of Uwe Ligges > Sent: 17 June 2011 12:04 > >> See the manual. It tells you > >> > >> R CMD INSTALL --build > >> > >> will generate a binary package. > > ... >> it seems I am not

[R] fitting an ARMA model

Dear all,   I would like to fit an ARMA model, but I'm not sure exactly how to fit it.   Here's an example of the problem.   This is my time variable, hourly data  t <- seq(as.POSIXct("2011-01-01 00:00:00"), as.POSIXct("2011-12-31 23:00:00"), by="hour")   my response y <- rnorm(length(t), 1000, 50

Re: [R] matching by gender and age

getting errors when running this file http://r.789695.n4.nabble.com/file/n3605311/il1ra_status.txt il1ra_status.txt > Age.co<-c(dat2b$Dr_Age) > Sex.co<-as.factor(dat2b$Sex) > casecont.co<-as.factor(dat2b$Self_T1D) > stat.co<-as.factor(dat2b$status) > > m <- matchControls(stat.co~Sex.co+Age.co,

[R] Factor labels

Please, can you help me with the following? I want to include a Greek letter as part of a factor label that I need to use as factor on a xyplot. Probably this was already solved so can I get the place where I can look for? Thanks, Hugo [[alternative HTML version deleted]] __

Re: [R] combining strings

You don't say what happens if both arrays have non-missing entries, but assuming that doesn't happen: > ifelse(is.na(xf),xg,xf) [1] "W" "k" "h" NA "g" "r" "j" NA "v" "d" NA "v" NA "z" "r" "r" "i" -- David L Carlson Associate Professor of Anthropolo

Re: [R] Size argument in sample function

Just use lapply(): set.seed(1) x <- rnorm(20) > lapply(1:3, function(i) sample(x, i)) [[1]] [1] -0.6264538 [[2]] [1] -0.3053884 -2.2146999 [[3]] [1] -2.21469989 -0.30538839 -0.04493361 HTH, Marc Schwartz On Jun 17, 2011, at 12:53 PM, Daniel Malter wrote: > I am sure there is a more elegan

Re: [R] Size argument in sample function

I am sure there is a more elegant version of doing this. But this works: x<-rnorm(20) y<-matrix(1:5) #number of points to sample f<-function(z){sample(x,z)} apply(y,1,f) Just adjust y to your liking. HTH, Daniel alfredo wrote: > > Hi All, > > I'd like to randomly sample a vector N times, whe

Re: [R] issue with strptime

You could loop over the list of reasonable formats, using the first one that works (does not return NA): f <- function (strings, formats) { times <- strptime(strings, formats[1]) i <- 1 while ((i < length(formats)) && any(isBad <- is.na(times))) { i <- i + 1 times[isBad

[R] issue with strptime

Hi everyone, I have bunch of date and time observations in the format %Y-%m-%d %I %M %S %p. I used strptime() to read this format. But the problem is some of the times are in the format of %I %M %p, so for those times, strptime is giving me NA values. For example, strptime(paste("2009-04-08","1

Re: [R] Discrete Test in R

I am looking for a discrete test for testing SNP data. Is ther any common ones on the market for the 2 x 3 SNP table? I don't want to use the Chi-square test. -- Thanks, Jim. [[alternative HTML version deleted]] __ R-help@r-project.org mailing

Re: [R] Multiply list objects

I still do not get the point for what task this expansion of data may be useful, by I guess you want (in this case probably very inefficient, but other can work out how to improve if interested) to insert after k <- lapply(h, function (x) x*0) the lines: for(i in seq_along(k)){ temp <- re

Re: [R] texi2dvi runs but produces no output

Thanks Uwe I have installed R in a path without spaces ( c:\R\R-2.13.0 ..) and ten installed patchDVI I Then rAn the command provided by Duncan on a file known to be correct from the Sweave examples: >example(Sweave) (lots of output) and then >patchDVI::SweavePDF("Sweave-test-1.tex") "C:

Re: [R] Porting "unmaintained" packages to post R 2.10.0 era

On 17.06.2011 14:39, Stephen Ellison wrote: -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Uwe Ligges Sent: 17 June 2011 12:04 See the manual. It tells you R CMD INSTALL --build will generate a binary package. ... it

Re: [R] Rcmdr Plugin error

Dear Erin, On Fri, Jun 17, 2011 at 1:45 AM, Erin Hodgess wrote: > Dear R People: (particularly those who have build Rcmdr Plugin packages): > > I'm building a new Plugin and keep getting the following error: > > > Error in if (is.null(where) || where >= n) rbind(object1, object2) > else if (where

Re: [R] question about split

try this: dat <- data.frame(x = runif(100), fac1 = rep(c("a", "b", "c", "d"), each = 25), fac2 = rep(c("A", "B"), 50)) splits <- c("fac1", "fac2") split(dat, dat[splits]) I hope it helps. Best, Dimitris On 6/17/2011 2:34 PM, carslaw wrote: Dear R-users I seem to be stumped on some

[R] question about split

Dear R-users I seem to be stumped on something simple. I want to split a data frame by factor levels given in one or more columns e.g. given dat <- data.frame(x = runif(100), fac1 = rep(c("a", "b", "c", "d"), each = 25), fac2 = rep(c("A", "B"), 50)) I know

Re: [R] Replacing values without looping

Thanks! na.locf was spot on Very fast approach to this type of problem. -- View this message in context: http://r.789695.n4.nabble.com/Replacing-values-without-looping-tp3602247p3605203.html Sent from the R help mailing list archive at Nabble.com. __

Re: [R] Multiply list objects

Sorry, forgot to quote: Hi, I am trying to use the objects from the list below to create more objects. For each year in h I am trying to create as many objects as there are B's keeping only the values of B. Example for 1999: $`1999`$`8025` B B 8025 8026 8027 8028 8029 802511

Re: [R] rgeos- installation help

It is very likely that your GEOS installation is too old. Are you trying to install rgeos 0.1-8, available since yesterday from CRAN? This should check that GEOS is >= 3.2.2, so perhaps you are trying to install an earlier version of rgeos? You have not provided your full configure report, so it is

[R] sink() and library() calls cause script to stop when called from php

Hello all, I have a php script that makes a call system(PATH_PROJECT_ROOT . '/sh/combineDays.sh ' . $dateString . " " . PATH_PROJECT_ROOT); and within combineDays.sh I run an R script: R --slave --vanilla --args $testDates $BASEDIR $defaultSearchVersion < $BASEDIR"R/combineDays.R"

[R] Hidden Markov Chains - mhsmm

Dear R-Users, I am currently trying my luck with Hidden Markov Chains and use the package mhsmm. I was very "shocked" to see how the initial values for the transition matrix and probability distributions affected the outcomes. The results were extremely different even for small changes in the init

Re: [R] optimization with Sparse matrices

Dube, Jean-Pierre chicagobooth.edu> writes: > > To whom it may concern, > > I am trying to maximize a log-likelihood function using optim. > This is a simple problem with only 18 > parameters. To conserve memory, I am using sparse matrices > (SLAM) for some of the data matrices used in the >

[R] the c implementation of loess

Hi All, I am trying to trace the origin of the current loess implementation in R. The reference mentions that Prof Ripley based it on the 1998 version of dloess. When I look at dloess in http://www.netlib.org/a, the file "changes" mentions dloess was made available in 1992 and that a memory

Re: [R] Rmpi fails to install

Juergen, > I try to install Rmpi as root with install.packages("Rmpi"). > It fails with: > ... I encountered a similar problem with openmpi-1.5.3 and Rmpi_0.5-9.tar.gz. The reason was an unchecked call to dlopen("libmpi.so.0"). At least openmpi-1.5.3 does not provide "libmpi.so.0" anymore. T

Re: [R] Porting "unmaintained" packages to post R 2.10.0 era

On 17.06.2011 12:04, Ben Rhelp wrote: - Original Message From: Uwe Ligges To: Ben Rhelp Cc: r-help@r-project.org Sent: Thu, 16 June, 2011 14:38:12 Subject: Re: [R] Porting "unmaintained" packages to post R 2.10.0 era [...] What about "--binary is deprecated"? What is the correct

Re: [R] Porting "unmaintained" packages to post R 2.10.0 era

Hi Prof Brian, Thank you for your email and for writing MASS. This book is brilliant. - Original Message > From: Prof Brian Ripley > To: Ben Rhelp > Cc: r-help@r-project.org > Sent: Thu, 16 June, 2011 14:48:00 > Subject: Re: [R] Porting "unmaintained" packages to post R 2.10.0 era >

Re: [R] Porting "unmaintained" packages to post R 2.10.0 era

- Original Message > From: Uwe Ligges > To: Ben Rhelp > Cc: r-help@r-project.org > Sent: Thu, 16 June, 2011 14:38:12 > Subject: Re: [R] Porting "unmaintained" packages to post R 2.10.0 era > [...] > > > > What about "--binary is deprecated"? What is the correct way now? > > See the m

Re: [R] Replacing values without looping

Hi > Hi, > > If you truly have an array, this is option that should be much faster > than a loop: > > index <- which(is.na(dat)) > dat[index] <- dat[index - 1] > > the only catch is that when there previous value is NA, you may have > to go through the process a few times to get them all. One

Re: [R] Polygon question

Well you need to recalculate the x values and need to interpolate for the position where you lines cross the m lines Uwe Ligges On 16.06.2011 23:35, Muhammad Rahiz wrote: Hi all, I have the following script which fills the values which are less than the mean of a given timeseries. If

Re: [R] texi2dvi runs but produces no output

On 17.06.2011 09:24, christiaan pauw wrote: Thanks for the advice. I tried to install the package as instructed. Command and results below: install.packages("patchDVI", repos="http://R-Forge.R-project.org";) Installing package(s) into ‘C:/Users/nova/Documents/R/win-library/2.13’ (as ‘lib’ is

Re: [R] Multiply list objects

Since this example is not reproducible (and you have not quuoted any former code) I can only give advice "in principle": 1. Never use 1:length(x) since this will seriously fail if x is a length 0 object. Instead, use seq_along(x) 2. If k is a list, then you probably want to use doubled brackets

Re: [R] problem with missing package

On 16.06.2011 16:19, Eric Elguero wrote: Hi everybody, I just tried to run R on one of my projects but it did not want to run: R version 2.12.1 (2010-12-16) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-pc-linux-gnu (64-bit) Loading req

[R] eRm question

Thank you for all the previous support. I have another question: I am using eRm to run a partial credit model. Is there a way to ask eRm to compute and present to me the probability of a person with a specific ability to achieve each one of the scores 0, 1, 2, 3, or 4 on a specific question with k

Re: [R] prediction intervals

Thank you for your post Greg. Do you have any useful references regarding this variability (papers etc)?   Many thanks. Dave From: Greg Snow To: Dave Evens ; "r-help@r-project.org" Sent: Thursday, 16 June 2011, 21:32 Subject: RE: [R] prediction intervals I don't think that this approach is app

Re: [R] NonParametric Anova

A function for non-parametric multivariate analysis of variance (should do univariate, too, I guess) allowing for interactions (as far as I can tell) is implemented in the anosim() function of the "vegan" package. See also: http://cc.oulu.fi/~jarioksa/opetus/metodi/vegantutor.pdf http://cran.r-

[R] rgeos- installation help

Hello Everybody. i am installing rgeo in R-2.13 on debian lenny. i am getting following error. Can anybody help me why i m not able to compile rgeos. rgeos_misc.c: In function ‘rgeos_hausdorffdistance’: rgeos_misc.c:55: error: ‘GEOSHausdorffDistance_r’ undeclared (first use in this function) rgeo

[R] Hidden Markov Chains - mhsmm

Dear R-Users, I am currently trying my luck with Hidden Markov Chains and use the package mhsmm. I was very "shocked" to see how the initial values for the transition matrix and probability distributions affected the outcomes. The results were extremely different even for small changes in the init

Re: [R] texi2dvi runs but produces no output

Thanks for the advice. I tried to install the package as instructed. Command and results below: > install.packages("patchDVI", repos="http://R-Forge.R-project.org";) Installing package(s) into ‘C:/Users/nova/Documents/R/win-library/2.13’ (as ‘lib’ is unspecified) Warning message: In getDependencie