Re: [R] elementwise matrix multiplication with a special structure

s, so i don't need to pass by everything, but still) which is quickly too heavy for R but still doable in C.  Thanks a lot and have  a great day ahead!Eugen From: David L Carlson To: eugen pircalabelu Sent: Friday, October 23, 2015 9:21 PM Subject: RE: [R] elementw

[R] elementwise matrix multiplication with a special structure

Hello R users, I have the following annoying matrix multiplication and I do not know how to avoid the nested loops, so maybe someone can help me with some ideas. I have searched the forum for past posts but nothing came up.Here it is:  aa=matrix(1:9,3,3) bb=matrix(seq(10,90,by=10),3,3) cc=matrix(

[R] expand.grid using a matrix and a vector as input

  Hello R-users,   I have the following question, for which my search did not really return any usable result. If I have a matrix a1, and a vector a2 like below   a1<-matrix(c(1:4),2,2) a2<-c(8,9)   is there any function like the expand.grid (or some clever calling of the function) such that it

Re: [R] Glmnet lambda value choice

 Checked the authors' paper on page 8 it explains how the \lambda sequence is chosen. I guess this is what you were looking for. http://www.stanford.edu/~hastie/Papers/glmnet.pdf   Eugen Pircalabelu (0032)471 842 140 (0040)727 839 293 From: Noah Silverma

Re: [R] SOLVED - export tables to excel files on multiple sheets with titles for each table

Hello R-users, I was able using RExcel and the VBA script offered by Eric Neuwirth (modifying it a little bit) to import dataframes on separate Excel sheets, each with its own particular title. As such the problem is solved. Thank you all for your suggestions. Have a nice day ahead! Eugen

[R] export tables to excel files on multiple sheets with titles for each table

hing for some time for a good solution. Thank you very much and have a great day ahead! Eugen Pircalabelu (0032)471 842 140 (0040)727 839 293 __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] Re :help to aggregate data

levels that have observation # tapply(dummy.data[,3], dummy.data[,1:2], sum) Eugen Pircalabelu (0032)471 842 140 (0040)727 839 293 - Original Message From: Uwe Ligges To: Mohan L Cc: r-help@r-project.org Sent: Thu, June 10, 2010 10:47:20 AM Subject: Re: [R] Re :help to aggreg

Re: [R] Help with TukeyHSD

Hi, Is this what you want? ?TukeyHSD summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks)) tk.fm1<-TukeyHSD(fm1, "tension", ordered = TRUE) tk.fm1 (tk.fm1[[1]][,4]) Eugen Pircalabelu (0032)471 842 140 (0040)727 839 293 - Original Message From: Amit Pat

Re: [R] Mosaic Plots

", shade = TRUE, legend = TRUE) Eugen Pircalabelu (0032)471 842 140 (0040)727 839 293 - Original Message From: Sunitap22 To: r-help@r-project.org Sent: Mon, March 22, 2010 12:21:14 PM Subject: [R] Mosaic Plots Hello Everyone I want to plot Moasic Plots, I have tried them using i

[R] ascending permutation within treatment

! Eugen Pircalabelu (0032)471 842 140 (0040)727 839 293 __ 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

[R] Changes in Predicted Probabilities for multinomial regression as in STATA software

c/polytomous-effect-displays/polytomous-effect-displays.pdf but could not obtain the changes in predicted probabilities. Can anyone give me a hint on how to obtain them? Thank you very much and have a great day ahead! Eugen Pircalabelu (0032)471 842 140 (0040)727 839 293 ___

Re: [R] ImageMagick not seen by movie3d function from rgl package --SOLVED

after the software installation, I should have loaded rgl package again in R. So now it works! Thank you once again! - Original Message From: Duncan Murdoch To: eugen pircalabelu Cc: R-help Sent: Thu, October 15, 2009 6:31:58 PM Subject: Re: [R] ImageMagick not seen by movie3d functio

[R] ImageMagick not seen by movie3d function from rgl package

Hello R-users, I have downloaded and installed a binary version for Windows OS of ImageMagick (ImageMagick6.5.6-10-q16-windows-dll.exe), I have installed the rgl package and i've tried running the following example from ?movie3d() library(rgl) open3d() plot3d( cube3d(col="green") ) M <- par3

Re: [R] crosstabulation and unlist function

Thank you! This is what i needed! I did not realize i could replicate my factor and then tabulate, but it makes sense! Thank you once again! - Original Message From: William Dunlap To: eugen pircalabelu ; r-help@r-project.org Sent: Mon, October 12, 2009 10:48:36 PM Subject: RE

