[R] ask for help about lme function

58.8%; 29.72%;AIC=462.4 ###2. lme1 lme1<-lme(data, y~x1+x2+x3,random = list(site=~1,year=~1)) r2(lme1)## NA(can't get the result) r.squaredGLMM(lme1)##92.4%; 30.89%;AIC=533 I will be appreciated for your help. Yours sincerely, Juan Liu Juan Liu PhD candidate Col

Re: [R] Learning Shinny

You could use the following rmarkdown book, with last chapter debited to shiny :) https://bookdown.org/yihui/rmarkdown/ El jueves, 9 de abril de 2020, Kumar t escribió: > Hello all , > > Very sorry to ask you question that might have been answered earlier . I > could not able to find right answ

Re: [R] k-nearest neighbours from distance matrix in spdep

Dear all, I got a message from the spdep package developer and suggested me to open a GitHub issue for the problem. Please follow the discussion at: https://github.com/r-spatial/spdep/issues/38 Cheers, J De: Juan Antonio Rodriguez Perez Enviado: venres 07

[R] k-nearest neighbours from distance matrix in spdep

How should I do this? Is this even the right way to proceed? Thanks all in advance!! Best wishes, Juan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/list

Re: [R] How to use Conda with R + RStudio Server

Solution: https://github.com/grst/rstudio-server-conda It works. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] How to use Conda with R + RStudio Server

uot; within my conda environment, but this is far more complicated, as I have not been able to install it in an easy way, so that, changing my conda environment, also changes my RStudio Server Environment: https://www.rstudio.com/products/rstudio/download-server/ Anyone has any gess on how to d

[R] Calibrate survey using rake

Hello, I need to calibrate a survey using raking. I already declared the complex design here: des_survey <- svydesign(id= ~idd, strata= ~strata, data = datos, weights = ~ft, check.strata=TRUE ) options(survey.lonely.psu = "remove") options(digits = 10) However, when I run the rake process, the t

[R] Error in unique() man page

Function unique() man page states that: "The array method calculates for each element of the dimension specified by MARGIN if the remaining dimensions are identical to those for an earlier element (in row-major order). " I think that the last precission: "(in row-major order)" is meaningless here,

[R] Fwd: error in duplicated() man page

-- Forwarded message - De: Juan Gomez Date: vie., 10 may. 2019 a las 11:12 Subject: Re: [R] error in duplicated() man page To: Bert Gunter Coming back to my example, which I think illustrates my point, we have 3x2 groups (as specified in MARGIN dimensions) with one cell in each

[R] error in duplicated() man page

I think there is an error in duplicated() help page when it states that: "The array method calculates for each element of the sub-array specified by MARGIN if the remaining dimensions are identical to those for an earlier (or later, when fromLast = TRUE) element (in row-major order). " Instead of:

Re: [R] which() function help page precision

Hi again (, I am the PO from my own email account) I agree that the word "basically" puts the NA issues aside. But my point is that R subsetting behavior when there are NAs in a logical index is quite tricky to say the less, and deserves the trouble of pointing it out in every place it is appro

Re: [R] ANOVA Permutation Test

Thank you all for your **very good** answers: Using aovp(..., perm="Exact") seems to be the way to go for small datasets, and also I should definitely try ?kruskal.test. Juan [[alternative HTML version deleted]] __ R-help@r-p

[R] ANOVA Permutation Test

