Displaying 20 results from an estimated 50000 matches similar to: "Problem of a function I wrote"
2017 Jun 01
0
Problem of a function I wrote
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
2017 Jun 01
1
Problem of a function I wrote
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
2017 Aug 16
4
{nlme} Question about modeling Level two heteroscedasticity in HLM
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
2010 Jun 15
3
Problem about zero
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
2009 Mar 09
1
How to optimize a matrix
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
2017 Aug 16
0
{nlme} Question about modeling Level two heteroscedasticity in HLM
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
2017 Aug 16
0
{nlme} Question about modeling Level two heteroscedasticity in HLM
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
2010 Feb 28
2
Calling SAS from R
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?
2010 Nov 05
1
Detect the Warning Message
Dear all,
I've written a function and repeated it for 5000 times with loops with
different value, and the messages returned are the output I set and 15
warnings.
I would like to trace the warnings by stopping the loop when warning came
out.
Does anyone know how to make it?
Thanks a lot for your help.
Yen
[[alternative HTML version deleted]]
2023 Feb 12
2
Removing variables from data frame with a wile card
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
2024 Feb 29
2
Initializing vector and matrices
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
2024 Feb 19
5
Looping
I need to read csv files repeatedly, named data1.csv, data2.csv,? data24.csv, 24 altogether. That is,
data<-read.csv(?data1.csv?)
?
data<-read.csv(?data24.csv?)
?
Is there a way to do this in a loop? Thank you.
Steven from iPhone
[[alternative HTML version deleted]]
2024 Aug 09
3
If loop
"Or use <<- assignment I think. (I usually return, but return can only
return one object and I think you want two or more"
You can return any number of objects by putting them in a list and
returning the list.
Use of "<<-" is rarely a good idea in R.
-- Bert
On Fri, Aug 9, 2024 at 1:53?AM CALUM POLWART <polc1410 at gmail.com> wrote:
>
> OK. The fact
2023 Aug 06
2
Stacking matrix columns
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,]
2023 Jan 14
2
Removing variables from data frame with a wile card
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:
>
2024 Aug 11
3
Printing
Thanks. Will try it.
Have not tried it but I think the following may work:
out$results<-NULL
out$results$ei<-ap
out$results$vi<-vap
All I need is printing by returning out (unless I turn it off). And,
retrieve ap and vap as needed as shown above. Guess I need to read more
about invisible.
On 8/11/2024 10:09 PM, Rui Barradas wrote:
> ?s 09:51 de 11/08/2024, Steven Yen escreveu:
2024 Feb 29
1
Initializing vector and matrices
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,
2023 Aug 06
1
Stacking matrix columns
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
2024 Aug 09
2
If loop
OK. The fact it's in a function is making things clearer.
Are you trying to update the values of an object from within the function,
and have them available outside the function. I don't speak functional
programming articulately enough but basically
v <- 1
funA <- function() {
v <- v+1
}
funA()
cat (v)
# 1
You either return the v from the function so
funB <- function() {
2018 Mar 14
1
Documenting R package with Rd file
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