[R] Help with script

there is only 1 ID and the values of each columns A...G are added. I n the example the new table would have only 3 IDs. a1, a2 and 3 and a2 has the values added by column: a2   0   2   1   3   2   4   2.. Thank you so much and have a wonderful year!. -- Pablo A. Ortiz-Pineda (Ph.D.) Molecular

[R] GLMM for Combined experiments and overdispersed data

I am analyzing data from 3 field experiments (farms=3) for a citrus flower disease: response variable is binomial because the flower can only be diseased or healthy. I have particular interest in comparing 5 fungicide spraying systems (trt=5). Each farm had 4 blocks (bk=4) including 2 trees as su

[R] Course: "Statistical Data Analysis with R"

the course. Below you find more information about this course. If you have any question don't hesitate to contact me. Best regards, Pablo ++ Course: Introduction to Statistical Data Analysis with R Where: Linux Hotel, Essen-Horst, Germ

[R] short course: Bayesian Data Analysis with R and WinBUGS

lides), an excursion to the Christmas' market of Essen and another excursion to one historical industrial area. If you have any question doesn’t hesitate to contact me. Best regards, Pablo ++ *Three days course in: Bayesian Data Anal

Re: [R] R lattice bwplot: Fill boxplots with specific color depending on factor level

Thank you both. I just wanted to point out that before assigning the order of colors in vector col, one should check that it corresponds with how levels are ordered in levels(mydata$Col3). Best, Pablo 2015-04-23 5:46 GMT+02:00 Richard M. Heiberger : > Pablo, > > I would do it simi

[R] R lattice bwplot: Fill boxplots with specific color depending on factor level

YOU REPRODUCE THE EXAMPLE YOU WILL SEE THAT THE COLORS ARE NOT RELATED TO THE LEVELS IN MY DATAFRAME AS YY BOX IS NOT ALWAYS GREEN. IS THERE A WAY TO ASSIGN YY:green, NN:red, YN:amber? You can see the resulting figure in: http://stackoverflow.com/questions/29802129/r-lattice-bwpl

[R] citation of a task view

Hi all, Which is a formal bibliography citation of an R's task view? For example if I want to make a citation of "MetaAnalysis" task view. Thanks in advance! Pablo __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

Re: [R] empty density plot for point pattern in spatstat

Hi Chris, Try the form: plot(dens, useRaster=FALSE, ribargs=list(useRaster=FALSE)) Pablo 2014-06-20 15:58 GMT-05:00 Christopher W Ryan : > I'm back in the office with the machine that was giving me trouble. > > # fresh start-up of R 3.1.0, installed on the my machine's h

Re: [R] empty density plot for point pattern in spatstat

Dear Adrian, Yes, effectively utilizanddo useRaster = FALSE, the plot is printed. Only the ribbon appears without color. I'm using R version 3.0.2, Windows 7 system, and spatstat 1.33-0 package version. Best wishes Pablo 2014-06-18 19:50 GMT-05:00 Adrian Baddeley : > Dear Pablo &g

Re: [R] empty density plot for point pattern in spatstat

and perform the plot, the graph is obtained without problem. Best wishes Pablo 2014-06-17 16:23 GMT-05:00 Rolf Turner : > On 18/06/14 03:49, Christopher W Ryan wrote: > >> I upgraded to version 3.1.0 of R, still on Win XP, but that did not >> solve the problem: >> >

[R] Course: Introduction to Bayesian Data Analysis using WinBUGS and R

Dear list members, Apologies for cross-posting. Please, find below the information of an introductory Bayesian Data Analysis course with R and WinBUGS. If you have any question don't hesitate to contact me. Best regards, Pablo ++ C

[R] RSA and lmer

. Thanks in advance! best, Pablo __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] multilevel sampling weight

Anyone is working with the upgrade of any multilevel package for apply sampling weights? Thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gu

[R] course: Bayesian Data Analysis with R and WinBUGS

Dear list members, Best greetings and apologies for cross-posting. There are available places for the course: "Bayesian Data Analysis with R and WinBUGS", please find the description of the course below. If you have any question don't hesitate to contact me. Best