sense? Or maybe this issue could be due to unbalanced classes? I also tried to weight observations based on proportions, but the function failed. Any alternative solution for performing a One-Way ANOVA Test over Non-Normal Data? Thank you. Juan [[alternative HTML version d

Re: [R] initiate elements in a dataframe with lists

By the way, this also works: dfl <- data.frame(x = 1:3, y = I(list(1:2, 1:3, 1:4))) As indicated in "Advanced R" book: http://adv-r.had.co.nz/Data-structures.html#data-frames __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://

Re: [R] initiate elements in a dataframe with lists

> Just for my understanding: > Is a data.frame with list columns still a data.frame? Isn't it then a list? * A data.frame (or tibble) is a list of columns. * In which each column must be from the same data type, in this case list(). __ R-help@r-project.

Re: [R] initiate elements in a dataframe with lists

Check tidyverse's purrr package: https://github.com/rstudio/cheatsheets/raw/master/purrr.pdf In the second page of the cheatsheet there is info on how to create list columns within a data.frame :) [[alternative HTML version deleted]] __ R-help

[R] Trouble with tibbles

Factor Data Type, indeed, is of typeof() Numeric. Try converting that column to character with as.character() Class transformation then must work :) Juan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To

Re: [R] parallel computing in r....

Maybe this R-bloggers post could also help you: https://www.google.es/amp/s/www.r-bloggers.com/implementing-parallel-processing-in-r/amp/ https://www.google.es/amp/s/www.r-bloggers.com/a-guide-to-parallelism-in-r/amp/ [[alternative HTML version deleted]]

Re: [R] Portable R in zip file for Windows

ricted machine) > but it might be safer to just ask the powers that be to download it > for you. You probably don't need a new version of R more than > once a year. > > > On Thu, Jan 25, 2018 at 3:04 PM, Juan Manuel Truppia > wrote: > > What is wrong with you gu

Re: [R] Portable R in zip file for Windows

t;https://maps.google.com/?q=3+Aranleigh+Park+Rathfarnham+Dublin&entry=gmail&source=g> > 14 > Ireland > www.tcd.ie/Economics/staff/frainj/home.html > mailto:fra...@tcd.ie > mailto:fra...@gmail.com > > On 25 January 2018 at 20:04, Juan Manuel Truppia > wrote: > &

Re: [R] Portable R in zip file for Windows

r Be sure to enter the desired > path in the Select Destination Location screen, choose Yes on the > Startup options screen and ensure that all boxes are unchecked on the > Select additional tasks screen. > > On Wed, Jan 24, 2018 at 10:11 PM, Juan Manuel Truppia > wrote: > >

Re: [R] Portable R in zip file for Windows

; have changed completely since then. > > > On Wednesday, January 24, 2018, 10:12:09 PM EST, Juan Manuel Truppia < > jmtrup...@gmail.com> wrote: > > > I read a message from 2009 or 2010 where it mentioned the availability of R > for Windows in a zip file, no installatio

Re: [R] Portable R in zip file for Windows

net/projects/rportable/ > -- > Sent from my phone. Please excuse my brevity. > > On January 24, 2018 7:11:49 PM PST, Juan Manuel Truppia < > jmtrup...@gmail.com> wrote: > >I read a message from 2009 or 2010 where it mentioned the availability > >of R > >for Window

[R] Portable R in zip file for Windows

I read a message from 2009 or 2010 where it mentioned the availability of R for Windows in a zip file, no installation required. It would be very useful for me. Is this still available somewhere? Thanks [[alternative HTML version deleted]] __ R

Re: [R] Best R GUIs

Thank you all, some of the best free IDEs are specified in the following article: http://www.linuxlinks.com/article/20110306113701179/GUIsforR.html Plus pluggable IDEs: A) ESS for Emacs B) IRkernel for Jupyter notebooks C) StatET for Eclipse D) Vim-R for vim Kind regards, Juan

[R] Best R GUIs

regards, Juan Telleria [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] Best way to study internals of R ( mix of C, C++, Fortran, and R itself)?

The R Community made a call for one person to be in charge of R Contributed Documentation, and I have done a request for being in charge of such duty. If assigned, my plan is to implement Atlassian's Confluence along the R Community (Accessed though R Project.org), in order to generate a Wiki and

[R] Print All Warnings that Occurr in All Parallel Nodes

ion)) Would this work? * And also, How could I check that a CSV can be opened before applying the function, and create an empty data.frame for those CSV. Thank you, Juan CODE ## DISPOIN Data Import Into Mari

[R] GLMM for Combined experiments and overdispersed data

levels... m3 = glmer(resp ~ trt * farm + (1|farm:bk) + (1|tree_id), family = binomial, data=df) I really appreciate your suggestions about my model specifications... *Juan​ Edwards- - - - - - - - - - - - - - - - - - - - - - - -# PhD student - ESALQ-USP/Brazil​* [[alternative

Re: [R] Display data by condition

