similar to: replication of latin squares

Dear all, this is a newbie's question. I have a 4x4 latin square replicated 3 times. That is: operators batches 1 2 3 4 1 A B C D 2 B C D A 3 C D A B 4 D A B C operators batches 1 2 3 4 5 A B C D 6 B C D A 7 C D A B 8 D A B C operators batches 1 2 3 4 9 A B C D

Hi all, The package crossdes could contruct a complete sets of mutually orthogonal latin squares. The construction works for prime powers only. I hope to know whether there is a way to construct a mutually orthogonal Lation square for 10 or other numbers that could not be prime powers. Thanks for any suggestions. Best wishes, Jinsong Zhao

Hello. I am working with a set of variables which are in columns. Three of them are of the same length while one has a different length. Typing '>data<-matrix(c(ca$value, mo$value,b2$value, y1), ncol = 3)' appears to read any three columns out of the four, though I can say exactly which of these three columns are read. However, typing '> data<-matrix(c(ca$value,

Dear R-gurus, being very new to R, (as well as lazy and not too smart !) I have some problems (and get lost in the docs) trying to write something to find the 9 values (A1,B1,C1,A2,B2,C2....C3) which are solutions of a 12 equations system of the form : > x1-(A1/(A1+B1+C1)) = 0 > y1-(B1/(A1+B1+C1))= 0 > z1-(C1/(A1+B1+C1)) = 0 > 3 same equations with subscript 2 > 3 same

Hi, I am running a simulation and have to perform ANOVA to determine the rank of factors. Used the aov() function and it works great for full factorial design. 1. For a massive set of data, I tried using biglm, while it can create the linear model, all the residuals (for assumption validation) are not recorded and the sum of squares are not there, just the estimated regression coefficient, 95%

Deal all: i want to do the probit with sample selection estimation, the following is my code: probit with sample selection can be done by stata :heckprob The heckprobll is the likelihood function shown in W.H. Greene 5th p714 ¡´ The question is the convergence is very slow compare with Stata using likellihood only. ¡´ Second i did the similar way in matlab using fminsearch , the estimated

Dear R-helpers, I have been trying to figure out how to plot a graph with an axis label consisting of a mixture of Latin, Greek and subscript characters. Specifically, I need to write A[beta]{1-42}, where A is Latin script A, [beta] is Greek lower case beta and {1-42} is subscript '1-42'. I can use xlab=expression(beta[1-42]) and obtain the [beta]{1-42} part of the label. But, I can't

I have a dataset that is basically structureless. Its dimension varies from row to row and sep(s) are a mixture of tab and semi colon (;) and example is HEADER1 HEADER2 HEADER3 HEADER3 A1 B1 C1 X11;X12;X13 A2 B2 C2 X21;X22;X23;X24;X25 A3 B3 C3 A4 B4 C4 X41;X42;X43 A5 B5 C5 X51 etc., say. Note that a blank

I think you should take this conversation private or seek local statistical expertise. This is about strategies for analysis, not about programming in R. Cheers, Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Wed, May 31, 2017 at

Hi there, I want to estimate simultaneous equation model with panel data. The model looks as follows Y1=a0+a1*X1+a2*X2 Y2=b0+b1*X2+b2*X1 X1=Z1-(Y1/a1) X2=Z2-(Y2/b1) I In this model Y1, Y2, X1 and X2 are endogenous variables; Z1, Z2 are exogenous variables and a0, a1, a2, b0, b1 and b2 are parameters. Could any one please help me how to estimate this model in R. Thanking you in anticipation

I want to measure the error in the estimation of a 2 dimensional density function that is calculated using an integral but run into problems trying to integrate with adaptIntegrate because the integrand also calls the function adaptIntegrate. In particular I want \int \hat{f}(x,y) - f(x,y) dx dy where \hat{f}(x,y) = \int K(a,b, x, y) da db and in this simulation study I know what the true value

Hello, I got one question on the Latin hypercube sampling. suppose there are three variables a, b, c, all of them follow the normal distribution. the mean value and standard deviation for each are  a(32, 2), b(35,5), c(37,3). I would like to use Latin hypercube sampling to random generate 1000 samples. but it needs to satisfy the condition that a<b<c. How can I implement this

Hi, I am having an issue when creating a postscript file from RGL window. It seems to cut off some of the axis labels. Here is the code I am using. I created a 3D plot using RGL_0.71 with R 2.5 on Windows XP. z1<-c(5,4,1,4.5,2,3,2,1,1) z2<-c(6,8,7,7.5,5,3.5,4,1,1) z3<-c(3,2,4,7,3,4.5,6,2,3) x1<-seq(1,9) x2<-seq(1,9) x3<-seq(10,18) y1<-seq(8,0) y2<--1*y1

Hi, In the FAQ I read about options to specify different fonts than the default ones for the console (in the file Rprofile) and for the graphical output (Rdevga). I would however like to replace Latin 1 with Latin 2 enconding for both (console and graphical) output in Windows and just graphical output in Linux. I guess it is possible but I did not find the way. How can I use the fonts

Hi all, I am writing just to keep a record for the solution for the problem we had with Notes Client. The problem is that latin characters (accentuation) appears strange at the emails using Notes Client. In order to the characters appears correctly you will have to *eliminate* the trebuchet font (either from msttcorefonts package or trebuc32.exe, installed via wine -

Hi all, I am a begginer in R and I have been trying to use the Conditioned Latin Hypercube to choose sample points only in areas close to roads due to the difficult thorough access in the study area. I could use a code to create the points throughout the area, but I need to create a code to make the R only comes up with points close to this road. I've been using the packages

Hi I want to du a sensitivity analysis using Latin Hypercubes. But my parameters have to fulfill two conditions: 1) ranging from 0 to 1 2) have to sum up to 1 So far I am using the lhs package and am doing the following: library(lhs) ws <- improvedLHS(1000, 7) wsSums <- rowSums(ws) wss <- ws / wsSums but I think I can't do that, as after the normalization > min(wss) [1]

Dear R experts I am looking for a package which gives me latin hyper cube samples from the grid of values produced from the command "expand.grid". Any pointers to this issue might be very useful. Basically, I am doing the following: > a<-(1:10) > b<-(20:30) > dataGrid<-expand.grid(a,b) Now, is there a way to use this "dataGrid" in the package

Hi all, I am attempting to use latin hypercube sampling to sample different variable functions in a series of simultaneous differential equations. There is very little code online about lhs or clhs, so from different other help threads I have seen, it seems I need to create a probability density function for each variable function, and then use latin hypercube sampling on this pdf. So far, I

Hi R-programmers ! I would like to perform a linear model regression month by month using the 'lm' function and i don't know how to do it. The data is organised as below: Month ExcessReturn Return STO 8 0.047595875 0.05274292 0.854352503 8 0.016134874 0.049226941 4.399372005 8 -0.000443869 0.004357305 -1.04980297 9 0.002206554 -0.089068828 0.544809429 9 0.021296551