Re: [R] foreign package problem with stata 13

dataset in stata. > > I don't think foreign library in R supports stata 13 yet. > > Richard > > On Wed, Oct 2, 2013 at 3:11 PM, Pablo Menese Camargo > wrote: > > foreign package does not support dataset saved at stata 13. > > anyone knows any wayt to make

[R] foreign package problem with stata 13

foreign package does not support dataset saved at stata 13. anyone knows any wayt to make it works? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] SEM polychoric lavaan

I have this mcfa<- "impulsivity=~imp1r+imp2+imp3+imp4 physical=~phys1+phys2+phys3+phys4 risky=~risk1+risk2+risk3+risk4 selfish=~self1+self2+self3+self4 simple=~simp1+simp2+simp3+simp4r temper=~temp1+temp2+temp3+temp4 control=~impulsivity+physical+risky+selfish+simple+temper" where: the model is a

Re: [R] SEM polychoric

I'm using lavaan indeed On Mon, Aug 12, 2013 at 6:21 AM, Jose Iparraguirre < jose.iparragui...@ageuk.org.uk> wrote: > Hi Pablo, > > Prof John Fox is a frequent contributor to this forum, and no-one is > better placed than him to respond on your particular point, but I

[R] SEM polychoric

I have this mcfa<- "impulsivity=~imp1r+imp2+imp3+imp4 physical=~phys1+phys2+phys3+phys4 risky=~risk1+risk2+risk3+risk4 selfish=~self1+self2+self3+self4 simple=~simp1+simp2+simp3+simp4r temper=~temp1+temp2+temp3+temp4 control=~impulsivity+physical+risky+selfish+simple+temper" where: the model is a

[R] eRm package, RSM: Error in solve.default(parets$hessian) :

I want to make an index using IRT, I use eRm package. I have a data frame (24 items, 4 categories, 427 cases) that works if I use PCM (partial credit model), but, when I run RSM (rating scale model) this appear: Error in solve.default(parest$hessian) : Lapack routine dgesv: system is exactly sin

[R] Courses: Statistical Analysis with R - Bayesian Data Analysis with R and WinBUGS

Dear list members, Apologies for cross-posting. Please, find below the information of two statistical courses with R: 1) Statistical Analysis with R 2) Bayesian Data Analysis with R and WinBUGS If you have any question don't hesitate to contact me. Best regards,

[R] Confirmatory Factor Analysis

