[R] Problems with foreign

Good day My name is José Oscar, I'm from Mexico and I have some questions about foreign in your write.foreig( ) function. We know that this function generates the inputs to be able to run them or execute them in another program like SPSS, SAS or Stata. In these cases, when creating an ex

Re: [R] Soft Question: Where to find this reference.

See book: Statistical Models in S (ISBN-10: 041283040X; ISBN-13: 978-0412830402) Chapter 5, p. 145. • Analysis of Variance; Designed Experiments, J. M. Chambers, A. E. Freeny and R. M. Heiberger, Statistical Models in S, J. M. Chambers and T. J. Hastie (editors), Wadsworth & Brooks/Cole, P

[R] R - Understanding output of the Knoxtest (Package (surveillance))

The function knox of the "Surveillance" package performs Knox test for space-time interaction. The output is supposed to give the numbers of events that occur in a specific distance in space and time define by the function arguments. This function also perform a Monte Carlo permutation test, which

[R] Krippendorff's alpha bootstrapping implementation

Hi Polina and Jim, Have you been able to implement the correct bootstrapping of Krippendorff's alpha? Kind Regards, Oscar [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

Re: [R] anova.coxph with subsets of data

Dear Dr. Therneau, Thanks a lot for your answer. I use your survival package often and I find it very useful. Cheers, Oscar Oscar M. Rueda, PhD. Postdoctoral Research Fellow, Caldas Lab, Breast Cancer Functional Genomics. University of Cambridge. Cancer Research UK Cambridge Institute. Li

Re: [R] anova.coxph with subsets of data

ecog.ps, data =sub.ovarian) > anova(fit) > fit2 <- coxph(Surv(futime, fustat) ~ resid.ds +rx + ecog.ps, data=sub.ovarian) > anova(fit2,fit) give the same answer? Thanks for your help, Oscar Oscar M. Rueda, PhD. Postdoctoral Research Fellow, Caldas Lab, Breast Cancer Functional Genomi

[R] anova.coxph with subsets of data

~ resid.ds *rx + ecog.ps, data = ovarian, > subset=ovarian$age>50) > anova(fit) > fit2 <- coxph(Surv(futime, fustat) ~ resid.ds +rx + ecog.ps, data=ovarian, > subset=ovarian$age>50) > anova(fit2,fit) The first p-value seems to be wrong. Would anybody please explain to

[R] Positive Smoothing in fda package

3.1061 215.9816 3.0551 314.4794 2.4754 410.2167 2.2024 58.462 1.8876 67.1204 1.8919 76.7124 1.8915 86.5114 1.9193 96.3748 1.864 Error: subscript out of bounds Would someone in the l

Re: [R] 回复： Bayesian Hidden Markov Models

Dear james, There is an argument to RJaCGH() named var.equal that sets equal variances to each hidden state. The default is TRUE, so you might want to set it as FALSE. Cheers, Oscar On 3/3/12 02:20, "monkeylan" wrote: > Dear Oscar, > > I have used the the following

Re: [R] Bayesian Hidden Markov Models

Dear James, The distances are normalized between zero and 1, so in your case all of them will be zero. You can check that with > res$Dist.for.model And do > Q.NH(summary(res)[[1]]$beta, x=0) To obtain the common transition matrix. Cheers, Oscar On 29/2/12 03:59, "monkey

Re: [R] Bayesian Hidden Markov Models

es a Reversible Jump algorithm and therefore jumps between models with different hidden states. I would suggest you take a look at the vignette that comes with the package or the paper that is referenced there for specific details of the model it fits. Hope it helps, Oscar On 28/2/12 04:5

Re: [R] Bayesian Hidden Markov Models

Dear James, Although designed for the analysis of copy number CGH microarrays, RJaCGH uses a Bayesian HMM model. Cheers, Oscar On 27/2/12 08:32, "monkeylan" wrote: > Dear R buddies, > > Recently, I attempt to model the US/RMB Exchange rate log-return time series >

[R] R not 'seeing' numbers?!

