[R] Code problem with R/Microeconomics

Dear Sir/Madam, I am working on a statistical evaluation with R at the Erasmus University and have a problem with one coding sentense. Who can help me with that? Thank you. Kind regards, Olga Pervushina [[alternative HTML version deleted

[R] ZINB multi-level model using MCMCglmm

Hi, I am Olga Viedma. I am running a Zero-inflated negative binomial (ZINB) multi-level model using MCMCglmm package. I have a doubt. Can I use the "Liab" outputs as fitted data, instead of the predicted values from "predict"? The liab outputs fit very well with the ob

[R] R in Excel

Hello! I need some help. I'd like to make a button in Excel which performes functions wrote in R. Is there is any macros to do this? I don't want to use RExcel. Thank you very much! Sincerely yours, Olga [[alternative HTML version deleted

[R] Making visible Rcommander

(see the pic.2) But in my work book there are no such an item (see the pic.3), but I've made a connection between R and Excel Please, help me to make visible Rcommander Thank you! Sincerely yours, Olga __ R-help@r-project.org mailing list https

[R] Complex Association Plots

I'm new to R plotting, so please be gentle. I want to make a large association plot, showing whether each combination of variables averages more or less than the whole for proportion in a certain category. So, for example: Dependent Variable: eye color (30% green, 45% brown, 25% blue on

Re: [R] Scaling Statistical

escreveu: Hello, See if the following does what you want. lapply(seq_len(obsv), function(i) adf.test(df[df$ID == i, 3])) Hope this helps, Rui Barradas Em 23-06-2013 19:12, Olga Musayev escreveu: Short question: Is it possible to use statistical tests, like the Augmented Dickey-Fuller

Re: [R] Scaling Statistical

. With a standard error like NA, maybe you don't have enough data points to run the tests. (Only 3 residuals, all zero, like seen above.) Rui Barradas Em 24-06-2013 20:12, Olga Musayev escreveu: Rui-- thanks so much for the help! I'm getting this error though, which is leaving me stumped: test

[R] Scaling Statistical

Short question: Is it possible to use statistical tests, like the Augmented Dickey-Fuller test, in functions with for-loops? If not, are there any alternative ways to scale measures? Detailed explanation: I am working with time-series, and I want to flag curves that are not stationary and which

Re: [R] permutation of vectors (1 or 0)

On 12/20/2012 02:13 PM, arun wrote: library(gtools) permutations(2,10,0:1,repeats.allowed=TRUE) Thanks, this is exactly what I need. Regards, Olga __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] permutation of vectors (1 or 0)

the maximum number of permutations that R can handle? Cheers, Olga __ 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] RV: problems for making grids from lmer models

, and my best regards, Olga Olga Viedma Sillero Profesora Ordenacion del Territorio Departamento Ciencias Ambientales Facultad Ciencias del Medio Ambiente Avd/ Carlos III, s/n 45071 Toledo (Spain) email: olga.vie...@uclm.es Tel: 925-268800 (ext. 5469) Tel movil: 675 880 297 Fax: 925 268844

Re: [R] rShowMessage Fatal error: unable to open the base package

*R_OpenLibraryFile(const char *file) { char buf[259]; FILE *fp; snprintf(buf, 259, %s/library/base/R/%s, C:/Rcpp_and_RInside/R/R-2.12.2, file); fp = R_fopen(buf, r); return fp; } Obviously, it returns NULL). I don't know, how it can be corrected. May be, you know? Olga

[R] vegan and sweave using xtable

100.52 -2.31224.084 4.2271 286.1372 293.7350 Many thanks, Olga __ 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

Re: [R] vegan and sweave using xtable

Thanks Ista, Very instructive and works like a charm! Cheers! Olga __ 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

[R] Linear mixed model: question about t-values

Dear all, I have a question about the output of linear mixed model fitted in R using nlme package. In particular, what are the t-values that are given in an output, how are they calculated and based on what test? I guess it cannot be a simple Student t-test, otherwise how can the simple

[R] Linear mixed model: question about t-values

Dear all, I have a question about the output of linear mixed model fitted in R using nlme package. In particular, what are the t-values that are given in an output, how are they calculated and based on what test? I guess it cannot be a simple Student t-test, otherwise how can the simple Student

[R] lower triangle of the correlation matrix with xtable

matrix? Using 'dist' doesn't seem to be useful as it doesnt allow to convert this table with xtable. print(xtable(cor(mydata),digits=3)) Any ideas? Cheers, Olga __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] lower triangle of the correlation matrix with xtable

)) Error in UseMethod(xtable) : no applicable method for 'xtable' applied to an object of class dist Jorge's solution with upper.tri works fine! Thanks everyone for your prompt response. Cheers Olga __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] standardize columns selectively within a dataframe