are trying to solve? > Tell me what you want to do, not how you want to do it. > > On Thu, Mar 16, 2017 at 11:16 AM, Juan Ceccarelli Arias > wrote: > >> Hello, >> I need to show the observations of a data set only if the earn more than >> $5000 (fact is its na

Re: [R] Display data by condition

data frame. Go re-read the > Introduction to R document section on indexing to figure out where the > comma goes. > -- > Sent from my phone. Please excuse my brevity. > > On March 16, 2017 8:16:29 AM PDT, Juan Ceccarelli Arias > wrote: > >Hello, > >I need to show the ob

[R] Display data by condition

Hello, I need to show the observations of a data set only if the earn more than $5000 (fact is its name in the date set). I use this: View(data[data$fact>5000]) The code above shows nothing. No error or message at all. What am i doing wrong? Thanks for your help and time. [[alternative H

[R] Seasonal package: "Error in x[[2]] : subscript out of bounds"

Hi Im running X13 with the seasonal package. However, in some series i get this error after specifying the x11 options: o_1=ts(o1,frequency=12,start=c(2010,2)) so_1=seas(o_1, x11="",transform.function ="log",arima.model="([3] 1 1)(1 1 0)", x11.seasonalma="S3x5",x11.trendma=13) Error in x[[2]] : su

Re: [R] svytable: How do i create a table informing a third variable?

Your help was everything i needed it. Please, declare this topic as solved. And thanks again. On Sat, Sep 3, 2016 at 10:06 PM, David Winsemius wrote: > > > On Sep 2, 2016, at 6:08 PM, Juan Ceccarelli Arias > wrote: > > > > Thanks a lot. Your code does the trick. > &

Re: [R] svytable: How do i create a table informing a third variable?

, svymean ) > > # median > svyquantile( ~ income_variable , NN ) > svyby( ~ income_variable , ~ age + sex , NN , svyquantile , 0.5 ) > > > > > On Fri, Sep 2, 2016 at 3:04 PM, Juan Ceccarelli Arias > wrote: > >> Hello >> Im analyzing a survey and i need to obt

[R] svytable: How do i create a table informing a third variable?

Hello Im analyzing a survey and i need to obtain some statistics per groups. Im able to create a table with sex and age. However, if i want to know how much income earns the population by sex and age, i can't. Im loading the dataset as describe the line below NN <- svydesign(ids = ~1, data = encues

Re: [R] Importint stata file and using value labels

I solved this problem using the nonint.factors=F,generate.factors=F) when i imported the dta files. Thanks to all. My issue can be declared as solved. On Sat, Aug 27, 2016 at 12:55 PM, Michael Friendly wrote: > On 8/26/2016 11:05 AM, Juan Ceccarelli Arias wrote: > >> Yep. Im a

Re: [R] Loop over rda list files and using the attach function

Hi I want to comment something. When i added the detach(get(yyz)) the RAM consumption was considerable reduced. So, i want to declare this issue as solved and thank you all for your assistance. Good luck to all. On Tue, Aug 30, 2016 at 6:24 PM, Juan Ceccarelli Arias wrote: > The attach(get(

Re: [R] Loop over rda list files and using the attach function

g 30, 2016 at 3:42 PM, wrote: > Hello, > > Try > > attach(get(yyz)) > > Hope this helps, > > Rui Barradas > > > > > > > > Citando Juan Ceccarelli Arias : > > Hi. >> I need to loop over rda files. >> I generated the list of t

[R] Loop over rda list files and using the attach function

Hi. I need to loop over rda files. I generated the list of them. That's ok. The problem is that the name of the files are as _mm (eg 2010_01 is january or 2010, 2016_03 is march of 2016). So, when i try to use the attach function to create a simple table(age, sex) it fails. The only way to atta

Re: [R] Importint stata file and using value labels

t; > > -Original Message----- > > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Juan > > Ceccarelli Arias > > Sent: Thursday, August 25, 2016 9:15 PM > > To: ruipbarra...@sapo.pt > > Cc: r-help@r-project.org > > Subject: Re: [R] Impor

Re: [R] Importint stata file and using value labels

"Hombre"],type[sex=="Hombre"] > > Rui Barradas > > > Citando Juan Ceccarelli Arias : > > str(sex) >>> >> Factor w/ 2 levels "Hombre","Mujer": 1 2 2 1 2 2 1 1 1 1 ... >> On Thu, Aug 25, 2016 at 1:42 PM, wrote: >&

Re: [R] Importint stata file and using value labels

> str(sex) Factor w/ 2 levels "Hombre","Mujer": 1 2 2 1 2 2 1 1 1 1 ... On Thu, Aug 25, 2016 at 1:42 PM, wrote: > Maybe sex is a factor and Man its label. Factors are coded internally as > integers, to see it use > > str(sex) > > Rui Barradas > &g

Re: [R] Importint stata file and using value labels

instead > > table(region[sex==1],type[sex==1]) > > To test for equality use == not =. > > Hope this helps, > > Rui Barradas > > > Citando Juan Ceccarelli Arias : > > Hi > Im a bit lost. > Ive imported a stata13 file. When i try to make a simple >

[R] Importint stata file and using value labels

Hi Im a bit lost. Ive imported a stata13 file. When i try to make a simple table(region[sex=1],type[sex=1]) i get something as type1 type2 region1 0 0 region2 0 0 region3 0 0 I don't know how to operate with the value labels. Sex is defined as 1 as "Man" and 2 as "Female". And region1 i

Re: [R] Loop over folder files

Ok. Please, declare this issue as solved. And thanks again for your help. On Wed, Aug 24, 2016 at 2:18 PM, wrote: > Maybe it's better to open a new thread. > > Rui Barradas > > > Citando Juan Ceccarelli Arias : > > The error wasn't in the loop. It was in

Re: [R] Loop over folder files

setwd("path/to/files") to set the right place where the files can be found. > > Rui Barradas > > > Citando Juan Ceccarelli Arias : > > I just doesn't work... > Im loading the read,dta13 package already. > When i try to perform a simple table(sex), i recei

Re: [R] Loop over folder files

i doing wrong... On Tue, Aug 23, 2016 at 5:08 PM, wrote: > Or maybe a print() statement on the table() in the loop. > > print(table(...)) > > Rui Barradas > > > Citando David Winsemius : > > On Aug 23, 2016, at 10:01 AM, Juan Ceccarelli Arias > wrote: > >

Re: [R] Loop over folder files

> > > The first package is readstata13 but there are others. > > Hope this helps, > > Rui Barradas > > > > Citando Michael Dewey : > > Dear Juan > > If this is a Stata 13 file the package readstata13 available from CRAN may > be of assistance. > >

Re: [R] Loop over folder files

o a data.frame. > > Hope this helps, > > Rui Barradas > > > Citando Juan Ceccarelli Arias : > > Hi > I need to apply some code over some stata files that are in folder. > I've wrote this > > library(foreign) > > fuente=list.files("C:/Users/J

[R] Loop over folder files

Hi I need to apply some code over some stata files that are in folder. I've wrote this library(foreign) fuente=list.files("C:/Users/Jceccarelli/Bases/Stata", pattern="dta$", full.names=FALSE) for (i in 1:length(fuente)){ xxx=read.dta(fuente[i], to.data.frame=TRUE) } But i get this error Err

Re: [R] Adding loess lines subsetting to each panel in lattice plot

Breathed in his "Bloom County" comic strip ) On Fri, Aug 12, 2016 at 9:10 AM, Juan Perez via R-help wrote: > Hello, I've created an xyplot and I want to add a loess line for x (Age)  > <=40 and another for values >40. In a way it is similar to this > https://stat

[R] Adding loess lines subsetting to each panel in lattice plot

Hello, I've created an xyplot and I want to add a loess line for x (Age)  <=40 and another for values >40. In a way it is similar to this https://stat.ethz.ch/pipermail/r-help/2009-May/390502.html but still not succcessful. This is my try: xyplot(MOE~Age|Species, groups=Site,    panel = fu

[R] Something weird

Can anybody explain me this weird result? a=3 as.integer(a) 1] 3 a=(3/10)/0.1 a [1] 3 as.integer(a) [1] 2 Thank's in advance Juan A Hernández [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIB

