Re: [R] Large determinant problem

On Sun, 9 Dec 2007, [EMAIL PROTECTED] wrote: > I tried crossprod(S) but the results were identical. The term > -0.5*log(det(S)) is a complexity penalty meant to make it unattractive to > include too many components in a finite mixture model. This case was for a > 9-component mixture. At least up

Re: [R] Large determinant problem

I tried crossprod(S) but the results were identical. The term -0.5*log(det(S)) is a complexity penalty meant to make it unattractive to include too many components in a finite mixture model. This case was for a 9-component mixture. At least up to 6 components the determinant behaved as expected an

Re: [R] Large determinant problem

Hmm, S'S is numerically singular. crossprod(S) would be a better way to compute it than crossprod(S,S) (it does use a different algorithm), but look at the singular values of S, which I suspect will show that S is numerically singular. Looks like the answer is 0. On Sun, 9 Dec 2007, [EMAIL P

[R] Large determinant problem

I thought I would have another try at explaining my problem. I think that last time I may have buried it in irrelevant detail. This output should explain my dilemma: > dim(S) [1] 1455 269 > summary(as.vector(S)) Min.1st Qu. Median Mean3rd Qu. Max. -1.160e+04 0.000e

Re: [R] adjusting "levels" after subset a table

running my.df$my.var <- factor(my.df$my.var) seems to do the job P.J. On Sat, 8 Dec 2007, Milton Cezar Ribeiro wrote: > Dear all, > > I have a data.frame with a factor collumn with about 10 levels. > After extract a subset of this data.frame, by selecting 2 of my 10 levels, > the new data.fram

Re: [R] adjusting "levels" after subset a table

On Sat, 2007-12-08 at 19:26 -0800, Milton Cezar Ribeiro wrote: > Dear all, > > I have a data.frame with a factor collumn with about 10 levels. > After extract a subset of this data.frame, by selecting 2 of my 10 > levels, the new data.frame continue with original number of levels. > How can I adj

[R] adjusting "levels" after subset a table

Dear all, I have a data.frame with a factor collumn with about 10 levels. After extract a subset of this data.frame, by selecting 2 of my 10 levels, the new data.frame continue with original number of levels. How can I adjust it in a manner that when I try levels(my.df) I receive the actualyzed

Re: [R] lm: how to calculate rsquared of the predicted values?

Did you try to install using install.packages()? > install.packages("caret") Warning in install.packages("caret") : argument 'lib' is missing: using '/scratch5/hb/R/R_LIBS/linux/library/' trying URL 'http://cran.cnr.Berkeley.edu/src/contrib/caret_3.08.tar.gz' Content type 'application/x-gzip' l

Re: [R] [OT] A free (as in freedom) replacement for StatTransfer

On Sun, 2007-12-09 at 12:26 +1100, Tim Churches wrote: > Marc Schwartz wrote: > > Perhaps the most notable format that is lacking is the SAS proprietary > > format (not the Transport format), which is not openly published and to > > my knowledge, has not been independently reverse engineered. > >

Re: [R] [OT] A free (as in freedom) replacement for StatTransfer

Marc Schwartz wrote: > Perhaps the most notable format that is lacking is the SAS proprietary > format (not the Transport format), which is not openly published and to > my knowledge, has not been independently reverse engineered. The SAS proprietary dataset and format catalogue structures were su

Re: [R] R + LaTeX formula