trying to use subscripts but did not succeed so far. Any tips? Many thanks, Olga __ 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

Re: [R] standardize columns selectively within a dataframe

Thanks! It is exactly what I was looking for! Cheers Olga __ 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

Re: [R] standardize columns selectively within a dataframe

for this problem but it didn't work. Thanks for the references! Your solution with scale() is nice and neat, but for the sake of learning it is useful to persist. Cheers, Olga __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] frequency table for a list matching some conditions

) newlist [1] A B C D E F I want to match mylist[1], mylist[2], etc with newlist and I am looking for the following output: A B C D E F x1 3 1 1 NA NA NA x2 NA NA NA 2 2 1 x3 3 1 NA NA NA NA Any thoughts? Many thanks Olga

Re: [R] frequency table for a list matching some conditions

Thanks Dimitris, It works nicely! Regards, Olga On Thu, 2010-08-26 at 11:55 +0200, Dimitris Rizopoulos wrote: one way is the following: mylist - list(x1 = c(A,A,A,B,C,Z,Y), x2 = c(D,D,E,E,F,Z,X), x3 = c(A,A,A,B,Y,Z)) newlist - c(A,B,C,D,E,F) tab - t(sapply(mylist

[R] decision tree display

info on how to adjust the tree display and text size display? Please help. Thank you! Olga [[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

[R] decision tree finetune

=0 instead of 2) Info_G0.5. Thank you in advance! And thanks for Eric who helped with my previous question about starting rpart. Olga fit - rpart(Retention ~ Info_G+AOPD+Mail+Xref_Umbr+Ins_Age+Discount+Xref_A + Con6 + + Con5 +Con4 + + Con3 + Con2 + + Con1 , data=Home

Re: [R] decision tree finetune

I figured it out myself, here it is: control=rpart.control(cp=.001)) Thank you! On Fri, Aug 13, 2010 at 12:58 PM, Olga Shaganova olga.sha...@gmail.comwrote: My decision tree grows only with one split and based on what I see in E-Miner it should split on more variables. How can I adjust

[R] rpart package

Hi, I am a brand new user and may be my question is too simple. I have R on our (not Unix) server. I am trying to build a decision tree and the error message says couldn't find function rpart. Does it mean I have to ask our server guy to install an additional package? Thank you, Olga

Re: [R] try-error within for loop

Wu Gong wrote: So, try() inside loop should work. Like for () { if (class(try(...,silent=T))==try-error) result[[i]] - NA ... } Thanks a lot, it worked for me! Cheers Olga __ R-help@r-project.org mailing list https://stat.ethz.ch

[R] PageUp/Down in gnome-terminal

Olga -- R version 2.11.1 (2010-05-31) 2010 x86_64 GNU/Linux Ubuntu 10.04.1 LTS __ 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

Re: [R] PageUp/Down in gnome-terminal

Thanks Liviu, On Thu, 2010-08-05 at 09:51 +0100, Liviu Andronic wrote: On Thu, 05 Aug 2010 08:56:30 +0200 Olga Lyashevska o...@herenstraat.nl wrote: I am using gnome-terminal to run R, and I noticed that standard PageUp/Down do not work but they do work for other programs within the same

[R] try-error within for loop

missing 5 rows were omitted. Any suggestions? Thanks, Olga result-list() for(i in 1:1000)try({ hclass-mydata2$tclass[mydata2$tspecies%in%com[[i]][[5]][[1]]] horder-mydata2$torder[mydata2$tspecies%in%com[[i]][[5]][[1]]] hfamily-mydata2$tfamily[mydata2$tspecies%in%com[[i]][[5]][[1

Re: [R] PageUp/Down in gnome-terminal

On Thu, 2010-08-05 at 10:59 +0100, Liviu Andronic wrote: On Thu, 05 Aug 2010 11:13:00 +0200 Olga Lyashevska o...@herenstraat.nl wrote: I do not have any desktop environment (running R on a server), can I still try xfce? I am not very sure. If you are running gnome-terminal, I assume

Re: [R] multiple R sessions from one working directory using GNU screen

is reaching its limits quite often resulting in shutdown, so I wanted to split one task into pieces which indeed similar to cluster computing. I will check documentation on this topic. Kind regards, Olga On Tue, 2010-08-03 at 13:04 -0700, Steven McKinney wrote: Dear Olga An R session is conducted

Re: [R] multiple R sessions from one working directory using GNU screen

to manage massive matrices with shared memory and memory-mapped files. Now I will try to see how it can be implemented in my code. Many thanks Olga On Wed, 2010-08-04 at 13:51 +0100, david.jes...@ubs.com wrote: Olga There was a presentation at the London R user group about doing what you

[R] multiple R sessions from one working directory using GNU screen

to save.image(), q() and then restart both sessions, but it does not help. Any suggestions? Many thanks Olga R version 2.11.1 (2010-05-31) 2010 x86_64 GNU/Linux Ubuntu 10.04.1 LTS Intel(R) Xeon(R) CPU X3220 @ 2.40GHz MemTotal: 4050180 kB __ R-help@r

[R] jtest to test significance of R2 change

to generalize it? Would jtest be appropriate technique if I want to find out whether change in R2 across models is significant? Many thanks Olga __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] U de Theil

Hi I need to calculate the indicator U de Theil from my forecasts, so I want to know if there is a specific command or function that provides me the indicator U de Theil Thanks.. Olga Maria Formigoni C. Walter E-mail: mailto:olgaformig...@gmail.com olgaformig...@gmail.com MSN

[R] how to plot forecast together with historical series in OLS or special ARIMA model

= TRUE) : non-time series not of the correct length Does it exists the code how to plot it together with historical data or how to include the additional variable in the first ARIMA model? Thank you in advance, Olga Series: Year HCPI 1993 4.8 1994 2.91 1995 2.75 1996 1.03 1997 1.81 1998 1.03

Re: [R] barplot with a probability density curve

be a suitable tool to use. However, it also means that I cannot plot probability density curve using this type of data because I would need to have a histogram instead. Did I get it right? Many thanks, Olga On 03/15/2010 08:32 AM, Olga Lyashevska wrote: Dear all, I am making a barplot

[R] barplot with a probability density curve

. Any tips? Thanks Olga __ 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] looking into function body without package installation

to explore function syntax I end up installing loads of packages I am never going to use. I have also tried to use package 'sos' which does provide a good facility to find relevant information, but it seems it does not give any info on function body. Thanks, Olga

Re: [R] Conditional sampling?

a need to be dropped. col1col2 3 b f 4 b g 5 c h 6 c i s - sample(col1,2,replace=T) col2[col1 %in% s] Cheers, Olga __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

[R] Conditional sampling?

Hello, I have a dataframe consisting of two columns. col1-factor(c(a,a,b,b,c,c)) col2-factor(c(a,b,c,d,e,f)) somedf-data.frame(col1,col2) somedf col1 col2 1 a d 2 a e 3 b f 4 b g 5 c h 6 c i sample(col1,2,replace=T) [1] b c

Re: [R] Conditional sampling?

Thanks Thierry! On 04.02.2010, at 13:22, ONKELINX, Thierry wrote: selection - subset(somedf, col1 != a) sample(selection$col2, 2, replace = TRUE) I simplified it, but in fact it is not only 'a' to be eliminated, but thousands of rows. Can I use: col1sample-sample(col1,2,replace=T)

[R] working with taxonomic trees: sampling

time obtaining different number of elements at each taxonomic level. Sorry for this long-winded post, I hope it is clear what I am trying to do. I would appreciate any tips! Thanks, Olga __ R-help@r-project.org mailing list https://stat.ethz.ch

Re: [R] RMySQL - Bulk loading data and creating FK links

like? Cheers, Olga __ 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] Microsoft SQL Server and R

Hi David, I have a client running Microsoft SQL Server. I am interested in ways of accessing data from this server using R. Why don't you try RMySQL package? It works fine for me. [[alternative HTML version deleted]] __

Re: [R] Microsoft SQL Server and R

Yes, checking it now, I see you are absolutely right. Sorry for this confusion. Cheers, Olga On 27.01.2010, at 12:45, Orvalho Augusto wrote: RMySQL is for MySQL. MySQL is not MSSQL. There is not RMSSQL. So the best you have something to conneto MSSQL... RODBC or you may connect

[R] Nested Anova and type II sum of square

of the factors as I hope it does?) 2-      Why can’t I realize my nested anova with type II sum of square?   I apologize if my knowledge of statistics is so weak that I’ve make or ask something stupid… Many thanks in advance for time and consideration. Olga [[alternative HTML

Re: [R] read.csv from a remote machine

is solved. Cheers, Olga __ 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] read.csv from a remote machine

Dear all, I am trying to import data with read.csv and my file is on remote machine. I believe that I need to open a connection, not sure about syntax though. I would appreciate any suggestions, Thanks! Olga __ R-help@r-project.org mailing list

Re: [R] read.csv from a remote machine

, from this machine I establish ssh with a Linux machine. I run R on Linux, while all my data files are stored on Mac. So in this case although physically I am using Mac, it is in fact remote. I hope it answers your question. Thanks again, Olga