[R] non-linear fit and statistical significance

ck. Regards. Juan C. Munoz-García [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] reikon: A package to interact with Thomson Reuters Eikon from R

mine. It has documentation and a vignette. Please report any bugs directly in Bitbucket, and I would gladly accept help with the development. Regards Juan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBS

[R] two problems

Hi, list experts I am totally a new user of the R package. I just tried to use qcomhd, but it gives me an error alarm, as below. qcomhd(RT1, RT2, q = seq(.1, .9, by=.1)) Error: could not find function "qcomhd" Then I try to install it, however, it seems that I am not able to install it. Also, in

[R] rdde : A package for DDE connections

Hi all, I'm building a new package for DDE connections on R. It's called rdde and lives in https :// bitbucket.org / juancentro

[R] rdde : DDE connections on R

Hi all, I'm building a new package for DDE connections on R. It's called rdde and lives in https://bitbucket.org/juancentro/rdde. It's in alpha stage, but operational. It has a very simple vignette which explains the main reason you should try rdde vs tcltk2 (the only other option available) : perf

[R] error using the huge R package

p would be very appreciated Juan [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting

[R] How to F noncentrality parameter

Does anybody know how to calculate noncentrality parameter for a three way interaction? It is easy for a two way interaction, but I can't find the way to do it with a three way. I have been looking for a packages about this, but it seems that none manages to do this. Thank's in ad

[R] DOT PLOT help!!

splay this as categories on the x axis and all values for each directly above the applicable category, thus not a scatter graph but in a line above for each. I have attached a quick picture I found on the internet to demonstrate roughly what I would like. Thank you all for your time Juan Scheun

[R] How to extract table results from survival summary object

e or assign this table to an object. Anyone knows how to solve it. Thank's in advance Juan A. Hernández [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

[R] sqldf problems

ndows 8.1 Do I need any other installations o packages or something else? The code was working in older Windows XP and Linux Thanks -- ****** Juan Daniel García Villabrille, Ingeniero de Montes, colegiado nº

[R] How to extract convergence code from lmer object?

I get this important information in an automatic way? Thank's in advance Juan A. Hernández Spain [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gui

[R] How to delete designated sequences in alignment (DNAbin object)

Thanks in advance for your help. -- Dr. Juan A. Balbuena Cavanilles Institute of Biodiversity and Evolutionary Biology University of Valencia [1]http://www.uv.es/~balbuena P.O. Box 22085 [2]http://www.uv.es/cophylpaco 46071 Valencia, Spain e-mail: [3]j.a.balbu...@uv.este

[R] How to print something in the same location in console?

Hi I need to print the iteration number of a procedure but in the same location in console. Using cat with or without fill argument does not produce the desired outcome. Does anybody know how to get it?. for(i in 1:10) cat('Iteration:',i,fill=T) Thank's in advance Ju

Re: [R] R License

Thanks Mensaje original De: Barry Rowlingson Fecha:26/05/2014 03:11 AM (GMT-05:00) A: Juan Ulises Bohorquez Carvajal CC: John Kane ,"r-help@R-project.org" Asunto: Re: [R] R License On Sat, May 24, 2014 at 4:11 PM, Juan Ulises Bohorquez Carvajal wrote: > Thanks, c

Re: [R] R License

Thanks for your information Jb Mensaje original De: John Kane Fecha:24/05/2014 10:33 AM (GMT-05:00) A: "Timothy W. Cook" ,Juan Ulises Bohorquez Carvajal CC: "r-help@R-project.org" Asunto: Re: [R] R License Aha, someone knowledeable has shown up. T

Re: [R] R License

Thanks, could you send to me an official certificate in pdf with this information? Or a certified email to support the information license? Best regards, Jb Mensaje original De: John Kane Fecha:24/05/2014 10:04 AM (GMT-05:00) A: Juan Ulises Bohorquez Carvajal ,"r-h

Re: [R] R License

Thanks Mensaje original De: Jeff Newmiller Fecha:24/05/2014 10:37 AM (GMT-05:00) A: Juan Ulises Bohorquez Carvajal ,"r-help@R-project.org" Asunto: Re: [R] R License This is not an appropriate forum to be obtaining legal advice in. That said, the licenses [1] fo

[R] R License

Hello, I work in the area of ​​Geomatics in Ecopetrol, Ecopetrol is an Oil & GAS company with 10,000 employees. The licensing of R allows me to use it in production at my company or academic use only? Sincerely, JB. Este mensaje y sus anexos está dirigid

Re: [R] multicore and mclapply problem in calculation server

e, the systems gets quickly overflowed. So the solution is to split length(linH) into blocks equal to the number of cpus and the message to end-users is "be careful with parallel processing functions in R". Cheers Juan Antonio Balbuena El 19/12/2013 22:03, p_connolly escribió:

[R] multicore and mclapply problem in calculation server

na/1MH_30.tre", open="rt") # READS TEXT DATA FROM EXTERNAL FILE conP = file("/scratch/ba/balbuena/1MP_30.tre", open="rt") #"" "" " " replicate (N.runs, body(N.perm)) # LOOPING body NUMBER OF RUNS

[R] Calculate external validation

e it allows to use offset. What do you recommend me? Thank you member of R-project Juan Manuel Reyes [[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] Question about compatibility

connection Además: Mensajes de aviso perdidos In file(file, "rt") : It is not posible to open the file 'example.text': No such file or directory I am using office 365 and my computer has Windows 8. What can I do? Thank you members of R- project Juan

[R] Same code - error in one PC but not in other

uch appreciated. Thank you very much for your attention Juan A. Balbuena -- Dr. Juan A. Balbuena Marine Zoology Unit Cavanilles Institute of Biodiversity and Evolutionary Biology University of Valencia [1]http://www.uv.es/~balbuena P.O. Box 22085 [2]http://www.uv

[R] barplot to the right side

tarting in the right side of x axes and the horizontal bars starting in the right side too. Thank's in advance Juan A. Hernandez-Cabrera [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

Re: [R] transform dataframe with look-up table

t (pun intentional): End <- merge(merge(Start, transformer, by.x="Right", by.y="input", all.x=TRUE), transformer, by.x="Left", by.y="input", all.x=TRUE) Best wishes Juan Antonio El 25/07/2013 18:02, Brian Diggs escribió: On 7/25/20

[R] transform dataframe with look-up table

output for [5, ] should be NA NA. Thank you very much for your time. Juan Antonio Balbuena -- Dr. Juan A. Balbuena Marine Zoology Unit Cavanilles Institute of Biodiversity and Evolutionary Biology University of Valencia [1]http://www.uv.es/~balbuena P.O. Box 22085 [2

[R] question about read.ssd function

sas.exe") : SAS return code was 2 Why does R act differently for the different sas datasets? They are created at the same platform. Can anybody share some insights on this? Thanks Juan Wu Biostatistics Roche Diagnostics Professional Diagnostics Clinical Operations Phone: +1 317-521-239

[R] Hyperparameters in ARIMA models with dlm package

Hi, i'm beginner in Bayesian methods, I'm reading the documentation about dlm package and kalman filters, I'm looking for a example of transformation of ARIMA in a state space equivalent to use the dlm package and calcualte the hyperparameters. Someone can help me about it?. If it's possible with a

[R] hyper-parameters

I'm searching a method to estimate the hyper-parameters in arima models. I'm reading about r-inla package, but in the examples section only talk about the AR part of the arima, but i need help about the MA part too. I'm beginner in Bayesian methods, I'm reading the documentation about dlm package

[R] binary data

0 1 ... Pepz 1 1 1 1 1 0 0 1 I would like to know some R package to calculte similarity matrix or to make some statistical diferences between cases and controls Any help it would be very appreciated BW Juan --- Juan Fernandez

[R] use row and col names of a matrix to create colnames of a new one.

Hello I would appreciate any help with the following: Given a matrix with col and row names, such as > matrix1 <- matrix(rnorm(16),4) > rownames(matrix1) <- LETTERS(1:4) > colnames(matrix1) <- letters[1:4] > matrix1 a bc d A 1.6845882 -0.27809792 -0.87479

[R] MALDIquant

Dear list, My name is Juan Fernandez-Tajes and I´m trying to use MALDIquant for analysing MS data from ABI 4700 series, I´ve exported a spot set to T2Dextractor in order to convert in mzXML format. I´ve imported this into R with mqReadMzXml function and have followed the suggested

[R] procrustes (vegan) plot of residual differences

+/- SD. Thank you very much for your attention. Juan Antonio Balbuena -- Dr. Juan A. Balbuena Marine Zoology Unit Cavanilles Institute of Biodiversity and Evolutionary Biology University of Valencia [1]http://www.uv.es/~balbuena P.O. Box 22085 [2]http://www.uv.es/ca

[R] How to stop a loop for?

ata[,i]), main=paste('histogram',i)) #here something like waituntil command } Thank's in advance Juan A. Hernandez Spain. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

