[R] An iterative function

Dear all, The function f() below is a function of m1 and m2, both of which are matrices with 3 rows. The function works sequentially one row after another. So altogether there are three stages. I am trying to update the coding to write a generic function that will work for arbitrary k stages. I

[R] question with integrate function

Hi all, The function h below is a function of c and it should be a monotone increasing function since the integrand is nonnegative and integral is taken from c to infinity. However, as we can see from the plot, it is not shown to be monotone. Something wrong with the usage of integrate function?

Re: [R] question with integrate function

Sorry. I meant in the previous email that the function h() is a monotone decreasing function. Thanks very much. 2018-02-06 13:32 GMT-05:00 li li : > Hi all, > The function h below is a function of c and it should be a monotone > increasing function since the integrand is nonneg

Re: [R] question with integrate function

tely constant (in > particular, zero) over nearly all its range it is possible that the result > and error estimate may be seriously wrong." > > You could try to integrate over a finite interval, say (7, 12). > > Göran Broström > > > On 2018-02-06 19:40, li li wrote: > >&

[R] Looking for

Dear all, For a piecewise function F similar to the attached graph, I would like to find inf{x| F(x) >=0}. I tried to uniroot. It does not seem to work. Any suggestions? Thank you very much!! Hanna F.pdf Description: Adobe PDF document ___

Re: [R] Finding Infimum in R

al. > > 2. Can we assume monotonic increasing, as is shown? > > > -- 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 "Bl

Re: [R] Finding Infimum in R

Yes. If the function f takes the value zero at some discontinuity point, then the code gives the inf of the set I described. Otherwise, it is an approximation since we need to worry about numerical accuracy. 2017-04-10 14:08 GMT-04:00 Bert Gunter : > Well, I haven't checked carefully, but of cou

[R] creat contingency tables with fixed row and column margins

Hi all, Is there an R function that can be used to enumerate all the contingency tables with fixed row and column margins. For example, can we list all 3 by 3 tables with row margins 3,6,6 and column margins 5,5,5. Thanks very much! Hanna [[alternative HTML version deleted]] _

Re: [R] creat contingency tables with fixed row and column margins

Hi Dave, Thanks for your reply but the post in your link above was certainly asked by me. Hanna 2017-05-27 16:04 GMT-04:00 David Winsemius : > > > On May 27, 2017, at 12:49 PM, li li wrote: > > > > Hi all, > > Is there an R function that can be used to en

Re: [R] creat contingency tables with fixed row and column margins

I meant that the post in your link above was NOT asked by me. 2017-05-27 16:53 GMT-04:00 li li : > Hi Dave, > Thanks for your reply but the post in your link above was certainly > asked by me. >Hanna > > 2017-05-27 16:04 GMT-04:00 David Winsemius : > >> >&g

[R] Question on function "scatterplot3d"

Hi all, I have a question with regard to making plots using function "scatterplot3d". Please see the example below. It looks like, for y axis, the tickmark text was cutoff. The number "10" does not show up completely. I tried to work with par(mpg). It does not seem to work. Hope to get some advic

Re: [R] [FORGED] Re: Question on function "scatterplot3d"

Thanks so much for the help!! Hanna 2017-06-01 11:02 GMT-04:00 Uwe Ligges : > > > On 01.06.2017 10:03, Rolf Turner wrote: > >> On 01/06/17 19:54, Uwe Ligges wrote: >> >>> A design flaw, whether the labels are cut depends somewhat on the sizce >>> of the device, hence there is the argument >>>

[R] subletting an array according to dimnames

