Edit - Stack Overflow

Skip to main content

1. Home
2. Questions
3. AI Assist
4. Tags
5. Challenges
6. Chat
7. Articles
8. Users
9. Companies
11. Communities for your favorite technologies. Explore all Collectives
Stack Internal

Stack Overflow for Teams is now called Stack Internal. Bring the best of human thought and AI automation together at your work.
Try for free Learn more
Bring the best of human thought and AI automation together at your work. Learn more

You are not logged in. Your edit will be placed in a queue until it is peer reviewed.

We welcome edits that make the post easier to understand and more valuable for readers. Because community members review edits, please try to make the post substantially better than how you found it, for example, by fixing grammar or adding additional resources and hyperlinks.

Title*

Body*

I am writing to implement a Simplex regression in this file (`simplex_reg.cpp`):

```
// Simplex Regression
#include <TMB.hpp>
template<class Type>
Type objective_function<Type>::operator() ()
{
 // Data provided from user
 DATA_VECTOR(Y);
 DATA_MATRIX(X);
 DATA_STRING(log);
 // Parameter of the model
 PARAMETER_VECTOR(Beta);
 // Variables used in the programming;
 int n = Y.size();
 vector<Type> eta(n); 
 vector<Type> d0(n); 
 vector<Type> mu(n); 
 Type sd = 0.0;
 Type pi = 3.141593;
 Type nll = 0.0;
 
 // Linear predictor
 eta = X*Beta;
 // logit is link function, and inv.logit is its inverse and is being applied
 mu = 1/(1+exp(-eta));
 // Density
 d0 = pow((Y-mu),2)/(Y*(1-Y)*pow(mu,2)*pow((1-mu),2));
 sd = sum(d0)/n;
 d0 = -0.5*log(2*pi) -0.5*log(sd) - (3/2)*log(Y*(1-Y)) - (1/(2*sd))*d0;
 
 if(log == "false"){
 d0 = exp(d0);
 }
 nll = -sum(d0);
 return nll;
}
```
However, it return an errors when I do `compile('simplex_reg.cpp', "-O0 -g")`. When I use `gdbsource()` as below it returns `Program returned without errors`:

```
library(TMB)
setwd("~/Google Drive/Mestrado/dissertacao/TMB")
compile('simplex_reg.cpp', "-O0 -g")
gdbsource("Teste.R")
```

I have already read the [Compilation and run time errors](https://kaskr.github.io/adcomp/_book/Errors.html), but I couldn't understand how `gdbsource`really works.

So, How do I run the debug/gdbsource into "simplex_reg.cpp" file to find my mistake?

Thanks in advance

Draft saved

Draft discarded

Tags

Edit Summary*

Cancel

6

building in terminal ubuntu 18.04 I get the following: -I/home/chris/R-361/R-3.6.1/src/include simplex_reg.cpp: In instantiation of ‘Type objective_function<Type>::operator()() [with Type = CppAD::AD<double>]’: tmb_core.hpp:1135:7: required from here simplex_reg.cpp:28:16: error: no match for call to ‘(std::__cxx11::string {aka std::__cxx11::basic_string<char>}) (CppAD::AD<double>)’ d0 = -0.5*log(2*pi) -0.5*log(sd) - (3/2)*log(Y*(1-Y)) - (1/(2*sd))*d0; as first error. And same within R using your > compile('simplex_reg.cpp', '-O0 -ggdb'). Is as far as I've gotten.

Chris
– Chris

2020年04月13日 21:46:07 +00:00
Commented Apr 13, 2020 at 21:46
1

Tks Chris! Did you look through all the compilation messages to find this message? (I already solved this problem, but I want to know for the next day). When you said "building in terminal ubuntu", which commands did you use? If you reply in the answer section I will give you the credits, thanks in advance!

Guilherme Parreira
– Guilherme Parreira

2020年04月14日 12:24:29 +00:00
Commented Apr 14, 2020 at 12:24
I haven't gotten near solving it, you should post your solution to guide the blind.

Chris
– Chris

2020年04月14日 15:47:46 +00:00
Commented Apr 14, 2020 at 15:47
I got the solution via commenting the code and running smaller parts of it. But I didn't understand the workflow of gdbsource(). For example, I didn't manage to find the error message that you found, this is the path that I need.

Guilherme Parreira
– Guilherme Parreira

2020年04月14日 17:18:16 +00:00
Commented Apr 14, 2020 at 17:18
gdbsource() only comes into play after simplex_reg.o or .so has been fully compiled, like working with with an R script in debug. You can str() objects and values and watch execution, but without the .o, nothing happens. This is what is misleading about the Program exited without error after running >gdbsource(simplex_reg.cpp), because the message refers to gdb, not simplex_reg.cpp. While nonsensical, I ran it many times. My approach to c++ is akin to 1000 monkeys with typewriters attempting Rilke poems. And I'd lost my path to your Teste.R, alas.

Chris
– Chris

2020年04月14日 21:28:52 +00:00
Commented Apr 14, 2020 at 21:28

| Show 1 more comment

How to Edit

Correct minor typos or mistakes
Clarify meaning without changing it
Add related resources or links
Always respect the author’s intent
Don’t use edits to reply to the author

How to Format

create code fences with backticks ` or tildes ~
```
like so
```
add language identifier to highlight code
```python
def function(foo):
print(foo)
```
put returns between paragraphs
for linebreak add 2 spaces at end
_italic_ or **bold**
indent code by 4 spaces
backtick escapes `like _so_`
quote by placing > at start of line
to make links (use https whenever possible)

<https://example.com>

[example](https://example.com)

<a href="https://example.com">example</a>

formatting help »
asking help »

How to Tag

A tag is a keyword or label that categorizes your question with other, similar questions. Choose one or more (up to 5) tags that will help answerers to find and interpret your question.

complete the sentence: my question is about...
use tags that describe things or concepts that are essential, not incidental to your question
favor using existing popular tags
read the descriptions that appear below the tag

If your question is primarily about a topic for which you can't find a tag:

combine multiple words into single-words with hyphens (e.g. python-3.x), up to a maximum of 35 characters
creating new tags is a privilege; if you can't yet create a tag you need, then post this question without it, then ask the community to create it for you

popular tags »

lang-r