[R] Error in ifelse(append, "a", "w") : , (list) object cannot be coerced to type 'logical'

imul3_2012.txt", sep ="\t ", append =T). I have read that including "\" in such statements can be problematic, because it is a scape character in R, but trying sep=" " instead of "\t" did not solve the problem. Any help will be much appreciate

[R] any DCCA function in R?

Dear members, I am performing multivariate analysis on marine benthic populations using R. At first glance I found ca and VEGANO packages to be the suitable for the task, but neither has incorporated Detrended Canonical Correspondence Analysis (DCCA), which is just the method I want to apply

[R] Discarding a matrix based on the rowSums value

Hello I would appreciate your help on the followig. I want to generate random binary matrices but I need to discard those with all-1 rows. That is, for a 10x10 matrix with five 1's [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,]000000000 1 [2

[R] permutation within rows of a matrix

Hello This is probably a basic question but I am quite new to R. I need to permute elements within rows of a binary matrix, such as [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,]000010000 0 [2,]001100011

[R] hazard rates in conditional frailty models

. Plot #1 would provide the rates when “ecogr” is below its mean and “elderly” is, respectively, one standard deviation above and one below its mean.Plot #2 would provide the rates when “ecogr” is above its mean and “elderly” is, respectively, one standard deviation above and b