Hi all, I have a three dimensional array with the corresponding dimension names. I would like to subset the array according to the dimension names. For example, suppose I want to extract the values corresponding to A=20, B=10, C=0. I know I can do: P2[dimnames(P2)$A==20, dimnames(P2)$B==10, dim

Re: [R] subletting an array according to dimnames

That works. Thank you! 2017-06-02 8:00 GMT-04:00 Adams, Jean : > Have you tried P2["20", "10", "0"] ? > > Jean > > On Thu, Jun 1, 2017 at 3:10 PM, li li wrote: > >> Hi all, >> I have a three dimensional array with the corres

[R] An R question

Hi all, In checking my R codes, I encountered the following problem. Is there a way to fix this? I tried to specify options(digits=). I did not fix the problem. Thanks so much for your help! Hanna > cdf(pmass)[2,2]==pcum[2,2][1] FALSE> cdf(pmass)[2,2][1] 0.758> > pcum[2,2][1] 0.75

Re: [R] An R question

https://www.researchgate.net/profile/Ivan_Calandra > > > On 07/06/2017 16:32, li li wrote: > >> Hi all, >>In checking my R codes, I encountered the following problem. Is there a >> way to fix this? >> I tried to specify options(digits=). I did not fix the pro

[R] Adding zeros in each dimension of an array

For a data frame, we can add an additional row or column easily. For example, we can add an additional row of zero and an additional row of column as follows. Is there an easy and similar way to add zeros in each dimension? For example, for array(1:12, dim=c(2,2,3))? Thanks for your help!! Han

Re: [R] Adding zeros in each dimension of an array

g that last expression without hardcoding the sequences > may not be obvious: > > expandArray <- function( a ) { > b <- array( 0, dim = dim( a ) + 1 ) > do.call( `[<-`, c( list( b ), lapply( dim( a ), seq.int ), list( a ) ) ) > } > expandArray( a ) > > > On

[R] problem with "unique" function

I have the joint distribution of three discrete random variables z1, z2 and z3 which is captured by "z" and "prob" as described below. For example, the probability for z1=0.46667, z2=-1 and z3=-1 is 2.752e-13. Also, the probability adds up to 1. > head(z) z1 z2 z3 [1,] -0.4666

Re: [R] problem with "unique" function

case, it seems like all numbers > are multiples of 1/30, so another idea could be to multiply by 30, round, > and divide by 30. > > -pd > > > On 28 Jul 2017, at 17:17 , li li wrote: > > > > I have the joint distribution of three discrete random variables z1, z2 >

[R] read multiple sheets of excel data into R

Hi all, I tried to use the package "XLConnect" to read excel data into R. I got the following error message: Error : .onLoad failed in loadNamespace() for 'rJava', details: call: fun(libname, pkgname) error: No CurrentVersion entry in Software/JavaSoft registry! Try re-installing Java and m

[R] model specification using lme

Hi all, For the following data, I consider the following random intercept and random slope model. Denote as y_ijk the response value from *j*th individual within *i*th method at time point *k*. Assume the following model for y_ijk: y_ijk= (alpha_0+ tau_i +a_j(i))+(beta_i+b_j(i)) T_k + e_ij

[R] Fwd: model specification using lme

kelinx Date: 2016-05-30 4:40 GMT-04:00 Subject: Re: [R] model specification using lme To: li li Cc: r-help Dear Hanna, None of the models are correct is you want the same random intercept for the different methods but different random slope per method. You can random = ~ 1 + time:

[R] multcomp package

Hi all, After fitting a random slope and random intercept model using lme function, I want to test whether each of the fixed slopes is equal to zero (The output of model is below). Can this be done (testing each individual slope) using multcomp package? Thanks much for the help. Hanna > sum

Re: [R] [R-sig-ME] multcomp package

Min Q1 Med Q3 Max -2.67575293 -0.51633192 0.06742723 0.59706762 2.81061874 Number of Observations: 69 Number of Groups: 7 2016-06-06 11:21 GMT-04:00 Gabriel Baud-Bovy : > On 06/06/2016 4:57 PM, li li wrote: > > Hi all, > After fitti

[R] Fixed Effects in lme function

Dear all, For the data below, I would like to fit a model with common random slope and common random intercept as shown below. I am interested in obtaining separate fixed effect estimates (intercept and slope and corresponding hypothesis test) for each method. Instead of performing the analys

[R] Help on improving an algorithm

Hi all, I have four matrix tmp_a, tmp_b, tmp_c and tmp_d whose dimensions are shown as below. I want to take one row from each of these matrices and then put the four selected rows into a matrix. I want to count the number of such matrices for which the vector of row sum is less than or equal to (

Re: [R] Help on improving an algorithm

t2) > print(torig) > > > A couple of minor points: > dim(tmp_a)[1] can be replaced by nrow(tmp_a) > in the original code, apply() can be replaced by rowSums() > It might be faster if rowSums() was replaced by .rowSums(); I haven't > tried that. > > -Don &

[R] plot.drm in "drc" package

Hi all, I am trying to use the drc package to fit 4L curve. I am so confused about the plot.drm function in the drc package. Particularly, I am confused about the scale of the xaxis in the plots generated using the plot.drm function. See the example below: ## generate data and fit the model dose

[R] plot with different symbols and colors according to the factor levels

Hi all, I have the following data. I want to plot the data (y ~ conc) with different symbols and colors corresponding to different levels of the factor sample. I could create a column with color and pch and then do the plot, but I am sure there are much better ways. Can anyone make suggestions?

Re: [R] plot with different symbols and colors according to the factor levels

..@math.utah.edu > Department of Ecology VOICE: (360) 407-6815 > PO Box 47600FAX:(360) 407-7534 > Olympia, WA 98504-7600 > > USPS: PO Box 47600, Olympia, WA 98504-7600 > Parcels:300 Desmond Drive, Lace

Re: [R] plot.drm in "drc" package

1 GMT-04:00 PIKAL Petr : > Hi > > Thanks for code. > see in line > > > -Original Message- > > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of li li > > Sent: Tuesday, August 30, 2016 5:07 PM > > To: r-help > > Subject: [R] plot.drm i

Re: [R] Question with uniroot function

-- > Sent from my phone. Please excuse my brevity. > > On April 15, 2015 7:20:04 PM PDT, li li wrote: > >Hi all, > >In the following code, I am trying to use uniroot function to solve for > >the root (a and b in code below) for function f

Re: [R] Question with uniroot function

Thank you. 2015-04-16 12:33 GMT-04:00 William Dunlap : > Use optimize() to find the minimum and feed that value into uniroot(). > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Thu, Apr 16, 2015 at 7:47 AM, li li wrote: > >> Hi Jeff, >> T

[R] Random effect in ancova model

Hi all, I have the following data and would like to consider the following model: value = type + batch (type) + beta1* month +beta2*type*month+ error Here type is fixed effect and batch is a random effect neste within type, and month is continuous. Since there is a random effect in this model.

[R] predict function in regression analysis

Hi all, I have the following data in which there is one factor lot with six levels and one continuous convariate time. I want to fit an Ancova model with common slope and different intercept. So the six lots will have seperate paralell regression lines.I wanted to find the upper 95% confidence li

[R] dotplot

Hi all, I wrote the following code and have two questions: (1) As you can see, I would like different colors for different types. It does not come out that way in the graph from this code. Anyone know how to fix this? (2) How do I made the lots number on x axis eligible to read? (3) How d

[R] lme function to obtain pvalue for fixed effect

Hi all, I am using the lme function to run a random coefficient model. Please see output (mod1) as below. I need to obtain the pvalue for the fixed effect. As you can see, the pvalues given using the summary function is different from the resutls given in anova function. Why should they be diff

Re: [R] [R-sig-ME] lme function to obtain pvalue for fixed effect

viado desde mi iPhone > >> El 26/5/2015, a las 13:18, li li escribió: >> >> Hi all, >> I am using the lme function to run a random coefficient model. Please >> see >> output (mod1) as below. >> I need to obtain the pvalue for the fixed effect.

[R] different results from lme and lmer function

Hi all, I am fitting a random slope and random intercept model using R. I used both lme and lmer funciton for the same model. However I got different results as shown below (different variance component estimates and so on). I think that is really confusing. They should produce close results. Any

[R] How to add legend to a dotplot

Hi all, I wanted to add the legend to a dotplot using legend funciton. If does not seem to be working? Anyone have any suggestions? Thanks! Hanna __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/list

Re: [R] How to add legend to a dotplot

omy and Soil Science > University of New England > Armidale NSW 2351 > Email: home: mac...@northnet.com.au > > > > > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of li li > Sent: Wednesday, 3 June 2015 13:07 > To: r-

[R] Confidence interval for the mean in the random intercept and random slople model

Hi all, I am fitting a random slope and random intercept model usign lme fucntion as shown below. Type is factor with two levels. I would like to to find a confidence interval for mean of this model. Note that the variance we use in finding the confidence interval should include the variariance c

[R] Warning message when using lmer function

I got the following warning message when using the lmer function. Does anyone know what is the implication? Thanks! Warning message: In anova(model, ddf = "lme4") : bytecode version mismatch; using eval __ R-help@r-project.org mailing list -- To UNSUBSC

[R] Different random intercepts but same random slope for groups

Hi all, I'd like to fit a random intercept and random slope model. In my data, there are three groups. I want to have different random intercept for each group but the same random slope effect for all three groups. I used the following R command. However, there seems to be some problem. Any sugge

Re: [R] Different random intercepts but same random slope for groups

is done may be no more >> than asking him to perform a post-mortem examination: he may be able to >> say >> what the experiment died of. ~ Sir Ronald Aylmer Fisher >> The plural of anecdote is not data. ~ Roger Brinner >> The combination of some data and an aching desire

[R] title of r plots

Hi all, I would like to add "-70°C ± 10°C/Ambient" as the title of my plot. Could anyone give some help on this? Thanks. Hanna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see htt

Re: [R] title of r plots

-70*degree*C%+-%10*degree*C/Ambient)) > > Bill Dunlap > TIBCO Software > wdunlap tibco.com > > On Fri, Feb 27, 2015 at 7:27 PM, li li wrote: > >> Hi all, >> I would like to add "-70°C ± 10°C/Ambient" as the title of my plot. >> Could any