Hi, I'm trying to perform a hierachical, second order CFA. That's the thing that I need to leave AMOS. I found some sim.hierarchical and omega stuff, but nothing clear enough. Can anyone help me? I just need a simple and clear manual. Best, Pablo. [[alternative HTML versi

Re: [R] issue creating a subset

you create the matrix using > cbind, as in my matrix3 example. > > -Don > > > -- > Don MacQueen > > Lawrence Livermore National Laboratory > 7000 East Ave., L-627 > Livermore, CA 94550 > 925-423-1062 > > > > > > On 2/28/13 5:53 PM, "Pablo M

[R] issue creating a subset

I'm performing item response theory with eRm packages I am excluding the persons that doesn't fit in the infit/outfit persons. for that I created a condition. then I have to create a new subset or matrix but with the condition. So: ORIGINAL matrix<-cbind(item1, item2, item3, item4) IF I PERFORM

Re: [R] Issue attaching a dataset

thanks a lot. works!!! On Tue, Feb 26, 2013 at 7:58 AM, Prof Brian Ripley wrote: > On 26/02/2013 09:55, Milan Bouchet-Valat wrote: > >> Le lundi 25 février 2013 à 20:26 -0300, Pablo Menese a écrit : >> >>> I use to work whit stata dataframe, so, when I use R I type re

[R] Issue attaching a dataset

I use to work whit stata dataframe, so, when I use R I type read.dta Until today I do that without any problem, after type: mydata<-read.dta("C:/dropbox/dataframe.dta") attach(mydata) Everything works great... but today, when I typed: mydata<-read.dta("C:/dropbox/dataframe.dta") attach(mydata)

Re: [R] problem with svyglm

81 .0410122 .4743798 .6359363 So, I thing that the problem is in the survey design... On Tue, Nov 27, 2012 at 11:49 PM, David Winsemius wrote: > > On Nov 27, 2012, at 2:31 PM, Pablo Menese wrote: > > Sorry, it send it alone... >> >> When I use it: >> >

Re: [R] Problems with weight

: > Le mardi 27 novembre 2012 à 18:33 -0300, Pablo Menese a écrit : > > I can't ... I don't know why but I can't > > > > When I use it: > > > > logit <- glm(bach ~ egp4 + programa, weight=wst7, > > family=quasibinomial(link"logit")) >

Re: [R] Problems with weight

> > http://journal.r-project.org/archive/2009-2/RJournal_2009-2_Damico.pdf > > once you're ready to run the regression, use svyglm() instead of glm() and > drop the weights argument (since it will already be part of the survey > design) :) > > > > On Fri, Nov 23, 2012

Re: [R] problem with svyglm

Sorry, it send it alone... When I use it: logit <- glm(bach ~ egp4 + programa, weight=wst7, family=quasibinomial(link"logit")) I reach the same betas that in STATA, but the hypothesis test, the t value, and the std. error is different. I think that the solution can't be so far from this...

Re: [R] problem with svyglm

I colud not, even without attach the dataset. The thing is, when I use this On Fri, Nov 23, 2012 at 5:56 PM, David Winsemius wrote: > > On Nov 23, 2012, at 12:08 PM, Pablo Menese wrote: > > I have this problem. >> >> test <- svydesign(id=~1,weights=~peso) >&g

[R] Problems with weight

Until a weeks ago I used stata for everything. Now I'm learning R and trying to move. But, in this stage I'm testing R trying to do the same things than I used to do in stata whit the same outputs. I have a problem with the logit, applying weights. in stata I have this output . svy: logit bach job

[R] problem with svyglm

be in design= argument I don't know what this mean... Please help. Pablo. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-

[R] Introductory R course: 6-8.12.2012

Dear list members, Apologies for cross-posting, there are some places available in an introductory R course. Please, find below the information of this R training course. If you have any question, don't hesitate to contact me. Best regards,

Re: [R] crosstable and regression for survey data (weighted)

gn=dstrat) is correct, using only the var names, not dstrat[smth]~dstrat[smth]+ dstrat[smth] If you write the names correctly it should work regards pablo [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.et

Re: [R] svyglm question (haps)

It seems the var names you've put are not the same as in the design object: "all variables must be in design= argument ": that means the object you've assigned in mydesign <- svydesign(ids=~vill_neigh_code+clust, strata=~strat, weights=~sweight, data=mydata) Check the spelling. Note that the "myd

Re: [R] crosstable and regression for survey data (weighted)

Regarding regression models, there's a bit of discussion on whether or not it is necessary to take the sample design into account (for instance, SPSS doesn't), so you can run them just normally without much remorse. Or get your life complicated (see below). Your xtabs call seems OK to me. However,

[R] COURSE: Introduction to Bayesian Data Analysis using WinBUGS and R, 21-23 Juni 2012

Dear list members, Apologies for cross-posting. Please, find below the information of an introductory Bayesian Data Analysis course with WinBUGS and R. If you have any question don't hesitate to contact me. Best regards, Pablo ++ C

[R] Stepwise in lme

I'm manually doing a form of stepwise regression in a mixed model but with many variables, it is time consuming. I thought I'd try to use an automated approach. stepAIC gave me false convergence when I used it with my model, so I thought it can't be hard to set up a basic program to do it based o

Re: [R] Summing x1 to x6

;- laply(1:6, function(a) get(paste('x', a, sep = ''))) rowMeans(x1.6) It is nice because it's easily extendable to matrices too which I want to do later - change rowMeans(x1.6) to colMeans(x1.6). Thank you all for your help, Pablo -- View this message in context: http

[R] Summing x1 to x6

ector(rnorm(10)) I would like the mean of x1 to x6 for each vector position. I would do something else with x7-x10. These vectors are not currently in a dataframe. I tried to be clever by trying get(paste(paste("x", 1:6, sep=""), collapse="+")) but it didn't

Re: [R] References for book "R In Action" by Kabacoff

The references are here: http://manning.com/kabacoff/excerpt_references.pdf (they will be included on the next printing too, got omitted by mistake) Regards, Pablo [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] References for book "R In Action" by Kabacoff (Uwe Ligges)

The references are here: http://manning.com/kabacoff/excerpt_references.pdf (they will be included on the next printing too, got omitted by mistake) regards, Pablo [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] factors in probit regression

stical package, it is possible to change the design so there is no "baseline" and the effects can be tested. I don't need a economy book, I'm an ecologist. Please, think before post anything like your mail. Pablo PS. David, many thanks for your help. __

[R] Changing x-axis when dealing with time

0:00:00") t.start t.vector <- t.start + (0:24)*3600 t.vector z.vector <- 1:length(t.vector) z.vector ##vector of decimal hours dh.vector <- decimal.hour(t.vector) dh.vector plot(x=dh.vector, y=z.vector, xlim=c(0,24), xaxt="n", xlab="Hour of day", ylab="Some

Re: [R] selecting data from table with timestamp

which has basically nothing to do with extracting a subset of rows from a > dataframe. > > -Don > > -- > Don MacQueen > > Lawrence Livermore National Laboratory > 7000 East Ave., L-627 > Livermore, CA 94550 > 925-423-1062 > > > > > > -Original

[R] selecting data from table with timestamp

13.35 17.99 4 2011-03-15 00:00:30 1 13.35 18.00 5 2011-03-15 00:01:00 2 13.35 17.98 6 2011-03-15 00:01:30 3 13.35 17.99 7 2011-03-15 00:02:00 4 13.35

[R] course: Introduction to Bayesian Data Analysis using WinBUGS and R; Summer School 2011 - Portugal

Dear Colleges, Apologies for cross-posting. Places are available in the following course: "Introduction to Bayesian Data Analysis using WinBUGS and R" Best Regards, Pablo ++ Course: Introduction to Bayesian Data Analysis usi

Re: [R] 1 continuous non-normal variable ~ 4 factors + 1 continuous covariate (with interactions)

in the package survival. Cheers, Pablo Thanks, Al -- View this message in context: http://r.789695.n4.nabble.com/1-continuous-non-normal-variable-4-factors-1-continuous-covariate-with-interactions-tp3444378p3444378.html Sent from the R help mailing list archiv

Re: [R] Converting a categorical variable to multiple dichotemousvariables

Hi Shane, An alternative is: cat<-as.factor(c(1,1,3,2,4)) outer(cat, levels(cat), "==")+0 Cheers, Pablo - Original Message - From: "andrija djurovic" To: "Shane Phillips" Cc: Sent: Tuesday, April 12, 2011 2:32 PM Subject: Re: [R] Converting a

Re: [R] Testing equality of coefficients in coxph model

;- survreg(Surv(Time, Status) ~ -1 + tr, data = dat[Ind,]) b <- coef(mod1) delta <- b[1] - b[2] # statistic of interest delta } # Bootstraping ... library(boot) boot.res <- boot(test1, boot.fun, R = 1999) boot.res boot.ci(boot.res) plot(boot.res, jack = TRUE) # end ... Cheers, Pablo

[R] Variance of random effects: survreg()

I-likelihood = -95.7 ... In this case I expected that the variance of random effects is close to 4. Thanks! Pablo [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo

[R] Course: "Introduction to Bayesian Data Analysis using WinBUGS and R"

Apologies for cross-posting. Places are available in the following course: "Introduction to Bayesian Data Analysis using WinBUGS and R" Best Regards, Pablo ++ Course: Introduction to Bayesian Data Analysis using WinBUGS an

[R] How to convert ARMA process to infinite AR?

Hi, I need a function to convert an ARMA process to an "infinite" AR process. I know that in the Stats package exist one to convert an ARMA to MA, but i can't find one function to convert to AR. There exists one? THANKS __ R-help@r-project.org maili

Re: [R] Nestad ANOVA with random Factors

Dear Dennis, I sampled five different plots inside six predetermined areas in three diferent days. Am I doing the analysis right? Thanks, Pablo -- View this message in context: http://r.789695.n4.nabble.com/Nestad-ANOVA-with-random-Factors-tp2340725p2341484.html Sent from the R help

[R] Nestad ANOVA with random Factors

y, I am sure its not right because of these "NaN" produced and because of the area beeing assigned as Random in the output. I would appreciate any help and thanks for your attention in advance, Pablo -- View this message in context: http://r.789695.n4.nabble.com/Nestad-ANOVA-with-ran

Re: [R] Plot in cartesian plane

Thanks! That was exactly what I was looking for. best regards *pablo de camargo cerdeira* pablo.cerde...@gmail.com | pa...@fgv.br Phone: +55-(21)-3799-6065 [image: Facebook] <http://www.facebook.com/pablo.cerdeira>[image: LinkedIn]<http://br.linkedin.com/in/pablocerdeira>[image: G

[R] Plot in cartesian plane

ne help me with this? If possible, I'd like to plot this curves in a perfect cartesian plane. f = function(x) { x^2 } f2 = function(x) { -x^2 } plot(f,-10,10) abline(h=0, v=0, col = "gray60") curve(f2,col="orange", add=T) Thanks in advanced! *pablo de camargo cerdeira*

Re: [R] The condition has length > 1 and only the first element will be used

e, Aug 3, 2010 at 8:03 AM, Pablo Cerdeira wrote: > Hi All, > > I'm trying to run the following script in R, but I'm getting a warning saying: > > Warning message: > In if (z < 0) { : >  the condition has length > 1 and only the first element will be used > >

[R] The condition has length > 1 and only the first element will be used

} else { mod <- x } } f = function(x) { f <- mod(x)/x } x <- seq(-1,1,0.01) x y <- f(x) y plot(f,xlim = c(-1,1)) remove(x,y,f,mod) best regards, -- Pablo de Camargo Cerdeira pa...@fgv.br pablo.cerde...@gmail.com +55 (21) 3799-6065 __

Re: [R] Pattern recognition

This is a good example of what I'm looking for: [image: dendrogram.jpg] Best On Thu, Jul 29, 2010 at 12:01 AM, Pablo Cerdeira wrote: > > Dear all, > > I'm trying to use some technic to do a pattern recognition over a large > dataset. I really don't have any

[R] Pattern recognition

dispersion graph, identifying the possible clusters. Does anybody have any idea on how to do something like this? Many thanks, in advanced, -- *Pablo de Camargo Cerdeira* pa...@fgv.br pablo.cerde...@gmail.com +55 (21) 3799-6065 [[alternative HTML version deleted]] ___

Re: [R] Cluster analysis

dataset to you. Best regards and thank you very much, On Tue, Jul 27, 2010 at 10:16 AM, Jim Porzak wrote: > Pablo, we've had success using > http://mephisto.unige.ch/traminer/preview.shtml to look at marketing > paths. Question would be how many distinct case step discript

Re: [R] Cluster analysis

8-10-12 00:00:00","MIN. PAULO BROSSARD",6677 1480183,2,20,3,"ADIADO O JULGAMENTO","1991-07-01 00:00:00",32,"FÓRUM DA COMARCA DE RANCHARIA","1988-10-12 00:00:00","MIN. PAULO BROSSARD",6677 1480183,2,21,158,"DECISAO PUBLICADA, DJ:&quo

[R] Cluster analysis

ave no idea on how to do it using R. Can someone help me? Thanks in advanced -- *Pablo de Camargo Cerdeira* pa...@fgv.br pablo.cerde...@gmail.com +55 (21) 3799-6065 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEAS

[R] Plot map Brazil and documentation

x27;d like to read more about this amazing tool.* * * *Thanks in advanced, again,* * * *best regards* -- Pablo de Camargo Cerdeira pablo.cerde...@gmail.com +55 (21) 3799-6065 [[alternative HTML version deleted]] __ R-help@r-project.org mail

Re: [R] Multiline and grouping in R

gt; ~ Roger Brinner > > The combination of some data and an aching desire for an answer does not > ensure that a reasonable answer can be extracted from a given body of > data. > ~ John Tukey > > > > -Oorspronkelijk bericht- > > Van: r-help-boun...@r-project.org &g

Re: [R] Multiline and grouping in R

he auto.key list. > > There are **many** options in xyplot() and other Lattice graphics > functions, so you have the capability of fine tuning a graph to meet your > specifications. > > HTH, > Dennis > > > On Tue, Jun 29, 2010 at 5:42 PM, Pablo Cerdeira > wrote: &

[R] Multiline and grouping in R

ar so good... But the problem is: I'd like to create a multiline plot, one line per AREA, showing the evolution of the number of CASES per YEAR. I know how to do it in Excel, using a Pivot Table. But I'm trying hard to do the same with R but I have no idea on how to do it. Can someone

[R] combining dataframes with different row lenght

Hello list, I want to combine dataframes from tsv files which have different row lenght. I tryed cbind but doesn't work. Is there an easy way to fill with NAs till a common point of rows for all the dataframes and then merge them? Any other idea? Thanks a lot, Juan [[alternative HTML ve

Re: [R] How to report or correct translations typos?

Dear Duncan, Thank you for clarifying this issue. Cheers, Pablo - Original Message - From: "Duncan Murdoch" To: "Kenneth Roy Cabrera Torres" Cc: ; Sent: Thursday, January 07, 2010 3:12 PM Subject: Re: [R] How to report or correct translations typos? >

[R] Help with RGDAL

y much for your help, Pablo and Margarita PS. We have already installed: GDAL (for QGIS) and "sp" (for R) _ [[alternative HTML version deleted]] _

[R] how to partitioning a clustered heatmap

ng is something like this: heatmap.2(x, scale="none", cexRow=0.1, keysize=0.8, hclustfun = hclust, sepcolor="none",Colv=FALSE, dendrogram="row", trace=c("none"),density.info=c("none")) If you can help me, will be great. Thanks in advance, j

[R] matrix merge problem

Dear r-help contributors, I have two questions: first: I have a matrix A and a vector B. I want to make a new matrix C, which is made of the rows of A having a value included in B. Second: I have two matrixes A and B, of different dimensions. B has unique values in column 2 and A has not uniqu

[R] matrix filtering and reordering

to the file with same [x,1] value of blue, adding a cero in the case that there is no [x,1] corresponding value in the red matrixes (case for row 2 of blue). Thanks in advance, juan Juan Pablo Fededa [[alternative HTML version deleted]] __ R

Re: [R] Interpreting GLM coefficients

Thank you very much! It helps a lot. You are right about the NA´s in the coefficients, the model needs some simplification. Thank you again. --- Pablo Pita Orduna Grupo de

[R] Interpreting GLM coefficients

38e-02 NA NA I would like to know if it is possible to extract the intercepts and the slopes for Year and Year^2 for this Model2. Thank you very much for your help. ____ Pablo Pita Orduna Gr

[R] filtering a simple matrix

Hi all, I have the next matrix: ab c 1223 2342 3054 4523 5022 6721 7120 8919 I want to filter the rows with the values of b higher than 1 in a way tha

Re: [R] How to create a chromosome location map by locus ID

hanx Sake -- *Pablo G Goicoechea* Bioteknología Saila / /Dpto Biotecnología/ NEIKER-Tecnalia Apdo 46 01080 Vitoria-Gasteiz (SPAIN) Phone: +34 902 540 546 Fax: +34 902 540 547 pgoikoet...@neiker.net <mailto:pgoikoet...@neiker.net> __ R-help@r-pro

Re: [R] How to create a chromosome location map by locus ID

__ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] variance components

