similar to: Detect the Warning Message

Hello everyone, I have been working on a code which simply repeatedly appends a number into a vector and write a file. However, it could not be properly implemented when I use it. It works when I run it line by line. I wonder what is the problem and I appreciate anyone who is willing to help. The function and the example code is attached. Any advice is appreciated! Best, Yen

Hello everyone, I know where is wrong. I forget to specify the parameters in my function. Thank you for anyone who was trying to help me! Best, Yen ?? b88207001 at ntu.edu.tw: > Hello everyone, > > It seems that I was not successfully attached the code. Here is the > code. I appreciate any help! > > Best, > Yen > > ?? b88207001 at ntu.edu.tw: > >> Hello

Hello everyone, It seems that I was not successfully attached the code. Here is the code. I appreciate any help! Best, Yen ?? b88207001 at ntu.edu.tw: > Hello everyone, > > I have been working on a code which simply repeatedly appends a > number into a vector and write a file. However, it could not be > properly implemented when I use it. It works when I run it line by

Hello dear uesRs, I am working on modeling both level one and level two heteroscedasticity in HLM. In my model, both error variance and variance of random intercept / random slope are affected by some level two variables. I found that nlme is able to model heteroscedasticity. I learned how to use it for level one heteroscedasticity but don't know how to use it to model the level

Hello, everyone, There's a problem about zero in R and I really need your help. I have a vector shown as x=c(0.1819711,0.4811463,0.1935151,0.1433675), The sum of this vector is shown as 1 in R, but when I type 1-sum(x), the value is not zero, but -2.220446e-16. I can accept that this value is quite small and could be seen as zero, but there would be a problem when it's not really

I would like to estimate the sigma matrix of multinormal distribution through ML. But I don't know how to optimize the parameter sigma. Could any one help me? Thank you so much~ Yen Lee

If you don't get a response it is because you did not read the Posting Guide which indicates that the R-sig-ME mailing list is where this question would have been on-topic. -- Sent from my phone. Please excuse my brevity. On August 16, 2017 6:17:03 AM PDT, b88207001 at ntu.edu.tw wrote: >Hello dear uesRs, > >I am working on modeling both level one and level two

A better place for this post would be on R's mixed models list: r-sig-mixed-models . 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, Aug 16, 2017 at 6:17 AM, <b88207001 at ntu.edu.tw> wrote: > Hello dear

I'm new to post in R-help and my native language is not English. I apologize if my sentence is not fluent to read. I am doing a simulation study and I need to execute SAS and read a SAS code in R. I try the following code but it doesn't work. system('"c:\\program files\\SAS\\SAS 9.1\\sas.exe" "c:\\syntax.sas"') can anyone give me some help with this?

x["V2"] is more efficient than using drop=FALSE, and perfectly normal syntax (data frames are lists of columns). I would ignore the naysayers, or put a comment in if you want to accelerate their uptake. As I understand it, one of the main reasons tibbles exist is because of drop=TRUE. List-slice (single-dimension) indexing works equally well with both standard and tibble types of data

You could declare a matrix much larger than you intend to use. This works with a few megabytes of data. It is not very efficient, so scaling up may become a problem. m22 <- matrix(NA, 1:600000, ncol=6) It does not work to add a new column to the matrix, as in you get an error if you try m22[ , 7] but convert to data frame and add a column m23 <- data.frame(m22) m23$x7 <- 12 The only

Your desire is not unusual among novices... but it is really not a good idea for your function to be making those decisions. Look at how R does things: The lm function prints nothing... it returns an object containing the result of a linear regression. If you happen to call it directly from the R command prompt and don't assign it to a variable, then the command interpreter notices that

I just like the subroutine to spit out results (Mean, Std.dev, etc.) and also be able to access the results for further processing, i.e., v$Mean v$Std.dev On 3/26/2024 11:24 AM, Richard O'Keefe wrote: > Not clear what you mean by "saved". > If you call a function and the result is printed, the result is > remembered for a wee while in > the variable .Last.value, so

Just FYI, the R interpreter typically saves the last value returned briefly in a variable called .Last.value that can be accessed before you do anything else. > sin(.5) [1] 0.4794255 > temp <- .Last.value > print(temp) [1] 0.4794255 > sin(.666) [1] 0.6178457 > .Last.value [1] 0.6178457 > temp [1] 0.4794255 > invisible(sin(0.2)) > .Last.value [1] 0.1986693 So perhaps if

You could also do dim(x) <- c(length(x), 1) On Sat, Aug 5, 2023, 20:12 Steven Yen <styen at ntu.edu.tw> wrote: > I wish to stack columns of a matrix into one column. The following > matrix command does it. Any other ways? Thanks. > > > x<-matrix(1:20,5,4) > > x > [,1] [,2] [,3] [,4] > [1,] 1 6 11 16 > [2,] 2 7 12 17 > [3,]

Thanks to all. Very helpful. Steven from iPhone > On Jan 14, 2023, at 3:08 PM, Andrew Simmons <akwsimmo at gmail.com> wrote: > > ?You'll want to use grep() or grepl(). By default, grep() uses extended > regular expressions to find matches, but you can also use perl regular > expressions and globbing (after converting to a regular expression). > For example: >

I have trouble documenting an R package. In my .Rd file (sixth line below), I have uhat<-m%*%y but when the package is built (successfully), the matrix multiplication part does not show up in the documentation. The line become (missing %*% y) uhat<-m === \examples{ x<-c(1,2,3,4,5) y<-c(1,1,2,2,4) x<-cbind(1,x) m<-mmat(x) uhat<-m%*%y dstat(uhat) } ?-- styen at ntu.edu.tw

Eric, I am not sure your solution is particularly economical albeit it works for arbitrary arrays of any dimension, presumably. But it seems to involve converting a matrix to a tensor just to undo it back to a vector. Other solutions offered here, simply manipulate the dim attribute of the data structure. Of course, the OP may have uses in mind which the package might make easier. We often get

Hi all, I tried to do a variance components using nlme, but I got the following warning mesage ##################################################### not meaningful for factors in: Ops.factor(y[revOrder], Fitted) ###################################################### Can someone point out what is the meaning of this warning message? I tried to look at Ops.factor, but I don't

x <- numeric(0) for (...) { x[length(x)+1] <- ... } works. You can build a matrix by building a vector one element at a time this way, and then reshaping it at the end. That only works if you don't need it to be a matrix at all times. Another approach is to build a list of rows. It's not a matrix, but a list of rows can be a *ragged* matrix with rows of varying length. On Wed,