[R] xlims of box() function

Hi all, It looks like the limits for x axis for box function in R is from -1 to 1. How can I change this limit? Thanks for your help. Hanna [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and

[R] Wrong results from anova

Hi all, I have following data. When I perform an anova, the residual sum of square returns to be zero. But this is wrong, since we can hand calculate the RSS to be around 0.0308. Did anyone come across the same problem before? Any suggestions? Thanks. Hanna > ydata Y sample.Y 1

Re: [R] Wrong results from anova

0.003725088 is the mean square error. 2015-04-14 13:44 GMT-04:00 Michael Dewey : > See in-line > > > On 14/04/2015 18:23, li li wrote: > >> Hi all, >>I have following data. When I perform an anova, the residual sum of >> square returns to be zero. >&g

[R] Question with uniroot function

Hi all, In the following code, I am trying to use uniroot function to solve for the root (a and b in code below) for function f1. I am not sure why uniroot function does not give the answer since when we look the graph, the function does cross 0 twice. Any suggestion? Thanks. Hanna u1

[R] Identify row indices corresponding to each distinct row of a matrix

Hi all, I use the following example to illustrate my question. As you can see, in matrix C some rows are repeated and I would like to find the indices of the rows corresponding to each of the distinct rows. For example, for the row c(1,9), I have used the "which" function to identify the row i

Re: [R] Identify row indices corresponding to each distinct row of a matrix

6-14" "7-15" "8-16" "1-9" >> >"2-10" >> >"3-11" "4-12" "5-13" "6-14" >> >[15] "7-15" "8-16" >> >> ## Then maybe: >> >> tapply(seq_along(vec),vec, I)

Re: [R] Identify row indices corresponding to each distinct row of a matrix

t; Cm <- as.matrix( C ) > D <- which( !duplicated( Cm, MARGIN=1 ) ) > nCm <- nrow( Cm ) > F <- lapply( D, function(d) { >idxrep <- rep( d, nCm ) >which( 0 == unname( rowSums( Cm[idxrep,] != Cm ) ) ) > } ) > > > On November 8, 2018 1:42:40 PM PST, l

[R] list contingency tables

Hi all, I am trying to list all the 4 by 2 tables with some fixed margins. For example, consider 4 by 2 tables with row margins 1,2,2,1 and column margins 3,3. I was able to do it using the code below. However, as seen below, I had to first count the total number of tables with the specific row

[R] Finding MLE

Hi all, I am trying to use the mle function in R to find the maximum likelihood estimator. The ll function below is the negative of the log likelihood. Suppose x0 is the observed values, I want to find the maximum likelihood for a and b. After running the code below, I get the error message "Erro

Re: [R] Finding MLE

:00 Bert Gunter : > You do not appear to provide initial values for a and b , i.e. the > "start" argument for mle. > > Cheers, > Bert > > Bert Gunter > Genentech Nonclinical Biostatistics > (650) 467-7374 > > "Data is not information. Information is not k

Re: [R] Finding MLE

22:04 GMT-05:00 li li : > Thanks Bert for the reply. I still get a message when adding the start > argument. > > > > > n <- 8 > > x0 <- c(2,3) > > > > ll<- function(a,b,x=x0,size=n){ > + > -sum(log(gamma((n-1)/2+a-1)/(gamma((n-1)/2)

[R] Adding a second Y axis on a dotplot

Hi all, I plotted a dotplot based on the data below and code below. I would like to add another yaxis on the right with a different col, different tickmarks and a different label. Can anyone give some help?Thanks very much!! Hanna > tmp1 result lot trt trtsymb trtcol 1 98 lot1 tr

[R] Order of boxplots

Hi all, I have the following data "tmp" and I want to plot boxplots for each level of the factor "type" and the order the factor should be c, b ,a. In other words, the boxplot corresponding to the level "c" should be the first and so on. Any suggestions? Li > tmp result type 1 101

Re: [R] Order of boxplots

gt; Research Engineer (Solar/BatteriesO.O#. #.O#. with > /Software/Embedded Controllers) .OO#. .OO#. rocks...1k > ------- > Sent from my phone. Please excuse my brevity. > > On S

Re: [R] Order of boxplots

Thank you! That worked. 2015-09-15 2:50 GMT-04:00 peter dalgaard : > > > On 15 Sep 2015, at 04:31 , li li wrote: > > > > Hi Jeff, > > Thanks for replying. I actually tried "ordered(tmp$type, levels=c("c", > > "b", "a")."

[R] cumulative distribtuion function for multinomial distribution in R

Hi all, In R, is there a function for the cumulative distribution function for multinomial distribution? I only see pmultinom and rmultinom which are the prabability mass function and the function for generating multinomial random variables respectively. Thanks! Hanna _

[R] correlated binomial random variables

__ 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 and provide commented, minimal, self-contained, reproducible co

Re: [R] correlated binomial random variables

. Say if X is binomial (n, p1) and Y is binomial (n, p2) and the correlation between X and Y is rho and we want to calculate P(X <= c, Y <= c). Thanks so much! Hanna 2015-10-05 11:27 GMT-04:00, li li : > > __ R-help@r-project.

Re: [R] [FORGED] Re: correlated binomial random variables

, and X_2 are both Binomial(n, p). X_0, X_1, and X_2 are all independent. Then X, Y are correlated and P(X <= t, Y <= t) can be exactly calculated. Thanks! Hanna 2015-10-05 18:00 GMT-04:00, Rolf Turner : > On 06/10/15 04:43, li li wrote: >> Hi all, >> Using the &qu

[R] question with drm function in "drc package"

Hi all, I am trying to use the drm function in drc package to fit a 4 PL or 3PL curve for an assay response. Please see the listed data below. When I do the curve fitting, it returns the following error message. Anyone who familiar with this have any input on what went wrong? Thanks so much in

Re: [R] question with drm function in "drc package"

not corrupt your code. > > > http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example > -- > Sent from my phone. Please excuse my brevity. > > On December 7, 2015 1:05:32 PM PST, li li wrote: >> >> Hi Jeff, >> I am not sure where the lo

[R] Plot in drc package

Hi all, I have the following data and I fit a log logistic model using the drm function in DRC package. I saved the fitted model in an object called "mod". (See below) I understand that we can use the plot function to plot the regression curve with the original data points using the code "plot

[R] Error in loading the drc package

Hi all, When trying to load the drc package. I got the following error. Any suggestions? Thanks. Hanna > install.packages("drc", dependencies=TRUE) --- Please select a CRAN mirror for use in this session --- trying URL 'https://cran.fhcrc.org/bin/windows/contrib/3.1/drc_2.5-12.zip' Conte

Re: [R] [FORGED] Error in loading the drc package

Thanks. 2015-12-08 15:20 GMT-05:00 Rolf Turner : > > Please keep communications on-list. Others may have relevant comments and > suggestions to make. > > > On 09/12/15 09:00, li li wrote: > >> Thanks for the reply. So a newer version of R can solve the problem? But &g

[R] change the x axis tickmarks when using plot function in drc package

Hi all, When plotting the dose response curve using plot function as in the example codes below, the scale of the axis should really be on the log scale of the dose given the shape of the graph. But as you can see, the tickmarks of the returned graph represent the original scale. How can I change

Re: [R] change the x axis tickmarks when using plot function in drc package

uot;,log="xy") > > Jim > > > On Wed, Dec 9, 2015 at 1:57 PM, li li wrote: > >> Hi all, >> When plotting the dose response curve using plot function as in the >> example codes below, the scale of the axis should really be on the log >> scale of the dose

Re: [R] change the x axis tickmarks when using plot function in drc package

Thanks for your help Giorgio! Both options worked. 2015-12-09 11:46 GMT-05:00 Giorgio Garziano : > Looking at the source code of the package drc, there is something that may > somehow explain what > you are experiencing: > > file: plot.drc.R, function addAxes(), lines 543-626 > > ceilingxTicks <-

[R] Restricted 4-PL curves using drc package

Hi all, In drc package, is there a function which can be used to fit restricted 4PL curves? For example, we restrict two 4PL curves have the same lower and upper asymptotes? Thanks for the help in advance! Hanna [[alternative HTML version deleted]]

Re: [R] Restricted 4-PL curves using drc package

n mind is that people keep coming along > and sticking things into it." > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) > > > On Wed, Dec 9, 2015 at 12:27 PM, li li wrote: > > Hi all, > > In drc package, is there a function which can

[R] exact trend test

Hi all, Is there an R function that does exact randomization trend test? For example, consider the 2 by 5 contingency table below: dose0dose 0.15dose 0.5dose 1.5dose 5 row margin Yes 43 4 5 8

Re: [R] exact trend test (enumerate all possible contingency tables with fixed row and column margins)

y : > You received a number of suggestions about where to look and packages that > might be suitable. Did you do that? If you did which ones did you look at > and why did you reject them? > > > On 07/01/2016 16:29, li li wrote: > >> Thanks for all the reply. Below is the data in

Re: [R] Fwd: exact trend test (enumerate all possible contingency tables with fixed row and column margins)

ble contingency > tables with fixed row and column margins) > To: li li > > > I do not know whether there is any package to do what you want. > > I **do** know that the algorithms required to do this are very > sophisticated and that with more than a few groups, all possibl

Re: [R] Fwd: exact trend test (enumerate all possible contingency tables with fixed row and column margins)

As a follow up, I found out the proc freq in SAS can perform the exact permutation trend test. Hanna 2016-01-08 9:30 GMT-05:00 li li : > Thanks Bert. > > > 2016-01-07 13:39 GMT-05:00 Bert Gunter : > >> Sorry -- neglected to reply to the list. -- Bert >> >> &

[R] permTREND function in package "perm"

Hi all, I am trying to figure out how to use the function permTREND correctly in package "perm". There does not seem to be any examples given for this function. In the help file, it says the following: ## Default S3 method: permTREND(x, y, alternative = c("two.sided", "less", "greater"), exact =

[R] Logical operator in R

Hi all, I encountered the following strange phenomenon. For some reason, the obs_p[1] and res1$st_p[89] have the same value but when I run "==", it returns FALSE. Can anyone help give some explanation on this? Thanks very much! Hanna > obs_p[1] [1] 0.002201438 > res1$st_p[89] [1] 0.0022014

Re: [R] Logical operator in R

I see. Thanks! 2016-01-22 10:57 GMT-05:00 Rmh : > FAQ 7.31 > > in this case subtract the two numbers and see that > they differ by about 1e-16 > > Sent from my iPhone > > > On Jan 22, 2016, at 10:46, li li wrote: > > > > Hi all, > > I encountered t

[R] Error because of large dimension

Hi all, I am doing some calculation with very large dimension. I need to create a matrix with three columns and a very large number of rows (3195*1290*495*35*35*35*15=1.312083e+15) i n order to allocate calculation result from a for loop. R does not allow me to create such a matrix because of the

Re: [R] Error because of large dimension

GiB = >> 28640 TiB = 28 PiB. Storing such a large matrix even on file is not >> possible. >> >> In other words, you need to figure out how to approach your original >> problem in a different way. >> >> /Henrik >> >> On Sun, Jan 24, 2016 at 8

[R] drm function in drc package for fitting dose response curve

Hi all, When using drm function in drc package, we can fit several does response curves simultaneously or fitting each curve separately. If we fit simultaneous curves without any constraining condition, for example, parallelism condition, will we get the same results as the results obtained by fi

[R] ancova for dependent response

Hi all, Is there an R function which can handles dependent response in Analysis of covariance model. The dependence structure is known and is there to account for it in ANCOVA analysis in R? Thanks. Hanna [[alternative HTML version deleted]] _

[R] weighted sum of independent chi square random variables

Hi all, Is there a function in R that can calculate the quantiles or percentiles for the weighted sum of independent chi square random variables. I found a few functions for calculating probability distribution function for such random variables (e.g. pchisqsum..), but can not find any function

[R] nested structure for Ancova

Hi All, I have a dataset which contains 4 variables: area, group, time, y, Area is a factor that has two levels A and B, group is a factor that is nested within area. There are four groups within each area. y is the response variable, and time refers to different days. Below is the how data lo

Re: [R] nested structure for Ancova

nd area:group and Residuals add up. > > > On Sun, Mar 12, 2017 at 10:39 li li wrote: > >> Hi All, >> I have a dataset which contains 4 variables: area, group, time, y, >> Area is a factor that has two levels A and B, group is a factor that is >> nested w

[R] Test individual slope for each factor level in ANCOVA

Hi all, Consider the data set where there are a continuous response variable, a continuous predictor "weeks" and a categorical variable "region" with five levels "a", "b", "c", "d", "e". I fit the ANCOVA model as follows. Here the reference level is region "a" and there are 4 dummy variables.

Re: [R] Test individual slope for each factor level in ANCOVA

- > Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 > > - - - snip - - - > > In this case, there are just two levels for partner.status, but for a > multi-level factor you can simply perform more than one test. > > > I hope this helps, > > John > >

[R] standard error for regression coefficients corresponding to factor levels

Hi all, I have the following data called "data1". After fitting the ancova model with different slopes and intercepts for each region, I calculated the regression coefficients and the corresponding standard error. The standard error (for intercept or for slope) are all the same for different regi

Re: [R] Is R's fast fourier transform function different from "fft2" in Matlab?

Thanks for both replies. Then I found the "ifft2" from Matlab gives different result from "fft( , inverse=T)" from R. An example: in R: > temp <- matrix(c(1,4,2, 20), nrow=2) > fft(temp) [,1] [,2] [1,] 27+0i -17+0i [2,] -21+0i 15+0i > fft(temp,inverse=T) [,1] [,2] [1,] 27+0i -

[R] Simplify the code

Hi all, I wrote the following code for the function that is attached. I want to simply the code. Can some one give me some help? The function is attached and the following is the code: ## f1 is the function calculate the valu

[R] histogam plots

Hi all, I want to draw a histgram for each row of a matrix and compare them. However the plot I got does not have the same y range and x range, which makes it difficult to make the comparison. Is there a easy way to fix the x range and y range in a xy plot for several plots, instead of specifyi

Re: [R] histogam plots

ttice) > histogram(~w|r, W) > > Identical axis systems will be used for all panels. > > > Bill Venables > CSIRO/CMIS Cleveland Laboratories > > > -Original Message- > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > On Behalf O

Re: [R] Simplify the code

Thank you for your reply. I was just trying to make the code less complicated. I got your point. Thanks! 2010/2/11 Petr PIKAL > Hi > > I did not see the answer yet, probably nobody is much interested in such > function, so I will try. > > r-help-boun...@r-project.org napsal dne 10.02.2010 21:57:

Re: [R] histogam plots

Many thanks for your kind help!! Li 2010/2/11 Dieter Menne > > > Dennis Murphy wrote: > > > > > > ## Histogram + density plots in lattice and ggplot2 > > > > > > There was a typo in your example (should be dd <- data.frame), but anyway. > ggplot2 needs 10 times as long. While Bill Ven

  1   2   3   >