both in the standard R console, and in RStudio (v. 0.94.110), and all have the same problem. Any ideas?! I'm stumped... Thanks, Oscar. Code (with comments): tmp=seq(20,90,0.005) # I want to find the position of 35.52: which(tmp==35.52) # integer(0) -- It's not there?! tmp[3100:31

[R] Ipad on R

It is possible to install R on Ipad 2? -- Oscar Ramírez A. Universidad Nacional, Escuela de Ciencias Biológicas. M.Sc. en Conservación y Manejo de Vida Silvestre osorami...@gmail.com [[alternative HTML version deleted]] __ R-help@r

Re: [R] Choropleth in R

, Oscar. --- Oscar Perpiñán Lamigueiro Dpto. Ingeniería Eléctrica EUITI-UPM http://procomun.wordpress.com - >Hi all, >I would like to use R in order to generate a Heatmap/Choropleth >overlaid on the UK map. Can someone please help wit

[R] [R-pkgs] rasterVis 0.10-3

Dear useRs, I'd like to announce the release of version 0.10-3 of rasterVis. This package provides a set of methods for enhanced visualization and interaction with the Raster* objects from the raster package. You will find some examples at: http://rastervis.r-forge.r-project.org/. Best,

[R] [R-pkgs] pxR 0.20

e(my.px.object) will create the data.frame my.px.data with the data in the corresponding PC-Axis file. Best, Oscar. - Oscar Perpiñán Lamigueiro Dpto. de Ingeniería Eléctrica EUITI-UPM http://procomun.wordpress.com ___ R-packages mail

[R] Warning in ansari.test

PS: I'm running 2.13.1 Cheers, Oscar Oscar M. Rueda, PhD. Postdoctoral Research Fellow, Breast Cancer Functional Genomics. Cancer Research UK Cambridge Research Institute. Li Ka Shing Centre, Robinson Way. Cambridge CB2 0RE England NOTICE AND DISCLAIMER This e-mail (including any attachmen

Re: [R] adding text to spplot

) l4 = list("sp.text", c(180500,329900), "0") l5 = list("sp.text", c(181000,329900), "500 m") spplot(meuse, c("ffreq"), sp.layout=list(l2,l3,l4,l5), col.regions= "black", pch=c(1,2,3), key.space=list(x=0.1,y=.

Re: [R] aggregating data

Hi, You can get it with "by": foo <- function(x)c(length(x$probe), mean(x$exp)) res <- by(df[c('exp', 'probe')], df['gene'], FUN=foo) do.call(rbind, res) Bye, Oscar. -- Oscar Perpiñán Lamigueiro Dpto. Ingeniería Eléctrica EUITI-UPM http://pro

Re: [R] RJMCMC.

additional MCMC techniques such as delayed rejection or coupled parallel chains, so probably it is a good start for learning RJMCMC. Hope it helps, Oscar On 17/4/11 01:31, "Marcus Vinicius" wrote: Dear R users, I´m studying about Bayesian Statistics. In this context, please, anyone

Re: [R] Adding margin text to lattice graphics

', 'MoreText2')) There are lots of examples in help(c.trellis). Cheers. Oscar. - Oscar Perpiñán Lamigueiro Dpto. de Ingeniería Eléctrica EUITI-UPM http://procomun.wordpress.com --- En Sat, 9 Apr 2011 18:33:42 -0700 Dennis Fisher

Re: [R] Time series of spatial data

Hi, You could try spacetime: http://cran.r-project.org/web/packages/spacetime/ Cheers. Oscar. - Oscar Perpiñán Lamigueiro Dpto. de Ingeniería Eléctrica EUITI-UPM http://procomun.wordpress.com --- En Thu, 7 Apr 2011 03:38:12 -0700 (PDT) idham

Re: [R] Venn Diagram corresponding to size in R

Hi Shira, you can also try with an R-forge package called Vennerable. It has been in alpha status for a while now, but so far it has been very useful for all our needs regarding Venn diagrams. https://r-forge.r-project.org/projects/vennerable/ https://r-forge.r-project.org/projects/vennerable/

[R] Introducing lags in my GAM

To whom it may concern, This is Oscar Saenz. I am working on my thesis and wanted to ask for your help. I am using "R" to estimate a Generalized Additive Model where the concentration of a specific pollutant is explained by a series of variables. I need to introduce a lag of the reg

[R] Estimating a differential equation parameters based on a data set

Hello R community: Here is another question, How can I estimate a differential equation parameters based on a data set? If you know where I can find some examples will be better, Thanks again, and cheers, Oscar Rodriguez Gonzalez Mobile: 519.823.3409 PhD Student Canadian

[R] Calculating weibull's distribution mean, standard deviation, and variance

Dear R community: Sorry if this question has a simple answer, but I am a new user of R. Do you know a command, or package that can estimate the weibull distribution's mean, standard deviation and variance? or can direct me to where to find it? Thanks in advance, Oscar Rodr

[R] conditional dataframe search and find

248 23 22 1.712378 24 23 1.391501 25 24 1.130752 I need to find the time when conc < 25. I can read it off the table but I am looking for a programmatic solution. Thanks. Oscar -- Oscar A. Linares, MD Clinical Assistant Professor of Medicine Department of Me

Re: [R] t-test problem

Hi Worik, You can try 2*pt(abs(t1$statistic), t1$parameter, lower.tail=FALSE) If the test is two sided. Cheers, Oscar Oscar M. Rueda, PhD Postdoc, Breast Cancer Functional Genomics Cancer Research UK Cambridge Research Institute Li Ka Shing Centre Robinson Way Cambridge CB2 0RE England

[R] Plotting the explanatory against the dependent in a GAM

  To whoever it may correspond, My name is Oscar Saenz and I am working on my thesis in Spain. I am using GAMs in "R" and, now that I have estimated my models, I need to plot the predicted smooth functions against the dependent variable (just as in Carlslaw et al. 2007, fig. 1*).

[R] Calculation error

46.5 6.0 279.0 8.0 4.6 2.4 Cmaxss (mg/L) Cminss (mg/L) Cfluct (mg/L) C(t) (mg/L) DL wVd (mg) DL wDm (mg) 18.810.8 8.014.3 658.0 658.0 But, there is an error from R? AR

[R] removing fields of the same group from a data frame

there some simple way to do this in R? Thanks in advance. /Oscar [[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/po

[R] wavethresh: length must be power of two

; wdS3<-wd(data3,type='station') Is there any other solution? By the way, does anybody knows why the wavethresh package at CRAN is not being updated? Thank you very much! Oscar Perpiñán Lamigueiro Department of Electrical Engineering Universidad Politécnica de Madrid ___

[R] [R-pkgs] New solaR package

rg/web/packages/solaR/vignettes/solaR.pdf. A handbook of photovoltaic systems (in Spanish) is available at http://procomun.wordpress.com/documentos/libroesf/. Best regards. Oscar Perpiñán Lamigueiro. ___ R-packages mailing list r-packa...@r-pr

[R] R programmer Ann Arbor Area

Dear Rs Our group is in need of an R programmer to work on Clinical Pharmacokinetics Project in Ann Arbor, Michigan Area. Interested R programmers send e-mail to: olina...@umich.edu Thanks:-) -- Oscar Oscar A. Linares, MD Translational Medicine Unit LaPlaisance Bay, Bolles Harbor Monroe

[R] boxplot help

E110.356 E16.313 E121.708 E112.651 E115.464 E237.2 E212.639 E216.791 E236.996 E222.808 CC1.975 CC3.125 CC4.433 CC6.154 CC4.175 NC9.301 NC13.531 NC12.84 NC14.336

[R] anova/factor

69666 51.920969 53.18234 3 4.32e-05 4.32e+00 0.000490530 0.0 16.030226 16.03023 Oscar -- Oscar Oscar A. Linares, MD Translational Medicine Unit LaPlaisance Bay, Bolles Harbor Monroe, Michigan 48161 Department of Medicine, University of Toledo College of Medicine Toledo, OH 43606

[R] data manipulation

Dear Wiza[R]ds, I have a data.frame header that looks like this: v2FfaPre15v2FfaPre10v2FfaPre5v2Ffa2v2Ffa3v2Ffa4 I need it to look like this, 1510523 4 i.e., with v2FfaPre and v2Ffa stripped off Any suggestions, Thanks in advance! -- Oscar Oscar A

[R] Plotting observed vs. Predicted values, change of symbols

it.wt.up, times), lty = "dashed") plot(c(datasu$concsu, dataup$concup) ~ c(datasu$time, dataup$time), xlab = "Time (minutes)", ylab = "3H-NE (dpm/ml)", log = "y", ylim=c(1,1000),cex=1.25, pch=16,col=1) lines(times, predictionssu(fit.wt.su, times)) lines

[R] Plotting observed vs. fitted values

(concsu ~ time, data=datasu) lines(times, predictionssu(fit.wt.su, times)) plot(concsu ~ time, data=datasu, log="y") lines(times, predictionssu(fit.wt.su, times)) # Need upright in same plots (e.g., dashed line) Many thanks in advance... -- Oscar Oscar A. Linares, MD Translational Me

[R] Plot fitted vs observed values

log(y) plots (i.e., linear and semi-log plots) # Help appreciated...I spent over 2 hours searching archives. # Thanks in advance. -- Oscar Oscar A. Linares, MD Translational Medicine Unit LaPlaisance Bay, Bolles Harbor Monroe, Michigan 48161 [[alternative HT

Re: [R] Switch Help

Thanks for the Rxpert advice! I am up and running again... Oscar On Wed, Nov 18, 2009 at 5:57 AM, Colin Millar wrote: > I think you just missed some commas out... > > aar <- > function(command = c("scrn", "dx", "df")) > { > command <-

[R] Switch Help

r("dx") # function does cat("dx:Diagnosis","\n") aar(c("dx","df")) # function does cat("dx:Diagnosis","\n") # function does df = cat("df:Don't Forget","\n") BUT IT IS NOT

[R] Rcmdr

Is there a UseR! manual in the works for the outstanding Rcmdr? -- Oscar Oscar A. Linares, MD Translational Medicine Unit Bolles Harbor Monroe, Michigan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

Re: [R] Two way joining vs heatmap

Hi all How can I un subscrive this forum? 2009/8/31 Schalk Heunis > Hi > > STATISTICA has a function called "Two-way joining" (see > http://www.statsoft.com/TEXTBOOK/stcluan.html#twotwo) and the > reference material states that this is based on the method as > published by Hartigan (found this p

[R] matplot in two axis

Hell, I want to plot two column vectors in a plot, line type, but they have different scale, how can I plot both series on the same graph but one on the left axis and the other scale reflecting on the right axis? Thanks in advance. [[alternative HTML version deleted]] __

[R] lm

Hi all, I want to make multiple least squares estimation on two matrix (without intercept) imagine matrix "a" and matrix "b", I want to "regress" each colum on matrix "a" (dependent variable) on each column of matrix b, I mean, regress first colum on a to first column on b. Second column on a to

[R] Hi

Hi I have a simple question. I want to run a "n times" a simple linear regession and save beta in a matrix but I´m not able. Imagine: Data.txt is a 10*5 file and want to run 4 different stimations always regressing first column on the rest. So I try this: First I run Data on memory This is my

[R] STACK

! -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[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

Re: [R] Question about datatypes/plotting issue

make them specific dates (the first of the month). Cheers, -Oscar On Mar 10, 2009, at 7:58 PM, David Winsemius wrote: You need to convert W$Date into a real date variable. At the moment it is just a character variable. > str(W) 'data.frame': 265 obs. of 23 variables: $ Date

[R] data management

| |--| Any suggestions for doing this in R? Many thanks in advance for your help. -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[alternative HTML version deleted

[R] ?grep

help! -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[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

Re: [R] read a xls file

uide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[alternative HTML version deleted]] __

[R] data management

The data is in a tab delimited text file and I want to get it into a data.frame(). Many thanks in advance. OAL p1 p1 -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[alternative HTML version deleted]] __ R-help

Re: [R] install.views()

Thanks! Your suggestion did it:-) OAL On Sat, Jan 10, 2009 at 3:00 PM, oscar linares wrote: > > Dear Rxperts, > > Using R 2.8.1 and trying > > install.views("Cluster") > > getting error > > Error: could not find function "install.views" > &

[R] install.views()

Dear Rxperts, Using R 2.8.1 and trying install.views("Cluster") getting error Error: could not find function "install.views" Please help:-( -- Oscar Oscar A. Linares Molecular Medicine Unit Bolles Harbor Monroe, Michigan [[alternative

[R] Converting english words to numeric equivalents

Hello, I am trying to convert english words to numeric equivalents, e.g., abc to 123. Is there a function or library or package for doing this in R? If not, can it be done easily in R? Thanks, Oscar [[alternative HTML version deleted

Re: [R] Manipulate Data (with regular expressions)

Check out sedit() in the Hmisc package Cheers! --- On Tue, 7/8/08, Kunzler, Andreas <[EMAIL PROTECTED]> wrote: > From: Kunzler, Andreas <[EMAIL PROTECTED]> > Subject: [R] Manipulate Data (with regular expressions) > To: r-help@r-project.org > Date: Tuesday, July 8, 2008, 7:11 AM > Dear Everyon

[R] Vector Size

Hello everybody!! I'm from Colombia (South America) and I'm new on R. I've been trying to generate all of the possible combinations for a 6 number combination with numbers that ranges from 1 to 53. I've used the following commands: datos<-c(1:53) M<-matrix(data=(combn(datos,6,FUN=NULL,simplify=

[R] install problem 2.6 mac os x 10.5

I am having the very SAME problem. I removed Version 2.5.1 before trying to install 2.6.0. Please help. TIA, Oscar __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R