Re: [R] crosstabulation and unlist function

erested in the NA occurence). Hopefully, this time my question was a lot more clear. Thank you very much ! - Original Message From: David Winsemius To: David Winsemius Cc: eugen pircalabelu ; R-help Sent: Mon, October 12, 2009 9:36:39 PM Subject: Re: [R] crosstabulation and un

[R] crosstabulation and unlist function

e do a crosstabulation between table(unlist(df[,1:3])) and df$dd? I want to find out when dd==A (or B, or C) how many times do the values 1, 2 ,3,.. appear in df[,1:3]? Thank you very much! Eugen Pircalabelu __ R-help@r-project.org mailing list h

[R] Can not replicate estimates with rScreen function from ROSSI "Bayesian Statistics and Marketing"

Hi R-users, I have the following problem: I am trying to learn something about bayes methodology and started paying around bayesm package, but could not replicate the Conjunctive model's estimates as they appear in Rossi et al "Bayesian Statistics and Marketing", 2005, JWS, pages 264-265, Tabl

Re: [R] joining "one-to-many"

A not so elegant way of obtaining your result (compared to merge()) would be: > t1$ val3<-rep(t2$val3, table(t1$loc)) > t1$ val4<-rep(t2$val4, table(t1$loc)) Eugen. --- On Tue, 2/17/09, Monica Pisica wrote: > From: Monica Pisica > Subject: Re: [R] joining "one-to-many" > To: ggrothendi...

[R] Beginner question about an Error message when loading SWinTypeLibs package