78934812 Hope this helps Pablo John Sorkin escribió: Please forgive the cross posting. I sent this to Medstats and have not received any response. I am looking for a good written explanation of components of variance (a.k.a va riance components), in particular as they are used i

Re: [R] Parallel R

Just out of curiosity, what system do you have? These are the results in my machine: > system.time(exp(m), gcFirst=TRUE) user system elapsed 0.520.040.56 > library(pnmath) > system.time(exp(m), gcFirst=TRUE) user system elapsed 0.660 0.016 0.175 J

Re: [R] Parallel R

a/R-271/lib/R/library/pnmath/libs/pnmath.so': libgomp.so.1: shared object cannot be dlopen()ed Error: package/namespace load failed for 'pnmath' I find it odd, because libgomp.so.1 is in /usr/lib, so R should find it. Juan Pablo On Sun, Jun 29, 2008 at 1:36 AM, Martin Morg

[R] Parallel R

is very simple, it doesn't require any communication between parallel tasks; only that it divides simetricaly the task between the available cores. Also, I don't want to run the code in a cluster, just my multicore machine (4 cores). What solution would you propose, given your experience

Re: [R] Equivalent of Mathematica's Apply

om: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On > Behalf Of Juan Pablo Romero Méndez > Sent: Tuesday, June 24, 2008 11:51 AM > To: [EMAIL PROTECTED] > Subject: [R] Equivalent of Mathematica's Apply > > Hello, > > Is there some equivalent function in R of Mathema