[R] ksvm and predict

g columns (I did it because I get better results that if I use just WA) but I whant to make a prediction of KL using just the first of them (WA). Any suggestions? Thanks a lot for your help. Regards, Juan -- View this message in context: http://r.789695.n4.nabble.com/ksvm-and-predict-tp3827634p3

Re: [R] dead code removal

Ideally you'd have the next two items available: - tests that ensure that your code carries out what it should and as it should. - a coverage analysis tool that reports what parts of your code have been and have not been executed by your tests above. Neither of those are mandatory though, but they

Re: [R] tm: Read a single text file into a corpus as single document?

Some hints: list.files() will return the list of files in a directory readLines() will allow you to load text files as vectors of lines strsplit() will allow you to break lines into words c(x,y) concatenates vectors x and y ; x <- c(x,y) appends vector y to x unique() will allow you to get rid of r

Re: [R] Gaussian low-pass filter

gfcoeffs <- function(s, n) { t <- seq(-n,n,1) ## assuming 2*n+1 taps return ( exp(-(t^2/(2*s^2)))/sqrt(2*pi*s^2) ) } 2011/6/29 Martin Wilkes : > I want to filter my time series with a low-pass filter using a Gaussian > smoothing function defined as: > > w(t) = (2πσ^2)^0.5  exp(-t^2/2σ^2) > > I

Re: [R] A masked function is a masked function by any other name

Ops! Thank-you Duncan for clarifying the 2 vs. 3 colon difference and a couple of other things. Working like a charm now. Cheers, jcb! > If you are using ::: (three colons), then you may be looking into the > unexported functions in log4r.  The only normal way to see unexported > functions is to u

[R] A masked function is a masked function by any other name

Dear all, It looks like I do not grasp the concept of masked functions enough as to solve this trivial problem. The code that replicates the problem (a source code tree that realizes a R package actually) is under github so one can call it clone it easily from the command line (though more experien

[R] How to export to pdf in landscape orientation?

Does anybody know how to get a pdf file with landscape orientation?. pdf(file= 'my_file.pdf' ,onefile=T,paper='A4') plot(sin, -pi, 2*pi) dev.off() Thank's in advance Juan A. Hernandez Spain [[alternative HTML version deleted]] ___

  1   2   >