On Sun, 2007-12-09 at 01:04 +, Jonas Stein wrote: > Hi, > > what is actually the best method to include R-plots into LaTeX > documents? > At the moment i use > > postscript("myplot.eps", width = 12.0, height = 9.0, horizontal = > FALSE, > onefile = TRUE, paper = "special",encoding = "TeXte

Re: [R] package "growth" ... where is it ?

On Sat, 2007-12-08 at 18:42 -0600, Maura E Monville wrote: > I would like to install the package "growth" as it contains the function > "corgram" and some other presumably useful stuff for time series analysis. > I can see it is in R standard library list: > http://hosho.ees.hokudai.ac.jp/~kubo/Rdo

[R] R + LaTeX formula

Hi, what is actually the best method to include R-plots into LaTeX documents? At the moment i use postscript("myplot.eps", width = 12.0, height = 9.0, horizontal = FALSE, onefile = TRUE, paper = "special",encoding = "TeXtext.enc") plot(foo,bar) dev.off() But it is a bit unhandy to scale late

Re: [R] Adding text outside lattice plot

On 12/8/07, Judith Flores <[EMAIL PROTECTED]> wrote: > Hello, > >I need to add some text in the upper left position, > outside a lattice plot. > > xyplot(x~y) > ltext(locator(1), label='My text') > > doesn't work. library(grid) ltext(grid.locator(), label='My text') should. -Deepayan __

[R] package "growth" ... where is it ?

I would like to install the package "growth" as it contains the function "corgram" and some other presumably useful stuff for time series analysis. I can see it is in R standard library list: http://hosho.ees.hokudai.ac.jp/~kubo/Rdoc/doc/html/packages.html

[R] Adding text outside lattice plot

Hello, I need to add some text in the upper left position, outside a lattice plot. xyplot(x~y) ltext(locator(1), label='My text') doesn't work. I would appreciate any help. Tahnk you, Judith Be

Re: [R] spliting windows 5 frames.

?layout try: layout(rbind(c(1, 2, 3), c(1, 4, 5)), widths=c(2, 1, 1)) layout.show(5) and see if this is what you want. On Dec 8, 2007 3:38 PM, Milton Cezar Ribeiro <[EMAIL PROTECTED]> wrote: > Dear all, > > I need split the output windows into two "frames" of equal size. > On firs

[R] spliting windows 5 frames.

Dear all, I need split the output windows into two "frames" of equal size. On first frame I will print only a image, and on second frame I will print out for graphs (like those automaticaly generated by plot of a glm model. Below follow a visual which I looking for. x11(width=10,height=8) oldp

Re: [R] OT: 3d surfaces with transparency

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Deepayan Sarkar wrote: > On 12/8/07, Ben Bolker <[EMAIL PROTECTED]> wrote: >> -BEGIN PGP SIGNED MESSAGE- >> Hash: SHA1 >> >> >> I would be grateful if anyone had suggestions >> about software that could (1) create 3D surface >> plots, (2) han

[R] time series tests

Citando [EMAIL PROTECTED]: > > Hi all, > > Can anyone clear my doubts about what conclusions to take with the following > outputs of some time series tests: > > > adf.test(melbmax) > > Augmented Dickey-Fuller Test > > data: melbmax > Dickey-Fuller = -5.4075, Lag order = 15, p-value =

Re: [R] OT: 3d surfaces with transparency

On 12/8/07, Ben Bolker <[EMAIL PROTECTED]> wrote: > -BEGIN PGP SIGNED MESSAGE- > Hash: SHA1 > > > I would be grateful if anyone had suggestions > about software that could (1) create 3D surface > plots, (2) handle transparency/alpha blending, > (3) generate output in some vector graphics

[R] level significance

Hi all! I am fitting a (mixed) model with a factor (F) and continuous response and predictor: y~F+F:x (How) can I check the significance of the model at each factor level (i.e. the model could be significant only at one of the levels)? Thank you! ___

[R] time series tests

Hi all, Can anyone clear my doubts about what conclusions to take with the following what puts of some time series tests: > adf.test(melbmax) Augmented Dickey-Fuller Test data: melbmax Dickey-Fuller = -5.4075, Lag order = 15, p-value = 0.01 alternative hypothesis: stationary Warni

[R] OT: 3d surfaces with transparency

-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 I would be grateful if anyone had suggestions about software that could (1) create 3D surface plots, (2) handle transparency/alpha blending, (3) generate output in some vector graphics format that preserved the transparency. I could also live with

Re: [R] Make natural splines constant outside boundary

On Sat, 8 Dec 2007, Charles C. Berry wrote: > Gad Abraham ms.unimelb.edu.au> writes: > >> >> Hi, >> >> I'm using natural cubic splines from splines::ns() in survival >> regression (regressing inter-arrival times of patients to a queue on >> queue size). The queue size fluctuates between 3600 and

Re: [R] lm: how to calculate rsquared of the predicted values?

Thank you Max Kuhn, So I decided to install caret (R version 2.6.0, gcc 4.1.2, caret 3.08, slackware 12 and kernel 2.6.17.11), and I've got a problem: R CMD INSTALL caret * Installing to library '/usr/local/lib/R/library' * Installing *source* package 'caret' ... ** libs //(it goes with no probl

Re: [R] MS Excel Data

On 08/12/2007 1:15 PM, diegol wrote: > > Prof Brian Ripley wrote: >> Why not read the 'R Data Import/Export' manual? It ships with R, or can >> be accessed from http://cran.r-project.org/manuals.html . >> > > Dear Professor, > > I'd like to make a suggestion. > > When I first read some artic

Re: [R] MS Excel Data

Prof Brian Ripley wrote: > > Why not read the 'R Data Import/Export' manual? It ships with R, or can > be accessed from http://cran.r-project.org/manuals.html . > Dear Professor, I'd like to make a suggestion. When I first read some articles ands posts about R's wonderful capabilities, th

Re: [R] FW: R memory management

The line: data. <- c(data., new.data) will eat both memory and time voraciously. You should change it by creating 'data.' to be the final size it will be and then subscript into it. If you don't know the final size, then you can grow it a lot a few times instead of growing it a little lots of

Re: [R] Make natural splines constant outside boundary

Gad Abraham ms.unimelb.edu.au> writes: > > Hi, > > I'm using natural cubic splines from splines::ns() in survival > regression (regressing inter-arrival times of patients to a queue on > queue size). The queue size fluctuates between 3600 and 3900. > > I would like to be able to run predict.

Re: [R] lm: how to calculate rsquared of the predicted values?

> > 1. can I measure de R-squared value between the predicted(by the model) and > real (observed) values.? > 2. Measure the RMSE error . > There is a function in the caret package called postResample that will do that (even if you aren't using resampling). As has been previously noted on this li

Re: [R] [OT] A free (as in freedom) replacement for StatTransfer

On 08/12/07 - 10:29, Marc Schwartz wrote: > > On Sat, 2007-12-08 at 12:24 +0100, Luca Braglia wrote: > > Hello * > > > > does anybody know a free (as in freedom) equivalent for Stat > > Transfer > > > > http://stattransfer.com/ > > > > Thank you > > Bye > > The answer is: it depends. > >

Re: [R] relationship between two factors

Thanks to everyone for their answers. What a helpful community! christopher snow wrote: > I have a dataset with two variables that are factors: > > 1) Decision Making Satisfaction (DMS), values = A - Completely, B - > Mostly, C - Partly, D - Not at all > 2) IT Satisfaction values (ITS), values =

[R] Constraint Partitioning

Dear all, I appreciate to let me know if there is any method in R to find all the possible solutions (non-negative integers) to the following constraint partitioning problem: X_1+X_2+...+X_(2k)=N X_1+X_2<=r_1 X_3+X_4<=r_2 X_5+X_6<=r_3 . X(2k-1)+X(2k)=r_k Thank you,

[R] help for segmented package

Hi, I am trying to find m breakpoints of a linear regression model. I used the segmented package. It works fine for small number of predicators and breakpoints.(3 r.v. 3 points). However, my model has 14 variables it even would not work even for just one breakpoints!. The error message

[R] lm: how to calculate rsquared of the predicted values?

Hi, I've built a linear model using multiple linear regression which leads me a R-squared value of 73.58%. After that, I used this model to predicted some values based on the test data. Now I'm wondering how: 1. can I measure de R-squared value between the predicted(by the model) and real (obse

Re: [R] [OT] A free (as in freedom) replacement for StatTransfer

On Sat, 2007-12-08 at 12:24 +0100, Luca Braglia wrote: > Hello * > > does anybody know a free (as in freedom) equivalent for Stat > Transfer > > http://stattransfer.com/ > > Thank you > Bye The answer is: it depends. If you are looking for a single free open source application that cove

Re: [R] Viewing a large data frame (was How to prevent fix() from converting Dates into numeric)

See ?View for your new question (and please change the subject line when you change the subject: see the R posting guide) You could also have used e.g. format() on the data frame before calling fix() if all you want to do was to view it. On Sat, 8 Dec 2007, Christian Gold wrote: > Thanks. So t

Re: [R] Junk or not Junk ???

Well thank you all, this is solved Safe domains do not work for lists News reader solves the problem but at the cost of needing to "go and look" Univ Wash Computing and Communications was no help And all you all said "make a filter" But I could not find custom filters in Mac Entourage But it was

Re: [R] seq_len

On Dec 8, 2007, at 1:02 AM, Joe W. Byers wrote: > In a post on R-devel, Prof Ripley add the following comment > | > BTW, 1:dim(names)[1] is dangerous: it could be 1:0. That was the > | > motivation for seq_len. > > I use the dim(names)[1] and dim(x)[2] along with length(x) with > varying > lev

Re: [R] How to extract numbers from ANOVA tables?

Dear Lars, In addition to looking at the relevant help pages, a nice thing about R is that the objects are there for you to examine. For example: > mod <- lm(mpg ~ ., data=mtcars) > av <- anova(mod) > names(av) [1] "Df" "Sum Sq" "Mean Sq" "F value" "Pr(>F)" > str(av) Classes 'anova' and

Re: [R] Using R function in Excel

Dear John, The R Workbench (biocep project) comes with a powerful Spreadsheet view fully connected to an R session. you can use it to do exactly what you want and it has much more features : from within the spreadsheet view you can : *import any R Data (numeric, integer, character, logic

[R] [OT] A free (as in freedom) replacement for StatTransfer

Hello * does anybody know a free (as in freedom) equivalent for Stat Transfer http://stattransfer.com/ Thank you Bye __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.

[R] FW: R memory management

Hi, I'm using R to collect data for a number of exchanges through a socket connection and constantly running into memory problems even though task I believe is not that memory consuming. I guess there is a miscommunication between R and WinXP about freeing up memory. So this is the code: f

Re: [R] How to prevent fix() from converting Dates into numeric

Thanks. So there is no solution, other than avoiding fix() and edit()? What would then be the recommended way to make visible and inspect large data.frames (i.e. that are to big for sensibly displaying on the console)? Would I need to write the data to a file and open in a spreadsheet programme?

Re: [R] How to extract numbers from ANOVA tables?

On Sat, 8 Dec 2007, Brinkmann, Lars wrote: > I need to extract numbers (eg. estimates, standard errors etc.) > summarized in an ANOVA table or summary() table for further > calculations. Can somebody help me or refer me to the needed resource on > the web? The Value section of the relevant hel

[R] How to extract numbers from ANOVA tables?

I need to extract numbers (eg. estimates, standard errors etc.) summarized in an ANOVA table or summary() table for further calculations. Can somebody help me or refer me to the needed resource on the web? [[alternative HTML version deleted]]

Re: [R] Junk or not Junk ???

On Thu, 06-Dec-2007 at 04:29PM -0800, Loren Engrav wrote: |> As for news readers |> I found R and R.mac and R.Bio on the sites you recommend, thank you very |> much, they would avoid the individual emails, but then I would have to go |> look at them, which might be Ok |> |> Deluge? Well, there ar