[R] Equivalent of Mathematica's Apply

5 10 and I'd like to extract sequences 1:6, 2:7, 3:8, ..., 5:10 from a vector x. What would be the most compact way of achieving this? I've tried sapply(1:5, function(i) x[ind[i,1] : ind[i,2]]) Although sapply(1:5, function(i) x[seq(ind[i,])]) would be better (if i

[R] value transformations in a vector

, Juan Pablo [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self

[R] problems installing rJava

Hello, I'm having trouble installing rJava on R-2.6.1 / Suse10.3. What could be wrong? Regards, Juan Pablo This is the output of > install.packages("rJava", dependencies=T) --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... do

Re: [R] tcltk issue remains

Hi Ingmar Have you opened X11? hope it helps Pablo Ingmar Visser escribió: > Dear R-help, > > I'm trying to load the fGarch package and keep running into problems > with tcltk: > > After succesfully instaling fGarch (and dependencies) I get: > > >library(fG

[R] Problems with quantreg and party package in 2.6.2 ver

I've problem with this pkg: - quantreg - party When i try to run this pkg under 2.6.2 (win 32). the system reports problems to reed the Rblas lib. (Not find the dynamic library links) Any suggestions to solve the problem ... tks, pablo. -- Pablo Fco. Fernández Al

[R] two histograms in the same graph

Dear Contributors: I have two vectors x and z, and I want to display the histograms of both vectors in the same graph, x in red bars, z in blue bars. If you have any clue on how to do that, I will be very glad to hear it!! Thanks in advance again, Juan Pablo Fededa [[alternative

[R] vector generation

Dear Contributors: I have the next vector: "Z" z = (526, 0.1, 110, 0.2, 34, 0.4, 614, 0.5, 1, 14, 15) I want to generate a vector containing the values higher than 1 corresponding to the z vector, wich in this case will mean to generate the next vector: (526, 110, 34, 614, 14, 15) If you ha

[R] vector generation

on, and so on. Is this doable in a simple way?? Thanks in advance again, Juan Pablo Fededa __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html a

Re: [R] Rating R Helpers

among old timers themselves. Packages reviews are another issue. But if anybody is going through all that work, why not to make the appropriate corrections to the packages? They are GPL, aren't they? Best Pablo Mark Kimpel escribió: I'll throw one more idea into the mi

[R] more outrageous plotting

Is this doable? Thanks in advance again, Juan Pablo __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] axis invertion

Hi all, I want to invert the y axis of a plot, how can I do that in an easy way? Thanks in advance, Juan Pablo __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org