) * (0 - (-Inf)) / (1000 - (-Inf))
= 1 + 0 * Inf / Inf
= 1 + NA / Inf
= 1 + NA
= NA
Zero times infinity is undefined.
[1] <https://en.m.wikipedia.org/wiki/Linear_equation>
On July 4, 2025 3:50:26 AM PDT, "Timothy Earl (Cefas) via R-help"
wrote:
>Hi,
>
>I
r Environment, Fisheries and
Aquaculture Science may be monitored in line with relevant UK legislation .
<https://www.gov.uk/government/organisations/centre-for-environment-fisheries-and-aquaculture-science/about/personal-information-charter>
______
R-hel
gt;> [4,]4 10
>> [5,]5 11
>> [6,]6 12
>> > b<-matrix(5:12,nrow=4); b
>> [,1] [,2]
>> [1,]59
>> [2,]6 10
>> [3,]7 11
>> [4,] 8 12
>> > cbind.fill(a,b)
>> Error in cbind.fill(
Not an expert, but from my read of the explanation you linked, the .svg
IS the source. You could think of it as the source for .svg render engines,
On 6/27/2025 12:16 PM, Ranjan Maitra via R-help wrote:
Hi,
Is there code for drawing the R logo? https://www.r-project.org/logo/ has the
logo
Hi,
Is there code for drawing the R logo? https://www.r-project.org/logo/ has the
logo in svg and png formats, and also the terms of the license, but I can not
find code to draw this.
Many thanks and best wishes,
Ranjan
__
R-help@r-project.org
gt;rare either.
>
>If R *did* guess from the first line, then the second line would not disappear.
>
>On Thu, 26 Jun 2025 at 04:58, Jeff Newmiller via R-help
> wrote:
>>
>> Indeed, the documentation appears to be incorrect... it is literally
>> impossible for the fun
rpreted as LF, LF, LF.
I believe that accounts for the extra new line on MacOS (and Linux) and in
binary mode on Windows, but leaves the disappearing line in text mode on Windows
a bit of a mystery still.
Best,
Mikko
[1]:
https://github.com/wch/r-source/blob/ff4be744953e2702f4194f01f669b7c93b44933e
o accommodate all of these possibilities, throwing an error if any
>> entry does not conform to any of the possibilities that you envisage.
>>
>> cheers,
>>
>> Rolf Turner
>>
>> --
>> Honorary Research Fellow
>> Department of Statistics
>> U
a major addition but often
> > useful.
> >
> > Dennis
> >
> > Dennis Fisher MD
> > P < (The "P Less Than" Company)
> > Phone: 1-415-307-4791
> > www.PLessThan.com
>
>
>
> [[alternative HTML version deleted]]
>
> ___
-Oorspronkelijk bericht-
Van: Enrico Schumann
Verzonden: woensdag 25 juni 2025 16:44
Aan: Heuvel, E.G. van den (Guido)
CC: 'r-help@R-project.org'
Onderwerp: Re: [R] Potential bug in readLines when reading empty lines
[Externe email]
Quoting "Heuvel, E.G. van de
, it appears to not be working
correctly.
-Oorspronkelijk bericht-
Van: Jeff Newmiller
Verzonden: woensdag 25 juni 2025 16:10
Aan: Heuvel, E.G. van den (Guido) ; Heuvel, E.G. van den
(Guido) via R-help ; 'r-help@R-project.org'
Onderwerp: Re: [R] Potential bug in readLines when re
pe of generic line end
>handling algorithm that you describe. However, it appears to not be working
>correctly.
>
>-Oorspronkelijk bericht-
>Van: Jeff Newmiller
>Verzonden: woensdag 25 juni 2025 16:10
>Aan: Heuvel, E.G. van den (Guido) ; Heuvel, E.G. van
>den (Guido
e this and submit a suggestion/patch to R and hope
someone agrees that such a change won't cause more havoc than it avoids for
other users. But that would be unlikely to happen in a timely fashion for your
current needs.
On June 24, 2025 11:59:58 PM PDT, "Heuvel, E.G. van den (Guido) via
ow should I read test.txt in such a way that the
empty 2nd line is left intact?
Best regards,
Guido van den Heuvel
Statistics Netherlands
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/
ice1 = paste0(var, '1_price'),
price2 = paste0(var, '2_price'),
volume1 = paste0(var, '1_volume'),
volume2 = paste0(var, '2_volume')
)
]
--
Best regards,
Ivan
______
R-help@r-project.org mailing list -- To U
fix or avoid it?
Most of our software's users (and I bet R's users as well) use Windows so this
has quite an impact.
Kind regards,
Don
[[alternative HTML version deleted]]
______
R-help@r-project.org mailing list -- To UNSUBSCRIBE and
be
> > done internally within the shell script that speed up many things. The UNIX
> > pipelines helped make fairly complex things seem almost routine to write
> > this way. A new language, perhaps somewhat based on R, might do well if it
> > did not do everything but a
Dear data analyst (or small investor?),
В Fri, 13 Jun 2025 11:49:37 + (UTC)
Small Investor via R-help пишет:
> 1. Version compatibility issues seem to be on the rise. Very
> often, you get the message that package x was built on R version y
> (and thus, won't work in your vers
r.
On June 13, 2025 7:49:37 AM EDT, Small Investor via R-help
wrote:
>Dear R community,
>I have been using R for over 15 years. I want to raise an issue which has been
>haunting me for some time now: It feels as if R is falling apart. I try to
>justify this feeling by providing thre
her people feel the same way, but I think I am
shifting away from R.
yours best,a data analyst dude
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-he
aes(x=x, y=y, colour=w, group=w)) +
> geom_line(linewidth=2) +
> ggtitle("A+B")
>ggplot(df[df$z=="A",], aes(x=x, y=y, colour=w, group=w)) +
> geom_line(linewidth=2) +
> ggtitle("A")
>ggplot(df[df$z=="B",], aes(x=x, y=y, colour=w, group=w
answers
against those given in the book for almost all the problems.
Let me know if you think this repository could be improved in any way.
Thank you to all the folks who helped me better understand statistics
and how to compute them with R.
-Kevin
______
M Brian Smith briansmith199...@gmail.com
>>> wrote:
>>>
>>> > Hi,
>>> >
>>> > Let say I have a range [0, 100]
>>> >
>>> > Now I need to simulate 1000 10 mid-points within the range with
>>> > accuracy upto
ence between any 2
> > consecutive mid-points shall be at-least 5.00.
> >
> > I wonder if there is any Statistical theory available to support this
> > kind of simulation.
> >
> > Alternately, is there any way in R to implement this?
> >
> > _
t;
>> I wonder if there is any Statistical theory available to support this
>> kind of simulation.
>>
>> Alternately, is there any way in R to implement this?
>>
>> __
>> R-help@r-project.org mailing list -- To UNSUBSCR
Could you please help why I am getting NA value?
>
>__
>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 https://www.R-project.org/posting-guide.html
>and provide
ould still have your x and y partition variables from when you
created the data frame, and you should be able to construct the matrix using
the matrix() constructor (not the as.matrix function) from your data frame z
column.
?matrix
On May 31, 2025 4:04:20 AM PDT, ravi via R-help wrote:
>H
phi = 20) :
At least one coordinate must be a matrix
How do I get z in the form of a matrix?
Thanks.
__
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
desperate to get it back.
>
>Hannah
>Sent from my iPhone
>__
>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 https://www.R-project.
gt; Paul Zachos, PhD
>> Director, Research and Evaluation
>> Association for the Cooperative Advancement of Science and Education
>> (ACASE)
>> 110 Spring Street Saratoga Springs, NY 12866 |
>> p...@acase.org | www.acase.org
>>
>>
>>
>&g
or 4
minor 5.0
year 2025
month 05
day 16
svn rev 88216
language R
version.string R version 4.5.0 Patched (2025-05-16 r88216)
nickname How About a Twenty-Six
----
On Tue May20'25 12:15:41AM, Ranjan Maitra via R-help wrote:
> From: Ranjan Maitra via R-help
> Date: Tue, 20 May 2025 00:15:41 -0500
> To: r-help@r-project.org
> Reply-To: Ranjan Maitra
> Subject: Re: [R] how to read a PSB file in R?
>
> On Tue May20'25 05:56:59A
On Tue May20'25 12:15:41AM, Ranjan Maitra via R-help wrote:
> From: Ranjan Maitra via R-help
> Date: Tue, 20 May 2025 00:15:41 -0500
> To: r-help@r-project.org
> Reply-To: Ranjan Maitra
> Subject: Re: [R] how to read a PSB file in R?
>
> On Tue May20'25 05:56:59A
t; 1-x
>
>> On 18.05.2025, at 19:40, paul zachos via R-help wrote:
>>
>> Dear R Community
>>
>> I am an R beginner
>>
>> I have a vector of ‘1’s and ‘0’s
>>
>> x
>> [1] 0 0 1 0 0 1 1 0 1 0 1 1 1 0 1 1 0 0 1 0 1 0 1 0 0 0 0
>> [28] 0
(In hindsight I should have used a `zero.print` value that looked less like a
typo, e.g. `zero.print = "--"` in my example. Sorry about the potential
confusion.)
-Original Message-----
From: R-help On Behalf Of Marttila Mikko via
R-help
Sent: Friday, 23 May 2025 13:21
To: How
place = TRUE)
[1] " "
Warning message:
In .format.zeros(x, zero.print, replace = replace.zero) :
'zero.print' is truncated to fit into formatted zeros; consider 'replace=TRUE'
When zero.print is only one character, it doesn't give a warning.
Tim
>
zero.print = " - ", replace = TRUE)
[1] " "
Warning message:
In .format.zeros(x, zero.print, replace = replace.zero) :
'zero.print' is truncated to fit into formatted zeros; consider 'replace=TRUE'
When zero.print is only one character, it doesn'
ribution, or reliance upon the contents of
this e-mail is strictly prohibited. If you have received this e-mail
transmission in error, please reply to the sender, so that they can arrange for
proper delivery, and then please delete the message from your computer systems.
Thank you.
mericans Independence Center;
>> Senior Statistician University of Maryland Center for Vascular Research;
>>
>> Division of Gerontology and Paliative Care,
>> 10 North Greene Street
>> GRECC (BT/18/GR)
>> Baltimore, MD 21201-1524
>> Cell phone 443-418
Surely doing
y <- 1 - x
would be the simplest way?
Best,
Tom
> -Original Message-
> From: R-help On Behalf Of Rui Barradas
> Sent: 19 May 2025 08:08
> To: paul zachos ; r-help@r-project.org
> Subject: Re: [R] Looking for a function or a set of steps
>
> Às 1
On Tue May20'25 05:56:59AM, Rui Barradas wrote:
> From: Rui Barradas
> Date: Tue, 20 May 2025 05:56:59 +0100
> To: Ranjan Maitra , r-help@r-project.org
> Subject: Re: [R] how to read a PSB file in R?
>
> Às 03:17 de 20/05/2025, Ranjan Maitra via R-help escreveu:
> >
file.
Many thanks and best wishes,
Ranjan
__
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 https://www.R-project.org/posting-guide.html
and provide commented
<https://www.coresystemtrust.org.uk/>
Other work web site: https://www.psyctc.org/psyctc/
Personal site: https://www.psyctc.org/pelerinage2016/
<https://ombook.psyctc.org/book>
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
http
and then in a browser I enter the URL http://localhost:8787/
HTH,
Eric
On Mon, May 19, 2025 at 6:01 PM Chris Evans via R-help
wrote:
I realise that shiny and Rstudio is often seen as OT here but
I have a
strong suspicion that the issue is R or R/OS related so
like to generate a new vector in which the ‘1’s in x become ‘0’s and
the ‘0’s in x become ‘1’s.
How should I go about this?
Thank you,
paz
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
lbuettel
wrote:
On 9 May 2025 at 16:41, Ivan Krylov via R-help wrote:
| This is not guaranteed to work, especially for (but not limited to)
| packages that compile source code and link to various interfaces inside
| R: https://stat.ethz.ch/pipermail/r-help/2025-April/480832.html
|
| This problem
on the
terminal?
What is a simple test to determine that I have the r2u up and running?
Thanks, Ravi
On Wednesday, 14 May 2025 at 13:53:26 CEST, Dirk Eddelbuettel
wrote:
On 9 May 2025 at 16:41, Ivan Krylov via R-help wrote:
| This is not guaranteed to work, especially for (but not limite
Date(0L, origin = "1970-01-01") |> typeof()
#> [1] "double"
```
Is it a bug or the intended behaviour? I don't see any mention of this in the
documentation.
I'm working on R 4.5
Thank you!
Tanguy BARTHELEMY
__
R-h
Thanks to Jeff and Ivan for your helpful comments.
Ravi
On Saturday, 10 May 2025 at 11:37:29 CEST, Ivan Krylov
wrote:
В Thu, 8 May 2025 19:39:40 + (UTC)
ravi via R-help пишет:
> In windows, I add the following line:
> .libPaths(c(“C:/Rownlib”,”C:/R/R-4.5.0/library”))
>
ative numbers into powers or log or sqrt... kind
of things. Read the documentation for your analysis packages carefully.
On May 9, 2025 2:31:29 AM PDT, Luigi Marongiu wrote:
>Dear R-Help,
>I am trying to determine the growth rate of bacteria under specific
>conditions using
>ipo
В Thu, 8 May 2025 19:39:40 + (UTC)
ravi via R-help пишет:
> In windows, I add the following line:
> .libPaths(c(“C:/Rownlib”,”C:/R/R-4.5.0/library”))
> in the file C:/R/R-4.5.0/etc/Rprofile.site
>
> I have my own list of packages in the Rownlib folder and the packages
&
ive())" and
the problem is gone. The reason the "$()" variant worked is probably due
to some subtle differences between the two.
<https://mywiki.wooledge.org/BashFAQ/082>
______
R-help@r-project.org mailing list -- To UNSUBSCRIBE and m
27;s
guidelines should address this?
__
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 https://www.R-project.org/posting-guide.html
and provide commented, mi
s the packages that are
"built-in" to R.
On May 8, 2025 12:39:40 PM PDT, ravi via R-help wrote:
>Hi,
>In Windows, I follow a method with customized library locations that has
>worked for me when upgrading to new R versions. I have not been able to
>follow the same method in U
d also appreciate explanations for why the library
system in ubuntu is structured in the way it is.
Thanks, Ravi
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the po
g does not corrupt your code example.
On May 6, 2025 11:53:42 PM PDT, "Deelen, M. (Mirjam) via R-help"
wrote:
>Hi guys,
>For my MSc. thesis i am using R studio. The goal is for me to merge a couple
>(6) of relatively large datasets (min of 200.000 and max of 2mil rows). I
set. Try merging that where it is easy to check. Then expand to using all
>the data. Your approach should work with a little modification. Alternatively
>make a couple of small fake data sets and play with those. Once you have the
>code correct there should be no problem expanding to the
(so I can drop the
additonal criteria and the final amount of cases, probably, increases).
* Or are there any tips as to figure out where in these steps cases are
lost ?
Thanks!
Mirjam
[[alternative HTML version deleted]]
__
R-help@r-proje
> > they’re more common in teaching. Meanwhile, classic Normal and t
> > methods aren’t going away – they’re fundamental and still useful.
> > Future students will continue to learn both, getting the best of both
> > worlds.
> >
> > Good luck in your studies!
>
; Thank you very much for your time and effort in answering my
> questions.
> I really appreciate the thoughts of the members of this group.
>
> -Kevin
>
>
>
>
______
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
h
what
> assumption or which test is problematic and why?
>
> Tim
__
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 https://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
till useful.
> Future students will continue to learn both, getting the best of both
> worlds.
>
> Good luck in your studies!
> gregg
______
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 https://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
the easiest and quickest part.
Thank you again, Bert, for replying. I always enjoy your contributions
to this group.
-Kevin
______
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 https://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
ago where this was discussed:
>
> https://www.stat.math.ethz.ch/pipermail/r-help/2023-August/477904.html
>
> I generally use:
>
> export OPENBLAS_NUM_THREADS=1
> export MKL_NUM_THREADS=1
>
> since in my experience the biggest performance gains come from switching to
going
away – they’re fundamental and still useful. Future students will continue to
learn both, getting the best of both worlds.
Good luck in your studies!
gregg
On Monday, May 5th, 2025 at 8:17 AM, Kevin Zembower via R-help
wrote:
>
>
> I marked this posting as Off Topic because i
much for your time and effort in answering my questions.
I really appreciate the thoughts of the members of this group.
-Kevin
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
if you have eigenvalues that are almost
> identical, though.)
>
> -pd
>
> > On 4 May 2025, at 17.57, smallepsilon via R-help r-help@r-project.org wrote:
> >
> > I am using MKL with R 4.5.0 on Linux, and eigen() is producing different
> > results with identical calls.
PM, Gregg Powell via R-help
wrote:
>
>
> Hello Christofer,
>
> That’s good information. Thanks for the precision.
> Think I can use this now to work on a script:
>
> Based on what you provided:
> 1. Two distinct sets of coefficients, β¹ and β², each associated with
et.seed(20250504)
Sigma <- rWishart(1, df = n, Sigma = diag(n))[,,1]
e1 <- eigen(Sigma)
e2 <- eigen(Sigma)
identical(e1, e2)
______
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
t; > > > > > nloptr::nloptr() with a manually coded log-likelihood function
> > > > > > > > for the cumulative logit model.
> >
> > > > > > > > Since you need:
> >
> > > > > > > > Coefficient bounds (e.g., lb ≤
t;Is there a vendor who sells a plug-and-play solution to this problem? I am
>guessing we are not unusual in a setup like this, though an upper price
>bound on the software here is of course just the cost of buying a giant
>homogeneous computer or using Amazon resources.
>
>Po
memory and CPU
core limits, but it might be unable to save you from inefficient job
plans. See 'future.batchtools' for an example of an R interface for
cluster job submission systems.
--
Best regards,
Ivan
[1] https://github.com/BOINC/boinc/wiki/BOINC-apps-(introduction)
__
alternate
return values to a local variable along different execution paths and require
all execution paths to reach the end of the function where that variable name
is alone on the last line.
On April 30, 2025 2:15:02 AM PDT, Ralf Goertz via R-help
wrote:
>I made a stupid error w
On Wed, 30 Apr 2025 11:15:02 +0200
Ralf Goertz via R-help wrote:
> If this is not an error what is its purpose?
>From the point of view of the R syntax, everything is an expression.
One of the uses of return() being an expression is base::callCC(). The
ability to use it inside an expr
ariant (which returns the value of
the numerator only) not a syntax error? I would have expected R to
report something like
> Error: unexpected '/' in "return(OR^2+6*OR+1)/"
If this is not an error what is its purpose?
______
R-help@r-p
; > > Equality constraint: sum(β) = C
> >
> > > > > > C. Likelihood
> > > > > > The probability of category k is given by:
> >
> > > > > > P(Y = k) = logistic(θ_k - Xβ) - logistic(θ_{k-1} - Xβ)
> >
> > > > > > Without intercepts, this becomes more like:
> >
> > > > > > P(Y ≤ k) = 1 / (1 + exp(-(c_k - Xβ)))
> >
> > > > > > …where c_k are fixed thresholds.
> >
> > > > > > Implementation using nloptr
> > > > > > This example shows a rough sketch using nloptr, which allows both
> > > > > > equality and bound constraints:
> >
> > > > > > > library(nloptr)
> >
> > > > > > > # Custom negative log-likelihood function
> > > > > > > negLL <- function(beta, X, y, K, cutpoints) {
> > > > > > > eta <- X %*% beta
> > > > > > > loglik <- 0
> > > > > > > for (k in 1:K) {
> > > > > > > pk <- plogis(cutpoints[k] - eta) - plogis(cutpoints[k - 1] - eta)
> > > > > > > loglik <- loglik + sum(log(pk[y == k]))
> > > > > > > }
> > > > > > > return(-loglik)
> > > > > > > }
> >
> > > > > > > # Constraint: sum(beta) = C
> > > > > > > sum_constraint <- function(beta, C) {
> > > > > > > return(sum(beta) - C)
> > > > > > > }
> >
> > > > > > > # Define objective and constraint wrapper
> > > > > > > objective <- function(beta) negLL(beta, X, y, K, cutpoints)
> > > > > > > eq_constraint <- function(beta) sum_constraint(beta, C = 2) #
> > > > > > > example >C
> >
> > > > > > > # Run optimization
> > > > > > > res <- nloptr(
> > > > > > > x0 = rep(0, ncol(X)),
> > > > > > > eval_f = objective,
> > > > > > > lb = lower_bounds,
> > > > > > > ub = upper_bounds,
> > > > > > > eval_g_eq = eq_constraint,
> > > > > > > opts = list(algorithm = "NLOPT_LD_SLSQP", xtol_rel = 1e-8)
> > > > > > > )
> >
> > > > > > The next step would be writing the actual log-likelihood for your
> > > > > > specific problem or verifying constraint implementation.
> >
> > > > > > Manually coding the log-likelihood for an ordinal model is
> > > > > > nontrivial... so a bit of a challenge :)
> >
> > > > > > All the best,
> > > > > > gregg powell
> > > > > > Sierra Vista, AZ
signature.asc
Description: OpenPGP digital signature
__
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 https://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
9.csv.gz") calls file(), which can transparently
switch to a gzfile() connection, which supports concatenated archives,
but gzcon() currently doesn't. Feature request submitted at
<https://bugs.r-project.org/show_bug.cgi?id=18887>.
--
Best regards,
Ivan
_
if there is only one folder for packages under 4.X... where
>any version of R $.X... would load from there the packages?
>
>Thanks
>
>Fer
>
>On 4/25/25 11:16, Peter Dalgaard via R-help wrote:
>> A couple of people have gotten themselves in trouble by copying the contents
.5.0. One of them has tools
in the 4.4 library and the other does not. I have no clue what the difference
might be
--
Peter Dalgaard, Professor,
Center for Statistics, Copenhagen Business School
Solbjerg Plads 3, 2000 Frederiksberg, Denmark
Phone: (+45)38153501
Office: A 4.23
Email: pd@c
nt"), lin)) :
> > > NA/NaN/Inf in foreign function call (arg 5)
> > > ```
> > > How can I properly set this regression model?
> > > Thank you
> >
> > > On Wed, Apr 16, 2025 at 7:08 AM Luigi Marongiu marongiu.lu...@gmail.com
> &g
- function(beta, X, y, K, cutpoints) {
> > > > > eta <- X %*% beta
> > > > > loglik <- 0
> > > > > for (k in 1:K) {
> > > > > pk <- plogis(cutpoints[k] - eta) - plogis(cutpoints[k - 1] - eta)
> > > > > loglik <
& Statistics
> > > University of Melbourne, VIC 3010 Australia
> > > Tel: (+61) 0403 138 955
> > > Email: a...@unimelb.edu.au
> > > Website: https://researchers.ms.unimelb.edu.au/~apro@unimelb/
> > >
> > > I acknowledge the Traditional Ow
; > > eq_constraint <- function(beta) sum_constraint(beta, C = 2) # example >C
> > >
> > > # Run optimization
> > > res <- nloptr(
> > > x0 = rep(0, ncol(X)),
> > > eval_f = objective,
> > > lb = lower_bounds,
> > >
gt;> >eq_constraint <- function(beta) sum_constraint(beta, C = 2) # example >C
>> >
>> ># Run optimization
>> >res <- nloptr(
>> > x0 = rep(0, ncol(X)),
>> > eval_f = objective,
>> > lb = lower_bounds,
>> > ub = uppe
implementation.
Manually coding the log-likelihood for an ordinal model is nontrivial... so a
bit of a challenge :)
All the best,
gregg powell
Sierra Vista, AZ
signature.asc
Description: OpenPGP digital signature
______
R-help@r-project.org mailing
the following error , at the following point in the code ..." )
On 2025-04-16 12:11 p.m., varin sacha via R-help wrote:
> R-experts,
>
> The R script here below.
>
> How to solve my problem ?
>
> Best,
>
>
> ##
> # Load packa
df_sample <- sample_n(base_df, n)
result <- gfci(data = df_sample, test = "kci", score = "sem-bic")
plot(result, main = paste("n =", n))
}
#####
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and mo
y of Melbourne, VIC 3010 Australia
Tel: (+61) 0403 138 955
Email: a...@unimelb.edu.au
Website: https://researchers.ms.unimelb.edu.au/~apro@unimelb/
I acknowledge the Traditional Owners of the land I inhabit, and pay my respects
to their Elders.
On 16 Apr 2025 at 1:01 AM +1000, Gregg Powell via R-
Norm[df$Norm<0] = 0 b_model = glm(Norm ~ Cycles, data=df,
> family=binomial) plot(Norm ~ Cycles, df, main="Normalized view",
> xlab=expression(bold("Time")), ylab=expression(bold("Signal (normalized)")),
> type="p", col="cyan") lines
ng execution
>time? Are you planning on looking at all t and T within a range?
>
>Tim
>
>-Original Message-
>From: R-help On Behalf Of Bert Gunter
>Sent: Monday, April 14, 2025 4:16 PM
>To: Duncan Murdoch
>Cc: r-help@r-project.org
>Subject: Re: [R] Drawing a sa
S3 or S4 structure, please?
>
>I know I will get lots of opinions, but that’s fine.
>
>Thanks,
>Erin
>
>
>Erin Hodgess, PhD
>mailto: erinm.hodg...@gmail.com
>
> [[alternative HTML version deleted]]
>
>__
&g
mats.
In addition: Warning message:
In unzip(file_loc, exdir = tmp) : error 1 in extracting from zip file
Error: Your geometry data download failed. Please try again later or check the
status of the Census Bureau website at https://www2.census.gov/geo/tiger/
_
pose certain constraints on each of the
> model parameters, based on some non-sample information.
>
> Is there any R function to estimate model coefficients with imposing
> some unser-defined constraints on the model parameters?
>
> Any pointer will be very helpful.
>
gards,
>
>[[alternative HTML version deleted]]
>
> __
> 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 https://www.R-project.org/pos
all been simulations, bootstraps, randomized distributions, etc.
Thank you, again, Rui and Michael, for your help for me.
-Kevin
On Sat, 2025-03-29 at 16:42 +, Rui Barradas wrote:
> Às 16:09 de 29/03/2025, Kevin Zembower via R-help escreveu:
> > Hello, all,
> >
> > We'r
hough. Tastes
> will
> vary.
>
> In dplyr style, add a column and use it in slice:
>
> df <- expand.grid(year=c(1980, 2010), DP=c("Approve", "Disapprove"))
> %>%
> mutate(freq = c(660, 640, 1000-660, 1000-640)) %>%
> slice(rep(1:n
> On Mar 29, 2025, at 9:59 AM, Kevin Zembower via R-help
> wrote:
>
> Hi, Rui and Michael, thank you both for replying.
>
> Yeah, I'm not supposed to know about Chi-squared yet. So far, all of
> our work with hypothesis tests has involved creating the sample dat
ble to do something like this:
=
(df <- data.frame(group1count = 660, #Or, group1prop = 0.66
group1samplesize = 1000,
group2count = 640, #Or, group2prop = 0.64
group2samplesize = 1000))
=
Am I overlooking a way to
quot;female","male","female","male","male")
color=c("red","blue","green","green","blue","red","blue","green","red","blue")
data<- data.frame(statut, sex,color)
rvations. I checked the
>vignettes but did not see this incorporated.
>
>Mariam.
>
> [[alternative HTML version deleted]]
>
>__
>R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>https://stat.ethz.ch/ma
1 - 100 of 3504 matches
Mail list logo