Hi List, I am trying to learn something about the exchange facilities between R and Windows Application with the RDCOMClient package from (http://www.omegahat.org/RDCOMClient/) and I have faced 2 problems: 1. I could not follow the link for the binary version of SWinTypeLibs package on the sit

[R] Question on Sweave-Latex and examples in the Sweave Manual

Hi List, I have a problem with using Latex and Sweave for creating a document. So I downloaded the Sweave manual from http://www.statistik.lmu.de/~leisch/Sweave/Sweave-manual.pdf and i have tried to replicate the example on pages 4-5, but i encounter the following problem: IT DOES NOT WORK. I

[R] Question On CrossTable function in gmodels package

Hi R-users, I have the following problem with CrossTable function within “gmodels” package: the output of the function (format “spss” and asresid=T) can not be stored within another object. For example: >library(gmodels) >data(infert, package = "datasets") > CrossTable(infert$education, infert$i

[R] newbie question on tcltk

Hi List, Can anyone tell me how could i put the "BACK" button in the following code, just under the "AAA" menu? I want this button to go back to the previous page, and since it has nothing to do with the "1" and "2" buttons, i want it somehow separated from these two buttons, but i don't know h

[R] question on lmList

led about the general error messages from different packages, somewhere? Thank you very much and have a great day ahead! Eugen Pircalabelu __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gu

[R] question on xyplot

(for eg SUBNUM 7) without picking them by inspecting the plot. Is there some way how I could access these coefficients for this apparently regression line? Thank you very much and have a great day ahead! Eugen Pircalabelu __ R-help@r-project.org

[R] Question on glm.nb vs zeroinfl vs hurdle models

Good afternoon, I’m in need of an advice regarding a proper use of glm.nb, zeroinfl or hurdle with my dataframe. I can not provide a self-contained example, since I need an advice on this current dataset and its “contradictory” results. So i have a dataset which contains 1309 cases

[R] Question regarding lmer vs glmmPQL vs glmm.admb model on a negative binomial distributed dependent variable

best suited for a highly skewed dependent variable (running these two commands on my data set the diffrence in SE estimation still remains, and I want to use a correct algorithm for my data set) .   Thank you very much and have a great day ahead!      Eugen Pircalabelu [[alterna

[R] R vs SAS and HLM on multilevel analysis- basic question

Hi R users! I am trying to learn some multilevel analysis, but unfortunately i am now very confused. The reason: http://www.ats.ucla.edu/stat/hlm/seminars/hlm_mlm/mlm_hlm_seminar.htm http://www.ats.ucla.edu/stat/sas/seminars/sas_mlm/mlm_sas_seminar.htm and MlmSoftRev. pdf from mlmRev package.

[R] multilevel basic lme question

Hi R users I want to use the lme package for a multilevel analysis on the following example: > math<-c(2, 3,2, 5, 6 ,7 , 7) > sex<-c(1, 2, 1, 2, 2, 2, 1) > school_A<-c(1,1,1,2,2,2,2) > school_B<-c(10,10,10,20,20,20,20) > mydata<-data.frame(math, sex, school_A, school_B) > mydata School_A and sc

[R] Is there a package assuming a DIRICHLET model for Marketing Consumer Panel

Hi List, I was wondering if there is a package in R that uses the DIRICHLET model for panel data like in Goodhart, Ehrenberg, Chatfield 1984. The model uses a combination of distributions and is specifically designed for marketing panel data, using input data such as penetration, brand share,

[R] factor analysis for binary data with FactoMineR package

Hi List, I'm following on the discussions regarding the use of exploratory factor analysis with binary data, from nabble: http://www.nabble.com/PCA--for-Binary-data-to11050371.html#a11092127 and i tried using FactorMineR package but unfortunately i have trouble understanding what are the typ

Re: [R] SOLVED plotting time series

Thank you very much for your replies , but what i don't understand from jim's reply is this command dates <- seq(tsp(x)[1]+5/12, tsp(x)[2], .5) Removing the +5/12 gives me the perfect solution, while keeping it makes the axis run from the second value plotted and the first value from my ts

[R] plotting time series

Hi List, I have the following time series and i want to be able to plot it while having the x-axis running from 1998 to 2006 but in a bi-annual format So here is my example: x<-sample (10:100, 17) x<-ts(x, start=c(1998,1), frequency=2) plot(x) When I plot the ts it goes from 1998 to 2006 by

Re: [R] question for aov and kruskal

Hi, My data was only a toy example that matched the real situation, with real data, but i could not have posted the entire data.set and so i gave a self contained example of what i thought was my problem. Of course you can see with the naked eye that the data is unbalanced, (this was done inten

[R] question for aov and kruskal

Hi R users! I have the following problem: how appropriate is my aov model under the violation of anova assumptions? Example: a<-c(1,1,1,1,1,1,1,1,1,1,2,2,2,3,3,3,3,3,3,3) b<-c(101,1010,200,300,400, 202, 121, 234, 55,555,66,76,88,34,239, 30, 40, 50,50,60) z<-data.frame(a, b) fligner.test(z$b, fa

Re: [R] question for crosstable

Yes, Mark! Thank you! I was repeatedly using the function with two vectors and i have forgotten that it could be used also with a matrix. Thank you once again! Marc Schwartz <[EMAIL PROTECTED]> wrote: eugen pircalabelu wrote: > Good evening R-users! > > I have the following

[R] question for crosstable

Good evening R-users! I have the following problem: i want to get a weighted crosstable along with the adjusted standardized residuals test Example: a<-c(1,1,1,1,1,2,2,2,3,3) b<-c(10,10,23,24,33,45,46,70,21,66) c<-c(3,3,2,3,4,1,1,1,3,3) d<-c("a", "b","b","c","a","a","a", "b", "c","c") z<-da

[R] question regarding using weights in the hierarchical/ kmeans clustering process

Hi R users! I have a bit of a problem with using an hierarchical clustering algorithm: a<-c(1:15) b<-rep(seq(1:3), 5) c<-rnorm(15, 0,1) d<-c(sample(1:100, 15, replace=T)) e<-c(sample(1:100, 15, replace=T)) f<-c(sample(1:100, 15, replace=T)) data<-data.frame(a,b,c,d,e,f) q<-data.frame(data

[R] multinomial confidence interval

Hi R users! Is there a function that extracts the simultaneous confidence intervals for multinomial proportions as described by Sison and Glaz 1995? or anyone else for that matter? I have seen that SAS has macro for this ( http://www.jstatsoft.org/v05/i06) and i was wondering if R had somet

[R] question regarding kpss tests from urca, uroot and tseries packages

Hi R users! I've come across using kpss tests for time series analysis and i have a question that troubles me since i don't have much experience with time series and the mathematical part underlining it. x<-c(253, 252, 275, 275, 272, 254, 272, 252, 249, 300, 244, 258, 255, 285, 301, 278, 279,

[R] question regarding hypothesis testing in contingency tables

Hi R-users! I have the following example: a<-data.frame(cat=c(5,10,15), dog=c(5,10, 15), mouse=c(10,10,20)) b<-data.frame(cat=c(15,10,5), dog=c(15, 10, 5), mouse=c(20,10,10)) rownames(b)<-c("scared", "happy", "sad") rownames(a)<-c("scared", "happy", "sad") Let's say that a and b are 2 contingency

[R] 4 questions regarding hypothesis testing, survey package, ts on samples, plotting

Good morning! I have 4 questions which trouble me: 1. I want to test the hypothesis that the 2 proportions (the mean of a binomial) which come from 2 different samples are equal. I want to use the following function z= (p1-p2)/ sqrt((p1(1-p1)/n1)+(p2(1-p2)/n2)) which is one of the standard f

[R] question regarding arima function and predicted values

Good evening! I have a question regarding forecast package and time series analysis. My syntax: x<-c(253, 252, 275, 275, 272, 254, 272, 252, 249, 300, 244, 258, 255, 285, 301, 278, 279, 304, 275, 276, 313, 292, 302, 322, 281, 298, 305, 295, 286, 327, 286, 270, 289, 293, 287, 267, 267, 288, 30

[R] modeling time series with ARIMA

Good afternoon! I'm trying to model a time series on the following data, which represent a monthly consumption of juices: >x<-scan() 1: 2859 3613 3930 5193 4523 3226 4280 3436 3235 3379 3517 6022 13: 4465 4604 5441 6575 6092 6607 6390 6150 6488 5912 6228 10196 25: 7612 72

[R] question regarding an error message from survey package

Hi R-users! Can anyone tell me what this error message refers to exactly > a<-c('m','f','m') > b<-c(1:3) > st<-c(1:3) > data.frame(a,b,st) a b st 1 m 1 1 2 f 2 2 3 m 3 3 > data.frame(a,b,st)->l > p<-svydesign(ids=~1, data=l) > pop<-c('(intercept)'=100, st1=10,st2=20) > calibrate(p,~st, pop)

[R] problem modeling time series

Good afternoon! I'm trying to model a ts but unfortunately i'.m very new to this kind of modeling so i 'll be very grateful if you have an advice. This was my syntax: t<-c(16115,17391,19011,20256,19034,18851,20016,18088,19166,21163,18463,19397,15800,16113,18879,20598,17252,19753,19110,19605,2

[R] question regarding time series packages

Good afternoon! I'm trying to learn time series but i have a bit of of a problem using R packages for this. 1. > LakeHuron > sample(500:600, 98) > sample(500:600, 98)->t > fit<-arima(LakeHuron, order=c(2,1,1), xreg=t) > fit > predict(fit, n.ahead=1, newxreg=t) Now, my problem is this: is it ok

[R] calibration question

Good afternoon! I have the following data: a<-data.frame (id_hh=c(1:5), strata=c(1,1,2,2,1), Nhstrata=c(100,100,200,200,100), Nrmemb=c(2,4,2,5,4)) a$ocmemb1<-c("wk","jl","st","jl","st") a$ocmemb2<-c("wk","jl","st","wk","wk") where id_hh is a code of identification for the household (my analys

[R] [non-R question] methodological problem

Sorry, I meant non-R question. Sorry for the inconvenience! __ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

[R] [non-statistics question]methodological problem

Good afternoon! As mentioned in the subject, my question regards more the methodological part that accompanies survey design and the statistical part that is involved. So, I have the following data: a<-data.frame (id_hh=c(1:5), strata=c(1,1,2,2,1), Nhstrata=c(100,100,200,200,100), Nrmemb=c(2,4,

[R] inference from dual frame surveys using R

Good afternoon! My question is more of a "is there a package for doing?" an inference regarding the mean, median, regression estimates, etc. by using the information from dual frame surveys? (this methodology is based on the work of BANKIER 1986, SKINNER 1991, LOHR and RAO 2000...)

[R] combining data from 2 independent surveys

Good afternoon! My question is more of a "is there a package for doing?" If one has 2 independent studies which both independently measure a trait of a population, but both of them are constructed from 2 separate frames (which TOGETHER refer to the population of interest),

[R] grossing-up weights in Survey package

Good afternoon! I can not find anywhere the formula by which the Survey package calculates the grossing-up weights for each element in a stratified sample which is then post-stratified by other variables. Can anyone help me? Thank you and have a great day! ---

[R] weighting question

Hi R-users, Can anyone tell me where can i find info about they way how post stratification weights are calculated when i have an already stratified survey design, especially in Survey Package (but any theoretical material would do me just fine) ? Thank you and have a nice day!

[R] The Survey package again

Good afternoon! 1 I have a bit of a problem with the Survey package, but one that is more theoretical than practical. If for example on wishes to use some additional information which comprises of 4 variables with 7, 5,5 and 5 categories (resulting in 875 post-strata if the whole crossin

[R] Survey package question

Good afternoon! 1 I have a bit of a problem with the Survey package, but one that is more theoretical than practical. If for example on wishes to use some additional information which comprises of 4 variables with 7, 5,5 and 5 categories (resulting in 875 post-strata if the whole crossin